Comparative analysis of four Zantedeschia chloroplast genomes: expansion and contraction of the IR region, phylogenetic analyses and SSR genetic diversity assessment

Sequence characteristics, genetic diversity and phylogenetic analysis of the Cucurbita ficifolia (Cucurbitaceae) chloroplasts genome

BACKGROUND

Curcubita ficifolia Bouché (Cucurbitaceae) has high value as a food crop and medicinal plant, and also has horticultural value as rootstock for other melon species. China is home to many different cultivars, but the genetic diversity of these resources and the evolutionary relationships among them, as well as the differences between C. ficifolia and other Cucurbita species, remain unclear.

RESULTS

We investigated the chloroplast (cp) genomes of 160 C. ficifolia individuals from 31 populations in Yunnan, a major C. ficifolia production area in China. We found that the cp genome of C. ficifolia is ~151 kb and contains 128 genes, of which 86 are protein coding genes, 34 encode tRNA, and eight encode rRNAs. We also identified 64 SSRs, mainly AT repeats. The cp genome was found to contain a total of 204 SNP and 57 indels, and a total of 21 haplotypes were found in the 160 study individuals. The reverse repeat (IR) region of C. ficifolia contained a few differences compared with this region in the six other Cucurbita species. Sequence difference analysis demonstrated that most of the variable regions were concentrated in the single copy (SC) region. Moreover, the sequences of the coding regions were found to be more similar among species than those of the non-coding regions. The phylogenies reconstructed from the cp genomes of 61 representative species of Cucurbitaceae reflected the currently accepted classification, in which C. ficifolia is sister to the other Cucurbita species, however, different interspecific relationships were found between Cucurbita species.

CONCLUSIONS

These results will be valuable in the classification of C. ficifolia genetic resources and will contribute to our understanding of evolutionary relationships within the genus Cucurbita.

Comparative analysis of plastid genomes reveals rearrangements, repetitive sequence features, and phylogeny in the Annonaceae

The Annonaceae stands as the most species rich family in the Magnoliales, a basal group of angiosperms. Widely distributed in tropical and subtropical regions, it holds significant ecological and economic value. The plastid genome (plastome) is often employed in studies related to plant phylogenetics, comparative genomics, evolutionary biology, and genetic engineering. Nonetheless, research progress on plastid genomics in the Annonaceae has been relatively slow. In this study, we analyzed the structure and repetitive sequence features of plastomes from 28 Annonaceae species. Among them, Mitrephora tomentosa and Desmos chinensis were newly sequenced, with sizes of 160,157 bp and 192,167 bp, and GC contents of 38.3% and 38.4%, respectively. The plastome size in the Annonaceae ranged from 158,837 bp to 202,703 bp, with inverted repeat (IR) region sizes ranging from 64,621 bp to 25,861 bp. Species exhibiting expansion in the IR region showed an increase in plastome size and gene number, frequent boundary changes, different expansion modes (bidirectional or unidirectional), and an increase in repetitive sequences. Specifically, a large number of dispersed repetitive sequences lead to an increase in the size of the LSC region in Goniothalamus tamirensis. Phylogenetic analysis revealed Annonoideae and Malmeoideae as monophyletic groups and sister clades, with Cananga odorata outside of them, followed by Anaxagorea javanica. This research uncovers the structural variation characteristics of plastomes in the Annonaceae, providing valuable information for understanding the phylogeny and plastome evolution of Annonaceae.

Molecular Characterization and Phylogenetic Analysis of Centipedegrass [Eremochloa ophiuroides (Munro) Hack.] Based on the Complete Chloroplast Genome Sequence

Centipedegrass (Eremochloa ophiuroides) is an important warm-season grass plant used as a turfgrass as well as pasture grass in tropical and subtropical regions, with wide application in land surface greening and soil conservation in South China and southern United States. In this study, the complete cp genome of E. ophiuroides was assembled using high-throughput Illumina sequencing technology. The circle pseudomolecule for E. ophiuroides cp genome is 139,107 bp in length, with a quadripartite structure consisting of a large single copyregion of 82,081 bp and a small single copy region of 12,566 bp separated by a pair of inverted repeat regions of 22,230 bp each. The overall A + T content of the whole genome is 61.60%, showing an asymmetric nucleotide composition. The genome encodes a total of 131 gene species, composed of 20 duplicated genes within the IR regions and 111 unique genes comprising 77 protein-coding genes, 30 transfer RNA genes, and 4 ribosome RNA genes. The complete cp genome sequence contains 51 long repeats and 197 simple sequence repeats, and a high degree of collinearity among E. ophiuroide and other Gramineae plants was disclosed. Phylogenetic analysis showed E. ophiuroides, together with the other two Eremochloa species, is closely related to Mnesithea helferi within the subtribe Rottboelliinae. These findings will be beneficial for the classification and identification of the Eremochloa taxa, phylogenetic resolution, novel gene discovery, and functional genomic studies for the genus Eremochloa.

The complete chloroplast genome sequence of Amorphophallus konjac (Araceae) from Yunnan, China and its phylogenetic analysis in the family Araceae

This work determined and analyzed the complete chloroplast genome sequence of Amorphophallus konjac K. Koch ex N.E.Br 1858 from Yunnan, China. The genome size was 167,470 bp, of which contains a large single-copy region (LSC 93,443 bp), a small single-copy region (SSC 21,575 bp), and a pair of inverted repeat regions (IR 26,226 bp). The chloroplast genome has 131 genes, including 86 protein-coding genes, 37 tRNAs, and eight rRNAs. A previous study reported deletion of accD, psbE, and trnG-GCC genes in the A. konjac chloroplast genome. Our study supports the conservative structure of A. konjac and does not support the gene deletion mentioned above. Phylogenetic analysis indicated that A. konjac shares a close relationship with another A. konjac (collected from Guizhou) and A. titanium by forming a clade in the genus Amorphophallus. Our results provide some useful information to the evolution of the family Araceae.

Phylogenetics Study to Compare Chloroplast Genomes in Four Magnoliaceae Species

Magnoliaceae, a family of perennial woody plants, contains several endangered species whose taxonomic status remains ambiguous. The study of chloroplast genome information can help in the protection of Magnoliaceae plants and confirmation of their phylogenetic relationships. In this study, the chloroplast genomes were sequenced, assembled, and annotated in Woonyoungia septentrionalis and three Michelia species (Michelia champaca, Michelia figo, and Michelia macclurei). Comparative analyses of genomic characteristics, repetitive sequences, and sequence differences were performed among the four Magnoliaceae plants, and phylogenetic relationships were constructed with twenty different magnolia species. The length of the chloroplast genomes varied among the four studied species ranging from 159,838 bp (Woonyoungia septentrionalis) to 160,127 bp (Michelia macclurei). Four distinct hotspot regions were identified based on nucleotide polymorphism analysis. They were petA-psbJ, psbJ-psbE, ndhD-ndhE, and rps15-ycf1. These gene fragments may be developed and utilized as new molecular marker primers. By using Liriodendron tulipifera and Liriodendron chinense as outgroups reference, a phylogenetic tree of the four Magnoliaceae species and eighteen other Magnoliaceae species was constructed with the method of Shared Coding Sequences (CDS). Results showed that the endangered species, W. septentrionalis, is relatively genetically distinct from the other three species, indicating the different phylogenetic processes among Magnoliaceae plants. Therefore, further genetic information is required to determine the relationships within Magnoliaceae. Overall, complete chloroplast genome sequences for four Magnoliaceae species reported in this paper have shed more light on phylogenetic relationships within the botanical group.

Development of Chloroplast Microsatellite Markers and Evaluation of Genetic Diversity and Population Structure of Cutleaf Groundcherry (Physalis angulata L.) in China

Cutleaf groundcherry (Physalis angulata L.), an annual plant containing a variety of active ingredients, has great medicinal value. However, studies on the genetic diversity and population structure of P. angulata are limited. In this study, we developed chloroplast microsatellite (cpSSR) markers and applied them to evaluate the genetic diversity and population structure of P. angulata. A total of 57 cpSSRs were identified from the chloroplast genome of P. angulata. Among all cpSSR loci, mononucleotide markers were the most abundant (68.24%), followed by tetranucleotide (12.28%), dinucleotide (10.53%), and trinucleotide (8.77%) markers. In total, 30 newly developed cpSSR markers with rich polymorphism and good stability were selected for further genetic diversity and population structure analyses. These cpSSRs amplified a total of 156 alleles, 132 (84.62%) of which were polymorphic. The percentage of polymorphic alleles and the average polymorphic information content (PIC) value of the cpSSRs were 81.29% and 0.830, respectively. Population genetic diversity analysis indicated that the average observed number of alleles (Na), number of effective alleles (He), Nei's gene diversity (h), and Shannon information indices (I) of 16 P. angulata populations were 1.3161, 1.1754, 0.1023, and 0.1538, respectively. Moreover, unweighted group arithmetic mean, neighbor-joining, principal coordinate, and STRUCTURE analyses indicated that 203 P. angulata individuals from 16 populations were grouped into four clusters. A molecular variance analysis (AMOVA) illustrated the considerable genetic variation among populations, while the gene flow (Nm) value (0.2324) indicated a low level of gene flow among populations. Our study not only provided a batch of efficient genetic markers for research on P. angulata but also laid an important foundation for the protection and genetic breeding of P. angulata resources.

Multiple domestication events explain the origin of Gossypium hirsutum landraces in Mexico

Several Mesoamerican crops constitute wild-to-domesticated complexes generated by multiple initial domestication events, and continuous gene flow among crop populations and between these populations and their wild relatives. It has been suggested that the domestication of cotton (Gossypium hirsutum) started in the northwest of the Yucatán Peninsula, from where it spread to other regions inside and outside of Mexico. We tested this hypothesis by assembling chloroplast genomes of 23 wild, landraces, and breeding lines (transgene-introgressed and conventional). The phylogenetic analysis showed that the evolutionary history of cotton in Mexico involves multiple events of introgression and genetic divergence. From this, we conclude that Mexican landraces arose from multiple wild populations. Our results also revealed that their structural and functional chloroplast organizations had been preserved. However, genetic diversity decreases as a consequence of domestication, mainly in transgene-introgressed (TI) individuals (π = 0.00020, 0.00001, 0.00016, 0, and 0, of wild, TI-wild, landraces, TI-landraces, and breeding lines, respectively). We identified homologous regions that differentiate wild from domesticated plants and indicate a relationship among the samples. A decrease in genetic diversity associated with transgene introgression in cotton was identified for the first time, and our outcomes are therefore relevant to both biosecurity and agrobiodiversity conservation.

CPJSdraw: analysis and visualization of junction sites of chloroplast genomes

Background

Chloroplast genomes are usually circular molecules, and most of them are tetrad structures with two inverted repeat (IR) regions, a large single-copy region, and a small single-copy region. IR contraction and expansion are among the genetic diversities during the evolution of plant chloroplast genomes. The only previously released tool for the visualization of junction sites of the regions does not consider the diversity of the starting point of genomes, which leads to incorrect results or even no results for the examination of IR contraction and expansion.

Results

In this work, a new tool named CPJSdraw was developed for visualizing the junction sites of chloroplast genomes. CPJSdraw can format the starting point of the irregular linearized genome, correct the junction sites of IR and single-copy regions, display the tetrad structure, visualize the junction sites of any number (≥1) of chloroplast genomes, show the transcription direction of genes adjacent to junction sites, and indicate the IR expansion or contraction of chloroplast genomes.

Conclusions

CPJSdraw is a software that is universal and reliable in analysis and visualization of IR expansion or contraction of chloroplast genomes. CPJSdraw has more accurate analysis and more complete functions when compared with previously released tool. CPJSdraw as a perl package and tested data are available at http://dx.doi.org/10.5281/zenodo.7669480 for English users. In addition, an online version with a Chinese interface is available at http://cloud.genepioneer.com:9929/#/tool/alltool/detail/335.

Comparative chloroplast genomes and phylogenetic relationships of Aglaonema modestum and five variegated cultivars of Aglaonema

Aglaonema, commonly called Chinese evergreens, are widely used for ornamental purposes. However, attempts to identify Aglaonema species and cultivars based on leaf morphology have been challenging. In the present study, chloroplast sequences were used to elucidate the phylogenetic relationships of cultivated Aglaonema in South China. The chloroplast genomes of one green species and five variegated cultivars of Aglaonema, Aglaonema modestum, ‘Red Valentine’, ‘Lady Valentine’, ‘Hong Yan’, ‘Hong Jian’, and ‘Red Vein’, were sequenced for comparative and phylogenetic analyses. The six chloroplast genomes of Aglaonema had typical quadripartite structures, comprising a large single copy (LSC) region (91,092–91,769 bp), a small single copy (SSC) region (20,816–26,501 bp), and a pair of inverted repeat (IR) regions (21,703–26,732 bp). The genomes contained 112 different genes, including 79–80 protein coding genes, 28–29 tRNAs and 4 rRNAs. The molecular structure, gene order, content, codon usage, long repeats, and simple sequence repeats (SSRs) were generally conserved among the six sequenced genomes, but the IR-SSC boundary regions were significantly different, and ‘Red Vein’ had a distinct long repeat number and type frequency. For comparative and phylogenetic analyses, Aglaonema costatum was included; it was obtained from the GenBank database. Single-nucleotide polymorphisms (SNPs) and insertions/deletions (indels) were determined among the seven Aglaonema genomes studied. Nine divergent hotspots were identified: trnH-GUG-CDS1_psbA, trnS-GCU_trnS-CGA-CDS1, rps4-trnT-UGU, trnF-GAA-ndhJ, petD-CDS2-rpoA, ycf1-ndhF, rps15-ycf1-D2, ccsA-ndhD, and trnY-GUA-trnE-UUC. Additionally, positive selection was found for rpl2, rps2, rps3, ycf1 and ycf2 based on the analyses of Ka/Ks ratios among 16 Araceae chloroplast genomes. The phylogenetic tree based on whole chloroplast genomes strongly supported monophyletic Aglaonema and clear relationships among Aroideae, Lasioideae, Lemnoideae, Monsteroideae, Orontioideae, Pothoideae and Zamioculcadoideae in the family Araceae. By contrast, protein coding gene phylogenies were poorly to strongly supported and incongruent with the whole chloroplast genome phylogenetic tree. This study provided valuable genome resources and helped identify Aglaonema species and cultivars.

The entire chloroplast genome sequence of Asparagus cochinchinensis and genetic comparison to Asparagus species

Asparagus cochinchinensis is a traditional Chinese medicinal plant. The chloroplast (cp) genome study on A. cochinchinensis is poorly understood. In this research, we collected the data from the cp genome assembly and gene annotation of A. cochinchinensis, followed by further comparative analysis with six species in the genus Asparagus. The cp genome of A. cochinchinensis showed a circular quadripartite structure in the size of 157,095 bp, comprising a large single-copy (LSC), a small single-copy (SSC), and two inverted repeat (IR) regions. A total of 137 genes were annotated, consisting of 86 protein-coding genes, 8 ribosomal RNAs, 38 transfer RNAs, and 5 pseudo-genes. Forty scattered repetitive sequences and 247 simple sequence repeats loci were marked out. In addition, A/T-ending codons were shown to have a basis in the codon analysis. A cp genome comparative analysis revealed that a similar gene composition was detected in the IR and LSC/SSC regions with Asparagus species. Based on the complete cp genome sequence in Asparagaceae, the result showed that A. cochinchinensis was closely related to A. racemosus by phylogenetic analysis. Therefore, our study providing A. cochinchinensis genomic resources could effectively contribute to the phylogenetic analysis and molecular identification of the genus Asparagus.

Comparison Analysis Based on Complete Chloroplast Genomes and Insights into Plastid Phylogenomic of Four Iris Species

Iris species, commonly known as rainbow flowers because of their attractive flowers, are extensively grown in landscape gardens. A few species, including Belamcanda chinensis, the synonym of I. domestica and I. tectorum, are known for their medicinal properties. However, research on the genomes and evolutionary relationships of Iris species is scarce. In the current study, the complete chloroplast (CP) genomes of I. tectorum, I. dichotoma, I. japonica, and I. domestica were sequenced and compared for their identification and relationship. The CP genomes of the four Iris species were circular quadripartite with similar lengths, GC contents, and codon usages. A total of 113 specific genes were annotated, including the ycf1 pseudogene in all species and rps19 in I. japonica alone. All the species had mononucleotide (A/T) simple sequence repeats (SSRs) and long forward and palindromic repeats in their genomes. A comparison of the CP genomes based on mVISTA and nucleotide diversity (Pi) identified three highly variable regions (ndhF-rpl32, rps15-ycf1, and rpl16). Phylogenetic analysis based on the complete CP genomes concluded that I. tectorum is a sister of I. japonica, and the subgenus Pardanthopsis with several I. domestica clustered into one branch is a sister of I. dichotoma. These findings confirm the feasibility of superbarcodes (complete CP genomes) for Iris species authentication and could serve as a resource for further research on Iris phylogeny.

Comparative and Phylogenetic Analysis of Complete Plastomes among Aristidoideae Species (Poaceae)

Aristidoideae is a subfamily in the PACMAD clade of family Poaceae, including three genera, Aristida, Stipagrostis, and Sartidia. In this study, the plastomes of Aristida adscensionis and Stipagrostis pennata were newly sequenced, and a total of 16 Aristidoideae plastomes were compared. All plastomes were conservative in genome size, gene number, structure, and IR boundary. Repeat sequence analysis showed that forward and palindrome repeats were the most common repeat types. The number of SSRs ranged from 30 (Sartidia isaloensis) to 54 (Aristida purpurea). Codon usage analysis showed that plastome genes preferred to use codons ending with A/T. A total of 12 highly variable regions were screened, including four protein coding sequences (matK, ndhF, infA, and rpl32) and eight non-coding sequences (rpl16-1-rpl16-2, ccsA-ndhD, trnY-GUA-trnD-GUC, ndhF-rpl32, petN-trnC-GCA, trnT-GGU-trnE-UUC, trnG-GCC-trnfM-CAU, and rpl32-trnL-UAG). Furthermore, the phylogenetic position of this subfamily and their intergeneric relationships need to be illuminated. All Maximum Likelihood and Bayesian Inference trees strongly support the monophyly of Aristidoideae and each of three genera, and the clade of Aristidoideae and Panicoideae was a sister to other subfamilies in the PACMAD clade. Within Aristidoideae, Aristida is a sister to the clade composed of Stipagrostis and Sartidia. The divergence between C₄ Stipagrostis and C₃ Sartidia was estimated at 11.04 Ma, which may be associated with the drought event in the Miocene period. Finally, the differences in carbon fixation patterns, geographical distributions, and ploidy may be related to the difference of species numbers among these three genera. This study provides insights into the phylogeny and evolution of the subfamily Aristidoideae.

Phylogenetic relationships and characterization of the complete chloroplast genome of Rosa sterilis

Rosa sterilis is an economically and important fruit that is extensively grown in Southwestern China. In this study, we determined the complete chloroplast genome of R. sterilis using high-throughput Illumina sequencing. The chloroplast genome of R. sterilis is 156,561 bp in size, containing a large single-copy region (LSC)(85,701 bp), a small single-copy region (SSC) (18,746 bp), and a pair of inverted repeat (IR) regions (each one of 26,057 bp). The overall GC content of the chloroplast genome is 37.23%, while the corresponding values of GC contents of the LSC, SSC, and IR regions are 35.20%, 31.37%, and 42.70%, respectively. The chloroplast genome of R. sterilis contains 130 genes, including 84 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. The phylogenetic maximum-likelihood tree revealed that Rosa chinensis or Rosa chinensis var. spontanea is the closest related to R. sterilis in the phylogenetic relationship. This complete chloroplast genome can be further used for genomic studies, evolutionary analyses, and genetic engineering studies of the family Rosaceae.

Variation in Plastome Sizes Accompanied by Evolutionary History in Monogenomic Triticeae (Poaceae: Triticeae)

To investigate the pattern of chloroplast genome variation in Triticeae, we comprehensively analyzed the indels in protein-coding genes and intergenic sequence, gene loss/pseudonization, intron variation, expansion/contraction in inverted repeat regions, and the relationship between sequence characteristics and chloroplast genome size in 34 monogenomic Triticeae plants. Ancestral genome reconstruction suggests that major length variations occurred in four-stem branches of monogenomic Triticeae followed by independent changes in each genus. It was shown that the chloroplast genome sizes of monogenomic Triticeae were highly variable. The chloroplast genome of Pseudoroegneria, Dasypyrum, Lophopyrum, Thinopyrum, Eremopyrum, Agropyron, Australopyrum, and Henradia in Triticeae had evolved toward size reduction largely because of pseudogenes elimination events and length deletion fragments in intergenic. The Aegilops/Triticum complex, Taeniatherum, Secale, Crithopsis, Herteranthelium, and Hordeum in Triticeae had a larger chloroplast genome size. The large size variation in major lineages and their subclades are most likely consequences of adaptive processes since these variations were significantly correlated with divergence time and historical climatic changes. We also found that several intergenic regions, such as petN-trnC and psbE-petL containing unique genetic information, which can be used as important tools to identify the maternal relationship among Triticeae species. Our results contribute to the novel knowledge of plastid genome evolution in Triticeae.