A review of the prevalence, utility, and caveats of using chloroplast simple sequence repeats for studies of plant biology

Phylogenomics and plastome evolution of Lithospermeae (Boraginaceae)

BACKGROUND

Lithospermeae is the largest tribe within Boraginaceae. The tribe has been the focus of multiple phylogenetic studies over the last 15 years, with most focused on one genus or a few genera. In the present study, we newly sequenced 69 species of Lithospermeae and relatives to analyze the phylogenomic relationships among its members as well as the evolution of the plastid genome.

RESULTS

The phylogeny of Lithospermeae resolved from the plastid genome and nrDNA cistron is generally congruent with prior studies, but is better resolved and supported. Increasing character sampling across the plastid genome results in gradually more similar trees to that from the entire plastid genome. Overall, plastid genome structure was quite consistent across Lithospermeae. Codon Usage Bias (CUB) analyses demonstrate that across Lithospermeae plastid genomes were rich in AT and poor in GC. Mutation may play a greater role than selection across the plastid genome of Lithospermeae. The present study is the first to highlight the CUB characteristics of Lithospermeae species, which can help elucidate the mechanisms underlying patterns of molecular evolution and improve the expression levels of exogenous genes by codon optimization.

CONCLUSIONS

This study provides a comprehensive phylogenomic analysis of Lithospermeae, significantly enhancing our understanding of the phylogenetic relationships and plastid genome evolution within this largest tribe of Boraginaceae. By utilizing an expanded genomic sampling approach, we have achieved increased resolution and support among the evolutionary relationships of the tribe, in line with but improving upon previous studies. The analyses of plastid genome structure revealed consistency across Lithospermeae, with a notable CUB. This study marks the first investigation into the CUB of Lithospermeae species and sets the stage for further research on the molecular evolution of plastid genomes across Boraginaceae.

The plastomes of Lepismium cruciforme (Vell.) Miq and Schlumbergera truncata (Haw.) Moran reveal tribe-specific rearrangements and the first loss of the trnT-GGU gene in Cactaceae

BACKGROUND

Recent studies have revealed atypical features in the plastomes of the family Cactaceae, the largest lineage of succulent species adapted to arid and semi-arid regions. Most plastomes sequenced to date are from short-globose and cylindrical cacti, while little is known about plastomes of epiphytic cacti. Published cactus plastomes reveal reduction and complete loss of IRs, loss of genes, pseudogenization, and even degeneration of tRNA structures. Aiming to contribute with new insights into the plastid evolution of Cactaceae, particularly within the tribe Rhipsalideae, we de novo assembled and analyzed the plastomes of Lepismium cruciforme and Schlumbergera truncata, two South American epiphytic cacti.

METHODS AND RESULTS

Our data reveal many gene losses in both plastomes and the first loss of functionality of the trnT-GGU gene in Cactaceae. The trnT-GGU is a pseudogene in L. cruciforme plastome and appears to be degenerating in the tribe Rhipsalideae. Although the plastome structure is conserved among the species of the tribe Rhipsalideae, with tribe-specific rearrangements, we mapped around 200 simple sequence repeats and identified nine nucleotide polymorphism hotspots, useful to improve the phylogenetic resolutions of the Rhipsalideae. Furthermore, our analysis indicated high gene divergence and rapid evolution of RNA editing sites in plastid protein-coding genes in Cactaceae.

CONCLUSIONS

Our findings show that some characteristics of the Rhipsalideae tribe are conserved, such as plastome structure with IRs containing only the ycf2 and two tRNA genes, structural degeneration of the trnT-GGU gene and ndh complex, and lastly, pseudogenization of rpl33 and rpl23 genes, both plastid translation-related genes.

Chloroplast and whole-genome sequencing shed light on the evolutionary history and phenotypic diversification of peanuts

Cultivated peanut (Arachis hypogaea L.) is a widely grown oilseed crop worldwide; however, the events leading to its origin and diversification are not fully understood. Here by combining chloroplast and whole-genome sequence data from a large germplasm collection, we show that the two subspecies of A. hypogaea (hypogaea and fastigiata) likely arose from distinct allopolyploidization and domestication events. Peanut genetic clusters were then differentiated in relation to dissemination routes and breeding efforts. A combination of linkage mapping and genome-wide association studies allowed us to characterize genes and genomic regions related to main peanut morpho-agronomic traits, namely flowering pattern, inner tegument color, growth habit, pod/seed weight and oil content. Together, our findings shed light on the evolutionary history and phenotypic diversification of peanuts and might be of broad interest to plant breeders.

The newly assembled chloroplast genome of Aeluropus littoralis: molecular feature characterization and phylogenetic analysis with related species

Aeluropus littoralis, a halophyte grass, is widely distributed from the Mediterranean to the Indian subcontinent through the Mongolian Gobi. This model halophyte has garnered increasing attention owing to its use as forage and its high tolerance to environmental stressors. The chloroplast genomes of many plants have been extensively examined for molecular, phylogenetic and transplastomic applications. However, no published research on the A. littoralis chloroplast (cp) genome was discovered. Here, the entire chloroplast genome of A. littoralis was assembled implementing accurate long-read sequences. The entire chloroplast genome, with an estimated length of 135,532 bp (GC content: 38.2%), has a quadripartite architecture and includes a pair of inverted repeat (IR) regions, IRa and IRb (21,012 bp each), separated by a large and a small single-copy regions (80,823 and 12,685 bp, respectively). The features of A. littoralis consist of 133 genes that synthesize 87 peptides, 38 transfer RNAs, and 8 ribosomal RNAs. Of these genes, 86 were unique, whereas 19 were duplicated in IR regions. Additionally, a total of forty-six simple sequence repeats, categorized into 32-mono, four-di, two-tri, and eight-tetranucleotides, were discovered. Furthermore, ten sets of repeats greater than 20 bp were located primarily in the LSC region. Evolutionary analysis based on chloroplast sequence data revealed that A. littoralis with A. lagopoides and A. sinensis belong to the Aeluropodinae subtribe, which is a sister to the Eleusininae in the tribe Cynodonteae and the subfamily Chloridoideae. This subfamily belongs to the PACMAD clade, which contains the majority of the C4 photosynthetic plants in the Poaceae. The newly constructed A. littoralis cp genome offers valuable knowledge for DNA barcoding, phylogenetic, transplastomic research, and other biological studies.

Complete chloroplast genomes of six neotropical palm species, structural comparison, and evolutionary dynamic patterns

The Arecaceae family has a worldwide distribution, especially in tropical and subtropical regions. We sequenced the chloroplast genomes of Acrocomia intumescens and A. totai, widely used in the food and energy industries; Bactris gasipaes, important for palm heart; Copernicia alba and C. prunifera, worldwide known for wax utilization; and Syagrus romanzoffiana, of great ornamental potential. Copernicia spp. showed the largest chloroplast genomes (C. prunifera: 157,323 bp and C. alba: 157,192 bp), while S. romanzoffiana and B. gasipaes var. gasipaes presented the smallest (155,078 bp and 155,604 bp). Structurally, great synteny was detected among palms. Conservation was also observed in the distribution of single sequence repeats (SSR). Copernicia spp. presented less dispersed repeats, without occurrence in the small single copy (SSC). All RNA editing sites were C (cytidine) to U (uridine) conversions. Overall, closely phylogenetically related species shared more sites. Almost all nodes of the phylogenetic analysis showed a posterior probability (PP) of 1.0, reaffirming the close relationship between Acrocomia species. These results elucidate the conservation among palm chloroplast genomes, but point to subtle structural changes, providing support for the evolutionary dynamics of the Arecaceae family.

New Insights Into The Evolution of Chloroplast Genomes in Ochna Species (Ochnaceae, Malpighiales)

Ochnaceae DC. includes more than 600 species that exhibit potential values for environmental ecology, ornamental, pharmaceutical, and timber industries. Although studies on phylogeny and phytochemicals have been intensively conducted, chloroplast genome data of Ochnaceae species have not been fully explored. In this study, the next-generation sequencing method was used to sequence the chloroplast genomes of Ochna integerrima and Ochna serrulata which were 157 329 and 157 835 bp in length, respectively. These chloroplast genomes had a quadripartite structure and contained 78 protein-coding genes, 30 tRNAs, and 4 rRNAs. Comparative analysis revealed 8 hypervariable regions, including trnK_UUU-trnQ_UUG, rpoB-psbM, trnS_GGA-rps4, accD-psaI, rpl33-rps18, rpl14-rpl16, ndhF-trnL_UAG, and rps15-ycf1 among 6 Ochnaceae taxa. Additionally, there were shared and unique repeats among 6 examined chloroplast genomes. The notable changes were the loss of rpl32 in Ochna species and the deletion of rps16 exon 2 in O. integerrima compared to other taxa. This study is the first comprehensive comparative genomic analysis of complete chloroplast genomes of Ochna species and related taxa in Ochnaceae. Consequently, the current study provides initial results for further research on genomic evolution, population genetics, and developing molecular markers in Ochnaceae and related taxa.

Chloroplast genome sequence of Nothapodytes nimmoniana (Icacinaceae), a vulnerable reservoir of camptothecin from Western Ghats

Nothapodytes nimmoniana belongs to family Icacinaceae and is a major source of compound Camptothecin. The global demand for Camptothecin has caused large-scale exploitation of N. nimmoniana from its wild habitat in Western Ghats of India, thereby making it vulnerable. The species is known to exhibit genetic diversity among the populations in Western Ghats. In this study, we report plastome sequence of N. nimmoniana, first for the genus. For the study, the species was collected from Western Ghats of Maharashtra. The plastome of N. nimmoniana was 150,726 bp in length and exhibited typical quadripartite structure with 83,771 bp LSC, 18,513 bp SSC and 24,221 IR region. The plastome was characterized by presence of 124 unique genes, 87 protein coding genes, 29 tRNA genes and four rRNA genes. Further, the plastome was compared with the available basal lamiid plastomes for gene order and composition. N. nimmoniana plastome exhibited SSC region in an inverted configuration. Phylogenomic study placed N. nimmoniana sister to Mappia mexicana. The SSR markers identified in this study, might help to distinguish genetically diverse populations, prioritizing the populations which need immediate conservation effects as well as for checking adulteration.

The chloroplast genome of black pepper (Piper nigrum L.) and its comparative analysis with related Piper species

Piper nigrum, also known as black pepper, is an economically and ecologically important crop of the genus Piper. It has been titled as the king of spices due to its wide consumption throughout the world. In the present investigation, the chloroplast genome of P. nigrum has been assembled from a whole genome sequence by integrating the short and long reads generated through Illumina and PacBio platforms, respectively. The chloroplast genome was observed to be 161,522 bp in size, having a quadripartite structure with a large single copy (LSC) region of 89,153 bp and a small single copy (SSC) region of 18,255 bp separated by a copy of inverted repeats (IRs), each 27,057 bp in length. Taking into consideration all the duplicated genes, a total of 131 genes were observed, which included 81 protein-coding genes, 37 tRNAs, 4 rRNAs, and 1 pseudogene. Individually, the LSC region consisted of 83 genes, the SSC region had 13 genes, and 18 genes were present in each IR region. Additionally, 216 SSRs were detected and 11 of these were validated through amplification in 12 species of Piper. The features of the chloroplast genome have been compared with those of the genus Piper. Our results provide useful insights into evolutionary and molecular studies of black pepper which will contribute to its further genetic improvement and breeding.

Analysis of Complete Chloroplast Genome: Structure, Phylogenetic Relationships of Galega orientalis and Evolutionary Inference of Galegeae

Galega orientalis, a leguminous herb in the Fabaceae family, is an ecologically and economically important species widely cultivated for its strong stress resistance and high protein content. However, genomic information of Galega orientalis has not been reported, which limiting its evolutionary analysis. The small genome size makes chloroplast relatively easy to obtain genomic sequence for phylogenetic studies and molecular marker development. Here, the chloroplast genome of Galega orientalis was sequenced and annotated. The results showed that the chloroplast genome of G. orientalis is 125,280 bp in length with GC content of 34.11%. A total of 107 genes were identified, including 74 protein-coding genes, 29 tRNAs and four rRNAs. One inverted repeat (IR) region was lost in the chloroplast genome of G. orientalis. In addition, five genes (rpl22, ycf2, rps16, trnE-UUC and pbf1) were lost compared with the chloroplast genome of its related species G. officinalis. A total of 84 long repeats and 68 simple sequence repeats were detected, which could be used as potential markers in the genetic studies of G. orientalis and related species. We found that the Ka/Ks values of three genes petL, rpl20, and ycf4 were higher than one in the pairwise comparation of G. officinalis and other three Galegeae species (Calophaca sinica, Caragana jubata, Caragana korshinskii), which indicated those three genes were under positive selection. A comparative genomic analysis of 15 Galegeae species showed that most conserved non-coding sequence regions and two genic regions (ycf1 and clpP) were highly divergent, which could be used as DNA barcodes for rapid and accurate species identification. Phylogenetic trees constructed based on the ycf1 and clpP genes confirmed the evolutionary relationships among Galegeae species. In addition, among the 15 Galegeae species analyzed, Galega orientalis had a unique 30-bp intron in the ycf1 gene and Tibetia liangshanensis lacked two introns in the clpP gene, which is contrary to existing conclusion that only Glycyrrhiza species in the IR lacking clade (IRLC) lack two introns. In conclusion, for the first time, the complete chloroplast genome of G. orientalis was determined and annotated, which could provide insights into the unsolved evolutionary relationships within the genus Galegeae.

The plastome of Melocactus glaucescens Buining & Brederoo reveals unique evolutionary features and loss of essential tRNA genes

The plastome of Melocactus glaucescens shows unique rearrangements, IR expansion, and unprecedented gene losses in Cactaceae. Our data indicate tRNA import from the cytosol to the plastids in this species. Cactaceae represents one of the richest families in keystone species of arid and semiarid biomes. This family shows various specific features comprehending morphology, anatomy, and metabolism, which allow them to grow under unfavorable environmental conditions. The subfamily Cactoideae contains the most divergence of species, which are highly variable in growth habit and morphology. This subfamily includes the endangered species Melocactus glaucescens (tribe Cereeae), which is a cactus endemic to the biome Caatinga in Brazil. Aiming to analyze the plastid evolution and develop molecular markers, we sequenced and analyzed in detail the plastome of M. glaucescens. Our analyses revealed that the M. glaucescens plastome is the most divergent among the species of the family Cactaceae sequenced so far. We characterized here unique rearrangements, expanded IRs containing an unusual set of genes, and several gene losses. Some genes related to the ndh complex were lost during the plastome evolution, while others have lost their functionality. Additionally, the loss of three tRNA genes (trnA-UGC, trnV-UAC, and trnV-GAC) suggests tRNA import from the cytosol to the plastids in M. glaucescens. Moreover, we identified high gene divergence, several putative positive signatures, and possible unique RNA-editing sites. Furthermore, we mapped 169 SSRs in the plastome of M. glaucescens, which are helpful to access the genetic diversity of natural populations and conservation strategies. Finally, our data provide new insights into the evolution of plastids in Cactaceae, which is an outstanding lineage adapted to extreme environmental conditions and a notorious example of the atypical evolution of plastomes.

Structural characterization of Platanthera ussuriensis chloroplast genome and comparative analyses with other species of Orchidaceae

BACKGROUND

The genus Tulotis has been classified into the genus Platanthera in the present taxonomic studies since the morphological characteristics of this genus is very similar to that of Platanthera. Platanthera ussuriensis, formerly named as Tulotis ussuriensis, is a small terrestrial orchid species and has been listed as wild plant under State protection (category II) in China. An improved understanding of the genomic information will enable future applications of conservation strategy as well as phylogenetic studies for this rare orchid species. The objective of this research was to characterize and compare the chloroplast genome of P. ussuriensis with other closely related species of Orchidaceae.

RESULTS

The chloroplast genome sequence of P. ussuriensis is 155,016 bp in length, which included a pair of inverted repeats (IRs) of 26,548 bp that separated a large single copy (LSC) region of 83,984 bp and a small single copy (SSC) region of 17,936 bp. The annotation contained a total of 132 genes, including 86 protein-coding genes, 38 tRNA genes and 8 rRNA genes. The simple sequence repeat (SSR) analysis showed that there were 104 SSRs in the chloroplast genome of P. ussuriensis. RNA editing sites recognition indicated 72 RNA editing events occurred, and all codon changes were C to T conversions. Comparative genomics showed that the chloroplast sequence of Platanthera related species were relatively conserved, while there were still some high variation regions that could be used as molecular markers. Moreover, Platanthera related species showed similar IR/SSC and IR/LSC borders. The phylogenetic analysis suggested that P. ussuriensis had a closer evolutionary relationship with P. japonica followed by the remaining Platanthera species.

CONCLUSION

Orchidaceae is a key group of biodiversity protection and also a hot spot group in the plant taxonomy and evolution studies due to their characteristics of high specialization and rapid evolution. This research determined the complete chloroplast genome of P. ussuriensis for the first time, and compared the sequence with other closely related orchid species. These results provide a foundation for future genomic and molecular evolution of the Orchidaceae species, and provide insights into the development of conservation strategy for Platanthera species.

Comparative and Phylogenetic Analysis Based on the Chloroplast Genome of Coleanthus subtilis (Tratt.) Seidel, a Protected Rare Species of Monotypic Genus

Coleanthus subtilis (Tratt.) Seidel (Poaceae) is an ephemeral grass from the monotypic genus Coleanthus Seidl, which grows on wet muddy areas such as fishponds or reservoirs. As a rare species with strict habitat requirements, it is protected at international and national levels. In this study, we sequenced its whole chloroplast genome for the first time using the next-generation sequencing (NGS) technology on the Illumina platform, and performed a comparative and phylogenetic analysis with the related species in Poaceae. The complete chloroplast genome of C. subtilis is 135,915 bp in length, with a quadripartite structure having two 21,529 bp inverted repeat regions (IRs) dividing the entire circular genome into a large single copy region (LSC) of 80,100 bp and a small single copy region (SSC) of 12,757 bp. The overall GC content is 38.3%, while the GC contents in LSC, SSC, and IR regions are 36.3%, 32.4%, and 43.9%, respectively. A total of 129 genes were annotated in the chloroplast genome, including 83 protein-coding genes, 38 tRNA genes, and 8 rRNA genes. The accD gene and the introns of both clpP and rpoC1 genes were missing. In addition, the ycf1, ycf2, ycf15, and ycf68 were pseudogenes. Although the chloroplast genome structure of C. subtilis was found to be conserved and stable in general, 26 SSRs and 13 highly variable loci were detected, these regions have the potential to be developed as important molecular markers for the subfamily Pooideae. Phylogenetic analysis with species in Poaceae indicated that Coleanthus and Phippsia were sister groups, and provided new insights into the relationship between Coleanthus, Zingeria, and Colpodium. This study presents the initial chloroplast genome report of C. subtilis, which provides an essential data reference for further research on its origin.

Comparing the population history of Neotropical annual species: The role of climate change and hybridization between Tithonia tubaeformis and T. rotundifolia (Asteraceae)

Historical climate changes and interspecific gene flow have played an important role in shaping the distribution and genetic diversity of the biota in the Neotropics. In this study, we explored the role of both the Pleistocene climate changes and hybridization on the contemporary geographic structure of two Neotropical weedy species: Tithonia tubaeformis and T. rotundifolia. Distribution shifts under past and current climate conditions were explored through ecological niche modelling (ENM). We then tested these hypotheses using chloroplast microsatellite (cpSSR) data in T. tubaeformis and compared the patterns deduced with those previously reported for T. rotundifolia using the same cpSSR loci. Lastly, we searched for shared haplotypes between species. Both species exhibited significant downwards altitudinal shifts during the last interglacial (LIG) and the Last Glacial Maximum (LGM). For instance, T. rotundifolia showed large suitable habitat areas since the LIG at xeric conditions in western Mesoamerica. Tithonia tubaeformis showed a ~five-fold range contraction during the LIG compared to current climate conditions, followed by a range expansion in the LGM. Despite the large shared refugial areas predicted through ENM, we found a low number of shared haplotypes, suggesting a minor role of hybridization in shaping the geographic structure of these species. Our results provide additional patterns of the population history of the northern Neotropics during the Quaternary, and we suggest that weedy widespread species are a well-suited group for the study of the effects of historic climatic changes on the biota of this region.

Complete Chloroplast Genome Variants Reveal Discrete Long-Distance Dispersal Routes of Rhizophora in the Western Indian Ocean

Historical processes of long-distance migration and ocean-wide expansion feature the global biogeographic pattern of Rhizophora species. Throughout the Indian Ocean, Rhizophora stylosa and Rhizophora mucronata seem to be a young phylogenetic group with an expansion of R. mucronata toward the Western Indian Ocean (WIO) driven by the South Equatorial Current (SEC). Nuclear microsatellites revealed genetic patterns and breaks; however, the estimation of propagule dispersal routes requires maternally inherited cytoplasmic markers. Here, we examine the phylogeography of 21 R. mucronata provenances across a >4,200 km coastal stretch in the WIO using R. stylosa as an outgroup. Full-length chloroplast genome (164,474 bp) and nuclear ribosomal RNA cistron (8,033 bp) sequences were assembled. The boundaries, junction point, sequence orientation, and stretch between LSC/IRb/SSC/IRa/LSC showed no differences with R. stylosa chloroplast genome. A total of 58 mutations in R. mucronata encompassing transitions/transversions, insertions-deletions, and mononucleotide repeats revealed three major haplogroups. Haplonetwork, Bayesian maximum likelihood (ML), and approximate Bayesian computation (ABC) analyses supported discrete historical migration events. An ancient haplogroup A in the Seychelles and eastern Madagascar was as different from other haplogroups as from R. stylosa. A star-like haplonetwork referred as the recent range expansion of haplogroup B from northern Madagascar toward the African mainland coastline, including a single variant spanning >1,800 km across the Mozambique Channel area (MCA). Populations in the south of Delagoa Bight contained haplogroup C and was originated from a unique bottleneck dispersal event. Divergence estimates of pre- and post-Last Glacial Maximum (LGM) illustrated the recent emergence of Rhizophora mangroves in the WIO compared to other oceans. Connectivity patterns could be aligned with the directionality of major ocean currents. Madagascar and the Seychelles each harbored haplogroups A and B, albeit among spatially separated populations, explained from a different migration era. Likewise, the Aldabra Atoll harbored spatially distinct haplotypes. Nuclear ribosomal cistron (8,033 bp) variants corresponded to haplogroups and confirmed admixtures in the Seychelles and Aldabra. These findings shed new light on the origins and dispersal routes of R. mucronata lineages that have shaped their contemporary populations in large regions of the WIO, which may be the important information for defining marine conservation units both at ocean scale and at the level of small islands.

Plastome evolution in the Caesalpinia group (Leguminosae) and its application in phylogenomics and populations genetics

The chloroplast genomes of Caesalpinia group species are structurally conserved, but sequence level variation is useful for both phylogenomic and population genetic analyses. Variation in chloroplast genomes (plastomes) has been an important source of information in plant biology. The Caesalpinia group has been used as a model in studies correlating ecological and genomic variables, yet its intergeneric and infrageneric relationships are not fully solved, despite densely sampled phylogenies including nuclear and plastid loci by Sanger sequencing. Here, we present the de novo assembly and characterization of plastomes from 13 species from the Caesalpinia group belonging to eight genera. A comparative analysis was carried out with 13 other plastomes previously available, totalizing 26 plastomes and representing 15 of the 26 known Caesalpinia group genera. All plastomes showed a conserved quadripartite structure and gene repertoire, except for the loss of four ndh genes in Erythrostemon gilliesii. Thirty polymorphic regions were identified for inter- or intrageneric analyses. The 26 aligned plastomes were used for phylogenetic reconstruction, revealing a well-resolved topology, and dividing the Caesalpinia group into two fully supported clades. Sixteen microsatellite (cpSSR) loci were selected from Cenostigma microphyllum for primer development and at least two were cross-amplified in different Leguminosae subfamilies by in vitro or in silico approaches. Four loci were used to assess the genetic diversity of C. microphyllum in the Brazilian Caatinga. Our results demonstrate the structural conservation of plastomes in the Caesalpinia group, offering insights into its systematics and evolution, and provides new genomic tools for future phylogenetic, population genetics, and phylogeographic studies.

Phylogeny and Comparative Analysis for the Plastid Genomes of Five Tulipa (Liliaceae)

Species of Tulipa (Liliaceae) are of great horticultural importance and are distributed across Europe, North Africa, and Asia. The Tien Shan Mountain is one of the primary diversity centres of Tulipa, but the molecular studies of Tulipa species from this location are lacking. In our study, we assembled four Tulipa plastid genomes from the Tien Shan Mountains, T. altaica, T. iliensis, T. patens, and T. thianschanica, combined with the plastid genome of T. sylvestris to compare against other Liliaceae plastid genomes. We focussed on the species diversity and evolution of their plastid genomes. The five Tulipa plastid genomes proved highly similar in overall size (151,691-152,088 bp), structure, gene order, and content. With comparative analysis, we chose 7 mononucleotide SSRs from the Tulipa species that could be used in further population studies. Phylogenetic analyses based on 24 plastid genomes robustly supported the monophyly of Tulipa and the sister relationship between Tulipa and Amana, Erythronium. T. iliensis, T. thianschanica, and T. altaica were clustered together, and T. patens was clustered with T. sylvestris, with our results clearly demonstrating the relationships between these five Tulipa species. Our results provide a more comprehensive understanding of the phylogenomics and comparative genomics of Tulipa.

Full-length SMRT transcriptome sequencing and microsatellite characterization in Paulownia catalpifolia

Paulownia catalpifolia is an important, fast-growing timber species known for its high density, color and texture. However, few transcriptomic and genetic studies have been conducted in P. catalpifolia. In this study, single-molecule real-time sequencing technology was applied to obtain the full-length transcriptome of P. catalpifolia leaves treated with varying degrees of drought stress. The sequencing data were then used to search for microsatellites, or simple sequence repeats (SSRs). A total of 28.83 Gb data were generated, 25,969 high-quality (HQ) transcripts with an average length of 1624 bp were acquired after removing the redundant reads, and 25,602 HQ transcripts (98.59%) were annotated using public databases. Among the HQ transcripts, 16,722 intact coding sequences, 149 long non-coding RNAs and 179 alternative splicing events were predicted, respectively. A total of 7367 SSR loci were distributed throughout 6293 HQ transcripts, of which 763 complex SSRs and 6604 complete SSRs. The SSR appearance frequency was 28.37%, and the average distribution distance was 5.59 kb. Among the 6604 complete SSR loci, 1-3 nucleotide repeats were dominant, occupying 97.85% of the total SSR loci, of which mono-, di- and tri-nucleotide repeats were 44.68%, 33.86% and 19.31%, respectively. We detected 112 repeat motifs, of which A/T (42.64%), AG/CT (12.22%), GA/TC (9.63%), GAA/TTC (1.57%) and CCA/TGG (1.54%) were most common in mono-, di- and tri-nucleotide repeats, respectively. The length of the repeat SSR motifs was 10-88 bp, and 4997 (75.67%) were ≤ 20 bp. This study provides a novel full-length transcriptome reference for P. catalpifolia and will facilitate the identification of germplasm resources and breeding of new drought-resistant P. catalpifolia varieties.

Complete chloroplast genome of a rare and endangered plant species Osteomeles subrotunda: genomic features and phylogenetic relationships with other Rosaceae plants

Osteomeles subrotunda is a rare and endangered plant species with extremely small populations. In our study, we sequenced the complete chloroplast (CP) genome of O. subrotunda and described its structural organization, and performed comparative genomic analyses with other Rosaceae CP genomes. The plastome of O. subrotunda was 159,902 bp in length with 36.6% GC content and contained a pair of inverted repeats of 26,367 bp which separated a large single-copy region of 87,933 bp and a small single-copy region of 19,235 bp. The CP genome included 130 genes, of which 85 were protein-coding genes, 37 were transfer RNAs, and eight were ribosomal RNAs. Two genes, rps19 and ycf1, which are located at the borders of IRB/SSC and IRB/LSC, were presumed to be pseudogenes. A total of 61 SSRs were detected, of which, 59 loci were mono-nucleotide repeats, and two were di-nucleotide repeats. The phylogenic analysis indicated that the 14 Rosaceae species were divided into three groups, among which O. subrotunda grouped with P. rupicola, E. japonica, P. pashia, C. japonica, S. torminalis, and M. florentina, and it was found to be a sister clade to C. japonica. Our newly sequenced CP genome of O. subrotunda will provide essential data for further studies on population genetics and biodiversity.

Morphological Diversity and Genetic Relationships in Pulque Production Agaves in Tlaxcala, Mexico, by Means of Unsupervised Learning and Gene Sequencing Analysis

Pulque is one of the oldest fermented beverages, with its origins dating back to pre-Hispanic Mexico. Recently, public consumption has increased. However, the majority of Agave plantations for pulque production have disappeared or been abandoned in recent decades. To create strategies for the conservation and production of pulque agaves, it is necessary to first determine their taxonomic identities and to better understand their genetic and morphological diversity. Despite the historical importance of pulque in Mexico, little attention has been placed on the study of Agave plants used for its production. Therefore, we analyzed the morphological diversity of vegetative characters of nine landraces of two Agave species (A. salmiana and A. mapisaga) which are widely cultivated for pulque production in Tlaxcala, Mexico. The analysis of morphological characters showed that the landraces largely clustered based on classic taxonomic relationships. One cluster of landraces associated with Agave mapisaga var. mapisaga and another with A. salmiana subsp. salmiana, but with the exception of A. salmiana subsp. salmiana "Ayoteco", which is more closely related with A. mapisaga var. mapisaga. Additionally, we analyzed the genetic relationships between 14 landraces and wild individuals using molecular markers (trnL and ITS). The identified genetic variants or haplotypes and genetic pools mainly corresponded with the species. In the case of "Ayoteco", incongruence between markers was observed. Low selection intensity, genetic flow events, and the plasticity of morphological traits may explain the high number of landraces without clear differences in their morphological diversity (vegetative characters) or genetic pools. The use of reproductive traits and massive sequencing might be useful for identifying possible morphological and genetic changes in the Agave landraces used for pulque production.

Complete Chloroplast Genomes and Comparative Analyses of L. chinensis, L. anhuiensis, and L. aurea (Amaryllidaceae)

The genus Lycoris (about 20 species) includes important medicinal and ornamental plants. Due to the similar morphological features and insufficient genomic resources, germplasm identification and molecular phylogeny analysis are very limited. Here, we sequenced the complete chloroplast genomes of L. chinensis, L. anhuiensis, and L. aurea; they have very similar morphological traits that make it difficult to identify. The full length of their cp genomes was nearly 158k bp with the same guanine-cytosine content of 37.8%. A total of 137 genes were annotated, including 87 protein-coding genes, 42 tRNAs, and eight rRNAs. A comparative analysis revealed the conservation in sequence size, GC content, and gene content. Some variations were observed in repeat structures, gene expansion on the IR-SC (Inverted Repeat-Single-Copy) boundary regions. Together with the cpSSR (chloroplast simple sequence repeats), these genetic variations are useful to develop molecular markers for germplasm identification. Phylogenetic analysis showed that seven Lycoris species were clustered into a monophyletic group, and closed to Narcissus in Amaryllidaceae. L. chinensis, L. anhuiensis, and L. longituba were clustered together, suggesting that they were very likely to be derived from one species, and had the same ancestor with L. squamigera. Our results provided information on the study of genetic diversity, origins or relatedness of native species, and the identification of cultivars.

Phylogenetic and evolutionary features of the plastome of Tropaeolum pentaphyllum Lam. (Tropaeolaceae)

Complete plastome sequence of Tropaeolum pentaphyllum revealed molecular markers, hotspots of nucleotide polymorphism, RNA editing sites and phylogenetic aspects Tropaeolaceae Juss. ex DC. comprises approximately 95 species across North and South Americas. Tropaeolum pentaphyllum Lam. is an unconventional and endangered species with occurrence in some countries of South America. Although this species presents nutritional, medicinal and ornamental uses, genetic studies involving natural populations or promising genotypes are practically non-existent. Here, we report the nucleotide sequence of T. pentaphyllum plastome. It represents the first complete plastome sequence of the family Tropaeolaceae to be fully sequenced and analyzed in detail. The sequencing data revealed that the T. pentaphyllum plastome is highly similar to the plastomes of other Brassicales. Notwithstanding, our analyses detected some specific features concerning events of IR expansion and structural changes in some genes such as matK, rpoA, and rpoC2. We also detected 251 SSR loci, nine hotspots of nucleotide polymorphism, and two specific RNA editing sites in the plastome of T. pentaphyllum. Moreover, plastid phylogenomic inference indicated a closed relationship between the families Tropaeolaceae and Akaniaceae, which formed a sister group to Moringaceae-Caricaceae. Finally, our data bring new molecular markers and evolutionary features to be applied in the natural population, germplasm collection, and genotype selection aiming conservation, genetic diversity evaluation, and exploitation of this endangered species.

Plastome comparative genomics in maples resolves the infrageneric backbone relationships

Maples (Acer) are among the most diverse and ecologically important tree genera of the north-temperate forests. They include species highly valued as ornamentals and as a source of timber and sugar products. Previous phylogenetic studies employing plastid markers have not provided sufficient resolution, particularly at deeper nodes, leaving the backbone of the maple plastid tree essentially unresolved. We provide the plastid genome sequences of 16 species of maples spanning the sectional diversity of the genus and explore the utility of these sequences as a source of information for genetic and phylogenetic studies in this group. We analyzed the distribution of different types of repeated sequences and the pattern of codon usage, and identified variable regions across the plastome. Maximum likelihood and Bayesian analyses using two partitioning strategies were performed with these and previously published sequences. The plastomes ranged in size from 155,212 to 157,023 bp and had structure and gene content except for Acer palmatum (sect. Palmata), which had longer inverted repeats and an additional copy of the rps19 gene. Two genes, rps2 and rpl22, were found to be truncated at different positions and might be non-functional in several species. Most dispersed repeats, SSRs, and overall variation were detected in the non-coding sequences of the LSC and SSC regions. Fifteen loci, most of which have not been used before in the genus, were identified as the most variable and potentially useful as molecular markers for barcoding and genetic studies. Both ML and Bayesian analyses produced similar results irrespective of the partitioning strategy used. The plastome-based tree largely supported the topology inferred in previous studies using cp markers while providing resolution to the backbone relationships but was highly incongruous with a recently published nuclear tree presenting an opportunity for further research to investigate the causes of discordance, and particularly the role of hybridization in the diversification of the genus. Plastome sequences are valuable tools to resolve deep-level relationships within Acer. The variable loci and SSRs identified in this study will facilitate the development of markers for ecological and evolutionary studies in the genus. This study underscores the potential of plastid genome sequences to improve our understanding of the evolution of maples.

Comparative analyses of 32 complete plastomes of Tef (Eragrostis tef ) accessions from Ethiopia: phylogenetic relationships and mutational hotspots

Eragrostis tef is an important cereal crop in Ethiopia with excellent storage properties, high–quality food, and the unique ability to thrive in extreme environmental conditions. However, the application of advanced molecular tools for breeding and conservation of these species is extremely limited. Therefore, developing chloroplast genome resources and high-resolution molecular markers are valuable to E. tef population and biogeographic studies. In the current study, we assembled and compared the complete plastomes of 32 E. tef accessions. The size of the plastomes ranged from 134,349 to 134,437 bp with similar GC content (∼38.3%). Genomes annotations revealed 112 individual genes, including 77 protein-coding, 31 tRNA, and 4 rRNA genes. Comparison of E. tef plastomes revealed a low degree of intraspecific sequence variations and no structural differentiations. Furthermore, we found 34 polymorphic sites (13 cpSSRs, 12 InDels, and 9 SNPs) that can be used as valuable DNA barcodes. Among them, the majority (88%) of the polymorphic sites were identified in the noncoding genomic regions. Nonsynonymous (ka) and synonymous (ks) substitution analysis showed that all PCGs were under purifying selection (ka/ks <1). The phylogenetic analyses of the whole plastomes and polymorphic region sequences were able to distinguish the accession from the southern population, indicating its potential to be used as a super-barcode. In conclusion, the newly generated plastomes and polymorphic markers developed here could be a useful genomic resource in molecular breeding, population genetics and the biogeographical study of E. tef.

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

The horticulturally important genus Zantedeschia (Araceae) comprises eight species of herbaceous perennials. We sequenced, assembled and analyzed the chloroplast (cp) genomes of four species of Zantedeschia (Z. aethiopica, Z. odorata, Z. elliottiana, and Z. rehmannii) to investigate the structure of the cp genome in the genus. According to our results, the cp genome of Zantedeschia ranges in size from 169,065 bp (Z. aethiopica) to 175,906 bp (Z. elliottiana). We identified a total of 112 unique genes, including 78 protein-coding genes, 30 transfer RNA (tRNA) genes and four ribosomal RNA (rRNA) genes. Comparison of our results with cp genomes from other species in the Araceae suggests that the relatively large sizes of the Zantedeschia cp genomes may result from inverted repeats (IR) region expansion. The sampled Zantedeschia species formed a monophylogenetic clade in our phylogenetic analysis. Furthermore, the long single copy (LSC) and short single copy (SSC) regions in Zantedeschia are more divergent than the IR regions in the same genus, and non-coding regions showed generally higher divergence than coding regions. We identified a total of 410 cpSSR sites from the four Zantedeschia species studied. Genetic diversity analyses based on four polymorphic SSR markers from 134 cultivars of Zantedeschia suggested that high genetic diversity (I = 0.934; Ne = 2.371) is present in the Zantedeschia cultivars. High genetic polymorphism from the cpSSR region suggests that cpSSR could be an effective tool for genetic diversity assessment and identification of Zantedeschia varieties.

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

Mangroves are a group of plant species that occupy the coastal intertidal zone and are major components of this ecologically important ecosystem. Mangroves belong to about twenty diverse families. Here, we sequenced and assembled chloroplast genomes of 14 mangrove species from eight families spanning five rosid orders and one asterid order: Fabales (Pongamia pinnata), Lamiales (Avicennia marina), Malpighiales (Excoecaria agallocha, Bruguiera sexangula, Kandelia obovata, Rhizophora stylosa, and Ceriops tagal), Malvales (Hibiscus tiliaceus, Heritiera littoralis, and Thespesia populnea), Myrtales (Laguncularia racemosa, Sonneratia ovata, and Pemphis acidula), and Sapindales (Xylocarpus moluccensis). These chloroplast genomes range from 149 kb to 168 kb in length. A conserved structure of two inverted repeats (IRa and IRb, ~25.8 kb), one large single-copy region (LSC, ~89.0 kb), and one short single-copy region (SSC, ~18.9 kb) as well as ~130 genes (85 protein-coding, 37 tRNAs, and 8 rRNAs) was observed. We found the lowest divergence in the IR regions among the four regions. We also identified simple sequence repeats (SSRs), which were found to be variable in numbers. Most chloroplast genes are highly conserved, with only four genes under positive selection or relaxed pressure. Combined with publicly available chloroplast genomes, we carried out phylogenetic analysis and confirmed the previously reported phylogeny within rosids, including the positioning of obscure families in Malpighiales. Our study reports 14 mangrove chloroplast genomes and illustrates their genome features and evolution.

Unraveling the genetic complexity of a cultivated breeding population of "yerba mate" (Ilex paraguariensis St. Hil.)

We report the first comprehensive multilocus molecular characterization of cultivated Ilex paraguariensis plants belonging to a breeding program. Using nuclear and homologous chloroplastidic molecular markers, we have genotyped 158 plants from four plantation sites. Analyses of the nuclear data (187 variable dominant loci) allowed detecting high diversity (0.569), the occurrence of four distinct genetic clusters, and a low but significant differentiation among sites. Additionally, 20 chloroplastidic alleles were identified applying five microsatellite polymorphic markers, and a high chloroplastidic diversity was recognized (0.505); two haplogroups were distinguished amongst the 63 haplotypes detected. Our results from both nuclear and plastidic markers indicate that most genetic variation reside within plantations sites (≥ 95%), and that these plantations were established on highly variable materials (either as seeds or plantlets) derived from, at least, 63 maternal lineages. Moreover, our study suggests that the genetic structure within each plantation site was most likely shaped by past admixture favored by farmers´ practices during the establishment of each plantation. Also, subsequent constraints in gene flow and/or a low level of shared polymorphism among plantations could have contributed to current structure.

The Complete Chloroplast Genome of Two Important Annual Clover Species, Trifolium alexandrinum and T. resupinatum: Genome Structure, Comparative Analyses and Phylogenetic Relationships with Relatives in Leguminosae

Trifolium L., which belongs to the IR lacking clade (IRLC), is one of the largest genera in the Leguminosae and contains several economically important fodder species. Here, we present whole chloroplast (cp) genome sequencing and annotation of two important annual grasses, Trifolium alexandrinum (Egyptian clover) and T. resupinatum (Persian clover). Abundant single nucleotide polymorphisms (SNPs) and insertions/deletions (In/Dels) were discovered between those two species. Global alignment of T. alexandrinum and T. resupinatum to a further thirteen Trifolium species revealed a large amount of rearrangement and repetitive events in these fifteen species. As hypothetical cp open reading frame (ORF) and RNA polymerase subunits, ycf1 and rpoC2 in the cp genomes both contain vast repetitive sequences and observed high Pi values (0.7008, 0.3982) between T. alexandrinum and T. resupinatum. Thus they could be considered as the candidate genes for phylogenetic analysis of Trifolium species. In addition, the divergence time of those IR lacking Trifolium species ranged from 84.8505 Mya to 4.7720 Mya. This study will provide insight into the evolution of Trifolium species.

Comparative Analyses of Chloroplast Genomes From 14 Zanthoxylum Species: Identification of Variable DNA Markers and Phylogenetic Relationships Within the Genus

Zanthoxylum L. is an economic crop with a long history of cultivation and domestication and has important economic, ecological, and medicinal value. To solve the classification problems caused by the similar morphological characteristics of Zanthoxylum and establish a credible phylogenetic relationship, we sequenced and annotated six Zanthoxylum chloroplast (cp) genomes (Z. piasezkii, Z. armatum, Z. motuoense, Z. oxyphyllum, Z. multijugum, and Z. calcicola) and combined them with previously published genomes for the Zanthoxylum species. We used bioinformatics methods to analyze the genomic characteristics, contraction, and expansion of inverted repeat (IR) regions; differences in simple sequence repeats (SSRs) and long repeat sequences; species pairwise Ka/Ks ratios; divergence hotspots; and phylogenetic relationships of the 14 Zanthoxylum species. The results revealed that cp genomes of Zanthoxylum range in size from 158,071 to 158,963 bp and contain 87 protein-coding, 37 tRNA, and 8 rRNA genes. Seven mutational hotspots were identified as candidate DNA barcode sequences to distinguish Zanthoxylum species. The phylogenetic analysis strongly supported the genus Fagara as a subgenus of Zanthoxylum and proposed the possibility of a new subgenus in Zanthoxylum. The availability of these genomes will provide valuable information for identifying species, molecular breeding, and evolutionary analysis of Zanthoxylum.

The plastomes of Astrocaryum aculeatum G. Mey. and A. murumuru Mart. show a flip-flop recombination between two short inverted repeats

The plastomes of Astrocaryum murumuru and A. aculeatum revealed a lineage-specific structural feature originated by flip-flop recombination, non-synonymous substitutions in conserved genes and several molecular markers. Astrocaryum murumuru Mart. and A. aculeatum G.Mey. are two palm species of Amazon forest that are economically important as source of food, oil and raw material for several applications. Genetic studies aiming to establish strategies for conservation and domestication of both species are still in the beginning given that the exploitation is mostly by extractive activity. The identification and characterization of molecular markers are essential to assess the genetic diversity of natural populations of both species. Therefore, we sequenced and characterized in detail the plastome of both species. We compared both species and identified 32 polymorphic SSR loci, 150 SNPs, 46 indels and eight hotspots of nucleotide diversity. Additionally, we reported a specific RNA editing site found in the ccsA gene, which is exclusive to A. murumuru. Moreover, the structural analysis in the plastomes of both species revealed a 4.6-kb inversion encompassing a set of genes involved in chlororespiration and plastid translation. This 4.6-kb inversion is a lineage-specific structural feature of the genus Astrocaryum originated by flip-flop recombination between two short inverted repeats. Furthermore, our phylogenetic analysis using whole plastomes of 39 Arecaceae species placed the Astrocaryum species sister to Acrocomia within the tribe Cocoseae. Finally, our data indicated substantial changes in the plastome structure and sequence of both species of the genus Astrocaryum, bringing new molecular markers, several structural and evolving features, which can be applied in several areas such as genetic, evolution, breeding, phylogeny and conservation strategies for both species.

Plastome sequencing of Myripnois dioica and comparison within Asteraceae

Myripnois is a monotypic shrub genus in the daisy family constricted to northern China. Although wild populations of Myripnois dioica are relatively rare, this plant may potentially be cultured as a fine ornamental. In the present study, we sequenced the complete plastome of M. dioica, generating the first plastome sequences of the subfamily Pertyoideae. The plastome of M. dioica has a typical quadripartite circular structure. A large ∼20-kb and a small ∼3-kb inversion were detected in the large single copy (LSC) region and shared by other Asteraceae species. Plastome phylogenomic analyses based on 78 Asteraceae species and three outgroups revealed four groups, corresponding to four Asteraceae subfamilies: Asteroideae, Cichorioideae, Pertyoideae and Carduoideae. Among these four subfamilies, Pertyoideae is sister to Asteroideae + Cichorioideae; Carduoideae is the most basal clade. In addition, we characterized 13 simple sequence repeats (SSRs) that may be useful in future studies on population genetics.

Sequencing and Structural Analysis of the Complete Chloroplast Genome of the Medicinal Plant Lycium chinense Mill

Lycium chinense Mill, an important Chinese herbal medicine, is widely used as a dietary supplement and food. Here the chloroplast (CP) genome of L. chinense was sequenced and analyzed, revealing a size of 155,756 bp and with a 37.8% GC content. The L. chinense CP genome comprises a large single copy region (LSC) of 86,595 bp and a small single copy region (SSC) of 18,209 bp, and two inverted repeat regions (IRa and IRb) of 25,476 bp separated by the single copy regions. The genome encodes 114 genes, 16 of which are duplicated. Most of the 85 protein-coding genes (CDS) had standard ATG start codons, while 3 genes including rps12, psbL and ndhD had abnormal start codons (ACT and ACG). In addition, a strong A/T bias was found in the majority of simple sequence repeats (SSRs) detected in the CP genome. Analysis of the phylogenetic relationships among 16 species revealed that L. chinense is a sister taxon to Lycium barbarum. Overall, the complete sequence and annotation of the L. chinense CP genome provides valuable genetic information to facilitate precise understanding of the taxonomy, species and phylogenetic evolution of the Solanaceae family.

Genetic, evolutionary and phylogenetic aspects of the plastome of annatto (Bixa orellana L.), the Amazonian commercial species of natural dyes

The plastome of B. orellana reveals specific evolutionary features, unique RNA editing sites, molecular markers and the position of Bixaceae within Malvales. Annatto (Bixa orellana L.) is a native species of tropical Americas with center of origin in Brazilian Amazonia. Its seeds accumulate the apocarotenoids, bixin and norbixin, which are only found in high content in this species. The seeds of B. orellana are commercially valued by the food industry because its dyes replace synthetic ones from the market due to potential carcinogenic risks. The increasing consumption of B. orellana seeds for dye extraction makes necessary the increase of productivity, which is possible accessing the genetic basis and searching for elite genotypes. The identification and characterization of molecular markers are essential to analyse the genetic diversity of natural populations and to establish suitable strategies for conservation, domestication, germplasm characterization and genetic breeding. Therefore, we sequenced and characterized in detail the plastome of B. orellana. The plastome of B. orellana is a circular DNA molecule of 159,708 bp with a typical quadripartite structure and 112 unique genes. Additionally, a total of 312 SSR loci were identified in the plastome of B. orellana. Moreover, we predicted in 23 genes a total of 57 RNA-editing sites of which 11 are unique for B. orellana. Furthermore, our plastid phylogenomic analyses, using the plastome sequences available in the plastid database belonging to species of order Malvales, indicate a closed relationship between Bixaceae and Malvaceae, which formed a sister group to Thymelaeaceae. Finally, our study provided useful data to be employed in several genetic and biotechnological approaches in B. orellana and related species of the family Bixaceae.

Comparative chloroplast genomes of Paris Sect. Marmorata: insights into repeat regions and evolutionary implications

Background

Species of Paris Sect. Marmorata are valuable medicinal plants to synthesize steroidal saponins with effective pharmacological therapy. However, the wild resources of the species are threatened by plundering exploitation before the molecular genetics studies uncover the genomes and evolutionary significance. Thus, the availability of complete chloroplast genome sequences of Sect. Marmorata is necessary and crucial to the understanding the plastome evolution of this section and facilitating future population genetics studies. Here, we determined chloroplast genomes of Sect. Marmorata, and conducted the whole chloroplast genome comparison.

Results

This study presented detailed sequences and structural variations of chloroplast genomes of Sect. Marmorata. Over 40 large repeats and approximately 130 simple sequence repeats as well as a group of genomic hotspots were detected. Inverted repeat contraction of this section was inferred via comparing the chloroplast genomes with the one of P. verticillata. Additionally, almost all the plastid protein coding genes were found to prefer ending with A/U. Mutation bias and selection pressure predominately shaped the codon bias of most genes. And most of the genes underwent purifying selection, whereas photosynthetic genes experienced a relatively relaxed purifying selection.

Conclusions

Repeat sequences and hotspot regions can be scanned to detect the intraspecific and interspecific variability, and selected to infer the phylogenetic relationships of Sect. Marmorata and other species in subgenus Daiswa. Mutation and natural selection were the main forces to drive the codon bias pattern of most plastid protein coding genes. Therefore, this study enhances the understanding about evolution of Sect. Marmorata from the chloroplast genome, and provide genomic insights into genetic analyses of Sect. Marmorata.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5281-x) contains supplementary material, which is available to authorized users.

Complete chloroplast genome sequence of Betula platyphylla: gene organization, RNA editing, and comparative and phylogenetic analyses

BACKGROUND

Betula platyphylla is a common tree species in northern China that has high economic and medicinal value. Our laboratory has been devoted to genome research on B. platyphylla for approximately 10 years. As primary organelle genomes, the complete genome sequences of chloroplasts are important to study the divergence of species, RNA editing and phylogeny. In this study, we sequenced and analyzed the complete chloroplast (cp) genome sequence of B. platyphylla.

RESULTS

The complete cp genome of B. platyphylla was 160,518 bp in length, which included a pair of inverted repeats (IRs) of 26,056 bp that separated a large single copy (LSC) region of 89,397 bp and a small single copy (SSC) region of 19,009 bp. The annotation contained a total of 129 genes, including 84 protein-coding genes, 37 tRNA genes and 8 rRNA genes. There were 3 genes using alternative initiation codons. Comparative genomics showed that the sequence of the Fagales species cp genome was relatively conserved, but there were still some high variation regions that could be used as molecular markers. The IR expansion event of B. platyphylla resulted in larger cp genomes and rps19 pseudogene formation. The simple sequence repeat (SSR) analysis showed that there were 105 SSRs in the cp genome of B. platyphylla. RNA editing sites recognition indicated that at least 80 RNA editing events occurred in the cp genome. Most of the substitutions were C to U, while a small proportion of them were not. In particular, three editing loci on the rRNA were converted to more than two other bases that had never been reported. For synonymous conversion, most of them increased the relative synonymous codon usage (RSCU) value of the codons. The phylogenetic analysis suggested that B. platyphylla had a closer evolutionary relationship with B. pendula than B. nana.

CONCLUSIONS

In this study, we not only obtained and annotated the complete cp genome sequence of B. platyphylla, but we also identified new RNA editing sites and predicted the phylogenetic relationships among Fagales species. These findings will facilitate genomic, genetic engineering and phylogenetic studies of this important species.

Genetic patterns in Neotropical Magnolias (Magnoliaceae) using de novo developed microsatellite markers

Conserving tree populations safeguards forests since they represent key elements of the ecosystem. The genetic characteristics underlying the evolutionary success of the tree growth form: high genetic diversity, extensive gene flow and strong species integrity, contribute to their survival in terms of adaptability. However, different biological and landscape contexts challenge these characteristics. This study employs 63 de novo developed microsatellite or SSR (Single Sequence Repeat) markers in different datasets of nine Neotropical Magnolia species. The genetic patterns of these protogynous, insect-pollinated tree species occurring in fragmented, highly-disturbed landscapes were investigated. Datasets containing a total of 340 individuals were tested for their genetic structure and degree of inbreeding. Analyses for genetic structure depicted structuring between species, i.e. strong species integrity. Within the species, all but one population pair were considered moderate to highly differentiated, i.e. no indication of extensive gene flow between populations. No overall correlation was observed between genetic and geographic distance of the pairwise species’ populations. In contrast to the pronounced genetic structure, there was no evidence of inbreeding within the populations, suggesting mechanisms favouring cross pollination and/or selection for more genetically diverse, heterozygous offspring. In conclusion, the data illustrate that the Neotropical Magnolias in the context of a fragmented landscape still have ample gene flow within populations, yet little gene flow between populations.

Population Structure and Genetic Diversity Within the Endangered Species Pityopsis ruthii (Asteraceae)

Pityopsis ruthii (Ruth's golden aster) is a federally endangered herbaceous perennial endemic to the Hiwassee and Ocoee Rivers in southeastern Tennessee, United States. Comprehensive genetic studies providing novel information to conservationists for preservation of the species are lacking. Genetic variation and gene flow were evaluated for 814 individuals from 33 discrete locations using polymorphic microsatellites: seven chloroplast and twelve nuclear. A total of 198 alleles were detected with the nuclear loci and 79 alleles with the chloroplast loci. Gene flow was estimated, with the Hiwassee River (N_m = 2.16; F_ST = 0.15) showing higher levels of gene flow and lower levels of population differentiation than the Ocoee River (N_m = 1.28; F_ST = 0.19). Population structure was examined using Bayesian cluster analyses. Nuclear and chloroplast analyses were incongruent. From the chloroplast microsatellites, three clusters were identified; all were present in sampling sites at both rivers, indicating a lack of allele fixation along rivers. Nuclear markers revealed two clusters and separated by river. When the Hiwassee River locations were analyzed, four clusters were identified for both the chloroplast and nuclear microsatellites, though the individuals clustered differently. Analysis of the Ocoee River revealed two clusters for the chloroplast microsatellites and three for the nuclear microsatellites. We recommend P. ruthii be managed as four populations for the Hiwassee River and three populations for the Ocoee River. Our results provide critical genetic information for P. ruthii that can be used for species management decisions to drive future population augmentation/reintroduction and ex situ conservation efforts.

The complete plastome of macaw palm [Acrocomia aculeata (Jacq.) Lodd. ex Mart.] and extensive molecular analyses of the evolution of plastid genes in Arecaceae

The plastome of macaw palm was sequenced allowing analyses of evolution and molecular markers. Additionally, we demonstrated that more than half of plastid protein-coding genes in Arecaceae underwent positive selection. Macaw palm is a native species from tropical and subtropical Americas. It shows high production of oil per hectare reaching up to 70% of oil content in fruits and an interesting plasticity to grow in different ecosystems. Its domestication and breeding are still in the beginning, which makes the development of molecular markers essential to assess natural populations and germplasm collections. Therefore, we sequenced and characterized in detail the plastome of macaw palm. A total of 221 SSR loci were identified in the plastome of macaw palm. Additionally, eight polymorphism hotspots were characterized at level of subfamily and tribe. Moreover, several events of gain and loss of RNA editing sites were found within the subfamily Arecoideae. Aiming to uncover evolutionary events in Arecaceae, we also analyzed extensively the evolution of plastid genes. The analyses show that highly divergent genes seem to evolve in a species-specific manner, suggesting that gene degeneration events may be occurring within Arecaceae at the level of genus or species. Unexpectedly, we found that more than half of plastid protein-coding genes are under positive selection, including genes for photosynthesis, gene expression machinery and other essential plastid functions. Furthermore, we performed a phylogenomic analysis using whole plastomes of 40 taxa, representing all subfamilies of Arecaceae, which placed the macaw palm within the tribe Cocoseae. Finally, the data showed here are important for genetic studies in macaw palm and provide new insights into the evolution of plastid genes and environmental adaptation in Arecaceae.

Chloroplast genomes of Lilium lancifolium, L. amabile, L. callosum, and L. philadelphicum: Molecular characterization and their use in phylogenetic analysis in the genus Lilium and other allied genera in the order Liliales

Chloroplast (cp) genomes of Lilium amabile, L. callosum, L. lancifolium, and L. philadelphicum were fully sequenced. Using these four novel cp genome sequences and five other previously sequenced cp genomes, features of the cp genomes were characterized in detail among species in the genus Lilium and other related genera in the order Liliales. The lengths and nucleotide composition showed little variation. No structural variation was found among the cp genomes in Liliales. Gene contents were conserved among four newly sequenced cp genome in Lilium species, the only differences being in two pseudogenes. We identified 112 genes in 13 functional categories, 18 of which carried introns that were conserved among the species in Liliales. There were 16-21 SSR loci (>12 bp, >3 repeats) in the cp genomes in Lilium and the genomic locations of these loci were highly variable among the species. Average mutations were 15 SNPs per 1kb and 5 indels per 1kb, respectively, in the cp genomes of the newly sequenced four Lilium species. Phylogenetic classifications revealed some discrepancies between trees based on the cp genomes and previous classifications based on the morphology and geographic distributions.

Chloroplast Genome Sequence of Clusterbean (Cyamopsis tetragonoloba L.): Genome Structure and Comparative Analysis

Clusterbean (Cyamopsis tetragonoloba L.), also known as guar, belongs to the family Leguminosae, and is an annual herbaceous legume. Guar is the main source of galactomannan for gas mining industries. In the present study, the draft chloroplast genome of clusterbean was generated and compared to some of the previously reported legume chloroplast genomes. The chloroplast genome of clusterbean is 152,530 bp in length, with a quadripartite structure consisting of large single copy (LSC) and small single copy (SSC) of 83,025 bp and 17,879 bp in size, respectively, and a pair of inverted repeats (IRs) of 25,790 bp in size. The chloroplast genome contains 114 unique genes, which includes 78 protein coding genes, 30 tRNAs, 4 rRNAs genes, and 2 pseudogenes. It also harbors a 50 kb inversion, typical of the Leguminosae family. The IR region of the clusterbean chloroplast genome has undergone an expansion, and hence, the whole rps19 gene is included in the IR, as compared to other legume plastid genomes. A total of 220 simple sequence repeats (SSRs) were detected in the clusterbean plastid genome. The analysis of the clusterbean plastid genome will provide useful insights for evolutionary, molecular and genetic engineering studies.

Appropriate homoplasy metrics in linked SSRs to predict an underestimation of demographic expansion times

BACKGROUND

Homoplasy affects demographic inference estimates. This effect has been recognized and corrective methods have been developed. However, no studies so far have defined what homoplasy metrics best describe the effects on demographic inference, or have attempted to estimate such metrics in real data. Here we study how homoplasy in chloroplast microsatellites (cpSSR) affects inference of population expansion time. cpSSRs are popular markers for inferring historical demography in plants due to their high mutation rate and limited recombination.

RESULTS

In cpSSRs, homoplasy is usually quantified as the probability that two markers or haplotypes that are identical by state are not identical by descent (Homoplasy index, P). Here we propose a new measure of multi-locus homoplasy in linked SSR called Distance Homoplasy (DH), which measures the proportion of pairwise differences not observed due to homoplasy, and we compare it to P and its per cpSSR locus average, which we call Mean Size Homoplasy (MSH). We use simulations and analytical derivations to show that, out of the three homoplasy metrics analyzed, MSH and DH are more correlated to changes in the population expansion time and to the underestimation of that demographic parameter using cpSSR. We perform simulations to show that Approximate Bayesian Computation (ABC) can be used to obtain reasonable estimates of MSH and DH. Finally, we use ABC to estimate the expansion time, MSH and DH from a chloroplast SSR dataset in Pinus caribaea. To our knowledge, this is the first time that homoplasy has been estimated in population genetic data.

CONCLUSIONS

We show that MSH and DH should be used to quantify how homoplasy affects estimates of population expansion time. We also demonstrate how ABC provides a methodology to estimate homoplasy in population genetic data.

Strategies for complete plastid genome sequencing

Plastid sequencing is an essential tool in the study of plant evolution. This high-copy organelle is one of the most technically accessible regions of the genome, and its sequence conservation makes it a valuable region for comparative genome evolution, phylogenetic analysis and population studies. Here, we discuss recent innovations and approaches for de novo plastid assembly that harness genomic tools. We focus on technical developments including low-cost sequence library preparation approaches for genome skimming, enrichment via hybrid baits and methylation-sensitive capture, sequence platforms with higher read outputs and longer read lengths, and automated tools for assembly. These developments allow for a much more streamlined assembly than via conventional short-range PCR. Although newer methods make complete plastid sequencing possible for any land plant or green alga, there are still challenges for producing finished plastomes particularly from herbarium material or from structurally divergent plastids such as those of parasitic plants.

Development and characterization of 24 chloroplast microsatellite markers for two species of Eranthis (Ranunculaceae)

Chloroplast microsatellites for two Korean endemic species, Eranthis byunsanensis and E. pungdoensis (Ranunculaceae), were isolated to address the questions of their distributional patterns and evolutionary relationships, using next-generation sequencing. Twenty-four polymorphic chloroplast microsatellite markers for these two species were developed, and then characterized in 65 individuals (55 individuals of E. byunsanensis and 10 individuals of E. pungdoensis). The number of alleles per locus ranged from 2 to 9; the average number of alleles across all the loci scored 4.792. The unbiased diversity per locus ranged from 0.089 to 0.880; the unbiased diversity averaged over all the loci was 0.646. The developed markers were successfully amplified for three congeneric species, E. stellata, E. pinnatifida, and E. longistipitata. The markers developed in this study can provide a valuable and important tool for understanding genetic variations, population structures, evolutionary histories and phylogeography of E. byunsanensis, E. pungdoensis, and related species.

Is biomass a reliable estimate of plant fitness?

The measurement of fitness is critical to biological research. Although the determination of fitness for some organisms may be relatively straightforward under controlled conditions, it is often a difficult or nearly impossible task in nature. Plants are no exception. The potential for long-distance pollen dispersal, likelihood of multiple reproductive events per inflorescence, varying degrees of reproductive growth in perennials, and asexual reproduction all confound accurate fitness measurements. For these reasons, biomass is frequently used as a proxy for plant fitness. However, the suitability of indirect fitness measurements such as plant size is rarely evaluated. This review outlines the important associations between plant performance, fecundity, and fitness. We make a case for the reliability of biomass as an estimate of fitness when comparing conspecifics of the same age class. We reviewed 170 studies on plant fitness and discuss the metrics commonly employed for fitness estimations. We find that biomass or growth rate are frequently used and often positively associated with fecundity, which in turn suggests greater overall fitness. Our results support the utility of biomass as an appropriate surrogate for fitness under many circumstances, and suggest that additional fitness measures should be reported along with biomass or growth rate whenever possible.

Complete Chloroplast Genomes of Erianthus arundinaceus and Miscanthus sinensis: Comparative Genomics and Evolution of the Saccharum Complex

The genera Erianthus and Miscanthus, both members of the Saccharum complex, are of interest as potential resources for sugarcane improvement and as bioenergy crops. Recent studies have mainly focused on the conservation and use of wild accessions of these genera as breeding materials. However, the sequence data are limited, which hampers the studies of phylogenetic relationships, population structure, and evolution of these grasses. Here, we determined the complete chloroplast genome sequences of Erianthus arundinaceus and Miscanthus sinensis by using 454 GS FLX pyrosequencing and Sanger sequencing. Alignment of the E. arundinaceus and M. sinensis chloroplast genome sequences with the known sequence of Saccharum officinarum demonstrated a high degree of conservation in gene content and order. Using the data sets of 76 chloroplast protein-coding genes, we performed phylogenetic analysis in 40 taxa including E. arundinaceus and M. sinensis. Our results show that S. officinarum is more closely related to M. sinensis than to E. arundinaceus. We estimated that E. arundinaceus diverged from the subtribe Sorghinae before the divergence of Sorghum bicolor and the common ancestor of S. officinarum and M. sinensis. This is the first report of the phylogenetic and evolutionary relationships inferred from maternally inherited variation in the Saccharum complex. Our study provides an important framework for understanding the phylogenetic relatedness of the economically important genera Erianthus, Miscanthus, and Saccharum.

Characterization and Phylogenetic Analysis of Ancient Italian Landraces of Pear

Pear is one of the oldest fruit tree crops and the third most important temperate fruit species. Its domestication took place independently in the Far East (China) and in the Caucasus region. While the origin of Eastern Asian cultivars is clear, that of European cultivars is still in doubt. Italy has a wealth of local varieties and genetic resources safeguarded by several public and private collections to face the erosion caused by the introduction of improved varieties in specialized orchards. The objectives of the present study were: (i) to characterize the existing germplasm through nuclear (SSR) and (ii) to clarify the genetic divergence between local and cultivated populations through chloroplast DNA (cpDNA) markers in order to provide insights into phylogenetic relationships of Pyrus spp. For this reason, 95 entries from five different germplasm collections, including nine European, Mediterranean and Eastern Asian species, were analyzed, and the intergenic accD-psaI sequences were compared to the worldwide distributed dataset encompassing a total of 298 sequences from 26 different Pyrus species. The nine nuclear SSRs were able to identify a total of 179 alleles, with a loci polymorphism P = 0.89. Most of the variation (97%) was found within groups. Five accessions from different sources were confirmed to be the same. Eight out of 20 accessions of unknown origin were identified, and six synonyms were detected. Locus NH030a was found to be monomorphic in all the cultivated accessions and in reference species interfertile with P. communis, leading to hypothesize selection pressures for adaptation to cultivation. The cpDNA sequences of the 95 accessions were represented by 14 haplotypes, six of which (derived from P. communis, P. cossonii and P. ussuriensis) are recorded here for the first time and may suggest the ancient origin of some local varieties. The network analysis of the 298 cpDNA sequences allowed two different haplogroups, Eastern and Western Eurasia, to be defined, supporting recent views of a clear division between Occidental and Oriental species. By combining the results from nuclear and uniparental markers, it was possible to better define many unknown accessions.

Genetic diversity and structure of wild populations of Carica papaya in Northern Mesoamerica inferred by nuclear microsatellites and chloroplast markers

BACKGROUND AND AIMS

Few studies have evaluated the genetic structure and evolutionary history of wild varieties of important crop species. The wild papaya (Carica papaya) is a key element of early successional tropical and sub-tropical forests in Mexico, and constitutes the genetic reservoir for evolutionary potential of the species. In this study we aimed to determine how diverse and structured is the genetic variability of wild populations of C. papaya in Northern Mesoamerica. Moreover, we assessed if genetic structure and evolutionary history coincide with hypothetized (1) pre-Pleistocene events (Isthmus of Tehuantepec sinking), (2) Pleistocene refugia or (3) recent patterns.

METHODS

We used six nuclear and two chloroplast (cp) DNA markers to assess the genetic diversity and phylogeographical structure of 19 wild populations of C. papaya in its natural distribution in Northern Mesoamerica.

KEY RESULTS

We found high genetic diversity (Ho = 0·681 for nuclear markers, and h = 0·701 for cpDNA markers) and gene flow between populations of C. papaya (migration r up to 420 km). A lack of phylogeographical structure was found with the cpDNA markers (NST < GST), whereas a recent population structure was inferred with the nuclear markers. Evidence indicates that pre-Pleistocene events or refugia did not play an important role in the genetic structuring of wild papaya.

CONCLUSIONS

Because of its life history characteristics and lack of an ancient phylogeographical structure found with the cpDNA markers, we suggest that C. papaya was dispersed throughout the lowland rain forests of Mexico (along the coastal plains and foothills of Sierras). This scenario supports the hypothesis that tropical forests in Northern Mesoamerica did not experience important climate fluctuations during the Pleistocene, and that the life history of C. papaya could have promoted long-distance dispersal and rapid colonization of lowland rainforests. Moreover, the results obtained with the nuclear markers suggest recent human disturbances. The fragmentation of tropical habitats in Northern Mesoamerica appears to be the main driver of genetic structuring, and the major threat to the dispersion and survival of the species in the wild.

Complete Chloroplast Genome Sequence of Corroborates Structural Heterogeneity of Inverted Repeats in Wild Progenitors of Cultivated Bananas and Plantains

Complete genome sequencing of cytoplasmically inherited chloroplast DNA provides novel insights into the origins of clonally propagated crops such as banana and plantain ( spp.). This study describes the structural organization of the chloroplast genome of Colla and its phylogenetic relationship with other wild progenitors of the domesticated banana cultivars. The chloroplast genome was sequenced using Illumina HiSeq 2000 platform, followed by a combination of de novo short-read assembly and reference-guided mapping of contigs to generate complete plastome sequence. The chloroplast genome is 169,503 bp in length, exhibits a typical quadripartite structural organization with a large single-copy (LSC; 87,828 bp) region and a small single-copy (SSC; 11,547 bp) region interspersed between inverted repeat (IRa/b; 35,064 bp) regions. Overall, its gene content, size, and gene order were identical to that of Colla with extensive expansion of the inverted repeat-small single-copy (IR-SSC) junctions. Comparative analyses revealed the conserved IRa-SSC expansion in three wild species and members of the order Zingiberales. In contrast, IRb-SSC expansion was conspicuously absent in the sister taxon Nee and related species of Zingiberales. Interestingly, phylogenomic assessment based on whole-plastome and protein-coding gene sets have provided robust support for the association of and as a sister group, despite the variation in IRb-SSC expansion. Although the current study substantiates the infrageneric IRb-SSC fluctuations in Musaceae, extensive taxon sampling is necessary to confirm whether the accessions of section have undergone independent IRb-SSC expansion relative to section .

Comparative assessment of genetic diversity in cytoplasmic and nuclear genome of upland cotton

The importance of the cytoplasmic genome for many economically important traits is well documented in several crop species, including cotton. There is no report on application of cotton chloroplast specific SSR markers as a diagnostic tool to study genetic diversity among improved Upland cotton lines. The complete plastome sequence information in GenBank provided us an opportunity to report on 17 chloroplast specific SSR markers using a cost-effective data mining strategy. Here we report the comparative analysis of genetic diversity among a set of 42 improved Upland cotton lines using SSR markers specific to chloroplast and nuclear genome, respectively. Our results revealed that low to moderate level of genetic diversity existed in both nuclear and cytoplasm genome among this set of cotton lines. However, the specific estimation suggested that genetic diversity is lower in cytoplasmic genome compared to the nuclear genome among this set of Upland cotton lines. In summary, this research is important from several perspectives. We detected a set of cytoplasm genome specific SSR primer pairs by using a cost-effective data mining strategy. We reported for the first time the genetic diversity in the cytoplasmic genome within a set of improved Upland cotton accessions. Results revealed that the genetic diversity in cytoplasmic genome is narrow, compared to the nuclear genome within this set of Upland cotton accessions. Our results suggested that most of these polymorphic chloroplast SSRs would be a valuable complementary tool in addition to the nuclear SSR in the study of evolution, gene flow and genetic diversity in Upland cotton.

Microsatellites in Pursuit of Microbial Genome Evolution

Microsatellites or short sequence repeats are widespread genetic markers which are hypermutable 1-6 bp long short nucleotide motifs. Significantly, their applications in genetics are extensive due to their ceaseless mutational degree, widespread length variations and hypermutability skills. These features make them useful in determining the driving forces of evolution by using powerful molecular techniques. Consequently, revealing important questions, for example, what is the significance of these abundant sequences in DNA, what are their roles in genomic evolution? The answers of these important questions are hidden in the ways these short motifs contributed in altering the microbial genomes since the origin of life. Even though their size ranges from 1 -to- 6 bases, these repeats are becoming one of the most popular genetic probes in determining their associations and phylogenetic relationships in closely related genomes. Currently, they have been widely used in molecular genetics, biotechnology and evolutionary biology. However, due to limited knowledge; there is a significant gap in research and lack of information concerning hypermutational mechanisms. These mechanisms play a key role in microsatellite loci point mutations and phase variations. This review will extend the understandings of impacts and contributions of microsatellite in genomic evolution and their universal applications in microbiology.