Comparative Genomics and Phylogenetic Analysis Revealed the Chloroplast Genome Variation and Interspecific Relationships of Corylus (Betulaceae) Species

Real-time expression and in silico characterization of pea genes involved in salt and water-deficit stress

BACKGROUND

To tolerate salt and water-deficit stress, the plant adapts to the adverse environment by regulating its metabolism and expressing certain stress-induced metabolic pathways. This research analyzed the relative expression of four pea genes (P5CR, PAL1, SOD, and POX) in three pea varieties (Climax, Green grass, and Meteor) under different levels of salt and water-deficit stress.

METHODS AND RESULTS

The experiments on salt stress and water-deficit stress were carried out within greenhouse settings under controlled environment. The saturation percentage was employed to create artificial salinity conditions: Control without NaCl treatment, Treatment 1: 50 mM NaCl treatment, Treatment 2: 75 mM NaCl treatment, and Treatment 3: 100 mM NaCl treatment. Field capacity (FC) was used for the development of artificial water-deficit treatments in the pots, i.e., Treatment 1 (Control; water application 100% of FC), Treatment 2 (water application 75% of FC), and Treatment 3 (water application 50% of FC). Pea genes involved in biosynthetic pathways of proline, flavonoids, and enzymatic antioxidant enzymes including P5CR, PAL1, SOD, and POX were selected based on literature. Quantitative real-time PCR using cDNA as a template was used to analyze the gene expression. Pea genes were analyzed for phylogenetic analysis in closely related crops having similarity percent identity 80% and above. In silico characterization of selected proteins including the family classification was done by the NCBI CDD and INTERPRO online servers. Results from RT-qPCR analysis showed increased expression of P5CR, PAL1, and POX genes, while SOD gene expression decreased under both stresses. Climax exhibited superior stress tolerance with elevated expression of P5CR and PAL1, while Meteor showed better tolerance through increased POX expression. Phylogenetic analysis revealed common ancestry with other species like chickpea, red clover, mung bean, and barrel clover, suggesting the cross relationship among these plant species. Conserved domain analysis of respective proteins revealed that these proteins contain PLNO 2688, PLN02457, Cu-Zn Superoxide dismutase, and secretory peroxidase conserved domains. Furthermore, protein family classification indicated that the oxidation-reduction process is the most common chemical process involved in these stresses given to pea plant which validates the relationship of these proteins.

CONCLUSIONS

Salt and water-deficit stresses trigger distinct metabolic pathways, leading to the up-regulation of specific genes and the synthesis of corresponding proteins. These findings further emphasize the conservation of stress-tolerance-related genes and proteins across various plant species. This knowledge enhances our understanding of plant adaptation to stress and offers opportunities for developing strategies to improve stress resilience in crops, thereby addressing global food security challenges.

Complete plastid genome structure of 13 Asian Justicia (Acanthaceae) species: comparative genomics and phylogenetic analyses

BACKGROUND

Justicia L. is the largest genus in Acanthaceae Juss. and widely distributed in tropical and subtropical regions of the world. Previous phylogenetic studies have proposed a general phylogenetic framework for Justicia based on several molecular markers. However, their studies were mainly focused on resolution of phylogenetic issues of Justicia in Africa, Australia and South America due to limited sampling from Asia. Additionally, although Justicia plants are of high medical and ornamental values, little research on its genetics was reported. Therefore, to improve the understanding of its genomic structure and relationships among Asian Justicia plants, we sequenced complete chloroplast (cp.) genomes of 12 Asian plants and combined with the previously published cp. genome of Justicia leptostachya Hemsl. for further comparative genomics and phylogenetic analyses.

RESULTS

All the cp. genomes exhibit a typical quadripartite structure without genomic rearrangement and gene loss. Their sizes range from 148,374 to 151,739 bp, including a large single copy (LSC, 81,434-83,676 bp), a small single copy (SSC, 16,833-17,507 bp) and two inverted repeats (IR, 24,947-25,549 bp). GC contents range from 38.1 to 38.4%. All the plastomes contain 114 genes, including 80 protein-coding genes, 30 tRNAs and 4 rRNAs. IR variation and repetitive sequences analyses both indicated that Justicia grossa C. B. Clarke is different from other Justicia species because its lengths of ndhF and ycf1 in IRs are shorter than others and it is richest in SSRs and dispersed repeats. The ycf1 gene was identified as the candidate DNA barcode for the genus Justicia. Our phylogenetic results showed that Justicia is a polyphyletic group, which is consistent with previous studies. Among them, J. grossa belongs to subtribe Tetramerinae of tribe Justicieae while the other Justicia members belong to subtribe Justiciinae. Therefore, based on morphological and molecular evidence, J. grossa should be undoubtedly recognized as a new genus. Interestingly, the evolutionary history of Justicia was discovered to be congruent with the morphology evolution.

CONCLUSION

Our study not only elucidates basic features of Justicia whole plastomes, but also sheds light on interspecific relationships of Asian Justicia plants for the first time.

Extensive characterization of 28 complete chloroplast genomes of Hydrangea species: A perspective view of their organization and phylogenetic and evolutionary relationships

The tribe Hydrangeeae displays a unique, distinctive disjunct distribution encompassing East Asia, North America and Hawaii. Despite its complex trait variations and polyphyletic nature, comprehensive phylogenomic and biogeographical studies on this tribe have been lacking. To address this gap, we sequenced and characterized 28 plastomes of Hydrangeeae. Our study highlights the highly conserved nature of Hydrangeaceae chloroplast (cp) genomes in terms of gene content and arrangement. Notably, synapomorphic characteristics of tandem repeats in the conserved domain of accD were observed in the Macrophyllae, Chinenses, and Dichroa sections within the Hydrangeeae tribe. Additionally, we found lower expression of accD in these sections using structure prediction and quantitative real-time PCR analysis. Phylogenomic analyses revealed the subdivision of the Hydrangeeae tribe into two clades with robust support values. Consistent with polyphyletic relationships, sect. Broussaisia was identified as the basal group in the tribe Hydrangeeae. Our study also provides insights into the phylogenetic relationships of Hydrangea petiolaris in the Jeju and Ulleung Island populations, suggesting the need for further studies with more samples and molecular data. Divergence time estimation and biogeographical analyses suggested that the common ancestors of the tribe Hydrangeeae likely originated from North America and East Asia during the Paleocene period via the Bering Land Bridge, potentially facilitating migration within the tribe between these regions. In conclusion, this study enhances our understanding of the evolutionary history and biogeography of the tribe Hydrangeeae, shedding light on the dispersal patterns and origins of this intriguing plant group with its unique disjunct distribution.

Dig up tall fescue plastid genomes for the identification of morphotype-specific DNA variants

BACKGROUND

Tall fescue (Festuca arundinacea Schreb.) is an important cool-season perennial grass species. Hexaploid tall fescue has three distinct morphotypes used either as forage or turf purposes. Its chloroplast genome is conserved due to it being maternally inherited to the next generation progenies. To identify morphotype-specific DNA markers and the genetic variations, plastid genomes of all three tall fescue morphotypes, i.e., Continental cv. Texoma MaxQ II, Rhizomatous cv. Torpedo, and Mediterranean cv. Resolute, have been sequenced using Illumina MiSeq sequencing platform.

RESULTS

The plastid genomes of Continental-, Rhizomatous-, and Mediterranean tall fescue were assembled into circular master molecules of 135,283 bp, 135,336 bp, and 135,324 bp, respectively. The tall fescue plastid genome of all morphotypes contained 77 protein-coding, 20 tRNAs, four rRNAs, two pseudo protein-coding, and three hypothetical protein-coding genes. We identified 630 SNPs and 124 InDels between Continental and Mediterranean, 62 SNPs and 20 InDels between Continental and Rhizomatous, and 635 SNPs and 123 InDels between Rhizomatous and Mediterranean tall fescue. Only four InDels in four genes (ccsA, rps18, accD, and ndhH-p) were identified, which discriminated Continental and Rhizomatous plastid genomes from the Mediterranean plastid genome. Here, we identified and reported eight InDel markers (NRITCHL18, NRITCHL35, NRITCHL43, NRITCHL65, NRITCHL72, NRITCHL101, NRITCHL104, and NRITCHL110) from the intergenic regions that can successfully discriminate tall fescue morphotypes. Divergence time estimation revealed that Mediterranean tall fescue evolved approximately 7.09 Mya, whereas the divergence between Continental- and Rhizomatous tall fescue occurred about 0.6 Mya.

CONCLUSIONS

To our knowledge, this is the first report of the assembled plastid genomes of Rhizomatous and Mediterranean tall fescue. Our results will help to identify tall fescue morphotypes at the time of pre-breeding and will contribute to the development of lawn and forage types of commercial varieties.

Comparative and phylogenetic analysis of complete chloroplast genomes from seven Neocinnamomum taxa (Lauraceae)

The genus Neocinnamomum is considered to be one of the most enigmatic groups in Lauraceae, mainly distributed in tropical and subtropical regions of Southeast Asia. The genus contains valuable oilseed and medicinal tree species. However, there are few studies on the genus Neocinnamomum at present, and its interspecific relationship is still unclear. In order to explore the genetic structure and evolutionary characteristics of the Neocinnamomum chloroplast genome and to resolve the species relationships within the genus, comparative genomic and phylogenetic analyses were performed on the whole chloroplast genome sequences of 51 samples representing seven Neocinnamomum taxa. The whole Neocinnamomum chloroplast genome size ranged from 150,753-150,956 bp, with a GC content of 38.8%-38.9%. A total of 128 genes were annotated within the Neocinnamomum chloroplast genome, including 84 protein coding genes, 8 rRNA genes, and 36 tRNA genes. Between 71-82 SSRs were detected, among which A/T base repeats were the most common. The chloroplast genome contained a total of 31 preferred codons. Three highly variable regions, trnN-GUU-ndhF, petA-psbJ, and ccsA-ndhD, were identified with Pi values > 0.004. Based on the whole chloroplast genome phylogenetic tree, the phylogenetic relationships among the seven Neocinnamomum taxa were determined. N. delavayi and N. fargesii were the most closely related species, and N. lecomtei was identified as the most basal taxon. In this study, the characteristics and sequence variation of the chloroplast genomes of seven Neocinnamomum taxa were revealed, and the genetic relationship among the species was clarified. The results of this study will provide a reference for subsequent molecular marker development and phylogenetic research of Neocinnamomum.

The complete chloroplast genome of Ulmus mianzhuensis with insights into structural variations, adaptive evolution, and phylogenetic relationships of Ulmus (Ulmaceae)

BACKGROUND

Ulmus mianzhuensis is an endemic tree species in China with high ornamental and economic value. Currently, little is known regarding its genomic architecture, phylogenetic position, or adaptive evolution. Here, we sequenced the complete chloroplast genome (cp genome) of U. mianzhuensis and further compared the variations in gene organization and structure within Ulmus species to define their genomic evolution, then reconstructed the phylogenomic relationship of 31 related Ulmus species to explore the systematic position of U. mianzhuensis and the utility of cp genome for resolving phylogenetics among Ulmus species.

RESULTS

Our results revealed that all the Ulmus species exhibited a typical quadripartite structure, with a large single copy (LSC) region of 87,170 - 88,408 bp, a small single copy (SSC) region of 18,650 - 19,038 bp and an inverted repeat (IR) region of 26,288 - 26,546 bp. Within Ulmus species, gene structure and content of cp genomes were highly conserved, although slight variations were found in the boundary of SC/IR regions. Moreover, genome-wide sliding window analysis uncovered the variability of ndhC-trnV-UAC, ndhF-rpl32, and psbI-trnS-GCU were higher among 31 Ulmus that may be useful for the population genetics and potential DNA barcodes. Two genes (rps15 and atpF) were further detected under a positive selection of Ulmus species. Comparative phylogenetic analysis based on the cp genome and protein-coding genes revealed consistent topology that U. mianzhuensis is a sister group to U. parvifolia (sect. Microptelea) with a relatively low-level nucleotide variation of the cp genome. Additionally, our analyses also found that the traditional taxonomic system of five sections in Ulmus is not supported by the current phylogenomic topology with a nested evolutionary relationship between sections.

CONCLUSIONS

Features of the cp genome length, GC content, organization, and gene order were highly conserved within Ulmus. Furthermore, molecular evidence from the low variation of the cp genome suggested that U. mianzhuensis should be merged into U. parvifolia and regarded as a subspecies of U. parvifolia. Overall, we demonstrated that the cp genome provides valuable information for understanding the genetic variation and phylogenetic relationship in Ulmus.

Analysis of the Chloroplast Genome of Ficus simplicissima Lour Collected in Vietnam and Proposed Barcodes for Identifying Ficus Plants

Ficus simplicissima Lour. is an Asian species of fig tree in the family Moraceae. The chloroplast (cp) genome of F. simplicissima m3 was sequenced using the Pacbio sequel platform. The F. simplicissima cpDNA has a size of 160,321 bp in length, of which GC content accounts for 36.13%. The cp genome of F. simplicissima consists of a single large copy (LSC) with a size of 91,346 bp, a single small copy (SSC) with a size of 20,131 bp, and a pair of inverted repeats with a size of 24,421 to 24,423 bp. The cp genome of F. simplicissima has 127 genes, including 85 protein-coding genes, eight rRNA genes, and 34 tRNA genes; 92 simple sequence repeats and 39 long repeats were detected in the cpDNA of F. simplicissim. A comparative cp genome analysis among six species in the Ficus genus indicated that the genome structure and gene content were highly conserved. The non-coding regions show more differentiation than the coding regions, and the LSC and SSC regions show more differences than the inverted repeat regions. Phylogenetic analysis supported that F. simplicissima m3 had a close relationship with F. hirta. The complete cp genome of F. simplicissima was proposed as a chloroplast DNA barcoding for genus-level in the Moraceae family and the psbA-trnH gene region for species-level identification.

Complete Plastome of Physalis angulata var. villosa, Gene Organization, Comparative Genomics and Phylogenetic Relationships among Solanaceae

Physalis angulata var. villosa, rich in withanolides, has been used as a traditional Chinese medicine for many years. To date, few extensive molecular studies of this plant have been conducted. In the present study, the plastome of P. angulata var. villosa was sequenced, characterized and compared with that of other Physalis species, and a phylogenetic analysis was conducted in the family Solanaceae. The plastome of P. angulata var. villosa was 156,898 bp in length with a GC content of 37.52%, and exhibited a quadripartite structure typical of land plants, consisting of a large single-copy (LSC, 87,108 bp) region, a small single-copy (SSC, 18,462 bp) region and a pair of inverted repeats (IR: IRA and IRB, 25,664 bp each). The plastome contained 131 genes, of which 114 were unique and 17 were duplicated in IR regions. The genome consisted of 85 protein-coding genes, eight rRNA genes and 38 tRNA genes. A total of 38 long, repeat sequences of three types were identified in the plastome, of which forward repeats had the highest frequency. Simple sequence repeats (SSRs) analysis revealed a total of 57 SSRs, of which the T mononucleotide constituted the majority, with most of SSRs being located in the intergenic spacer regions. Comparative genomic analysis among nine Physalis species revealed that the single-copy regions were less conserved than the pair of inverted repeats, with most of the variation being found in the intergenic spacer regions rather than in the coding regions. Phylogenetic analysis indicated a close relationship between Physalis and Withania. In addition, Iochroma, Dunalia, Saracha and Eriolarynx were paraphyletic, and clustered together in the phylogenetic tree. Our study published the first sequence and assembly of the plastome of P. angulata var. villosa, reported its basic resources for evolutionary studies and provided an important tool for evaluating the phylogenetic relationship within the family Solanaceae.

Comparative and phylogenetic analysis of complete chloroplast genomes from five Artemisia species

Artemisia Linn. is a large genus within the family Asteraceae that includes several important medicinal plants. Because of their similar morphology and chemical composition, traditional identification methods often fail to distinguish them. Therefore, developing an effective identification method for Artemisia species is an urgent requirement. In this study, we analyzed 15 chloroplast (cp) genomes, including 12 newly sequenced genomes, from 5 Artemisia species. The cp genomes from the five Artemisia species had a typical quadripartite structure and were highly conserved across species. They had varying lengths of 151,132-151,178 bp, and their gene content and codon preferences were similar. Mutation hotspot analysis identified four highly variable regions, which can potentially be used as molecular markers to identify Artemisia species. Phylogenetic analysis showed that the five Artemisia species investigated in this study were sister branches to each other, and individuals of each species formed a monophyletic clade. This study shows that the cp genome can provide distinguishing features to help identify closely related Artemisia species and has the potential to serve as a universal super barcode for plant identification.

Comparative Genomic Analysis Uncovers the Chloroplast Genome Variation and Phylogenetic Relationships of Camellia Species

Camellia is the largest genus in the family Theaceae. Due to phenotypic diversity, frequent hybridization, and polyploidization, an understanding of the phylogenetic relationships between Camellia species remains challenging. Comparative chloroplast (cp) genomics provides an informative resource for phylogenetic analyses of Camellia. In this study, 12 chloroplast genome sequences from nine Camellia species were determined using Illumina sequencing technology via de novo assembly. The cp genome sizes ranged from 156,545 to 157,021 bp and were organized into quadripartite regions with the typical angiosperm cp genomes. Each genome harbored 87 protein-coding, 37 transfer RNA, and 8 ribosomal RNA genes in the same order and orientation. Differences in long and short sequence repeats, SNPs, and InDels were detected across the 12 cp genomes. Combining with the complete cp sequences of seven other species in the genus Camellia, a total of nine intergenic sequence divergent hotspots and 14 protein-coding genes with high sequence polymorphism were identified. These hotspots, especially the InDel (~400 bp) located in atpH-atpI region, had sufficient potential to be used as barcode markers for further phylogenetic analysis and species identification. Principal component and phylogenetic analysis suggested that regional constraints, rather than functional constraints, strongly affected the sequence evolution of the cp genomes in this study. These cp genomes could facilitate the development of new molecular markers, accurate species identification, and investigations of the phylogenomic relationships of the genus Camellia.

Comparative chloroplast genome analysis of Ficus (Moraceae): Insight into adaptive evolution and mutational hotspot regions

As the largest genus in Moraceae, Ficus is widely distributed across tropical and subtropical regions and exhibits a high degree of adaptability to different environments. At present, however, the phylogenetic relationships of this genus are not well resolved, and chloroplast evolution in Ficus remains poorly understood. Here, we sequenced, assembled, and annotated the chloroplast genomes of 10 species of Ficus, downloaded and assembled 13 additional species based on next-generation sequencing data, and compared them to 46 previously published chloroplast genomes. We found a highly conserved genomic structure across the genus, with plastid genome sizes ranging from 159,929 bp (Ficus langkokensis) to 160,657 bp (Ficus religiosa). Most chloroplasts encoded 113 unique genes, including a set of 78 protein-coding genes, 30 transfer RNA (tRNA) genes, four ribosomal RNA (rRNA) genes, and one pseudogene (infA). The number of simple sequence repeats (SSRs) ranged from 67 (Ficus sagittata) to 89 (Ficus microdictya) and generally increased linearly with plastid size. Among the plastomes, comparative analysis revealed eight intergenic spacers that were hotspot regions for divergence. Additionally, the clpP, rbcL, and ccsA genes showed evidence of positive selection. Phylogenetic analysis indicated that none of the six traditionally recognized subgenera of Ficus were monophyletic. Divergence time analysis based on the complete chloroplast genome sequences showed that Ficus species diverged rapidly during the early to middle Miocene. This research provides basic resources for further evolutionary studies of Ficus.

Comparative Analysis of Complete Chloroplast Genomes of Nine Species of Litsea (Lauraceae): Hypervariable Regions, Positive Selection, and Phylogenetic Relationships

Litsea is a group of evergreen trees or shrubs in the laurel family, Lauraceae. Species of the genus are widely used for a wide range of medicinal and industrial aspects. At present, most studies related to the gene resources of Litsea are restricted to morphological analyses or features of individual genomes, and currently available studies of select molecular markers are insufficient. In this study, we assembled and annotated the complete chloroplast genomes of nine species in Litsea, carried out a series of comparative analyses, and reconstructed phylogenetic relationships within the genus. The genome length ranged from 152,051 to 152,747 bp and a total of 128 genes were identified. High consistency patterns of codon bias, repeats, divergent analysis, single nucleotide polymorphisms (SNP) and insertions and deletions (InDels) were discovered across the genus. Variations in gene length and the presence of the pseudogene ycf1^Ψ, resulting from IR contraction and expansion, are reported. The hyper-variable gene rpl16 was identified for its exceptionally high Ka/Ks and Pi values, implying that those frequent mutations occurred as a result of positive selection. Phylogenetic relationships were recovered for the genus based on analyses of full chloroplast genomes and protein-coding genes. Overall, both genome sequences and potential molecular markers provided in this study enrich the available genomic resources for species of Litsea. Valuable genomic resources and divergent analysis are also provided for further research of the evolutionary patterns, molecular markers, and deeper phylogenetic relationships of Litsea.

Fifteen complete chloroplast genomes of Trapa species (Trapaceae): insight into genome structure, comparative analysis and phylogenetic relationships

BACKGROUND

Trapa L. is a floating-leaved aquatic plant with important economic and ecological values. However, the species identification and phylogenetic relationship within Trapa are still controversial, which necessitates the need for plastid genome information of Trapa. In this study, complete chloroplast genomes of 13 Trapa species/taxa were sequenced and annotated. Combined with released sequences, comparative analyses of chloroplast genomes were performed on the 15 Trapa species/taxa for the first time.

RESULTS

The Trapa chloroplast genomes exhibited typical quadripartite structures with lengths from 155,453 to 155,559 bp. The gene orders and contents within Trapa were conservative, but several changes were found in the microstructure. The intron loss of rpl2, also detected in Lythraceae, was found in all Trapa species/taxa, suggesting close genetic relationship between Lythraceae and Trapaceae. Notably, two small-seed species (T. incisa and T. maximowiczii) showed the smallest genome size with 155,453 and 155,477 bp, respectively. Each cp genome contained the same 130 genes consisting of 85 protein-coding genes, 37 tRNA genes and 8 rRNA genes. Trapa species/taxa showed 37 (T. incisa and T. maximowiczii) to 41 (T. sibirica) long repeats, including forward, palindromic, reversed and complementary repeats. There were 110 (T. quadrispinosa) to 123 (T. incisa and T. maximowiczii) SSR (simple sequence repeat) loci in Trapa chloroplast genomes. Comparative analyses revealed that two hotspot regions (atpA-atpF and rps2-rpoC2) in Trapa chloroplast genomes could be served as potential molecular markers. Three phylogenetic analyses (ML, MP and BI) consistently showed that there were two clusters within Trapa, including large- and small-seed species/taxa, respectively; for the large-seed Trapa, they clustered according to their geographical origin and tubercle morphology on the surface of seeds.

CONCLUSION

In summary, we have acquired the sequences of 13 Trapa chloroplast genomes, and performed the comparative analyses within Trapa for the first time. The results have helped us better identify the Trapa species/taxa and deepen the understanding of genetic basis and phylogenetic relationship of Trapa, which will facilitate the effective management and utilization of the important genetic resources in the future.

Comparative Analysis of Chloroplast Genome Structure and Phylogenetic Relationships Among Six Taxa Within the Genus Catalpa (Bignoniaceae)

Species within the Genus Catalpa are mostly semievergreen or deciduous trees with opposite or whorled leaves. C. bungei, C. fargesii f. duclouxii and C. fargesii are sources of traditional precious wood in China, known as the “kings of wood”. Due to a lack of phenotypic and molecular studies and insufficient sequence information, intraspecific morphological differences, common DNA barcodes and partial sequence fragments cannot clearly reveal the phylogenetic or intraspecific relationships within Catalpa. Therefore, we sequenced the complete chloroplast genomes of six taxa of the genus Catalpa and analyzed their basic structure and evolutionary relationships. The chloroplast genome of Catalpa shows a typical tetrad structure with a total length ranging from 157,765 bp (C. fargesii) to 158,355 bp (C. ovata). The length of the large single-copy (LSC) region ranges from 84,599 bp (C. fargesii) to 85,004 bp (C. ovata), that of the small single-copy (SSC) region ranges from 12,662 bp (C. fargesii) to 12,675 bp (C. ovata), and that of the inverted repeat (IR) regions ranges from 30,252 bp (C. fargesii) to 30,338 bp (C. ovata). The GC content of the six chloroplast genomes were 38.1%. In total, 113 unique genes were detected, and there were 19 genes in IR regions. The 113 genes included 79 protein-coding genes, 30 tRNA genes and four rRNA genes. Five hypervariable regions (trnH-psbA, rps2-rpoC2, rpl22, ycf15-trnl-CAA and rps15) were identified by analyzing chloroplast nucleotide polymorphisms, which might be serve as potential DNA barcodes for the species. Comparative analysis showed that single nucleotide polymorphisms (SNPs) and simple sequence repeats (SSRs) were highly diverse in the six species. Codon usage patterns were highly similar among the taxa included in the present study. In addition to the stop codons, all codons showed a preference for ending in A or T. Phylogenetic analysis of the entire chloroplast genome showed that all taxa within the genus Catalpa formed a monophyletic group, clearly reflecting the relationships within the genus. This study provides information on the chloroplast genome sequence, structural variation, codon bias and phylogeny of Catalpa, which will facilitate future research efforts.

Supporting in situ conservation of the genetic diversity of crop wild relatives using genomic technologies

The last decade has witnessed huge technological advances in genomics, particularly in DNA sequencing. Here, we review the actual and potential application of genomics in supporting in situ conservation of crop wild relatives (CWRs). In addition to helping in prioritization of protection of CWR taxa and in situ conservation sites, genome analysis is allowing the identification of novel alleles that need to be prioritized for conservation. Genomics is enabling the identification of potential sources of important adaptive traits that can guide the establishment or enrichment of in situ genetic reserves. Genomic tools also have the potential for developing a robust framework for monitoring and reporting genome‐based indicators of genetic diversity changes associated with factors such as land use or climate change. These tools have been demonstrated to have an important role in managing the conservation of populations, supporting sustainable access and utilization of CWR diversity, enhancing accelerated domestication of new crops and forensic genomics thus preventing misappropriation of genetic resources. Despite this great potential, many policy makers and conservation managers have failed to recognize and appreciate the need to accelerate the application of genomics to support the conservation and management of biodiversity in CWRs to underpin global food security. Funding and inadequate genomic expertise among conservation practitioners also remain major hindrances to the widespread application of genomics in conservation.

The complete chloroplast genomes of three Hamamelidaceae species: Comparative and phylogenetic analyses

Hamamelidaceae is an important group that represents the origin and early evolution of angiosperms. Its plants have many uses, such as timber, medical, spice, and ornamental uses. In this study, the complete chloroplast genomes of Loropetalum chinense (R. Br.) Oliver, Corylopsis glandulifera Hemsl., and Corylopsis velutina Hand.‐Mazz. were sequenced using the Illumina NovaSeq 6000 platform. The sizes of the three chloroplast genomes were 159,402 bp (C. glandulifera), 159,414 bp (C. velutina), and 159,444 bp (L. chinense), respectively. These chloroplast genomes contained typical quadripartite structures with a pair of inverted repeat (IR) regions (26,283, 26,283, and 26,257 bp), a large single‐copy (LSC) region (88,134, 88,146, and 88,160 bp), and a small single‐copy (SSC) region (18,702, 18,702, and 18,770 bp). The chloroplast genomes encoded 132–133 genes, including 85–87 protein‐coding genes, 37–38 tRNA genes, and 8 rRNA genes. The coding regions were composed of 26,797, 26,574, and 26,415 codons, respectively, most of which ended in A/U. A total of 37–43 long repeats and 175–178 simple sequence repeats (SSRs) were identified, and the SSRs contained a higher number of A + T than G + C bases. The genome comparison showed that the IR regions were more conserved than the LSC or SSC regions, while the noncoding regions contained higher variability than the gene coding regions. Phylogenetic analyses revealed that species in the same genus tended to cluster together. Chunia Hung T. Chang, Mytilaria Lecomte, and Disanthus Maxim. may have diverged early and Corylopsis Siebold & Zucc. was closely related to Loropetalum R. Br. This study provides valuable information for further species identification, evolution, and phylogenetic studies of Hamamelidaceae plants.

Comparative Analysis the Complete Chloroplast Genomes of Nine Musa Species: Genomic Features, Comparative Analysis, and Phylogenetic Implications

Musa (family Musaceae) is monocotyledonous plants in order Zingiberales, which grows in tropical and subtropical regions. It is one of the most important tropical fruit trees in the world. Herein, we used next-generation sequencing technology to assemble and perform in-depth analysis of the chloroplast genome of nine new Musa plants for the first time, including genome structure, GC content, repeat structure, codon usage, nucleotide diversity and etc. The entire length of the Musa chloroplast genome ranged from 167,975 to 172,653 bp, including 113 distinct genes comprising 79 protein-coding genes, 30 transfer RNA (tRNA) genes and four ribosomal RNA (rRNA) genes. In comparative analysis, we found that the contraction and expansion of the inverted repeat (IR) regions resulted in the doubling of the rps19 gene. The several non-coding sites (psbI-atpA, atpH-atpI, rpoB-petN, psbM-psbD, ndhf-rpl32, and ndhG-ndhI) and three genes (ycf1, ycf2, and accD) showed significant variation, indicating that they have the potential of molecular markers. Phylogenetic analysis based on the complete chloroplast genome and coding sequences of 77 protein-coding genes confirmed that Musa can be mainly divided into two groups. These genomic sequences provide molecular foundation for the development and utilization of Musa plants resources. This result may contribute to the understanding of the evolution pattern, phylogenetic relationships as well as classification of Musa plants.

Comparative plastomes of Carya species provide new insights into the plastomes evolution and maternal phylogeny of the genus

Carya, in the Juglandiodeae subfamily, is to a typical temperate-subtropical forest-tree genus for studying the phylogenetic evolution and intercontinental disjunction between eastern Asia (EA) and North America (NA). Species of the genus have high economic values worldwide for their high-quality wood and the rich healthy factors of their nuts. Although previous efforts based on multiple molecular markers or genome-wide SNPs supported the monophyly of Carya and its two EA and NA major subclades, the maternal phylogeny of Carya still need to be comprehensively evaluated. The variation of Carya plastome has never been thoroughly characterized. Here, we novelly present 19 newly generated plastomes of congeneric Carya species, including the recently rediscovered critically endangered C. poilanei. The overall assessment of plastomes revealed highly conservative in the general structures. Our results indicated that remarkable differences in several plastome features are highly consistent with the EA-NA disjunction and showed the relatively diverse matrilineal sources among EA Carya compared to NA Carya. The maternal phylogenies were conducted with different plastome regions and full-length plastome datasets from 30 plastomes, representing 26 species in six genera of Juglandoideae and Myrica rubra (as root). Six out of seven phylogenetic topologies strongly supported the previously reported relationships among genera of Juglandoideae and the two subclades of EA and NA Carya, but displayed significant incongruencies between species within the EA and NA subclades. The phylogenetic tree generated from full-length plastomes demonstrated the optimal topology and revealed significant geographical maternal relationships among Carya species, especially for EA Carya within overlapping distribution areas. The full-length plastome-based phylogenetic topology also strongly supported the taxonomic status of five controversial species as separate species of Carya. Historical and recent introgressive hybridization and plastid captures might contribute to plastome geographic patterns and inconsistencies between topologies built from different datasets, while incomplete lineage sorting could account for the discordance between maternal topology and the previous nuclear genome data-based phylogeny. Our findings highlight full-length plastomes as an ideal tool for exploring maternal relationships among the subclades of Carya, and potentially in other outcrossing perennial woody plants, for resolving plastome phylogenetic relationships.

Comparative chloroplast genome analysis of medicinally important Veratrum (Melanthiaceae) in China: Insights into genomic characterization and phylogenetic relationships

Members of Veratrum are perennial herbs widely used in traditional Chinese medicine to induce vomiting, resolve blood stasis and relieve pain. However, the intrageneric classification and phylogenetic relationships within Veratrum have long been controversial due to the complexity of morphological variations and lack of high-resolution molecular markers. In this study, we reevaluated the infrageneric relationships with the genus Veratrum using complete chloroplast genome sequence data. Herein, the complete cp genomes of ten species of Veratrum were newly sequenced and characterized. The complete cp genomes of ten species of Veratrum had the typical quadripartite structure, ranging from 151,597 bp to 153,711 bp in size and comprising a total of 135 genes. The structure of Veratrum cp genomes (i.e., gene order, content, and genome components) was highly similar across species. The number of simple sequence repeats (SSRs) ranged from 63 to 78, and of long repeats ranged from 31 to 35. Eight highly divergent regions (ndhF, psbC-psbZ, psbK-psbI, rpoB-trnC_GCA, trnK_UUU-trnQ_UUG, trnS_GCU-trnG_UCC, trnT_UGU-trnL_UAA and ycf1) were identified and are potentially useful for the DNA barcoding of Veratrum. Phylogenetic analysis among 29 taxa based on cp genomes, total genes, protein-coding genes and intergenic regions strongly supported the monophyly of Veratrum. The circumscription and relationships of the infrageneric taxa of Veratrum were well-presented with great resolution. These results will facilitate the identification, taxonomy, and utilization of Veratrum plants as well as the evolutionary studies of Melanthiaceae.

Complete Chloroplast Genome Sequence of Sonchus brachyotus Helps to Elucidate Evolutionary Relationships with Related Species of Asteraceae

Sonchus brachyotus DC. possesses both edible and medicinal properties and is widely distributed throughout China. In this study, the complete cp genome of S. brachyotus was sequenced and assembled. The total length of the complete S. brachyotus cp genome was 151,977 bp, including an LSC region of 84,553 bp, SSC region of 18,138 bp, and IR region of 24,643 bp. Sequence analyses revealed that the cp genome encoded 132 genes, including 87 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. The GC content was 37.6%. One hundred mononucleotide microsatellites, 4 dinucleotide microsatellites, 67 trinucleotide microsatellites, 4 tetranucleotide microsatellites, and 1 long repeat were identified. The SSR frequency of the LSC region was significantly greater than that of the IR and SSC regions. In total, 175 SSRs and highly variable regions were recognized as potential cp markers. By analyzing the IR/LSC and IR/SSC boundaries, structural differences between S. brachyotus and 6 other species were detected. According to phylogenetic analyses, S. brachyotus was most closely related to S. arvensis and S. oleraceus. Overall, this study provides complete cp genome resources for S. brachyotus that will be beneficial for identifying potential molecular markers and evolutionary patterns of S. brachyotus and its closely related species.

New Insight into the Phylogeny and Taxonomy of Cultivated and Related Species of Crataegus in China, Based on Complete Chloroplast Genome Sequencing

Hawthorns (Crataegus L.) are one of the most important processing and table fruits in China, due to their medicinal properties and health benefits. However, the interspecific relationships and evolution history of cultivated Crataegus in China remain unclear. Our previously published data showed C. bretschneideri may be derived from the hybridization of C. pinnatifida with C. maximowiczii, and that introgression occurs between C. hupehensis, C. pinnatifida, and C. pinnatifida var. major. In the present study, chloroplast sequences were used to further elucidate the phylogenetic relationships of cultivated Crataegus native to China. The chloroplast genomes of three cultivated species and one related species of Crataegus were sequenced for comparative and phylogenetic analyses. The four chloroplast genomes of Crataegus exhibited typical quadripartite structures and ranged from 159,607 bp (C. bretschneideri) to 159,875 bp (C. maximowiczii) in length. The plastomes of the four species contained 113 genes consisting of 79 protein-coding genes, 30 tRNA genes, and 4 rRNA genes. Six hypervariable regions (ndhC-trnV(UAC)-trnM(CAU), ndhA, atpH-atpI, ndhF, trnR(UCU)-atpA, and ndhF-rpl32), 196 repeats, and a total of 386 simple sequence repeats were detected as potential variability makers for species identification and population genetic studies. In the phylogenomic analyses, we also compared the entire chloroplast genomes of three published Crataegus species: C. hupehensis (MW201730.1), C. pinnatifida (MN102356.1), and C. marshallii (MK920293.1). Our phylogenetic analyses grouped the seven Crataegus taxa into two main clusters. One cluster included C. bretschneideri, C. maximowiczii, and C. marshallii, whereas the other included C. hupehensis, C. pinnatifida, and C. pinnatifida var. major. Taken together, our findings indicate that C. maximowiczii is the maternal origin of C. bretschneideri. This work provides further evidence of introgression between C. hupehensis, C. pinnatifida, and C. pinnatifida var. major, and suggests that C. pinnatifida var. major might have been artificially selected and domesticated from hybrid populations, rather than evolved from C. pinnatifida.

Comparative chloroplast genome analysis of Impatiens species (Balsaminaceae) in the karst area of China: insights into genome evolution and phylogenomic implications

BACKGROUND

Impatiens L. is a genus of complex taxonomy that belongs to the family Balsaminaceae (Ericales) and contains approximately 1000 species. The genus is well known for its economic, medicinal, ornamental, and horticultural value. However, knowledge about its germplasm identification, molecular phylogeny, and chloroplast genomics is limited, and taxonomic uncertainties still exist due to overlapping morphological features and insufficient genomic resources.

RESULTS

We sequenced the chloroplast genomes of six different species (Impatiens chlorosepala, Impatiens fanjingshanica, Impatiens guizhouensis, Impatiens linearisepala, Impatiens loulanensis, and Impatiens stenosepala) in the karst area of China and compared them with those of six previously published Balsaminaceae species. We contrasted genomic features and repeat sequences, assessed sequence divergence and constructed phylogenetic relationships. Except for those of I. alpicola, I. pritzelii and I. glandulifera, the complete chloroplast genomes ranging in size from 151,366 bp (I. alpicola) to 154,189 bp (Hydrocera triflora) encoded 115 distinct genes [81 protein-coding, 30 transfer RNA (tRNA), and 4 ribosomal RNA (rRNA) genes]. Moreover, the characteristics of the long repeat sequences and simple sequence repeats (SSRs) were determined. psbK-psbI, trnT-GGU-psbD, rpl36-rps8, rpoB-trnC-GCA, trnK-UUU-rps16, trnQ-UUG, trnP-UGG-psaJ, trnT-UGU-trnL-UAA, and ycf4-cemA were identified as divergence hotspot regions and thus might be suitable for species identification and phylogenetic studies. Additionally, the phylogenetic relationships based on Maximum likelihood (ML) and Bayesian inference (BI) of the whole chloroplast genomes showed that the chloroplast genome structure of I. guizhouensis represents the ancestral state of the Balsaminaceae family.

CONCLUSION

Our study provided detailed information about nucleotide diversity hotspots and the types of repeats, which can be used to develop molecular markers applicable to Balsaminaceae species. We also reconstructed and analyzed the relationships of some Impatiens species and assessed their taxonomic statuses based on the complete chloroplast genomes. Together, the findings of the current study might provide valuable genomic resources for systematic evolution of the Balsaminaceae species.

Plant super-barcode: a case study on genome-based identification for closely related species of Fritillaria

Background

Although molecular analysis offers a wide range of options for species identification, a universal methodology for classifying and distinguishing closely related species remains elusive. This study validated the effectiveness of utilizing the entire chloroplast (cp) genome as a super-barcode to help identify and classify closely related species.

Methods

We here compared 26 complete cp genomes of ten Fritillaria species including 18 new sequences sequenced in this study. Each species had repeats and the cp genomes were used as a whole DNA barcode to test whether they can distinguish Fritillaria species.

Results

The cp genomes of Fritillaria medicinal plants were conserved in genome structure, gene type, and gene content. Comparison analysis of the Fritillaria cp genomes revealed that the intergenic spacer regions were highly divergent compared with other regions. By constructing the phylogenetic tree by the maximum likelihood and maximum parsimony methods, we found that the entire cp genome showed a high discrimination power for Fritillaria species with individuals of each species in a monophyletic clade. These results indicate that cp genome can be used to effectively differentiate medicinal plants from the genus Fritillaria at the species level.

Conclusions

This study implies that cp genome can provide distinguishing differences to help identify closely related Fritillaria species, and has the potential to be served as a universal super-barcode for plant identification.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13020-021-00460-z.

New Insight into Taxonomy of European Mountain Pines, Pinus mugo Complex, Based on Complete Chloroplast Genomes Sequencing

The Pinus mugo complex is a large group of closely related mountain pines, which are an important component of the ecosystems of the most important mountain ranges, such as the Alps, Carpathians and Pyrenees. The phylogenetic relationships between taxa in this complex have been under discussion for many years. Despite the use of many different approaches, they still need to be clarified and supplemented with new data, especially those obtained with high-throughput methods. Therefore, in this study, the complete sequences of the chloroplast genomes of the three most recognized members of the Pinus mugo complex, i.e., Pinus mugo, Pinus rotundata and Pinus uncinata, were sequenced and analyzed to gain new insight into their phylogenetic relationships. Comparative analysis of their complete chloroplast genome sequences revealed several mutational hotspots potentially useful for the genetic identification of taxa from the Pinus mugo complex. Phylogenetic inference based on sixteen complete chloroplast genomes of different coniferous representatives showed that pines from the Pinus mugo complex form one distinct monophyletic group. The results obtained in this study provide new and valuable omics data for further research within the European mountain pine complex. They also indicate which regions may be useful in the search for diagnostic DNA markers for the members of Pinus mugo complex and set the baseline in the conservation of genetic resources of its endangered taxa.

Holocene chloroplast genetic variation of shrubs (Alnus alnobetula, Betula nana, Salix sp.) at the siberian tundra-taiga ecotone inferred from modern chloroplast genome assembly and sedimentary ancient DNA analyses

Climate warming alters plant composition and population dynamics of arctic ecosystems. In particular, an increase in relative abundance and cover of deciduous shrub species (shrubification) has been recorded. We inferred genetic variation of common shrub species (Alnus alnobetula, Betula nana, Salix sp.) through time. Chloroplast genomes were assembled from modern plants (n = 15) from the Siberian forest-tundra ecotone. Sedimentary ancient DNA (sedaDNA; n = 4) was retrieved from a lake on the southern Taymyr Peninsula and analyzed by metagenomics shotgun sequencing and a hybridization capture approach. For A. alnobetula, analyses of modern DNA showed low intraspecies genetic variability and a clear geographical structure in haplotype distribution. In contrast, B. nana showed high intraspecies genetic diversity and weak geographical structure. Analyses of sedaDNA revealed a decreasing relative abundance of Alnus since 5,400 cal yr BP, whereas Betula and Salix increased. A comparison between genetic variations identified in modern DNA and sedaDNA showed that Alnus variants were maintained over the last 6,700 years in the Taymyr region. In accordance with modern individuals, the variants retrieved from Betula and Salix sedaDNA showed higher genetic diversity. The success of the hybridization capture in retrieving diverged sequences demonstrates the high potential for future studies of plant biodiversity as well as specific genetic variation on ancient DNA from lake sediments. Overall, our results suggest that shrubification has species-specific trajectories. The low genetic diversity in A. alnobetula suggests a local population recruitment and growth response of the already present communities, whereas the higher genetic variability and lack of geographical structure in B. nana may indicate a recruitment from different populations due to more efficient seed dispersal, increasing the genetic connectivity over long distances.

Chloroplast genome variation and phylogenetic relationships of Atractylodes species

Background

Atractylodes DC is the basic original plant of the widely used herbal medicines “Baizhu” and “Cangzhu” and an endemic genus in East Asia. Species within the genus have minor morphological differences, and the universal DNA barcodes cannot clearly distinguish the systemic relationship or identify the species of the genus. In order to solve these question, we sequenced the chloroplast genomes of all species of Atractylodes using high-throughput sequencing.

Results

The results indicate that the chloroplast genome of Atractylodes has a typical quadripartite structure and ranges from 152,294 bp (A. carlinoides) to 153,261 bp (A. macrocephala) in size. The genome of all species contains 113 genes, including 79 protein-coding genes, 30 transfer RNA genes and four ribosomal RNA genes. Four hotspots, rpl22-rps19-rpl2, psbM-trnD, trnR-trnT^(GGU), and trnT^(UGU)-trnL, and a total of 42–47 simple sequence repeats (SSR) were identified as the most promising potentially variable makers for species delimitation and population genetic studies. Phylogenetic analyses of the whole chloroplast genomes indicate that Atractylodes is a clade within the tribe Cynareae; Atractylodes species form a monophyly that clearly reflects the relationship within the genus.

Conclusions

Our study included investigations of the sequences and structural genomic variations, phylogenetics and mutation dynamics of Atractylodes chloroplast genomes and will facilitate future studies in population genetics, taxonomy and species identification.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07394-8.

Identification of plastid genomic regions inferring species identity from de novo plastid genome assembly of 14 Korean-native Iris species (Iridaceae)

Iris is one of the largest genera in the family Iridaceae, comprising hundreds of species, including numerous economically important horticultural plants used in landscape gardening and herbal medicine. Improved taxonomic classification of Iris species, particularly the endangered Korean-native Iris, is needed for correct species delineation. To this end, identification of diverse genetic markers from Iris genomes would facilitate molecular identification and resolve ambiguous classifications from molecular analyses; however, only two Iris plastid genomes, from Iris gatesii and Iris sanguinea, have been sequenced. Here, we used high-throughput next-generation sequencing, combined with Sanger sequencing, to construct the plastid genomes of 14 Korean-native Iris species with one outgroup and predict their gene content. Using these data, combined with previously published plastid genomes from Iris and one outgroup (Sisyrinchium angustifolium), we constructed a Bayesian phylogenetic tree showing clear speciation among the samples. We further identified sub-genomic regions that have undergone neutral evolution and accurately recapitulate Bayesian-inferred speciation. These contain key markers that could be used to identify and classify Iris samples into taxonomic clades. Our results confirm previously reported speciation patterns and resolve questionable relationships within the Iris genus. These data also provide a valuable resource for studying genetic diversity and refining phylogenetic relationships between Iris species.

Insights into phylogeny, age and evolution of Allium (Amaryllidaceae) based on the whole plastome sequences

BACKGROUND AND AIMS

The genus Allium L., one of the largest monocotyledonous genera and one that includes many economically important crops with nutritional and medicinal value, has been the focus of classification or phylogeny studies for centuries. Recent studies suggested that the genus can be divided into 15 subgenera and 72 sections, which were further classified into three evolutionary lineages. However, the phylogenetic relationships reconstructed by one or two loci showed weaker support, especially for the third evolutionary lineage, which might not show the species relationships very clearly and could hinder further adaptive and evolutionary study.

METHODS

In this study, a total of 39 complete chloroplast genomes of Allium (covering 12 Allium subgenera) were collected, and combining these with 125 species of plastomes from 19 other families of monocots, we reconstructed the phylogeny of the genus Allium, estimated the origin and divergence time of the three evolutionary lineages and investigated the adaptive evolution in this genus and related families.

RESULTS

Our phylogenetic analysis confirmed the monophyly and three evolutionary lineages of Allium, while new species relationships were detected within the third evolutionary lineage. The divergence time of the three evolutionary lineages was estimated to be in the early Eocene to the middle Miocene, and numerous positive selected genes (PSGs) and PSGs with high average Ka/Ks values were found in Allium species.

CONCLUSIONS

Our results detected a well-supported phylogenetic relationship of Allium. The PSGs and PSGs with high Ka/Ks values, as well as diversified morphologies, complicated chromosome characteristics and unique reproductive modes may play important roles in the adaptation and evolution of Allium species. This is the first study that conducted phylogenetic and evolutionary analyses on the genus Allium combined with the plastome and morphological and cytological data. We hope that this study can contribute to further analysis of Allium for other researchers.

Complete chloroplast genomes of four Physalis species (Solanaceae): lights into genome structure, comparative analysis, and phylogenetic relationships

BACKGROUND

Physalis L. is a genus of herbaceous plants of the family Solanaceae, which has important medicinal, edible, and ornamental values. The morphological characteristics of Physalis species are similar, and it is difficult to rapidly and accurately distinguish them based only on morphological characteristics. At present, the species classification and phylogeny of Physalis are still controversial. In this study, the complete chloroplast (cp) genomes of four Physalis species (Physalis angulata, P. alkekengi var. franchetii, P. minima and P. pubescens) were sequenced, and the first comprehensive cp genome analysis of Physalis was performed, which included the previously published cp genome sequence of Physalis peruviana.

RESULTS

The Physalis cp genomes exhibited typical quadripartite and circular structures, and were relatively conserved in their structure and gene synteny. However, the Physalis cp genomes showed obvious variations at four regional boundaries, especially those of the inverted repeat and the large single-copy regions. The cp genomes' lengths ranged from 156,578 bp to 157,007 bp. A total of 114 different genes, 80 protein-coding genes, 30 tRNA genes, and 4 rRNA genes, were observed in four new sequenced Physalis cp genomes. Differences in repeat sequences and simple sequence repeats were detected among the Physalis cp genomes. Phylogenetic relationships among 36 species of 11 genera of Solanaceae based on their cp genomes placed Physalis in the middle and upper part of the phylogenetic tree, with a monophyletic evolution having a 100% bootstrap value.

CONCLUSION

Our results enrich the data on the cp genomes of the genus Physalis. The availability of these cp genomes will provide abundant information for further species identification, increase the taxonomic and phylogenetic resolution of Physalis, and assist in the investigation and utilization of Physalis plants.

Comparative Analysis of the Complete Chloroplast Genome of Mainland Aster spathulifolius and Other Aster Species

Aster spathulifolius, a common ornamental and medicinal plant, is widely distributed in Korea and Japan, and is genetically classified into mainland and island types. Here, we sequenced the whole chloroplast genome of mainland A. spathulifolius and compared it with those of the island type and other Aster species. The chloroplast genome of mainland A. spathulifolius is 152,732 bp with a conserved quadripartite structure, has 37.28% guanine-cytosine (GC) content, and contains 114 non-redundant genes. Comparison of the chloroplast genomes between the two A. spathulifolius lines and the other Aster species revealed that their sequences, GC contents, gene contents and orders, and exon-intron structure were well conserved; however, differences were observed in their lengths, repeat sequences, and the contraction and expansion of the inverted repeats. The variations were mostly in the single-copy regions and non-coding regions, which, together with the detected simple sequence repeats, could be used for the development of molecular markers to distinguish between these plants. All Aster species clustered into a monophyletic group, but the chloroplast genome of mainland A. spathulifolius was more similar to the other Aster species than to that of the island A. spathulifolius. The accD and ndhF genes were detected to be under positive selection within the Aster lineage compared to other related taxa. The complete chloroplast genome of mainland A. spathulifolius presented in this study will be helpful for species identification and the analysis of the genetic diversity, evolution, and phylogenetic relationships in the Aster genus and the Asteraceae.

Chloroplast genomes of seven species of Coryloideae (Betulaceae): structures and comparative analysis

Coryloideae is a subfamily in the family Betulaceae consisting of four extant genera: Carpinus, Corylus, Ostrya, and Ostryopsis. We sequenced the plastomes of six species of Corylus and one species of Ostryopsis for comparative and phylogenetic analyses. The plastomes are 159-160 kb long and possess typical quadripartite cp architecture. The plastomes show moderate divergence and conserved arrangement. Five mutational hotspots were identified by comparing the plastomes of seven species of Coryloideae: trnG-atpA, trnF-ndhJ, accD-psaI, ndhF-ccsA, and ycf1. We assembled the most complete phylogenomic tree for the family Betulaceae using 68 plastomes. Our cp genomic sequence phylogenetic analyses placed Carpinus, Ostrya, and Ostryopsis in a clade together and left Corylus in a separate clade. Within the genus Corylus, these analyses indicate the existence of five subclades reflecting the phylogeographical relationships among the species. The data offer significant genetic information for the identification of species of the Coryloideae, taxonomic and phylogenetic studies, and molecular breeding.

Comparison of different annotation tools for characterization of the complete chloroplast genome of Corylus avellana cv Tombul

BACKGROUND

Several bioinformatics tools have been designed for assembly and annotation of chloroplast (cp) genomes, making it difficult to decide which is most useful and applicable to a specific case. The increasing number of plant genomes provide an opportunity to accurately obtain cp genomes from whole genome shotgun (WGS) sequences. Due to the limited genetic information available for European hazelnut (Corylus avellana L.) and as part of a genome sequencing project, we analyzed the complete chloroplast genome of the cultivar 'Tombul' with multiple annotation tools.

RESULTS

Three different annotation strategies were tested, and the complete cp genome of C. avellana cv Tombul was constructed, which was 161,667 bp in length, and had a typical quadripartite structure. A large single copy (LSC) region of 90,198 bp and a small single copy (SSC) region of 18,733 bp were separated by a pair of inverted repeat (IR) regions of 26,368 bp. In total, 125 predicted functional genes were annotated, including 76 protein-coding, 25 tRNA, and 4 rRNA unique genes. Comparative genomics indicated that the cp genome sequences were relatively highly conserved in species belonging to the same order. However, there were still some variations, especially in intergenic regions, that could be used as molecular markers for analyses of phylogeny and plant identification. Simple sequence repeat (SSR) analysis showed that there were 83 SSRs in the cp genome of cv Tombul. Phylogenetic analysis suggested that C. avellana cv Tombul had a close affinity to the sister group of C. fargesii and C. chinensis, and then a closer evolutionary relationship with Betulaceae family than other species of Fagales.

CONCLUSION

In this study, the complete cp genome of Corylus avellana cv Tombul, the most widely cultivated variety in Turkey, was obtained and annotated, and additionally phylogenetic relationships were predicted among Fagales species. Our results suggest a very accurate assembly of chloroplast genome from next generation whole genome shotgun (WGS) sequences. Enhancement of taxon sampling in Corylus species provide genomic insights into phylogenetic analyses. The nucleotide sequences of cv Tombul cp genomes can provide comprehensive genetic insight into the evolution of genus Corylus.

Complete Chloroplast Genome Sequence and Phylogenetic Analysis of Quercus bawanglingensis Huang, Li et Xing, a Vulnerable Oak Tree in China

Quercus bawanglingensis Huang, Li et Xing, an endemic evergreen oak of the genus Quercus (Fagaceae) in China, is currently listed in the Red List of Chinese Plants as a vulnerable (VU) plant. No chloroplast (cp) genome information is currently available for Q. bawanglingensis, which would be essential for the establishment of guidelines for its conservation and breeding. In the present study, the cp genome of Q. bawanglingensis was sequenced and assembled into double-stranded circular DNA with a length of 161,394 bp. Two inverted repeats (IRs) with a total of 51,730 bp were identified, and the rest of the sequence was separated into two single-copy regions, namely, a large single-copy (LSC) region (90,628 bp) and a small single-copy (SSC) region (19,036 bp). The genome of Q. bawanglingensis contains 134 genes (86 protein-coding genes, 40 tRNAs and eight rRNAs). More forward (29) than inverted long repeats (21) are distributed in the cp genome. A simple sequence repeat (SSR) analysis showed that the genome contains 82 SSR loci, involving 84.15% A/T mononucleotides. Sequence comparisons among the nine complete cp genomes, including the genomes of Q. bawanglingensis, Q. tarokoensis Hayata (NC036370), Q. aliena var. acutiserrata Maxim. ex Wenz. (KU240009), Q. baronii Skan (KT963087), Q. aquifolioides Rehd. et Wils. (KX911971), Q. variabilis Bl. (KU240009), Fagus engleriana Seem. (KX852398), Lithocarpus balansae (Drake) A. Camus (KP299291) and Castanea mollissima Bl. (HQ336406), demonstrated that the diversity of SC regions was higher than that of IR regions, which might facilitate identification of the relationships within this extremely complex family. A phylogenetic analysis showed that Fagus engleriana and Trigonobalanus doichangensis form the basis of the produced evolutionary tree. Q. bawanglingensis and Q. tarokoensis, which belong to the group Ilex, share the closest relationship. The analysis of the cp genome of Q. bawanglingensis provides crucial genetic information for further studies of this vulnerable species and the taxonomy, phylogenetics and evolution of Quercus.

Development of a reliable Corylus sp. reference database through the implementation of a DNA fingerprinting test

This DNA fingerprinting test confirmed 195 unique Corylus sp. accessions that were used to build a reference database for identity verification of unknown hazelnut trees from three locations in Ontario. Hazelnut is one of the most profitable tree nuts worldwide. Development of a hazelnut industry in Ontario is urgently required, but economically important cultivars must be genetically verified first in order to meet industry standards. Traditional methods for cultivar identification are largely trait-based and unreliable. In this study, a multiplexed fingerprinting test was modified to allow for hazelnut cultivar discrimination at the DNA level. Fourteen highly polymorphic SSR markers covering the 11 linkage groups of Corylus genome were PCR amplified in multiplex using fluorescent-labelled primers. PCR conditions and primer physical properties were optimized to generate a clear signal for each locus. The 14 SSRs were used to fingerprint 195 unique Corylus accessions collected from the USDA-NCGR. Fragment sizes were subjected to a UPGMA clustering analysis which separated Corylus accessions based on species and geographic origin. For validation purposes, hazelnut leaves from three locations in Ontario were collected for identity verification using this DNA fingerprinting test. As a result, 33.3% of the unknown trees were duplicates of seven distinct genotypes and a small percentage (8.3%) of these were identical to reference Corylus hybrids. These results reflect common mislabelling issues and genotype duplications that can prevent a uniform plant propagation system. Implementation of this test together with the addition of more unique accessions to the reference database will help verification of trueness-to-type of economically important cultivars for the hazelnut industry.

The complete chloroplast genomes of three Betulaceae species: implications for molecular phylogeny and historical biogeography

Background

Previous phylogenetic conclusions on the family Betulaceae were based on either morphological characters or traditional single loci, which may indicate some limitations. The chloroplast genome contains rich polymorphism information, which is very suitable for phylogenetic studies. Thus, we sequenced the chloroplast genome sequences of three Betulaceae species and performed multiple analyses to investigate the genome variation, resolve the phylogenetic relationships, and clarify the divergence history.

Methods

Chloroplast genomes were sequenced using the high-throughput sequencing. A comparative genomic analysis was conducted to examine the global genome variation and screen the hotspots. Three chloroplast partitions were used to reconstruct the phylogenetic relationships using Maximum Likelihood and Bayesian Inference approaches. Then, molecular dating and biogeographic inferences were conducted based on the whole chloroplast genome data.

Results

Betulaceae chloroplast genomes consisted of a small single-copy region and a large single copy region, and two copies of inverted repeat regions. Nine hotspots can be used as potential DNA barcodes for species delimitation. Phylogenies strongly supported the division of Betulaceae into two subfamilies: Coryloideae and Betuloideae. The phylogenetic position of Ostryopsis davidiana was controversial among different datasets. The divergence time between subfamily Coryloideae and Betuloideae was about 70.49 Mya, and all six extant genera were inferred to have diverged fully by the middle Oligocene. Betulaceae ancestors were probably originated from the ancient Laurasia.

Discussions

This research elucidates the potential of chloroplast genome sequences in the application of developing molecular markers, studying evolutionary relationships and historical dynamic of Betulaceae.It also reveals the advantages of using chloroplast genome data to illuminate those phylogenies that have not been well solved yet by traditional approaches in other plants.