Sonia Ouadi
Nicolas Sierro
Felix Kessler
Nikolai V. Ivanov- 1Faculty of Sciences, Laboratory of Plant Physiology, University of Neuchâtel, Neuchâtel, Switzerland
- 2Philip Morris International R&D, Philip Morris Products S.A., Neuchâtel, Switzerland
Syzygium is a large and diverse tree genus in the Myrtaceae family. Genome assemblies for clove (Syzygium aromaticum, 370 Mb) and sea apple (Syzygium grande, 405 Mb) provided the first insights into the genomic features and evolution of the Syzygium genus. Here, we present additional de novo chromosome-scale genome assemblies for Syzygium malaccense, Syzygium aqueum, Syzygium jambos, and Syzygium syzygioides. Genome profiling analyses show that S. malaccense, like S. aromaticum and S. grande, is diploid (2n = 2x = 22), while the S. aqueum, S. jambos, and S. syzygioides specimens are autotetraploid (2n = 4x = 44). The genome assemblies of S. malaccense (430 Mb), S. aqueum (392 Mb), S. jambos (426 Mb), and S. syzygioides (431 Mb) are highly complete (BUSCO scores of 98%). Comparative genomics analyses showed conserved organization of the 11 chromosomes with S. aromaticum and S. grande, and revealed species-specific evolutionary dynamics of the long terminal repeat retrotransposon elements belonging to the Gypsy and Copia lineages. This set of Syzygium genomes is a valuable resource for future structural and functional comparative genomic studies on Myrtaceae species.
1 Introduction
Syzygium is the largest tree genus with about 1,200 species naturally occurring from the Old World tropics and subtropics to the Pacific (POWO, 2023; Craven and Biffin, 2010; Beech et al., 2017). In addition to their ecological importance, the genus includes several species grown for their edible fruit, medicinal properties, timber, and for the horticulture industry (e.g., S. malaccense, S. aqueum, S. jambos, and S. cumini), the most important economically being the clove tree (S. aromaticum) (Parnell et al., 2007; Nurdjannah and Bermawie, 2012; Nair, 2017; Cock and Cheesman, 2018).
The Syzygium genus belongs to the Myrtaceae family—the eighth largest family of flowering plants—and includes economically important species such as eucalyptus, myrtle, and guava (Grattapaglia et al., 2012; Christenhusz and Byng, 2016; Saber et al., 2023). Although the majority of species of the Myrtaceae family are diploids (2n = 22) with small to intermediate genome sizes (234–1785 Mb), occasional polyploids derived from the most conserved chromosome number x = 11 were also reported (e.g., within the Eugenia, Syzygium, and Psidium genera) (Wilson, 2010; Grattapaglia et al., 2012; Tuler et al., 2019; Pellicer and Leitch, 2020; Machado and Forni-Martins, 2022). The Eucalyptus grandis genome was released in 2014 as the first reference genome for the Myrtales order and the Myrtaceae family (Myburg et al., 2014). New chromosome-scale assemblies were subsequently published, enabling comparative genomics analyses within the family. Published chromosome-scale genome assemblies for the Myrtaceae currently represent major tribes of the family: Eucalypteae (Eucalyptus grandis, Corymbia citriodora, Eucalyptus urophylla × Eucalyptus grandis), Leptospermeae (Leptospermum scoparium), Myrteae (Psidium guajava, Rhodomyrtus tomentosa), Metrosidereae (Metrosideros polymorpha), Melaleuceae (Melaleuca alternifolia), and Syzygieae (S. aromaticum, Syzygium grande). These assemblies were generated from diploid specimens, and their size ranged from 297 Mb to 690 Mb (Myburg et al., 2014; Izuno et al., 2019; Thrimawithana et al., 2019; Feng et al., 2021; Healey et al., 2021; Low et al., 2022; Ouadi et al., 2022; Zheng et al., 2022; Li et al., 2023; Shen et al., 2023).
The clove (S. aromaticum (L.) Merr. & L.M. Perry) and sea apple (S. grande) genomes were constructed using a combination of Oxford Nanopore Technologies long-reads and Illumina short-reads and anchored on 11 chromosomes using Hi-C technologies (Low et al., 2022; Ouadi et al., 2022). The sea apple genome assembly (405 Mb), 182 re-sequenced Syzygium species and 58 re-sequenced unidentified taxa were used to generate whole genome-level phylogenies of the Syzygium genus, thus providing new insights into the infrageneric classification of Syzygium, as well as into the genus diversification patterns and their drivers. The clove genome assembly (370 Mb) was exploited to investigate the genetic basis of the biosynthesis of eugenol, the major biocompound of clove products (Kamatou et al., 2012; Otunola, 2022). To provide insights into the clove genome evolution, comparative genomics analyses were also performed between S. aromaticum and E. grandis. The synteny analysis performed between these two Myrtaceae species’ genomes assemblies revealed good genome structure conservation. The structures of chromosomes 1, 3, 5, and 7 were found to be highly conserved between E. grandis and S. aromaticum, and 10 intrachromosomal rearrangements occurring on the 7 other chromosomes were observed (chromosomes 2, 4, 6, 8, 9, 10, and 11). Interestingly, the intrachromosomal rearrangements detected between the two eucalypt species, E. grandis and C. citriodora, were located on the same seven chromosomes (Butler et al., 2017; Healey et al., 2021). Long terminal repeat retrotransposons (LTR-RTs) are transposable elements (TEs) that move through the genome via a copy-and-paste mechanism using an RNA intermediate. They are considered the most abundant TE component in plant genomes and important drivers of genome size variation and diversification (Wicker et al., 2007; Zhou et al., 2021). Comparing the LTR-RTs repertoires of S. aromaticum and E. grandis revealed a differential accumulation of the LTR-RTs belonging to the superfamilies Copia and Gypsy between the two species. In S. aromaticum genome assembly, the LTR-RTs belonging to the Gypsy superfamily were more abundant than those belonging to the Copia superfamily. In contrast, a higher number of LTR-RTs Copia versus Gypsy was found in the E. grandis genome assembly.
No infrageneric comparison of chromosome-scale assemblies has been performed for the Syzygium genus. To further investigate the evolution of the genome architecture of Syzygium species and verify whether the rearrangements found between S. aromaticum and E. grandis chromosomes were the consequences of evolutionary events or due to sequencing and assembly artifacts, we generated additional chromosome-scale genome assemblies for Syzygium malaccense (L.) Merr. & L.M. Perry, Syzygium aqueum (Burm.f.) Alston, Syzygium jambos (L.) Alston, and Syzygium syzygioides (Miq.) Merr. & L.M. Perry. Like S. aromaticum and S. grande, the four species belong to the subgenus Syzygium, the largest of the five Syzygium subgenera for which the crown age was estimated at 9.4 Mya by Low et al. (Low et al., 2022). Previous karyotype studies indicated that S. malaccense is a diploid with 2n = 22 chromosomes (Pedrosa et al., 1999) and that S. jambos is a tetraploid (2n = 44); however, different chromosome numbers were also reported in the literature for the species (2n = 28, 33, 46, ~54, 66) (Van Lingen, 1991; Oginuma et al., 1993). The chromosomal numbers reported in the literature indicate that S. aqueum is also a tetraploid (2n = 44) (Panggabean, 1991).
Here, we describe the de novo assembly and annotation for S. malaccense, S. aqueum, S. jambos, and S. syzygioides. To enable subsequent comparative genomic analyses, the four genomes consisting of monoploid consensus (11 chromosomes and unplaced sequences) were generated to achieve the same level of quality for the four species’ genome assemblies and comparable to those of published chromosome-scale assemblies of their Myrtaceae relatives. Then, we compared the genome architecture of the four newly Syzygium assembled genomes with those of S. aromaticum and S. grande and their genome features (gene sets and LTR-RTs repertoires) with those of S. aromaticum to investigate genomic evolution from their common ancestors.
2 Materials and methods
2.1 Biological materials
The S. malaccense, S. aqueum, S. jambos, and S. syzygioides genome assemblies were generated from trees growing in the Masoala Hall of the Zurich Zoo in Switzerland. Voucher specimens were deposited in the Zürich herbarium (S. malaccense (ZT-00170996), S. aqueum (ZT-00170994), S. jambos (ZT-00170999), and S. syzygioides (ZT-00170991)). Samples collected from the trees were stored at -80°C until nucleic acid extraction.
2.2 DNA and RNA isolation
High-molecular-weight genomic DNA was isolated from frozen leaves using the “ONT high-molecular-weight gDNA extraction from plant leaves” protocol (Oxford Nanopore Technologies, Oxford, UK). Following the extraction, we performed a size selection step using the Circulomics Nanobind Plant Nuclei Big DNA Kit from PacBio (Menlo Park, CA, USA). (NB-900-801-001).
Total RNA from S. malaccense, S. aqueum, S. jambos, and S. syzygioides were isolated in triplicate from whole leaves (young and mature), lamina of mature leaves, and stems. Total RNA was also isolated in triplicate from S. syzygioides’ buds (in the fruiting stage) and S. jambos’ buds (before and after flowering) and flowers.
Total RNA was extracted from frozen powder using Ambion PureLink Plant RNA Reagent (Ambion by Life Technologies, Carlsbad, CA, USA). The concentration and quality of the total RNA were assessed with an Agilent Bioanalyzer using the Agilent RNA 6000 Nano Kit (Agilent, Santa Clara, CA, USA).
2.3 Illumina sequencing library preparation and sequencing
DNAseq libraries were prepared from total gDNA using the Celero PCR workflow with an enzymatic fragmentation kit from Tecan (Männedorf, Switzerland). DNAseq libraries were loaded on an Illumina S2 flow cell and sequenced on the Illumina Novaseq 6000 instrument (Illumina, San Diego, CA, USA) as 2 x 151 bp paired-end reads.
Hi-C libraries were prepared from 0.2 g of frozen leaves using the Proximo Hi-C Kit following the manufacturer’s instructions (Phase Genomics, Seattle, WA, USA) and sequenced on an Illumina HiSeq 4000 instrument (Illumina) as 2 x 151 bp paired-end reads.
mRNA stranded libraries were prepared from 500 ng of total RNA using the Tecan Universal Plus mRNA-Seq library preparation kit with NuQuant® and sequenced on an Illumina HiSeq 4000 instrument as 2 x 151 bp paired-end reads.
Illumina raw reads generated from DNAseq libraries and Hi-C libraries were cleaned using fastp 0.23.2 (--length_required 75 --low_complexity_filter) (Chen et al., 2018).10.1038/s41597-021-00968-x
2.4 ONT sequencing library preparation and sequencing
Sequencing libraries were generated from high-molecular-weight gDNA and prepared for sequencing on PromethION flow cells (FLO-R0002) by using the ligation sequencing (SQK-LSK109) and flow cell priming (EXP-FLP002) kits (Oxford Nanopore Technologies, Oxford, UK). The base calling was performed by using Guppy 6.1.1 and the super accuracy plant model. Raw ONT reads were cleaned using seqkit 2.2.0 (--min-qual 9 --min-len 5000) (Shen et al., 2016) to discard reads shorter than 5,000 bp or with quality scores lower than 9.
2.5 Genome profiling
Cleaned Illumina paired-end reads from DNAseq libraries were analyzed by GenomeScope 2.0 and smudgeplot 0.2.4 to estimate the genome size, percentage of heterozygosity, and the ploidy level using a k-mer size equal to 21 (Ranallo-Benavidez et al., 2020).
2.6 De Novo genome assembly
ONT cleaned reads were corrected with fmlrc2 0.1.7 (--cache_size 13 –K 21 59 79) (Mak et al., 2023) using cleaned Illumina paired-end short-reads from DNAseq libraries. The corrected ONT reads were then assembled using flye 2.9 (--read-error 0.01 --nano-hq) (Kolmogorov et al., 2019) and iteratively polished with ntedit 1.3.5 (-m 2 -i 3 -d 3 -X 0.5 -Y 0.5) using kmer profiles created with nthits 0.0.1 (--solid --outbloom -b 36) for kmers of lengths 60, 50, 40 and 30 (Warren et al., 2019) using Illumina paired-end short reads from DNAseq libraries. Haplotigs were detected and removed from the polished contigs using purge_dups 1.2.5 (Guan et al., 2020) using cutoff of 10, 315 and 645 for S. malaccense, 70, 440 and 960 for S. aqueum, and 60, 410, 960 for S. jambos, 10, 410 and 960 for S. syzygioides.
Cleaned Illumina read pairs generated from Hi-C libraries were mapped to the genomes to remove reads with low mapping scores, duplicated reads, and paired-end reads. Illumina Hi-C read pairs were mapped to the haplotig-purged contigs using minimap2 2.24 (Li, 2018) rather than bwa (Li, 2013) since we noticed that it results in assemblies of equivalent qualities in a shorter time. The scaffolding to a chromosome-scale assembly was performed using yahs 1.1a2 (-r 1000,2000,5000,10000,20000,50000,100000,200000,500000,1000000,2000000,5000000) (Zhou et al., 2022). Hi-C map files were generated with PretextMap 0.1.9 (https://github.com/wtsi-hpag/PretextMap) and used to manually curate the assemblies using PretextView 0.2.5 (https://github.com/wtsi-hpag/PretextView).
The curated genome assemblies were mapped to the S. aromaticum genome (Ouadi et al., 2022) using minimap2 2.24, visualized using a custom R script, and the orientation and names of the chromosomes were set in accordance with those of S. aromaticum. Chromosome-scale assembly completeness was assessed by using the genome evaluation mode of BUSCO 5.4.4 and the eudicots_odb10 lineage dataset (Simão et al., 2015). The QVs of the final assemblies were estimated using yak 0.1 (qv -K 2000000000) with kmer profiles created using yak 0.1 (count -k 31 -K 2000000000 -b37) (Cheng et al., 2021).
2.7 Gene annotation
The Illumina RNAseq reads from S. malaccense, S. aqueum, S. jambos and S. syzygioides as well as those used for the clove genome annotation were cleaned, and overlapping paired-reads were merged using fastp 0.23.2 (--length_required 75 --low_complexity_filter --merge) (Chen et al., 2018) before being mapped as single cDNA reads to the assemblies using minimap2 2.24 (-ax splice:hq -G5K -N50) (Li, 2018). Gene models were then created for each RNASeq sample using scallop 0.10.5 (--min_transcript_coverage 5 --min_single_exon_coverage 50 --min_splice_bundary_hits 5 --min_mapping_quality 0) (Shao and Kingsford, 2017).This approach was used for the annotation of the clove genome, where it was observed to produce better gene models than by directly mapping paired-reads with a dedicated mapper.
To obtain models for genes that are not expressed in the RNAseq samples, the transcripts from S. aromaticum and E. grandis gene annotations were mapped to the assemblies using minimap2 2.24 (-ax splice:hq -I5G -G5K -N50 -uf) (Li, 2018), and gene models created using bedtools 2.30.0 (bamtobed -bed12) (Quinlan and Hall, 2010) and custom gawk scripts to convert the obtained bed file into a gtf file.
The final gene models were obtained by merging the RNAseq, S. aromaticum, and E. grandis gene models using taco 0.7.3 (--gtf-expr-attr TPM --filter-min-expr 10) (Niknafs et al., 2017) and adding coding sequences using Transdecoder 5.5.0 (LongOrfs -S -m 64; Predict --single_best_only --retain_blastp_hits dmd.tsv) (https://github.com/TransDecoder/TransDecoder/wiki), diamond 2.0.15 (blastp --query longest_orfs.pep --db uniref-malvids.dmnd --max-target-seqs 1 --outfmt 6 --evalue 1e-6) (Buchfink et al., 2015) and gffread 0.12.7 (Pertea and Pertea, 2020).
The eudicotyledons portion of UniProt filtered to remove proteins with poor descriptions was used to annotate the gene models with their best hit using diamond 2.0.15 (blastx --query tx.fa --db eudicotyledons.filtered.dmnd --top 10 --min-score 200 --ultra-sensitive --iterate). The illustration of the regions where genes encoding for putative eugenol synthase were predicted was generated using gggenes 0.4.0 (https://github.com/wilkox/gggenes).
2.8 Repeat annotation
Annotation of transposable elements was carried out using TE-greedy-nester 1.0.0 (--discovery_tool LTRharvest) (Lexa et al., 2020), genometools LTRharvest 1.6.2 (Ellinghaus et al., 2008) and TEsorter 1.3.0 (-db rexdb-plant --min-coverage 10 --max-evalue 0.01 --pass2-rule 70-30-80) (Zhang et al., 2022) with REXdb (Neumann et al., 2019). The insertion age of the predicted transposable elements was then calculated as previously reported (Marcon et al., 2015). In addition, Red 2.0 (Girgis, 2015), GRF 1.0 (Shi and Liang, 2019) and cd-hit 4.8.1 (grf-main -i genome.fa -c 1 -o genome.MITE --min_tr 10; cd-hit-est -i genome.MITE/candidate.fasta -o genome.MITE/clusteredCandidate.fasta -c 0.90 -n 5 -d 0 -aL 0.99 -s 0.8 -M 0; grf-mite-cluster -i genome.MITE/clusteredCandidate.fasta.clstr -g genome.fa -o genome.MITE) (Fu et al., 2012), EAHelitron (Hu et al., 2019), and tantan 39 (-f4) (Frith, 2011) were used to predict repeats, Miniature Inverted-repeat Transposable Elements (MITEs), helitron, and tandem repeats, respectively.
2.9 Synteny analyses
Synteny between the Syzygium species was done by pairwise mapping whole genomes using minimap2 2.24 (Li, 2018), identifying structural variants using syri 1.6 (Goel et al., 2019), and plotting syntenic blocks larger than 20 kb using plotsr 0.5.4 (Goel and Schneeberger, 2022).
2.10 Orthologue analyses
Orthologous genes were clustered into HOGs with OrthoFinder 2.5.4 (Emms and Kelly, 2019) using the set of predicted protein sequences from the five species assemblies.
3 Results
3.1 Genome profiling
Smudgeplot and GenomeScope 2.0 were used to perform a genome profiling step using Illumina PE short-reads from DNAseq libraries as input and a K-mer length of 21 bp (Ranallo-Benavidez et al., 2020) (Table 1; Supplementary Table 1; Supplementary Figures 1, 2). The ploidy level predicted by Smudgeplot was in accordance with previous karyotype studies for the studied S. malaccense and S. jambos specimens (Oginuma et al., 1993; Pedrosa et al., 1999). S. malaccense was predicted to be a diploid specimen (2n = 2x = 22) like S. aromaticum and S. grande. The S. aqueum, S. jambos, and S. syzygioides specimens were predicted as being autotetraploid (2n = 4x = 44). The estimated monoploid genome sizes were similar among the four Syzygium species (343–372 Mb), a size range consistent with the small genome assembly sizes of S. aromaticum (370 Mb) and S. grande (405 Mb) (Low et al., 2022; Ouadi et al., 2022). The heterozygosity rate estimated by the GenomeScope 2.0 ranged from 2.3% for the diploid specimen S. malaccense to 4.3% for the autotetraploid specimen S. aqueum. These heterozygosity rates appeared to be higher than for S. aromaticum (0.18%) (Ouadi et al., 2022) and the average reported by Ellestad et al., who performed a literature review of the genome-wide heterozygosity values estimated using the software GenomeScope and GenomeScope 2.0 (Ellestad et al., 2022). They found that the average value inferred for all plant species assessed was 1.59% (1.10% for diploid plants only) noting that over half of the plant species considered were cultivated for human usage, which could affect the average value accuracy.
Table 1 Genome profiling summary.
3.2 Genome De Novo assembly
The four de novo chromosome-scale assemblies were constructed using long-reads from Oxford Nanopore Technologies (ONT), short paired-end reads from Illumina DNAseq libraries, and Hi-C libraries generated for each Syzygium species (Supplementary Tables 1, 2).
To prevent assembly artifacts possibly caused by heterozygosity and polyploidy of the Syzygium specimens, haplotigs were detected and removed from the polished contigs. The effect of the haplotig removal step was assessed using BUSCO (Benchmarking Universal Single-Copy Orthologs) in genome mode (Simão et al., 2015). After the haplotig removal step, the number of complete and duplicated BUSCOs genes was considerably reduced in the haplotig-purged contigs (3.3% to 6.1%) when compared to the polished contigs (93.6% to 97.1%) (Figure 1A). Hi-C data enabled the scaffolding of contigs into 11 chromosomes. On the Hi-C contact matrices, a strong intra-chromosomal signal indicates efficient scaffolding, with the 11 chromosomes of each Syzygium assembly supported by a high number of their respective Hi-C reads (Figure 1B).
Figure 1 Assessment of the efficiency of the haplotig removal step and Hi-C scaffolding. (A) BUSCO completeness score comparison of the polished contigs before and after the haplotig removal step for S. malaccense (Smal), S. aqueum (Saqu), S. jambos (Sjam), and S. syzygioides (Ssyz) (BUSCO version 5.4.4 - dataset: eudicots_odb10 (n = 2326)). (B) Hi-C contact maps showing the Hi-C interactions among the 11 assembled chromosomes and unplaced scaffolds (un) for each species.
The final chromosome-scale assemblies for S. malaccense (430 Mb), S. aqueum (392 Mb), S. jambos (426 Mb), and S. syzygioides (431 Mb) consisted of monoploid consensus (11 chromosomes and unplaced sequences) with comparable quality metrics. A high level of quality at the base-scale (quality value [QV] between 44.006 and 45.114), of contiguity (97.5% to 99.8% of the assemblies length anchored on 11 chromosomes) and completeness (BUSCO complete genes scores of 98%) was reached for the four new assembled Syzygium genomes (Table 2; Figure 2; Supplementary Tables 3, 4).
Figure 2 BUSCO completeness assessment. Assessment of the final genome assembly, transcript set, and protein set of S. aromaticum (Saro), S. malaccense (Smal), S. aqueum (Saqu), S. jambos (Sjam), and S. syzygioides (Ssyz) (BUSCO version 5.4.4 - dataset: eudicots_odb10 (n = 2326)).
Table 2 Assembly and annotation statistics.
Despite their high heterozygosity rate, the quality metrics for the genome assemblies of the diploid specimen S. malaccense and the autotetraploids S. aqueum, S. jambos, and S. syzygioides were comparable to those reported for S. aromaticum assembly (370 Mb) (Ouadi et al., 2022). Nevertheless, BUSCO scores revealed a higher percentage of complete and duplicated BUSCOs in the four new assemblies compared to S. aromaticum (2.2%), principally in the genome assembly of the three autotetraploid specimens (3.3% to 5.5%) (Figure 2).
3.3 Genome annotation
The average number of protein-coding genes predicted for the four newly assembled genomes is 31,119, representing 26.52% of the genome assemblies’ size (Table 2).
The annotation completeness was assessed using the BUSCO method in transcriptome and protein modes and by selecting the whole set of predicted transcripts and proteins for each gene as inputs, respectively (Figure 2; Supplementary Table 3). BUSCO results indicated that the annotation completeness is comparable among the four newly assembled Syzygium species, with complete BUSCO scores ranging from 91.9% in S. aqueum assembly to 93.5% in S. malaccense assembly in transcript mode and from 89.3% in S. aqueum assembly to 90.9% in S. malaccense assembly in protein mode. BUSCO scores obtained for S. aromaticum by using the same assessment methods (95% in transcriptome mode and 93.7% in protein mode) were slightly superior to those of newly assembled genomes but still comparable. The loss of complete BUSCOs between the genome and protein mode assessments ranged from 7.2% in S. malaccense assembly to 8.7% in S. aqueum assembly, indicating acceptable quality of the predicted gene models and protein sets.
The genome assembly of S. aromaticum comprised multiple copies of a gene encoding for putative eugenol synthase (EGS), the enzyme that catalyzes the synthesis of eugenol from coniferyl acetate. In total, 15 copies split into 2 loci were reported: a first locus on chromosome 10 comprising 14 copies and a second locus on chromosome 11 with 1 copy (Ouadi et al., 2022). The functional annotation of the four newly assembled Syzygium species genomes revealed fewer genes encoding for putative EGS. One gene encoding for putative EGS was identified in the genome assembly of S. malaccense, two in the genome assembly of S. aqueum, and three copies were found in the genome assemblies of S. jambos and S. syzygioides. All putative EGS genes were located on chromosome 10 except for one of the three copies of S. syzygioides located on chromosome 11 (Figure 3; Supplementary Table 5).
Figure 3 Illustration of the regions of chromosomes 10 and 11 of S. aromaticum (Saro), S. malaccense (Smal), S. aqueum (Saqu), S. jambos (Sjam), and S. syzygioides (Ssyz) where genes encoding for EGS were predicted. The position (bp) and orientation of the predicted genes on the chromosomes are indicated by arrows colored according to the functional annotation. EGS, accelerated cell death (ACD1), Protochlorophyllide-dependent translocon component Tic52 (PTC52), leucine-rich repeat receptor-like protein kinase (LRR-RLK), Pre-mRNA-processing protein 40C-like (PRP40C), TATA-binding protein-associated factor 7 (TBP-associated factor 7), LTR-RTs.
Effective lengths of repeat elements, which are different from their genomic length, were calculated by removing the length of the nested elements they contained. The proportions of genome assembly length occupied by predicted genes (25.97% to 27.37%) and repeat sequences (41.34% to 43.02%) appear to be conserved among the four newly sequenced Syzygium genomes (Table 2). Using the same method, repeat elements in Syzygium aromaticum genome assembly represents 39.98%. The most abundant repeat elements identified in the four newly sequenced Syzygium genomes were the LTR-RTs spanning 16.97% of the assembly length for S. syzygioides to 22.35% for S. malaccense. As reported for S. aromaticum and S. grande, LTR-RTs belonging to the Gypsy superfamily were more abundant than elements belonging to the Copia superfamily in the four newly sequenced genomes (Table 2; Supplementary Tables 6–9) (Low et al., 2022; Ouadi et al., 2022).
3.4 Synteny analyses
To identify evolutionary structural changes among the Syzygium species chromosomes, we performed a synteny analysis on the four newly assembled genomes, S. aromaticum and S. grande. The alignment of the 11 chromosomes’ DNA sequences of the 6 Syzygium species revealed a high conservation of the chromosomal organization (Figure 4A).
Figure 4 Identification of syntenic and rearranged regions between the 11 chromosomes of S. aromaticum (Saro), S. malaccense (Smal), S. aqueum (Saqu), S. jambos (Sjam), S. syzygioides (Ssyz), and S. grande (Sgra). (A) Representation of the alignment of the chromosomal DNA sequences showing syntenic regions, interchromosomal, and intrachromosomal rearrangements larger than 20 kb (inversions, translocations, and duplications). Grey arrows indicate regions where rearrangements were reported between chromosomes of E grandis and S. aromaticum. (B) Pairwise comparison of the percentage of chromosome length occupied by syntenic regions and rearrangements between the chromosome-scale assembly (Chr01-Chr11) and 11 chromosomes (Chr01 to Chr11) of S. aromaticum with those of S. malaccense, S. aqueum, S. jambos, S. syzygioides, and S. grande.
No large interchromosomal rearrangements were detected between the chromosomes of the six Syzygium species. A high percentage of the five species’ chromosome lengths were syntenic with S. aromaticum, ranging from 68.45% between S. aromaticum and S. jambos to 73.02% between S. aromaticum and S. aqueum. Intrachromosomal rearrangements such as inversions, translocations, and duplications between the chromosomes of S. aromaticum and those of the other five Syzygium species represented 5% of their 11 chromosomes length on average. In terms of number, the most frequent rearrangements observed between S. aromaticum and the five other species were duplications and translocations with average numbers of 1348 and 1325, respectively, spanning an average of 0.85% to 1.43% of the 11 chromosome lengths. Inversions were found less frequently for all species but occupied a larger fraction of the genome assemblies’ length than duplications and translocations except for S. syzygioides. The percentage of assembly lengths comprising inversions between S. aromaticum and the five other Syzygium species ranged from 0.68% between S. aromaticum and S. syzygioides to 4.83% between S. aromaticum and S. grande. Overall, the size of the inversions was relatively small. For instance, 11 inversions were detected, between chromosome 5 of S. aromaticum and S. grande, representing 17.32% of the chromosome length of S. grande (41,797,999 bp) and 1.87% of its 11 chromosomes (387,620,547 bp) (Figure 4B; Supplementary Table 4). In contrast, the synteny analysis performed between S. aromaticum and E. grandis revealed 10 intrachromosomal rearrangements on chromosomes 2, 4, 6, 8, 9, and 10 that included large terminal inversions representing up to 40% of the chromosome length of S. aromaticum. The other four chromosomes (1, 3, 5, and 7) of the two Myrtaceae species were highly syntenic (Ouadi et al., 2022). To further investigate the chromosomal architecture evolution of the Syzygium species and verify that these rearrangements were due to biological events rather than assembly artifacts, we also performed DNA alignment of the chromosome sequences of E. grandis with the those of S. malaccense, S. aqueum, S. jambos, and S. syzygioides. Chromosomes 1, 3, 5, and 7 of E. grandis and those of the four newly assembled species were also highly syntenic, and we observed the same 10 rearrangements on chromosomes 2, 4, 6, 8, 9, 10 and 11 (Figure 4A; Supplementary Figure 3).
3.5 Gene orthology
To investigate the phylogenetic relationships among gene sequences of S. aromaticum, S. malaccense, S. aqueum, S. jambos, and S. syzygioides, the sets of predicted protein sequences from the five species assemblies were analyzed using OrthoFinder (Emms and Kelly, 2019).
A total of 49,269 hierarchical orthogroups (HOGs) were identified, including 93.7 to 95.2% of each species gene set (Figure 5A). Of these, 18,963 (38.5%) HOGs contained genes from all five species, and 4,928 (10%) were species specific. In more detail, 789 were specific to S. aromaticum, 950 were specific to S. malaccense, 940 HOGs were specific to S. aqueum, 1009 HOGs were specific to S. jambos, and 1240 HOGs were specific to S. syzygioides. Pairwise, S. aromaticum and S. aqueum appear to share the lowest number of orthogroups (625). The highest number of shared HOGs inferred between each pair of studied species was found between S. aqueum and S. syzygioides (1218), followed by S. aqueum and S. malaccense (1152), and S. jambos and S. malaccense (1027). The species tree resulting from the analysis of the HOGs divided the Syzygium species studied into two groups based on closer relationships: the first group comprising S. aromaticum and S. aqueum and a second group comprising S. jambos, S. malaccense, and S. syzygioides (Figure 5B).
Figure 5 Hierarchical orthogroups (HOGs) inferred by OrthoFinder between S. aromaticum (Saro), S. malaccense (Smal), S. aqueum (Saqu), S. jambos (Sjam), and S. syzygioides (Ssyz). (A) Number of HOGs inferred by OrthoFinder using the set of predicted proteins for the five Syzygium species. (B) Rooted species tree inferred by OrthoFinder.
3.6 Annotation and comparison of LTR-RTs Gypsy and Copia repertoires
To clarify the dynamic activity of full-length LTR-RTs belonging to the superfamilies Gypsy and Copia within the Syzygium genus, we identified the lineages belonging to each superfamily located on the chromosomes of S. malaccense (429 Mbp), S. aqueum (387 Mbp), S. jambos (416 Mbp), and S. syzygioides (425 Mbp) and estimated their insertion time. Then, we compared the repertoires’ compositions and repeat element insertion times of the four species with those of S. aromaticum (368 Mbp).
We found that S. malaccense and S. aromaticum, the largest and smallest chromosome-scale assemblies of this study, contained the highest (8427) and lowest number (6167) of LTR-RTs in Gypsy and Copia, respectively (Figure 6A; Supplementary Tables 6–9). In the five Syzygium species’ chromosomes, we identified a higher number of LTR-RTs for Gypsy than Copia, with a ratio of Gypsy to Copia content ranging from 1.09 for S. syzygioides to 1.45 for S. malaccense. The Gypsy superfamily comprised a higher proportion of nested elements (17.37% to 24.47%) compared to the Copia superfamily (7.01% to 9.44%), suggesting distinct accumulation and mobile activity of both superfamilies in all five species. Our results revealed little variation in the number of Copia elements on the chromosomes of S. aqueum (2705 elements) and S. aromaticum (2809), the two smallest chromosome-scale assemblies, and on the chromosomes of S. syzygioides (3290), S. jambos (3324), and S. malaccense (3433). In contrast, we found a notably higher accumulation of Gypsy elements (4994) in the chromosomes of S. malaccense compared to the four other species. The ratio of Gypsy content varied from 1.35 when comparing S. malaccense with S. jambos to 1.49 when comparing S. malaccense with S. aromaticum. It represented a difference in Gypsy effective length of 19,402,234 bp to 21,766,176bp, respectively. In the five Syzygium chromosome-scale assemblies, the most abundant lineage belonged to the Gypsy superfamily, but it varied according to the species. The Gypsy lineage Tekay was the most represented for S. aromaticum (1534 elements), S. jambos (1674 elements), and S. syzygioides (2090 elements). At the same time, for S. malaccense and S. aqueum, we found a higher abundance of the gypsy lineage Ogre (2382 and 1799 elements, respectively). Among the Gypsy superfamily, the most abundant lineages, Tekay and Ogre, were those with the highest proportion of nested elements (19.10% to 28.55% and 16.69% to 27.92%, respectively) in all five species. For S. aromaticum, S. malaccense, and S. syzygioides, the proportion of nested elements belonging to the Athila lineage was also among the highest identified (Figure 6B; Supplementary Figures 4, 5).
Figure 6 Composition of the full-length LTR-RTs Gypsy and Copia repertoires. (A) Number of elements belonging to the Gypsy and Copia lineages identified on the 11 chromosomes of S. aromaticum (Saro), S. malaccense (Smal), S. aqueum (Saqu), S. jambos (Sjam), and S. syzygioides (Ssyz). (B) Proportion of nested and non-nested elements. Gypsy (others) group comprises the lineages non-chromo-outgroup, Reina, Retand, tatIII, and elements Gypsy to which no lineages were assigned. Copia (others) group comprises the lineages Alesia, Bianca, Gymco-I, Gymco-IV, Gymco-II, and Osser.
Regarding the Copia superfamily, the most represented lineages on the chromosomes of the five Syzygium species were Ale (608 to 873 elements), followed by the lineage Tork (456 to 762 elements) for S. malaccense, S. aqueum, S. jambos and S. Syzygoides, and the lineage SIRE (502 elements) for S. aromaticum.
The insertion times of 97.13% of the full-length Gypsy and Copia elements identified in the five Syzygium species were estimated (33,861elements). Nearly all elements (97.33%) were inserted in the last 5 million years (32,958 elements) (Figure 7). During this time period, distinct insertion activities of the two superfamilies occurred in the five Syzygium species.
Figure 7 Distribution of insertion times of full-length LTR-RTs of S. aromaticum (Saro), S. malaccense (Smal), S. aqueum (Saqu), S. jambos (Sjam), and S. syzygioides (Ssyz). (A) LTR-RTs Gypsy. (B) LTR-RTs Copia.
Compared to the other four Syzygium species, the chromosomes of S. aromaticum underwent a more ancient wave of Gypsy insertions (peak at ~2.5 million years ago [Mya]), principally attributed to the Tekay elements, the most abundant lineage in this species (Figure 7A). We also found that a few recent insertions (18.02% of insertions) occurred in S. aromaticum chromosomes within the last one million years. In contrast, a recent burst of Gypsy insertions (~0–1 Mya) occurred in four other species chromosomes: most insertions of Gypsy in S. malaccense (44.53%), S. aqueum (44.43%), S. jambos (52.55%), and S. syzygioides (36.45%) were less than one million years old. We inferred that the high number of Gypsy LTR-RTs found in S. malaccense may be attributable to two successive waves of insertions: a peak of Tekay insertions at ~2 Mya and a more recent peak of Ogre at ~1 Mya.
Similar to what we observed for the Gypsy superfamily, the insertion of Copia elements occurred earlier in S. aromaticum compared to the four other species, with fewer recent insertions (Figure 7B). Compared to the Gypsy elements, a smaller proportion of recent Copia insertions (less than one million years old) were detected in S. aromaticum (10.66%), S. malaccense (24.16%), S. aqueum (26.49%), and S. jambos (36.90%) suggesting a distinct recent insertion pattern of the two superfamilies in the four species. However, we found a comparable proportion of Gypsy (36.45%) and Copia (32.22%) elements that were less than one million years old in S. syzygioides, the species for which we found the lowest ratio of Gypsy to Copia content (1.09).
4 Discussion
Plant genome size, ploidy level, and heterozygosity rates are challenges for genome assembly and annotation. However, lower sequencing costs and recent advances in long-read sequencing technologies, Hi-C technologies, and bioinformatics tools have facilitated the generation of assemblies with high contiguity up to the chromosome-scale also for non-model plants or non-major plant crops (Kyriakidou et al., 2018; Pucker et al., 2022). Newly assembled and annotated genomes from related species can then be used to perform comparative genomics analyses to investigate plant genome evolution and function. Third-generation long-reads from Oxford Nanopore Technologies and Illumina short-reads combined with the Hi-C technology enabled the de novo assembly of the chromosome-scale genome for S. malaccense, S. aqueum, S. jambos, and S. syzygioides. A high level of quality at the base level, contiguity, and completeness was reached for the four newly sequenced genomes. The quality of the newly assembled Syzygium species genomes were comparable to that of the S. aromaticum genome. The slight differences found between the species assemblies’ quality metrics may be linked to the combined impact of the ploidy level and high heterozygosity rates of the four newly sequenced species on the assembly process.
Previous infrageneric comparative genetic mapping analyses revealed high levels of synteny and collinearity among the Eucalyptus genus (Hudson et al., 2012; Li et al., 2015). In addition, genomic synteny analyses conducted between the de novo assembly of E. urophylla × E. grandis (EUC) and 30 Eucalyptus species revealed that the genome structure of EUC, E. grandis, and E. globulus showed the higher collinearity, and the absence of large-scale structural variation. Nevertheless, large structural variations among the different chromosomes of the EUC and other Eucalyptus species were also detected (Shen et al., 2023). We found that the six Syzygium genomes studied were highly syntenic. The intrachromosomal rearrangements (duplications, translocations, and inversions) observed between S. aromaticum and the five other Syzygium species represent a small percentage (~5% on average) of the 11 chromosomes’ length. These intrachromosomal rearrangements could result from contigs that were not well placed because of Hi-C signals that were not strong enough to correctly determine their position and orientation; however, they may also result from the six species’ distinct genome evolutions.
Organizational conservation of chromosomes 2, 4, 6, 8, 9, 10, and 11 among the six Syzygium species studied constitutes new evidence supporting the 10 intrachromosomal rearrangements previously reported on these chromosomes between S. aromaticum and E. grandis genomes (Ouadi et al., 2022). These 10 rearrangements were also observed when aligning the DNA sequences of the chromosomes of E. grandis with those of S. malaccense, S. aqueum, S. jambos, and S. syzygioides. Among the rearrangements reported between the chromosomes of S. aromaticum and E. grandis, similar large terminal inversions on chromosomes 4, 9, 10, and 11 were also reported in the two eucalypts E. grandis and C. citriodora suggesting that these terminal inversions occurred on E. grandis chromosomes (Butler et al., 2017). Two other large terminal inversions were detected between chromosomes 4 and 9 of S. aromaticum and E. grandis but not between C. citriodora and E. grandis. These inversions were also observed when comparing the chromosome sequences of E. grandis with those of the four newly assembled genomes, suggesting that these inversions resulted from an evolution of the chromosome organization rather than from sequencing and assembly artifacts. Further comparative genomics analyses will be needed with additional Syzygium and Myrtaceous species to determine if these inversions are specific to the Syzygium genus or subgenus, for which the crown ages were estimated at 51.2 Mya and 9.4 Mya, respectively, (Low et al., 2022).
The analyses of the phylogenetic relationships between gene sequences of S. aromaticum, S. malaccense, S. aqueum, S. jambos, and S. syzygioides and comparisons of their full-length LTR-RTs repertoires provided insights into the distinct genome evolution of each species following the divergence of the Syzygium subg. Syzygium species 9.4 Mya (Low et al., 2022). The species tree inferred by OrthoFinder indicated that pairwise S. aromaticum and S. aqueum and S. malaccense and S. jambos were closely related, which is consistent with the genome-level phylogenetic trees generated by Low et al. (Low et al., 2022). We observed older waves of LTR-RTs Gypsy and Copia insertions in S. aromaticum and fewer insertions less than 1 million years old in the S. aromaticum chromosomes compared to those of the four other species studied. In plants, the RNA Directed DNA Methylation (RdDM) pathway, a de novo DNA methylation mechanism involving small interfering RNA, plays an important role in TE repression (Wambui Mbichi et al., 2020). Further detailed analysis such as DNA methylation studies will be valuable to clarify the molecular causes of the recent low insertion number of LTR-RTs elements observed in S. aromaticum.
S. aromaticum is cultivated to produce clove bud (the dried, unopened flower bud), essential oil (EO), and oleoresins rich in eugenol (Nurdjannah and Bermawie, 2012). The EO of S. aromaticum contains ~72 to 96.6% of eugenol, while the EO of S. aqueum has 0.19% eugenol (Razafimamonjison et al., 2014; Sobeh et al., 2016). Eugenol is a phenylpropane with multiple pharmaceutical activities and is considered a promising alternative drug for human health (e.g., cancer and pathogenic microorganism resistance, diabetes, obesity, and autoimmune diseases) (Kamatou et al., 2012; Batiha et al., 2020; Otunola, 2022). The genome assembly of S. aromaticum was exploited to investigate the genetic basis of this important characteristic. The identification of gene families involved in eugenol biosynthesis revealed the presence of multiple copies of genes encoding EGS, which catalyzes the synthesis of eugenol from coniferyl acetate. A cluster of 14 copies was reported on chromosome 10, and additional copies were located on chromosome 11 of S. aromaticum. In the genome assembly of the four newly sequenced species, we found fewer gene copies on chromosome 10 (1 to 3 copies) and no copies on chromosome 11 of S. malaccense, S. aqueum, and S. jambos. The presence of this structural variation suggested that a gene-dosage effect may be associated with the high amount of eugenol. Further studies are needed to elucidate the biological functions of the EGS gene copies in S. aromaticum and the four other species (e.g., in vitro characterization).
S. malaccense, S. aqueum, and S. jambos are grown for their edible fruit. Like S. aromaticum and other Syzygium species, they are also used in traditional medicine. Research on their numerous pharmaceutical properties has been undertaken (e.g., analgesic, anti-inflammatory, antioxidant, hepatoprotective, antidiabetic, antifungal, antibacterial, antiviral, and anticancer activities) (Nair, 2017; Cock and Cheesman, 2018). For instance, S. jambos is traditionally used to treat hemorrhages, wounds, and ulcers; S. malaccense is used to treat mouth ulcers and diabetes; and S. aqueum to treat diabetes and childbirth pain (Uddin et al., 2022). The chromosome-scale assemblies for these species are new valuable resources for the Myrtaceae family. Combined with other comparative genomics and multi-omics studies, they can be used to further investigate the genomic evolution of the Myrtaceous species and to study the genetic basis of important agronomical traits and biosynthesis of secondary metabolites.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/, PRJNA962868 https://www.ncbi.nlm.nih.gov/, PRJNA962711 https://www.ncbi.nlm.nih.gov/, PRJNA962713 https://www.ncbi.nlm.nih.gov/, PRJNA962712 https://www.ncbi.nlm.nih.gov/genbank/, JASUUE000000000 https://www.ncbi.nlm.nih.gov/genbank/, JASUUB000000000 https://www.ncbi.nlm.nih.gov/genbank/, JASUUC000000000 https://www.ncbi.nlm.nih.gov/genbank/, JASUUD000000000 https://zenodo.org/, 7870328 https://zenodo.org/, 7870326 https://zenodo.org/, 7870330 https://zenodo.org/, 7870334.
Author contributions
SO performed the laboratory work, analyzed data, and wrote the manuscript. NS performed computational analysis of sequencing data, conceived, and supervised the study, and contributed to manuscript writing. FK, and NI conceived and supervised the study and contributed to manuscript writing. All authors contributed to the article and approved the submitted version.
Funding
The authors declare that this study received funding from the company Philip Morris International. The funder had the following involvement in the study: the study design, collection, analysis, interpretation of data, the writing of this article and the decision to submit it for publication.
Acknowledgments
We would like to thank Dr. Leyla Davis, curator of the Masoala Hall in the Zürich Zoo (Switzerland), for authorizing the sampling of the Syzygium trees for this project. We would like to also thank Remi Dulize for his technical contributions, and Lindsay Reese and Rebecca Higgins for manuscript revision.
Conflict of interest
Authors SO, NS, and NI were employed by the company Philip Morris International.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2023.1248780/full#supplementary-material
References
Batiha, G. E.-S., Alkazmi, L. M., Wasef, L. G., Beshbishy, A. M., Nadwa, E. H., Rashwan, E. K. (2020). Syzygium aromaticum L.(Myrtaceae): Traditional uses, bioactive chemical constituents, pharmacological and toxicological activities. Biomolecules 10 (2), 202. doi: 10.3390/biom10020202
Beech, E., Rivers, M., Oldfield, S., Smith, P. (2017). GlobalTreeSearch: The first complete global database of tree species and country distributions. J. Sustain. For. 36 (5), 454–489. doi: 10.1080/10549811.2017.1310049
Buchfink, B., Xie, C., Huson, D. H. (2015). Fast and sensitive protein alignment using DIAMOND. Nat. Methods 12 (1), 59–60. doi: 10.1038/nmeth.3176
Butler, J., Vaillancourt, R., Potts, B., Lee, D., King, G. J., Baten, A., et al. (2017). Comparative genomics of Eucalyptus and Corymbia reveals low rates of genome structural rearrangement. BMC Genom. 18 (1), 397. doi: 10.1186/s12864-017-3782-7
Chen, S., Zhou, Y., Chen, Y., Gu, J. (2018). fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34 (17), i884–i890. doi: 10.1093/bioinformatics/bty560
Cheng, H., Concepcion, G. T., Feng, X., Zhang, H., Li, H. (2021). Haplotype-resolved de novo assembly using phased assembly graphs with hifiasm. Nat. Methods 18 (2), 170–175. doi: 10.1038/s41592-020-01056-5
Christenhusz, M. J., Byng, J. W. (2016). The number of known plants species in the world and its annual increase. Phytotaxa 261 (3), 201–217-201–217. doi: 10.11646/phytotaxa.261.3.1
Cock, I. E., Cheesman, M. (2018). Plants of the genus Syzygium (Myrtaceae): A review on ethnobotany, medicinal properties and phytochemistry. Bioactive Compounds Medicinal Plants: Properties Potential Hum. Health 35–84. doi: 10.1201/b22426
Craven, L. A., Biffin, E. (2010). An infrageneric classification of Syzygium (Myrtaceae). Blumea-Biodiver. Evol. Biogeogr. Plants 55 (1), 94–99. doi: 10.3767/000651910X499303
Ellestad, P., Pérez-Farrera, M. A., Buerki, S. (2022). Genomic Insights into Cultivated Mexican Vanilla planifolia Reveal High Levels of Heterozygosity Stemming from Hybridization. Plants 11 (16), 2090. doi: 10.3390/plants11162090
Ellinghaus, D., Kurtz, S., Willhoeft, U. (2008). LTRharvest, an efficient and flexible software for de novo detection of LTR retrotransposons. BMC Bioinf. 9, 1–14. doi: 10.1186/1471-2105-9-18
Emms, D. M., Kelly, S. (2019). OrthoFinder: phylogenetic orthology inference for comparative genomics. Genome Biol. 20 (1), 1–14. doi: 10.1186/s13059-019-1832-y
Feng, C., Feng, C., Lin, X., Liu, S., Li, Y., Kang, M. (2021). A chromosome-level genome assembly provides insights into ascorbic acid accumulation and fruit softening in guava (Psidium guajava). Plant Biotechnol. J. 19 (4), 717–730. doi: 10.1111/pbi.13498
Frith, M. C. (2011). A new repeat-masking method enables specific detection of homologous sequences. Nucleic Acids Res. 39 (4), e23–e23. doi: 10.1093/nar/gkq1212
Fu, L., Niu, B., Zhu, Z., Wu, S., Li, W. (2012). CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics 28 (23), 3150–3152. doi: 10.1093/bioinformatics/bts565
Girgis, H. Z. (2015). Red: an intelligent, rapid, accurate tool for detecting repeats de-novo on the genomic scale. BMC Bioinf. 16 (1), 1–19. doi: 10.1186/s12859-015-0654-5
Goel, M., Schneeberger, K. (2022). plotsr: visualizing structural similarities and rearrangements between multiple genomes. Bioinformatics 38 (10), 2922–2926. doi: 10.1093/bioinformatics/btac196
Goel, M., Sun, H., Jiao, W.-B., Schneeberger, K. (2019). SyRI: finding genomic rearrangements and local sequence differences from whole-genome assemblies. Genome Biol. 20 (1), 1–13. doi: 10.1186/s13059-019-1911-0
Grattapaglia, D., Vaillancourt, R. E., Shepherd, M., Thumma, B. R., Foley, W., Külheim, C., et al. (2012). Progress in Myrtaceae genetics and genomics: Eucalyptus as the pivotal genus. Tree Genet. Genomes 8 (3), 463–508. doi: 10.1007/s11295-012-0491-x
Guan, D., McCarthy, S. A., Wood, J., Howe, K., Wang, Y., Durbin, R. (2020). Identifying and removing haplotypic duplication in primary genome assemblies. Bioinformatics 36 (9), 2896–2898. doi: 10.1093/bioinformatics/btaa025
Healey, A. L., Shepherd, M., King, G. J., Butler, J. B., Freeman, J. S., Lee, D. J., et al. (2021). Pests, diseases, and aridity have shaped the genome of Corymbia citriodora. Commun. Biol. 4 (1), 1–13. doi: 10.1038/s42003-021-02009-0
Hu, K., Xu, K., Wen, J., Yi, B., Shen, J., Ma, C., et al. (2019). Helitron distribution in Brassicaceae and whole Genome Helitron density as a character for distinguishing plant species. BMC Bioinf. 20 (1), 1–20. doi: 10.1186/s12859-019-2945-8
Hudson, C. J., Kullan, A. R., Freeman, J. S., Faria, D. A., Grattapaglia, D., Kilian, A., et al. (2012). High synteny and colinearity among Eucalyptus genomes revealed by high-density comparative genetic mapping. Tree Genet. Genomes 8 (2), 339–352. doi: 10.1007/s11295-011-0444-9
Izuno, A., Wicker, T., Hatakeyama, M., Copetti, D., Shimizu, K. K. (2019). Updated genome assembly and annotation for metrosideros polymorpha, an emerging model tree species of ecological divergence. G3-Genes Genom. Genet. 9 (11), 3513–3520. doi: 10.1534/g3.119.400643
Kamatou, G. P., Vermaak, I., Viljoen, A. M. (2012). Eugenol—from the remote Maluku Islands to the international market place: a review of a remarkable and versatile molecule. Molecules 17 (6), 6953–6981. doi: 10.3390/molecules17066953
Kolmogorov, M., Yuan, J., Lin, Y., Pevzner, P. A. (2019). Assembly of long, error-prone reads using repeat graphs. Nat. Biotechnol. 37 (5), 540–546. doi: 10.1038/s41587-019-0072-8
Kyriakidou, M., Tai, H. H., Anglin, N. L., Ellis, D., Strömvik, M. V. (2018). Current strategies of polyploid plant genome sequence assembly. Front. Plant Sci. 9, 1660. doi: 10.3389/fpls.2018.01660
Lexa, M., Jedlicka, P., Vanat, I., Cervenansky, M., Kejnovsky, E. (2020). TE-greedy-nester: structure-based detection of LTR retrotransposons and their nesting. Bioinformatics 36 (20), 4991–4999. doi: 10.1093/bioinformatics/btaa632
Li, H. (2013). Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv:1303.3997v2. doi: 10.48550/arXiv.1303.3997
Li, H. (2018). Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics 34 (18), 3094–3100. doi: 10.1093/bioinformatics/bty191
Li, F., Xu, S., Xiao, Z., Wang, J., Mei, Y., Hu, H., et al. (2023). Gap-free genome assembly and comparative analysis reveal the evolution and anthocyanin accumulation mechanism of Rhodomyrtus tomentosa. Hortic. Res 10 (3). doi: 10.1093/hr/uhad005
Li, F., Zhou, C., Weng, Q., Li, M., Yu, X., Guo, Y., et al. (2015). Comparative genomics analyses reveal extensive chromosome colinearity and novel quantitative trait loci in Eucalyptus. PloS One 10 (12), e0145144. doi: 10.1371/journal.pone.0145144
Low, Y. W., Rajaraman, S., Tomlin, C. M., Ahmad, J. A., Ardi, W. H., Armstrong, K., et al. (2022). Genomic insights into rapid speciation within the world’s largest tree genus Syzygium. Nat. Commun. 13 (1), 1–15. doi: 10.1038/s41467-022-32637-x
Machado, R. M., Forni-Martins, E. R. (2022). Psidium cattleyanum Sabine (Myrtaceae), a neotropical polyploid complex with wide geographic distribution: insights from cytogenetic and DNA content analysis. Braz. J. Bot. 45 (3), 943–955. doi: 10.1007/s40415-022-00829-w
Mak, Q. C., Wick, R. R., Holt, J. M., Wang, J. R. (2023). Polishing de novo nanopore assemblies of bacteria and eukaryotes with FMLRC2. Mol. Biol. Evol. 40 (3), msad048. doi: 10.1093/molbev/msad048
Marcon, H. S., Domingues, D. S., Silva, J. C., Borges, R. J., Matioli, F. F., de Mattos Fontes, M. R., et al. (2015). Transcriptionally active LTR retrotransposons in Eucalyptus genus are differentially expressed and insertionally polymorphic. BMC Plant Biol. 15 (1), 1–16. doi: 10.1186/s12870-015-0550-1
Myburg, A. A., Grattapaglia, D., Tuskan, G. A., Hellsten, U., Hayes, R. D., Grimwood, J., et al. (2014). The genome of Eucalyptus grandis. Nature 510 (7505), 356–362. doi: 10.1038/nature13308
Nair, K. N. (2017). The genus Syzygium: Syzygium Cumini and Other Underutilized Species (United States: CRC Press).
Neumann, P., Novák, P., Hoštáková, N., Macas, J. (2019). Systematic survey of plant LTR-retrotransposons elucidates phylogenetic relationships of their polyprotein domains and provides a reference for element classification. Mobile DNA 10, 1–17. doi: 10.1186/s13100-018-0144-1
Niknafs, Y. S., Pandian, B., Iyer, H. K., Chinnaiyan, A. M., Iyer, M. K. (2017). TACO produces robust multisample transcriptome assemblies from RNA-seq. Nat. Methods 14 (1), 68–70. doi: 10.1038/nmeth.4078
Nurdjannah, N., Bermawie, N. (2012). “Cloves,” in Handbook of herbs and spices (Amsterdam, Neatherlands: Elsevier), 197–215.
Oginuma, K., Kato, A., Tobe, H., Mathenge, S., Juma, F. (1993). Chromosomes of some woody plants in Kenya. Acta Phytotax. Geobot. 44 (1), 53–58.
Otunola, G. A. (2022). Culinary spices in food and medicine: an overview of Syzygium aromaticum (L.) Merr. and LM Perry [Myrtaceae]. Front. Pharmacol. 12, 3817. doi: 10.3389/fphar.2021.793200
Ouadi, S., Sierro, N., Goepfert, S., Bovet, L., Glauser, G., Vallat, A., et al. (2022). The clove (Syzygium aromaticum) genome provides insights into the eugenol biosynthesis pathway. Commun. Biol. 5 (1), 1–13. doi: 10.1038/s42003-022-03618-z
Panggabean, G. (1991). “Syzygium aqueum (Burm. f.) Alst., Syzygium malaccense (L.) M. & P, and Syzygium samarangense (Blume) M. & P. Plant Resources of South-East Asia 2,” in Edible fruits and nuts (Pudoc, Wageningen: Pudoc Scientific Publishers), 292–294.
Parnell, J. A., Craven, L. A., Biffin, E. (2007). “Matters of scale: dealing with one of the largest genera of angiosperms,” in Reconstructing the tree of life: taxonomy and systematics of species rich taxa (Boca Raton, FL: CRC Press LLC), 253–270.
Pedrosa, A., Gitaí, J., Silva, A. E. B., Felix, L. P., Guerra, M. (1999). Cytogenetics of angiosperms collected in the state of Pernambuco: V. Acta Bot. Bras. 13 (1), 49–60. doi: 10.1590/S0102-33061999000100006
Pellicer, J., Leitch, I. J. (2020). The Plant DNA C-values database (release 7.1): an updated online repository of plant genome size data for comparative studies. New Phytol. 226 (2), 301–305. doi: 10.1111/nph.16261
Pertea, G., Pertea, M. (2020). GFF Utilities: GffRead and GffCompare [version 2; peer review: 3 approved]. F1000Research 9 (304). doi: 10.12688/f1000research.23297.2
POWO (2023) Plants of the World Online. Facilitated by the Royal Botanic Gardens, Kew. Available at: http://www.plantsoftheworldonline.org/. Retrieved 11 April 2023.
Pucker, B., Irisarri, I., de Vries, J., Xu, B. (2022). Plant genome sequence assembly in the era of long reads: Progress, challenges and future directions. Quant. Plant Biol. 3 (5), e5. doi: 10.1017/qpb.2021.18
Quinlan, A. R., Hall, I. M. (2010). BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics 26 (6), 841–842. doi: 10.1093/bioinformatics/btq033
Ranallo-Benavidez, T. R., Jaron, K. S., Schatz, M. C. (2020). GenomeScope 2.0 and Smudgeplot for reference-free profiling of polyploid genomes. Nat. Commun. 11 (1), 1–10. doi: 10.1038/s41467-020-14998-3
Razafimamonjison, G., Jahiel, M., Duclos, T., Ramanoelina, P., Fawbush, F., Danthu, P. (2014). Bud, leaf and stem essential oil composition of Syzygium aromaticum from Madagascar, Indonesia and Zanzibar. Int. J. Basic Appl. Sci. 3 (3), 224. doi: 10.14419/ijbas.v3i3.2473
Saber, F. R., Munekata, P. E., Rizwan, K., El-Nashar, H. A., Fahmy, N. M., Aly, S. H., et al. (2023). Family Myrtaceae: The treasure hidden in the complex/diverse composition. Crit. Rev. Food Sci. Nutr., 1–19. doi: 10.1080/10408398.2023.2173720
Shao, M., Kingsford, C. (2017). Accurate assembly of transcripts through phase-preserving graph decomposition. Nat. Biotechnol. 35 (12), 1167–1169. doi: 10.1038/nbt.4020
Shen, W., Le, S., Li, Y., Hu, F. (2016). SeqKit: a cross-platform and ultrafast toolkit for FASTA/Q file manipulation. PloS One 11 (10), e0163962. doi: 10.1371/journal.pone.0163962
Shen, C., Li, L., Ouyang, L., Su, M., Guo, K. (2023). E. urophylla× E. grandis high-quality genome and comparative genomics provide insights on evolution and diversification of eucalyptus. BMC Genom. 24 (1), 1–10. doi: 10.1186/s12864-023-09318-0
Shi, J., Liang, C. (2019). Generic repeat finder: a high-sensitivity tool for genome-wide de novo repeat detection. Plant Physiol. 180 (4), 1803–1815. doi: 10.1104/pp.19.00386
Simão, F. A., Waterhouse, R. M., Ioannidis, P., Kriventseva, E. V., Zdobnov, E. M. (2015). BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31 (19), 3210–3212. doi: 10.1093/bioinformatics/btv351
Sobeh, M., Braun, M. S., Krstin, S., Youssef, F. S., Ashour, M. L., Wink, M. (2016). Chemical profiling of the essential oils of Syzygium aqueum, Syzygium samarangense and Eugenia uniflora and their discrimination using chemometric analysis. Chem. Biodivers. 13 (11), 1537–1550. doi: 10.1002/cbdv.201600089
Thrimawithana, A. H., Jones, D., Hilario, E., Grierson, E., Ngo, H. M., Liachko, I., et al. (2019). A whole genome assembly of Leptospermum scoparium (Myrtaceae) for mānuka research. N. Z. J. Crop Hortic. Sci. 47 (4), 233–260. doi: 10.1080/01140671.2019.1657911
Tuler, A. C., Carrijo, T. T., Peixoto, A. L., Garbin, M. L., da Silva Ferreira, M. F., Carvalho, C. R., et al. (2019). Diversification and geographical distribution of Psidium (Myrtaceae) species with distinct ploidy levels. Trees 33 (4), 1101–1110. doi: 10.1007/s00468-019-01845-2
Uddin, A. N., Hossain, F., Reza, A. A., Nasrin, M. S., Alam, A. K. (2022). Traditional uses, pharmacological activities, and phytochemical constituents of the genus Syzygium: A review. Food Sci. Nutr. 10 (6), 1789–1819. doi: 10.1002/fsn3.2797
Van Lingen, T. (1991). “Syzygium jambos (L.) Alston. Plant Resources of South-East Asia 2,” in Edible fruits and nuts (Pudoc, Wageningen: Pudoc Scientific Publishers), 296–298.
Wambui Mbichi, R., Wang, Q.-F., Wan, T. (2020). RNA directed DNA methylation and seed plant genome evolution. Plant Cell Rep. 39, 983–996. doi: 10.1007/s00299-020-02558-4
Warren, R. L., Coombe, L., Mohamadi, H., Zhang, J., Jaquish, B., Isabel, N., et al. (2019). ntEdit: scalable genome sequence polishing. Bioinformatics 35 (21), 4430–4432. doi: 10.1093/bioinformatics/btz400
Wicker, T., Sabot, F., Hua-Van, A., Bennetzen, J. L., Capy, P., Chalhoub, B., et al. (2007). A unified classification system for eukaryotic transposable elements. Nat. Rev. Genet. 8 (12), 973–982. doi: 10.1038/nrg2165
Wilson, P. G. (2010). “Myrtaceae,” in Flowering Plants. Eudicots (Berlin, Heidelberg: Springer), 212–271.
Zhang, R.-G., Li, G.-Y., Wang, X.-L., Dainat, J., Wang, Z.-X., Ou, S., et al. (2022). TEsorter: an accurate and fast method to classify LTR-retrotransposons in plant genomes. Hortic. Res. 9. doi: 10.1093/hr/uhac017
Zheng, X., Chen, X., Lin, G., Chen, J., Li, H., Xiao, Y., et al. (2022). The chromosome-level Melaleuca alternifolia genome provides insights into the molecular mechanisms underlying terpenoids biosynthesis. Ind. Crops Prod. 189, 115819. doi: 10.1016/j.indcrop.2022.115819
Zhou, C., McCarthy, S. A., Durbin, R. (2022). YaHS: yet another Hi-C scaffolding tool. Bioinformatics 39 (1). doi: 10.1093/bioinformatics/btac808
Keywords: Syzygium, Myrtaceae, de novo assembly, comparative genomics, synteny, long terminal repeat retrotransposons
Citation: Ouadi S, Sierro N, Kessler F and Ivanov NV (2023) Chromosome-scale assemblies of S. malaccense, S. aqueum, S. jambos, and S. syzygioides provide insights into the evolution of Syzygium genomes. Front. Plant Sci. 14:1248780. doi: 10.3389/fpls.2023.1248780
Received: 27 June 2023; Accepted: 28 August 2023;
Published: 06 October 2023.
Edited by:
Agnieszka Zmienko, Polish Academy of Sciences, PolandReviewed by:
Pawel Wojciechowski, Poznań University of Technology, PolandXiaojun Nie, Northwest A&F University, China
Jian-Feng Mao, Beijing Forestry University, China
Copyright © 2023 Ouadi, Sierro, Kessler and Ivanov. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Nikolai V. Ivanov, nikolai.ivanov@unine.ch