Complete Chloroplast Genomes of 14 Mangroves: Phylogenetic and Comparative Genomic Analyses

Genome assembly of the milky mangrove Excoecaria agallocha

The milky mangrove Excoecaria agallocha is a latex-secreting mangrove that are distributed in tropical and subtropical regions. While its poisonous latex is regarded as a potential source of phytochemicals for biomedical applications, the genomic resources of E. agallocha remains limited. Here, we present a chromosomal level genome of E. agallocha, assembled from the combination of PacBio long-read sequencing and Omni-C data. The resulting assembly size is 1,332.45 Mb and has high contiguity and completeness with a scaffold N50 of 58.9 Mb and a BUSCO score of 98.4%, with 86.08% of sequences anchored to 18 pseudomolecules. 73,740 protein-coding genes were also predicted. The milky mangrove genome provides a useful resource for further understanding the biosynthesis of phytochemical compounds in E. agallocha.

Phylogenomic and syntenic data demonstrate complex evolutionary processes in early radiation of the rosids

Some of the most vexing problems of deep level relationship that remain in angiosperms involve the superrosids. The superrosid clade contains a quarter of all angiosperm species, with 18 orders in three subclades (Vitales, Saxifragales and core rosids) exhibiting remarkable morphological and ecological diversity. To help resolve deep-level relationships, we constructed a high-quality chromosome-level genome assembly for Tiarella polyphylla (Saxifragaceae) thus providing broader genomic representation of Saxifragales. Whole genome microsynteny analysis of superrosids showed that Saxifragales shared more synteny clusters with core rosids than Vitales, further supporting Saxifragales as more closely related with core rosids. To resolve the ordinal phylogeny of superrosids, we screened 122 single copy nuclear genes from genomes of 36 species, representing all 18 superrosid orders. Vitales were recovered as sister to all other superrosids (Saxifragales + core rosids). Our data suggest dramatic differences in relationships compared to earlier studies within core rosids. Fabids should be restricted to the nitrogen-fixing clade, while Picramniales, the Celastrales-Malpighiales (CM) clade, Huerteales, Oxalidales, Sapindales, Malvales and Brassicales formed an "expanded" malvid clade. The Celastrales-Oxalidales-Malpighiales (COM) clade (sensu APG IV) was not monophyletic. Crossosomatales, Geraniales, Myrtales and Zygophyllales did not belong to either of our well-supported malvids or fabids. There is strong discordance between nuclear and plastid phylogenetic hypotheses for superrosid relationships; we show that this is best explained by a combination of incomplete lineage sorting and ancient reticulation.

Interspecific variation and phylogenetic relationship between mangrove and non-mangrove species of a same family (Meliaceae)-insights from comparative analysis of complete chloroplast genome

The mahogany family, Meliaceae, contains 58 genera with only one mangrove genus: Xylocarpus. Two of the three species of the genus Xylocarpus are true mangroves (X. granatum and X. moluccensis), and one is a non-mangrove (X. rumphii). In order to resolve the phylogenetic relationship between the mangrove and non-mangrove species, we sequenced chloroplast genomes of these Xylocarpus species along with two non-mangrove species of the Meliaceae family (Carapa guianensis and Swietenia macrophylla) and compared the genome features and variations across the five species. The five Meliaceae species shared 130 genes (85 protein-coding genes, 37 tRNA, and eight rRNA) with identical direction and order, with a few variations in genes and intergenic spacers. The repetitive sequences identified in the rpl22 gene region only occurred in Xylocarpus, while the repetitive sequences in accD were found in X. moluccensis and X. rumphii. The TrnH-GUG and rpl32 gene regions and four non-coding gene regions showed high variabilities between X. granatum and the two non-mangrove species (S. macrophylla and C. guianensis). In addition, among the Xylocarpus species, only two genes (accD and clpP) showed positive selection. Carapa guianensis and S. macrophylla owned unique RNA editing sites. The above genes played an important role in acclimation to different stress factors like heat, low temperature, high UV light, and high salinity. Phylogenetic analysis with 22 species in the order Sapindales supported previous studies, which revealed that the non-mangrove species X. rumphii is closer to X. moluccensis than X. granatum. Overall, our results provided important insights into the variation of genetic structure and adaptation mechanism at interspecific (three Xylocarpus species) and intergeneric (mangrove and non-mangrove genera) levels.

Comparative chloroplast genome analyses of Amomum: insights into evolutionary history and species identification

BACKGROUND

Species in genus Amomum always have important medicinal and economic values. Classification of Amomum using morphological characters has long been a challenge because they exhibit high similarity. The main goals of this study were to mine genetic markers from cp genomes for Amomum species identification and discover their evolutionary history through comparative analysis.

RESULTS

Three species Amomum villosum, Amomum maximum and Amomum longipetiolatum were sequenced and annotated for the complete chloroplast (cp) genomes, and the cp genomes of A. longipetiolatum and A. maximum were the first reported. Three cp genomes exhibited typical quadripartite structures with 163,269-163,591 bp in length. Each genome encodes 130 functional genes including 79 protein-coding, 26 tRNAs and 3 rRNAs genes. 113-152 SSRs and 99 long repeats were identified in the three cp genomes. By designing specific primers, we amplified the highly variable loci and the mined genetic marker ccsA exhibited a relatively high species identification resolution in Amomum. The nonsynonymous and synonymous substitution ratios (Ka/Ks) in Amomum and Alpinia showed that most genes were subjected to a purifying selection. Phylogenetic analysis revealed the evolutionary relationships of Amomum and Alpinia species and proved that Amomum is paraphyletic. In addition, the sequenced sample of A. villosum was found to be a hybrid, becoming the first report of natural hybridization of this genus. Meanwhile, the high-throughput sequencing-based ITS2 analysis was proved to be an efficient tool for interspecific hybrid identification and with the help of the chloroplast genome, the hybrid parents can be also be determined.

CONCLUSION

The comparative analysis and mined genetic markers of cp genomes were conducive to species identification and evolutionary relationships of Amomum.

Comparative chloroplast genome and phylogenetic analyses of Chinese Polyspora

Polyspora Sweet (Theaceae) are winter ornamental landscape plants native to southern and southeastern Asia, some of which have medicinal value. The chloroplast (cp) genome data of Polyspora are scarce, and the gene evolution and interspecific relationship are still unclear. In this study, we sequenced and annotated Polyspora chrysandra cp genome and combined it with previously published genomes for other Chinese Polyspora species. The results showed that cp genomes of six Chinese Polyspora varied in length between 156,452 bp (P. chrysandra) and 157,066 bp (P. speciosa), but all contained 132 genes, with GC content of 37.3%, and highly similar genes distribution and codon usage. A total of eleven intergenic spacer regions were found having the highest levels of divergence, and eight divergence hotspots were identified as molecular markers for Phylogeography and genetic diversity studies in Polyspora. Gene selection pressure suggested that five genes were subjected to positive selection. Phylogenetic relationships among Polyspora species based on the complete cp genomes were supported strongly, indicating that the cp genomes have the potential to be used as super barcodes for further analysis of the phylogeny of the entire genus. The cp genomes of Chinese Polyspora species will provide valuable information for species identification, molecular breeding and evolutionary analysis of genus Polyspora.

Comparative Analysis and Phylogenetic Relationships of Ceriops Species (Rhizophoraceae) and Avicennia lanata (Acanthaceae): Insight into the Chloroplast Genome Evolution between Middle and Seaward Zones of Mangrove Forests

Ceriops and Avicennia are true mangroves in the middle and seaward zones of mangrove forests, respectively. The chloroplast genomes of Ceriops decandra, Ceriops zippeliana, and Ceriops tagal were assembled into lengths of 166,650, 166,083 and 164,432 bp, respectively, whereas Avicennia lanata was 148,264 bp in length. The gene content and gene order are highly conserved among these species. The chloroplast genome contains 125 genes in A. lanata and 129 genes in Ceriops species. Three duplicate genes (rpl2, rpl23, and trnM-CAU) were found in the IR regions of the three Ceriops species, resulting in expansion of the IR regions. The rpl32 gene was lost in C. zippeliana, whereas the infA gene was present in A. lanata. Short repeats (<40 bp) and a lower number of SSRs were found in A. lanata but not in Ceriops species. The phylogenetic analysis supports that all Ceriops species are clustered in Rhizophoraceae and A. lanata is in Acanthaceae. In a search for genes under selective pressures of coastal environments, the rps7 gene was under positive selection compared with non-mangrove species. Finally, two specific primer sets were developed for species identification of the three Ceriops species. Thus, this finding provides insightful genetic information for evolutionary relationships and molecular markers in Ceriops and Avicennia species.

Genetic and molecular mechanisms underlying mangrove adaptations to intertidal environments

Mangroves are halophytic plants belonging to diverse angiosperm families that are adapted to highly stressful intertidal zones between land and sea. They are special, unique, and one of the most productive ecosystems that play enormous ecological roles and provide a large number of benefits to the coastal communities. To thrive under highly stressful conditions, mangroves have innovated several key morphological, anatomical, and physio-biochemical adaptations. The evolution of the unique adaptive modifications might have resulted from a host of genetic and molecular changes and to date we know little about the nature of these genetic and molecular changes. Although slow, new information has accumulated over the last few decades on the genetic and molecular regulation of the mangrove adaptations, a comprehensive review on it is not yet available. This review provides up-to-date consolidated information on the genetic, epigenetic, and molecular regulation of mangrove adaptive traits.

Comparative genomics and phylogenetic relationships of two endemic and endangered species (Handeliodendron bodinieri and Eurycorymbus cavaleriei) of two monotypic genera within Sapindales

BACKGROUND

Handeliodendron Rehder and Eurycorymbus Hand.-Mazz. are the monotypic genera in the Sapindaceae family. The phylogenetic relationship of these endangered species Handeliodendron bodinieri (Lévl.) Rehd. and Eurycorymbus cavaleriei (Lévl.) Rehd. et Hand.-Mazz. with other members of Sapindaceae s.l. is not well resolved. A previous study concluded that the genus Aesculus might be paraphyletic because Handeliodendron was nested within it based on small DNA fragments. Thus, their chloroplast genomic information and comparative genomic analysis with other Sapindaceae species are necessary and crucial to understand the circumscription and plastome evolution of this family.

RESULTS

The chloroplast genome sizes of Handeliodendron bodinieri and Eurycorymbus cavaleriei are 151,271 and 158,690 bp, respectively. Results showed that a total of 114 unique genes were annotated in H. bodinieri and E. cavaleriei, and the ycf1 gene contained abundant SSRs in both genomes. Comparative analysis revealed that gene content, PCGs, and total GC content were remarkably similar or identical within 13 genera from Sapindaceae, and the chloroplast genome size of four genera was generally smaller within the family, including Acer, Dipteronia, Aesculus, and Handeliodendron. IR boundaries of the H. bodinieri showed a significant contraction, whereas it presented a notable expansion in E. cavaleriei cp genome. Ycf1, ndhC-trnV-UAC, and rpl32-trnL-UAG-ccsA were remarkably divergent regions in the Sapindaceae species. Analysis of selection pressure showed that there are a few positively selected genes. Phylogenetic analysis based on different datasets, including whole chloroplast genome sequences, coding sequences, large single-copy, small single-copy, and inverted repeat regions, consistently demonstrated that H. bodinieri was sister to the clade consisting of Aesculus chinensis and A. wangii and strongly support Eurycorymbus cavaleriei as sister to Dodonaea viscosa.

CONCLUSION

This study revealed that the cp genome size of the Hippocastanoideae was generally smaller compared to the other subfamilies within Sapindaceae, and three highly divergent regions could be used as the specific DNA barcodes within Sapindaceae. Phylogenetic results strongly support that the subdivision of four subfamilies within Sapindaceae, and Handeliodendron is not nested within the genus Aesculus.

Complete chloroplast genome sequences of five Bruguiera species (Rhizophoraceae): comparative analysis and phylogenetic relationships

Bruguiera is a genus of true mangroves that are mostly distributed in the Indo-West Pacific region. However, the number of published whole chloroplast genome sequences of Bruguiera species are limited. Here, the complete chloroplast sequences of five Bruguiera species were sequenced and assembled using Illumina data. The chloroplast genomes of B. gymnorhiza, B. hainesii, B. cylindrica, B. parviflora and B. sexangula were assembled into 161,195, 164,295, 164,297, 163,228 and 164,170 bp, respectively. All chloroplast genomes contain 37 tRNA and eight rRNA genes, with either 84 or 85 protein-coding genes. A comparative analysis of these genomes revealed high similarity in gene structure, gene order and boundary position of the LSC, SSC and two IR regions. Interestingly, B. gymnorhiza lost a rpl32 gene in the SSC region. In addition, a ndhF gene in B. parviflora straddles both the SSC and IRB boundary regions. These genes reveal differences in chloroplast evolution among Bruguiera species. Repeats and SSRs in the chloroplast genome sequences were found to be highly conserved between B. cylindrica and B. hainesii as well as B. gymnorhiza and B. sexangula indicating close genetic relationships based on maternal inheritance. Notably, B. hainesii, which is considered a hybrid between B. gymnorhiza and B. cylindrica, appears to have inherited the chloroplast from B. cylindrica. Investigating the effects of selection events on shared protein-coding genes showed a positive selection in rps7 and rpl36 genes in all species compared to land-plant species. A phylogenetic analysis, based on 59 conserved chloroplast protein-coding genes, showed strong support that all Bruguiera species are in the clade Rhizophoraceae. This study provides valuable genetic information for the study of evolutionary relationships and population genetics in Bruguiera and other mangrove species.

Complete chloroplast genome of a semi-mangrove plant Hibiscus tiliaceus (Malvaceae)

Hibiscus tiliaceus is a semi-mangrove species that is widely distributed in tropical and subtropical coastal areas around the world. Here, the complete chloroplast (cp) genome sequence of H. tiliaceus was assembled and characterized. The cp genome was 161,748 bp in length, consisting of a large single copy (LSC) region of 89,190 bp and a small single copy (SSC) region of 19,616 bp, which were separated by a pair of 26,471 bp inverted repeat (IR) regions. The overall GC content was 36.88%. A total of 131 genes, including 85 protein-coding genes, 37 tRNA genes and 8 rRNA genes were identified. Phylogenetic tree reconstructed by 15 complete cp genomes revealed that H. tiliaceus was sister to the congeneric species H. cannabinus.

Characterization of the complete chloroplast genome of mangrove Rhizophora apiculata Blume (Rhizophoraceae)

The chloroplast (cp) genome sequence of Rhizophora apiculata was characterized. The cp genome length was 164,343 bp in length, containing a typical structure of a large single copy (LSC) of 93,155 bp, a small single copy (SSC) of 19,376 bp, and two inverted repeats (IRs) of 25,906 bp, with a GC content of 34.9%. There were 131 genes were annotated in the cp genome, including 85 protein-coding genes, 38 tRNA genes, and 8 rRNA genes. A phylogenetic analysis using cp genomes of mangroves and ecologically associated species resolved R. apiculata in Rhizophora with R. stylosa and R. x lamarckii. This complete chloroplast sequence offers a promising tool for further species identification and evolutionary studies of Rhizophora, as well as for mangroves.

Characterization of the complete chloroplast genome of mangrove Bruguiera gymnorrhiza (L.) Lam. ex Savigny

The chloroplast (cp) genome sequence of Bruguiera gymnorrhiza was characterized. The cp genome length was 163,795 bp in length, with a GC content of 35.3%, containing a large single copy (LSC) of 90,830 bp, a small single copy (SSC) of 20,207 bp, and a pair of inverted repeats (IRs) of 26,379 bp. The genome contained 121 genes, including 84 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. A phylogenetic analysis using cp genomes of mangroves and ecologically associated species resolved B. gymnorrhiza in Bruguiera with B. sexangula var. rhynchopetala. This complete chloroplast sequence offers a promising tool for further species identification and evolutionary studies of Bruguiera, as well as for mangroves.

Complete chloroplast genome sequence and phylogenetic analysis of winter oil rapeseed (Brassica rapa L.)

Winter oil rapeseed ‘18 R-1’ (Brassica rapa L.) is a new variety that can survive in northern China where the extreme low temperature is −20 °C to −32 °C. It is different from traditional B. rapa and Brassica napus. In this study, the complete chloroplast (cp) genome of ‘18 R-1’ was sequenced and analyzed to assess the genetic relationship. The size of cp genome is 153,494 bp, including one large single copy (LSC) region of 83,280 bp and one small single copy (SSC) region of 17,776 bp, separated by two inverted repeat (IR) regions of 26,219 bp. The GC content of the whole genome is 36.35%, while those of LSC, SSC, and IR are 34.12%, 29.20%, and 42.32%, respectively. The cp genome encodes 132 genes, including 87 protein-coding genes, eight rRNA genes, and 37 tRNA genes. In repeat structure analysis, 288 simple sequence repeats (SSRs) were identified. Cp genome of ‘18 R-1’ was closely related to Brassica chinensis, B. rapa and Brassica pekinesis.