The first complete plastid genomes of Melastomataceae are highly structurally conserved

Comparative chloroplast genomics and insights into the molecular evolution of Tanaecium (Bignonieae, Bignoniaceae)

Species of Tanaecium (Bignonieae, Bignoniaceae) are lianas distributed in the Neotropics and centered in the Amazon. Members of the genus exhibit exceptionally diverse flower morphology and pollination systems. Here, we sequenced, assembled, and annotated 12 complete and four partial chloroplast genomes representing 15 Tanaecium species and more than 70% of the known diversity in the genus. Gene content and order were similar in all species of Tanaecium studied, with genome sizes ranging between 158,470 and 160,935 bp. Tanaecium chloroplast genomes have 137 genes, including 80-81 protein-coding genes, 37 tRNA genes, and four rRNA genes. No rearrangements were found in Tanaecium plastomes, but two different patterns of boundaries between regions were recovered. Tanaecium plastomes show nucleotide variability, although only rpoA was hypervariable. Multiple SSRs and repeat regions were detected, and eight genes were found to have signatures of positive selection. Phylogeny reconstruction using 15 Tanaecium plastomes resulted in a strongly supported topology, elucidating several relationships not recovered previously and bringing new insights into the evolution of the genus.

A pipeline for assembling low copy nuclear markers from plant genome skimming data for phylogenetic use

Background

Genome skimming is a popular method in plant phylogenomics that do not include a biased enrichment step, relying on random shallow sequencing of total genomic DNA. From these data the plastome is usually readily assembled and constitutes the bulk of phylogenetic information generated in these studies. Despite a few attempts to use genome skims to recover low copy nuclear loci for direct phylogenetic use, such endeavor remains neglected. Causes might include the trade-off between libraries with few reads and species with large genomes (i.e., missing data caused by low coverage), but also might relate to the lack of pipelines for data assembling.

Methods

A pipeline and its companion R package designed to automate the recovery of low copy nuclear markers from genome skimming libraries are presented. Additionally, a series of analyses aiming to evaluate the impact of key assembling parameters, reference selection and missing data are presented.

Results

A substantial amount of putative low copy nuclear loci was assembled and proved useful to base phylogenetic inference across the libraries tested (4 to 11 times more data than previously assembled plastomes from the same libraries).

Discussion

Critical aspects of assembling low copy nuclear markers from genome skims include the minimum coverage and depth of a sequence to be used. More stringent values of these parameters reduces the amount of assembled data and increases the relative amount of missing data, which can compromise phylogenetic inference, in turn relaxing the same parameters might increase sequence error. These issues are discussed in the text, and parameter tuning through multiple comparisons tracking their effects on support and congruence is highly recommended when using this pipeline. The skimmingLoci pipeline (https://github.com/mreginato/skimmingLoci) might stimulate the use of genome skims to recover nuclear loci for direct phylogenetic use, increasing the power of genome skimming data to resolve phylogenetic relationships, while reducing the amount of sequenced DNA that is commonly wasted.

Out of chaos: Phylogenomics of Asian Sonerileae

Sonerileae is a diverse Melastomataceae lineage comprising ca. 1000 species in 44 genera, with >70% of genera and species distributed in Asia. Asian Sonerileae are taxonomically intractable with obscure generic circumscriptions. The backbone phylogeny of this group remains poorly resolved, possibly due to complexity caused by rapid species radiation in early and middle Miocene, which hampers further systematic study. Here, we used genome resequencing data to reconstruct the phylogeny of Asian Sonerileae. Three parallel datasets, viz. single-copy ortholog (SCO), genomic SNPs, and whole plastome, were assembled from genome resequencing data of 205 species for this purpose. Based on these genome-scale data, we provided the first well resolved phylogeny of Asian Sonerileae, with 34 major clades identified and 74% of the interclade relationships consistently resolved by both SCO and genomic data. Meanwhile, widespread phylogenetic discordance was detected among SCO gene trees as well as species trees reconstructed using different tree estimation methods (concatenation/site-based coalescent method/summary method) or different datasets (SCO/genomic/plastome). We explored sources of discordance using multiple approaches and found that the observed discordance in Asian Sonerileae was mainly caused by a combination of biased distribution of missing data, random noise from uninformative genes, incomplete lineage sorting, and hybridization/introgression. Exploration of these sources can enable us to generate hypotheses for future testing, which is the first step towards understanding the evolution of Asian Sonerileae. We also detected high levels of homoplasy for some characters traditionally used in taxonomy, which explains current chaotic generic delimitations. The backbone phylogeny of Asian Sonerileae revealed in this study offers a solid basis for future taxonomic revision at the generic level.

Comprehensive analysis of complete chloroplast genome and phylogenetic aspects of ten Ficus species

BACKGROUND

The large genus Ficus comprises approximately 800 species, most of which possess high ornamental and ecological values. However, its evolutionary history remains largely unknown. Plastome (chloroplast genome) analysis had become an essential tool for species identification and for unveiling evolutionary relationships between species, genus and other rank groups. In this work we present the plastomes of ten Ficus species.

RESULTS

The complete chloroplast (CP) genomes of eleven Ficus specimens belonging to ten species were determined and analysed. The full length of the Ficus plastome was nearly 160 kbp with a similar overall GC content, ranging from 35.88 to 36.02%. A total of 114 unique genes, distributed in 80 protein-coding genes, 30 tRNAs, and 4 rRNAs, were annotated in each of the Ficus CP genome. In addition, these CP genomes showed variation in their inverted repeat regions (IR). Tandem repeats and mononucleotide simple sequence repeat (SSR) are widely distributed across the Ficus CP genome. Comparative genome analysis showed low sequence variability. In addition, eight variable regions to be used as potential molecular markers were proposed for future Ficus species identification. According to the phylogenetic analysis, these ten Ficus species were clustered together and further divided into three clades based on different subgenera. Simultaneously, it also showed the relatedness between Ficus and Morus.

CONCLUSION

The chloroplast genome structure of 10 Ficus species was similar to that of other angiosperms, with a typical four-part structure. Chloroplast genome sizes vary slightly due to expansion and contraction of the IR region. And the variation of noncoding regions of the chloroplast genome is larger than that of coding regions. Phylogenetic analysis showed that these eleven sampled CP genomes were divided into three clades, clustered with species from subgenus Urostigma, Sycomorus, and Ficus, respectively. These results support the Berg classification system, in which the subgenus Ficus was further decomposed into the subgenus Sycomorus. In general, the sequencing and analysis of Ficus plastomes, especially the ones of species with no or limited sequences available yet, contribute to the study of genetic diversity and species evolution of Ficus, while providing useful information for taxonomic and phylogenetic studies of Ficus.

Reference-free discovery of nuclear SNPs permits accurate, sensitive identification of Carya (hickory) species and hybrids

PREMISE

DNA-based species identification is critical when morphological identification is restricted, but DNA-based identification pipelines typically rely on the ability to compare homologous sequence data across species. Because many clades lack robust genomic resources, we present here a bioinformatics pipeline capable of generating genome-wide single-nucleotide polymorphism (SNP) data while circumventing the need for any reference genome or annotation data.

METHODS

Using the SISRS bioinformatics pipeline, we generated de novo ortholog data for the genus Carya, isolating sites where genetic variation was restricted to a single Carya species (i.e., species-informative SNPs). We leveraged these SNPs to identify both full-species and hybrid Carya specimens, even at very low sequencing depths.

RESULTS

We identified between 46,000 and 476,000 species-identifying SNPs for each of eight diploid Carya species, and all species identifications were concordant with the species of record. For all putative F₁ hybrid specimens, both parental species were correctly identified in all cases, and more punctate patterns of introgression were detectable in more cryptic crosses.

DISCUSSION

Bioinformatics pipelines that use only short-read sequencing data provide vital new tools enabling rapid expansion of DNA identification assays for model and non-model clades alike.

A nuclear phylogenomic study of the angiosperm order Myrtales, exploring the potential and limitations of the universal Angiosperms353 probe set

PREMISE

To further advance the understanding of the species-rich, economically and ecologically important angiosperm order Myrtales in the rosid clade, comprising nine families, approximately 400 genera and almost 14,000 species occurring on all continents (except Antarctica), we tested the Angiosperms353 probe kit.

METHODS

We combined high-throughput sequencing and target enrichment with the Angiosperms353 probe kit to evaluate a sample of 485 species across 305 genera (76% of all genera in the order).

RESULTS

Results provide the most comprehensive phylogenetic hypothesis for the order to date. Relationships at all ranks, such as the relationship of the early-diverging families, often reflect previous studies, but gene conflict is evident, and relationships previously found to be uncertain often remain so. Technical considerations for processing HTS data are also discussed.

CONCLUSIONS

High-throughput sequencing and the Angiosperms353 probe kit are powerful tools for phylogenomic analysis, but better understanding of the genetic data available is required to identify genes and gene trees that account for likely incomplete lineage sorting and/or hybridization events.

Comparative analyses of Mikania (Asteraceae: Eupatorieae) plastomes and impact of data partitioning and inference methods on phylogenetic relationships

We assembled new plastomes of 19 species of Mikania and of Ageratina fastigiata, Litothamnus nitidus, and Stevia collina, all belonging to tribe Eupatorieae (Asteraceae). We analyzed the structure and content of the assembled plastomes and used the newly generated sequences to infer phylogenetic relationships and study the effects of different data partitions and inference methods on the topologies. Most phylogenetic studies with plastomes ignore that processes like recombination and biparental inheritance can occur in this organelle, using the whole genome as a single locus. Our study sought to compare this approach with multispecies coalescent methods that assume that different parts of the genome evolve at different rates. We found that the overall gene content, structure, and orientation are very conserved in all plastomes of the studied species. As observed in other Asteraceae, the 22 plastomes assembled here contain two nested inversions in the LSC region. The plastomes show similar length and the same gene content. The two most variable regions within Mikania are rpl32-ndhF and rpl16-rps3, while the three genes with the highest percentage of variable sites are ycf1, rpoA, and psbT. We generated six phylogenetic trees using concatenated maximum likelihood and multispecies coalescent methods and three data partitions: coding and non-coding sequences and both combined. All trees strongly support that the sampled Mikania species form a monophyletic group, which is further subdivided into three clades. The internal relationships within each clade are sensitive to the data partitioning and inference methods employed. The trees resulting from concatenated analysis are more similar among each other than to the correspondent tree generated with the same data partition but a different method. The multispecies coalescent analysis indicate a high level of incongruence between species and gene trees. The lack of resolution and congruence among trees can be explained by the sparse sampling (~ 0.45% of the currently accepted species) and by the low number of informative characters present in the sequences. Our study sheds light into the impact of data partitioning and methods over phylogenetic resolution and brings relevant information for the study of Mikania diversity and evolution, as well as for the Asteraceae family as a whole.

Comparative analysis of chloroplast genome structure and molecular dating in Myrtales

BACKGROUND

Myrtales is a species rich branch of Rosidae, with many species having important economic, medicinal, and ornamental value. At present, although there are reports on the chloroplast structure of Myrtales, a comprehensive analysis of the chloroplast structure of Myrtales is lacking. Phylogenetic and divergence time estimates of Myrtales are mostly constructed by using chloroplast gene fragments, and the support for relationships is low. A more reliable method to reconstruct the species divergence time and phylogenetic relationships is by using whole chloroplast genomes. In this study, we comprehensively analyzed the structural characteristics of Myrtales chloroplasts, compared variation hotspots, and reconstructed the species differentiation time of Myrtales with four fossils and one secondary calibration point.

RESULTS

A total of 92 chloroplast sequences of Myrtales, representing six families, 16 subfamilies and 78 genera, were obtained including nine newly sequenced chloroplasts by whole genome sequencing. Structural analyses showed that the chloroplasts range in size between 152,214-171,315 bp and exhibit a typical four part structure. The IR region is between 23,901-36,747 bp, with the large single copy region spanning 83,691-91,249 bp and the small single copy region spanning 11,150-19,703 bp. In total, 123-133 genes are present in the chloroplasts including 77-81 protein coding genes, four rRNA genes and 30-31 tRNA genes. The GC content was 36.9-38.9%, with the average GC content being 37%. The GC content in the LSC, SSC and IR regions was 34.7-37.3%, 30.6-36.8% and 39.7-43.5%, respectively. By analyzing nucleotide polymorphism of the chloroplast, we propose 21 hypervariable regions as potential DNA barcode regions for Myrtales. Phylogenetic analyses showed that Myrtales and its corresponding families are monophyletic, with Combretaceae and the clade of Onagraceae + Lythraceae (BS = 100%, PP = 1) being sister groups. The results of molecular dating showed that the crown of Myrtales was most likely to be 104.90 Ma (95% HPD = 87.88-114.18 Ma), and differentiated from the Geraniales around 111.59 Ma (95% HPD = 95.50-118.62 Ma).

CONCLUSIONS

The chloroplast genome structure of Myrtales is similar to other angiosperms and has a typical four part structure. Due to the expansion and contraction of the IR region, the chloroplast genome sizes in this group are slightly different. The variation of noncoding regions of the chloroplast genome is larger than those of coding regions. Phylogenetic analysis showed that Combretaceae and Onagraceae + Lythraceae were well supported as sister groups. Molecular dating indicates that the Myrtales crown most likely originated during the Albian age of the Lower Cretaceous. These chloroplast genomes contribute to the study of genetic diversity and species evolution of Myrtales, while providing useful information for taxonomic and phylogenetic studies of Myrtales.

Tigridiopalma exalata, a new and endangered species of Melastomataceae from China

A new species of the genus Tigridiopalma, formerly considered monotypic, is here described as T. exalata and illustrated based on molecular and morphological evidence. It is morphologically similar to T. magnifica in having a short stem, huge basal leaves, scorpioid cymes, and 5-merous flowers, but differs in having ribbed and pale yellow puberulent petioles, purple petals with a small white apical patch, connectives of longer stamens with a distinct dorsal short spur at their base, and wingless capsules. Due to the restricted distribution, small populations and horticultural potential of this new species, it should be categorized as an Endangered species (EN).

Extremely low nucleotide diversity among thirty-six new chloroplast genome sequences from Aldama (Heliantheae, Asteraceae) and comparative chloroplast genomics analyses with closely related genera

Aldama (Heliantheae, Asteraceae) is a diverse genus in the sunflower family. To date, nearly 200 Asteraceae chloroplast genomes have been sequenced, but the plastomes of Aldama remain undescribed. Plastomes in Asteraceae usually show little sequence divergence, consequently, our hypothesis is that species of Aldama will be overall conserved. In this study, we newly sequenced 36 plastomes of Aldama and of five species belonging to other Heliantheae genera selected as outgroups (i.e., Dimerostemma asperatum, Helianthus tuberosus, Iostephane heterophylla, Pappobolus lanatus var. lanatus, and Tithonia diversifolia). We analyzed the structure and gene content of the assembled plastomes and performed comparative analyses within Aldama and with other closely related genera. As expected, Aldama plastomes are very conserved, with the overall gene content and orientation being similar in all studied species. The length of the plastome is also consistent and the junction between regions usually contain the same genes and have similar lengths. A large ∼20 kb and a small ∼3 kb inversion were detected in the Large Single Copy (LSC) regions of all assembled plastomes, similarly to other Asteraceae species. The nucleotide diversity is very low, with only 1,509 variable sites in 127,466 bp (i.e., 1.18% of the sites in the alignment of 36 Aldama plastomes, with one of the IRs removed, is variable). Only one gene, rbcL, shows signatures of positive selection. The plastomes of the selected outgroups feature a similar gene content and structure compared to Aldama and also present the two inversions in the LSC region. Deletions of different lengths were observed in the gene ycf2. Multiple SSRs were identified for the sequenced Aldama and outgroups. The phylogenetic analysis shows that Aldama is not monophyletic due to the position of the Mexican species A. dentata. All Brazilian species form a strongly supported clade. Our results bring new understandings into the evolution and diversity of plastomes at the species level.

Phylogeny and biogeography of Myrcia sect. Aguava (Myrtaceae, Myrteae) based on phylogenomic and Sanger data provide evidence for a Cerrado origin and geographically structured clades

Myrcia is one of the largest exclusively Neotropical angiosperm genera, including ca. 800 species divided into nine sections. Myrcia sect. Aguava is one of most complex sections of Myrcia due to high morphological variation and wide distribution range of some species, including M. guianensis, with distribution throughout South America and a complex taxonomic history. We used complete plastid DNA sequences data generated using next-generation sequencing of 45 terminals, mostly from Myrcia sect. Aguava. These data were combined with five target DNA regions (ITS, psbA-trnH, trnL-trnF, trnQ-rps16, ndhF) of additional terminals to increase taxonomic coverage. Phylogenetic analyses were conducted using a maximum likelihood approach, and divergence times and ancestral range distributions were estimated. Myrcia sect. Aguava is monophyletic and exclusively comprises species with trilocular ovaries but has no relationship with other groups within Myrcia that possess trilocular ovaries. Three main lineages that correspond to geographical distribution are recognized within Myrcia sect. Aguava. Multiple accessions reveal a non-monophyletic Myrcia guianensis and stress the biogeographical structure inside the group. Myrcia sect. Aguava had a probable mid-Miocene origin in the Cerrado, but lineages that persisted there diversified only more recently, when the present-day vegetation started to stabilize. Posterior migrations to Atlantic Forest, Amazon and Caribbean occurred at the end of Miocene, evidencing transitions from open and dry to forested and more humid areas that are less frequent in the Neotropics. Overall, it is observed that related lineages remained in ecologically similar environments. Future perspectives on Myrcia and Myrteae in the phylogenomic era are also discussed.

Insights Into Chloroplast Genome Evolution Across Opuntioideae (Cactaceae) Reveals Robust Yet Sometimes Conflicting Phylogenetic Topologies

Chloroplast genomes (plastomes) are frequently treated as highly conserved among land plants. However, many lineages of vascular plants have experienced extensive structural rearrangements, including inversions and modifications to the size and content of genes. Cacti are one of these lineages, containing the smallest plastome known for an obligately photosynthetic angiosperm, including the loss of one copy of the inverted repeat (∼25 kb) and the ndh gene suite, but only a few cacti from the subfamily Cactoideae have been sufficiently characterized. Here, we investigated the variation of plastome sequences across the second-major lineage of the Cactaceae, the subfamily Opuntioideae, to address (1) how variable is the content and arrangement of chloroplast genome sequences across the subfamily, and (2) how phylogenetically informative are the plastome sequences for resolving major relationships among the clades of Opuntioideae. Our de novo assembly of the Opuntia quimilo plastome recovered an organelle of 150,347 bp in length with both copies of the inverted repeat and the presence of all the ndh gene suite. An expansion of the large single copy unit and a reduction of the small single copy unit was observed, including translocations and inversion of genes, as well as the putative pseudogenization of some loci. Comparative analyses among all clades within Opuntioideae suggested that plastome structure and content vary across taxa of this subfamily, with putative independent losses of the ndh gene suite and pseudogenization of genes across disparate lineages, further demonstrating the dynamic nature of plastomes in Cactaceae. Our plastome dataset was robust in resolving three tribes with high support within Opuntioideae: Cylindropuntieae, Tephrocacteae and Opuntieae. However, conflicting topologies were recovered among major clades when exploring different assemblies of markers. A plastome-wide survey for highly informative phylogenetic markers revealed previously unused regions for future use in Sanger-based studies, presenting a valuable dataset with primers designed for continued evolutionary studies across Cactaceae. These results bring new insights into the evolution of plastomes in cacti, suggesting that further analyses should be carried out to address how ecological drivers, physiological constraints and morphological traits of cacti may be related with the common rearrangements in plastomes that have been reported across the family.

Is dispersal mode a driver of diversification and geographical distribution in the tropical plant family Melastomataceae?

Species of plants with different life history strategies may differ in their seed dispersal mechanisms, impacting their distribution and diversification patterns. Shorter or longer distance dispersal is favored by different dispersal modes, facilitating (or constraining) population isolation, which can, in turn, impact speciation and species range sizes. While these associations are intuitive, few studies have explicitly tested these hypotheses for large clades of angiosperms. The plant family Melastomataceae is found on disparate habitats with different dispersal modes, representing a good model to address these questions. In this study, we reconstruct the phylogeny of Melastomataceae and gather data on their dispersal mode and range size to test the impact of dispersal mode on diversification and range size evolution. We found that abiotic dispersal is ancestral in the family, while biotic dispersal evolved multiple times. Species richness distribution is very similar across dispersal modes, although abiotically dispersed species tend to be relatively more diverse in seasonal environments. Range sizes across dispersal modes are not significantly different, although biotically dispersed species have slightly wider distributions. Model comparisons indicate that factors other than dispersal mode might have driven diversification heterogeneity. We did not find evidence for the role of dispersal mode driving diversification rates or range size in the Melastomataceae, suggesting a complex macroevolutionary scenario for this diverse angiosperm family. The bulk of changes to biotic dispersal coinciding with an increase in passerine diversification suggests a possible "past" key innovation in Melastomataceae. Future studies should investigate the role of other diversification drivers in the family and the relatively higher diversity of abiotically dispersed species in open habitats.

Characteristics analysis of the complete Wurfbainia villosa chloroplast genome

Wurfbainia villosa, which belongs to the huge family Zingiberaceae, is used in the clinic for the treatment of spleen and stomach diseases in southern China. The complete chloroplast genome of W. villosa was sequenced and analyzed using next-generation sequencing technology in the present work. The results showed that the W. villosa chloroplast genome is a circular molecule with 163,608 bp in length. It harbors a pair of inverted repeat regions (IRa and IRb) of 29,820 bp in length, which separate the large single copy (LSC, 88,680 bp) region and the small single copy (SSC, 15,288 bp) region. After annotation, 134 genes were identified in this plastome in total, comprising of 87 protein-coding genes, 38 transfer RNA genes, 8 ribosomal RNA genes and one pseudogene (ycf1). Codon usage, RNA editing sites and single/long sequence repeats were investigated to understand the structural characteristics of the W. villosa chloroplast genome. Furthermore, IR contraction and expansion were analyzed by comparison of complete chloroplast genomes of W. villosa and four other Zingiberaceae species. Finally, a phylogeny study based on the chloroplast genome of W. villosa, along with that of 15 different species, was conducted to further investigate the relationship among these lineages. Overally, our results represented the first insight into the chloroplast genome of W. villosa, and could serve as a significant reference for species identification, genetic diversity analysis and phylogenetic research between W. villosa and other species within Zingiberaceae.

Analyses of Plastome Sequences Improve Phylogenetic Resolution and Provide New Insight Into the Evolutionary History of Asian Sonerileae/Dissochaeteae

Sonerileae/Dissochaeteae (Melastomataceae) comprises ca. 50 genera, two thirds of which occur in Southeast Asia. Phylogenetic relationships within this clade remain largely unclear, which hampers our understanding of its origin, evolution, and biogeography. Here, we explored the use of chloroplast genomes in phylogenetic reconstruction of Sonerileae/Dissochaeteae, by sampling 138 species and 23 genera in this clade. A total of 151 complete plastid genomes were assembled for this study. Plastid genomic data provided better support for the backbone of the Sonerileae/Dissochaeteae phylogeny, and also for relationships among most closely related species, but failed to resolve the short internodes likely resulted from rapid radiation. Trees inferred from plastid genome and nrITS sequences were largely congruent regarding the major lineages of Sonerileae/Dissochaeteae. The present analyses recovered 15 major lineages well recognized in both nrITS and plastid phylogeny. Molecular dating and biogeographical analyses indicated a South American origin for Sonerileae/Dissochaeteae during late Eocene (stem age: 34.78 Mya). Two dispersal events from South America to the Old World were detected in late Eocene (33.96 Mya) and Mid Oligocene (28.33 Mya) respectively. The core Asian clade began to diversify around early Miocene in Indo-Burma and dispersed subsequently to Malesia and Sino-Japanese regions, possibly promoted by global temperature changes and East Asian monsoon activity. Our analyses supported previous hypothesis that Medinilla reached Madagascar by transoceanic dispersal in Miocene. In addition, generic limits of some genera concerned were discussed.

The complete chloroplast genome sequence of Blastus pauciflorus (Melastomataceae)

Blastus pauciflorus, a shrub endemic to Hong Kong and Guangdong, south China, growing on low-altitude hillsides, under the forest. The species is controversial in classification. Herein, we report the complete chloroplast genome sequence assembled from Illumina pair-end sequencing data, with aims to resolve its relationship with the related species. The complete chloroplast genome was 155,983 bp in length, includes two inverted repeat regions (IRs) of 26,716 bp each, which were separated by a large single copy region (LSC) 86,101 bp and a small single copy region (SSC) 16,450 bp. The chloroplast genome contained 129 genes, including 82 protein-coding genes, 2 pseudogenes, 37 tRNA genes and 8 rRNA genes. The overall GC content in the chloroplast genome of B. pauciflorus was 37.0%. Phylogenetic analysis showed that B. pauciflorus is closed to B. cochinchinensis.

The complete chloroplast genome sequence of monotypic Cyphotheca (Melastomataceae), an endemic genus in China

Cyphotheca, a shrub narrowly endemic to Yunnan, southwest China, is a monotypic genus of Sonerileae in Melastomataceae. The species, Cyphotheca montana, is dwindling and is listed as 'Near Threatened (NT)' by IUCN Red List of China. Herein, we report the complete chloroplast genome (cpDNA) sequence assembled from Illumina pair-end sequencing data. The complete chloroplast genome was 156,424 bp in length, which includes two inverted repeat regions (IRs) of 26,776 bp each, which were separated by a large single copy region (LSC) 85,900 bp and a small single copy region (SSC) 16,972 bp. The chloroplast genome contained 130 genes, including 85 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. The overall GC content in the chloroplast genome of C. montana was 37.0%. The phylogenomic analysis showed that C. montana is sister to Barthea barthei.

Characterization of the complete chloroplast genome sequence of Blastus cochinchinensis (Melastomataceae)

Blastus Lour. belongs to tribe Sonerileae (Melastomataceae), comprising about 18 species worldwide. Herein, we presented, assembled, and annotated the first complete chloroplast genome of Blastus (B. cochinchinensis). The complete chloroplast genome (cp) size of B. cochinchinensis was 156,005 bp in length, containing a pair of 26,812 bp inverted repeat (IR) regions, which were separated by a large single-copy region (LSC) 85,926 bp and a small single-copy (SSC) region 16,455 bp. The overall GC content of the genome was 37.0%. The whole genome contained 129 unique genes, including 81 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. The phylogenetic analysis showed that B. cochinchinens was sister to Barthea barthei. The present study provided B. cochinchinensis genomic information that may be found useful in conservation and molecular phylogenetic studies on Blastus.

Structure and features of the complete chloroplast genome of Melastoma dodecandrum

Melastoma dodecandrum, the only creeping species in the Melastoma genus, serves as a medicinal herb in southeast China. It belongs to the huge family Melastomataceae, which contains over 5000 species worldwide. In this study, we used next-generation sequencing to determine the complete chloroplast genome sequences of M. dodecandrum, which is a circular molecule of 156,611 bp in length. After annotation, we identified 131 putative genes in total, comprised of 85 protein-coding genes, 38 transfer RNA genes and 8 ribosomal RNA genes. Genome structure, GC content, repeat sequences and codon usage were investigated to gain a comprehensive understanding of this genome. Furthermore, we conducted comparative genome analyses between the M. dodecandrum genome and that of four other Melastomataceae species. Additionally, a phylogenetic analysis was performed based on available chloroplast genomes of Melastomataceae species and several Myrtaceae species, revealing the taxonomic relationships between M. dodecandrum and related species. In conclusion, our study represents the first look into the complete chloroplast genome of M. dodecandrum, providing abundant information for further studies such as species identification, taxonomy and phylogenetic resolution of Melastomataceae species.

Comparative analyses of chloroplast genomes from 22 Lythraceae species: inferences for phylogenetic relationships and genome evolution within Myrtales

BACKGROUND

Lythraceae belongs to the order Myrtales, which is part of Archichlamydeae. The family has 31 genera containing approximately 620 species of herbs, shrubs and trees. Of these 31 genera, five large genera each possess 35 or more species. They are Lythrum, with 35; Rotala, with 45; Nesaea, with 50; Lagerstroemia, with 56; and Cuphea, with 275 species.

RESULTS

We reported six newly sequenced chloroplast (cp) genomes (Duabanga grandiflora, Trapa natans, Lythrum salicaria, Lawsonia inermis, Woodfordia fruticosa and Rotala rotundifolia) and compared them with 16 other cp genomes of Lythraceae species. The cp genomes of the 22 Lythraceae species ranged in length from 152,049 bp to 160,769 bp. In each Lythraceae species, the cp genome contained 112 genes consisting of 78 protein coding genes, four ribosomal RNAs and 30 transfer RNAs. Furthermore, we detected 211-332 simple sequence repeats (SSRs) in six categories and 7-27 long repeats in four categories. We selected ten divergent hotspots (ndhF, matK, ycf1, rpl22, rpl32, trnK-rps16, trnR-atpA, rpl32-trnL, trnH-psbA and trnG-trnR) among the 22 Lythraceae species to be potential molecular markers. We constructed phylogenetic trees from 42 Myrtales plants with 8 Geraniales plants as out groups. The relationships among the Myrtales species were effectively distinguished by maximum likelihood (ML), maximum parsimony (MP) and Bayesian inference (BI) trees constructed using 66 protein coding genes. Generally, the 22 Lythraceae species gathered into one clade, which was resolved as sister to the three Onagraceae species. Compared with Melastomataceae and Myrtaceae, Lythraceae and Onagraceae differentiated later within Myrtales.

CONCLUSIONS

The study provided ten potential molecular markers as candidate DNA barcodes and contributed cp genome resources within Myrtales for further study.

De novo Assembly and Characterization of the Floral Transcriptomes of Two Varieties of Melastoma malabathricum

Melastoma malabathricum is an important medicinal and landscape plant that is globally distributed in temperate and subtropical regions. However, available genomic information for the entire Melastomataceae family is notably limited. In view of the application potential of floral parts in secondary metabolite extraction, we characterized for the first time the floral transcriptomes of two key M. malabathricum varieties, purple variety and white variety. Our transcriptome assembly generated 52,498 and 49,380 unigenes with an N50 of 1,906 and 1,929 bases for the purple and white varieties, respectively. Comparative analysis of two transcriptomes demonstrated that they are highly similar but also highlighted genes that are presumably lineage specific, which explains the phenotypes of each variety. Additionally, a shared transcriptional signature across the floral developmental stages was identified in both M. malabathricum varieties; this signature included pathways related to secondary metabolite synthesis, plant hormone signaling and production, energy homeostasis and nutrient assimilation pathways, and cellular proliferation. The expression levels of flavonoid accumulation and candidate flavonoid biosynthesis-related genes in M. malabathricum flower development stages validated the transcriptome findings. The transcriptome data presented in this study will serve as a valuable resource for future work on the exploitation of M. malabathricum and other related species. The gene expression dynamics during flower development will facilitate the discovery of lineage-specific genes associated with phenotypic characteristics and will elucidate the mechanism of the ontogeny of individual flower types.

Comprehensive Analysis of Rhodomyrtus tomentosa Chloroplast Genome

In the last decade, several studies have relied on a small number of plastid genomes to deduce deep phylogenetic relationships in the species-rich Myrtaceae. Nevertheless, the plastome of Rhodomyrtus tomentosa, an important representative plant of the Rhodomyrtus (DC.) genera, has not yet been reported yet. Here, we sequenced and analyzed the complete chloroplast (CP) genome of R. tomentosa, which is a 156,129-bp-long circular molecule with 37.1% GC content. This CP genome displays a typical quadripartite structure with two inverted repeats (IRa and IRb), of 25,824 bp each, that are separated by a small single copy region (SSC, 18,183 bp) and one large single copy region (LSC, 86,298 bp). The CP genome encodes 129 genes, including 84 protein-coding genes, 37 tRNA genes, eight rRNA genes and three pseudogenes (ycf1, rps19, ndhF). A considerable number of protein-coding genes have a universal ATG start codon, except for psbL and ndhD. Premature termination codons (PTCs) were found in one protein-coding gene, namely atpE, which is rarely reported in the CP genome of plants. Phylogenetic analysis revealed that R. tomentosa has a sister relationship with Eugenia uniflora and Psidium guajava. In conclusion, this study identified unique characteristics of the R. tomentosa CP genome providing valuable information for further investigations on species identification and the phylogenetic evolution between R. tomentosa and related species.

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.

Comparative Chloroplast Genomics at Low Taxonomic Levels: A Case Study Using Amphilophium (Bignonieae, Bignoniaceae)

Chloroplast (cp) genome organization, gene order, and content have long been considered conserved among land plants. Despite that, the generation of thousands of complete plastomes through next-generation sequencing (NGS) has challenged their conserved nature. In this study, we analyze 11 new complete plastomes of Amphilophium (Bignonieae, Bignoniaceae), a diverse genus of Neotropical lianas, and that of Anemopaegma prostratum. We explored the structure and content of the assembled plastomes and performed comparative analyses within Amphilophium and among other plastomes available for Bignoniaceae. The overall gene content and orientation of plastomes is similar in all species studied. Plastomes are not conserved among Amphilophium, showing significant differences in length (155,262-164,786 bp), number of genes duplicated in the IRs (eight, 18, or 19), and location of the SC/IR boundaries (i.e., LSC/IRa junction between rps19 and rpl2 genes, within petD, or within petB). Length differences reflect expansions of the IRs and contractions of the LSC regions. The plastome of A. prostratum is 168,172 bp, includes 19 duplicated genes, and has the LSC/IRa boundary located within the petB gene. Amphilophium plastomes show high nucleotide diversity, with many hypervariable regions, and 16 genes with signatures of positive selection. Multiple SSRs and repeat regions were identified for Amphilophium and Anemopaegma prostratum. The differences in structure detected within Amphilophium plastomes in terms of LSC/IR and IR/SSC boundaries, number of duplicated genes, and genome sizes are mostly shared between taxa that belong to the same clade. Our results bring new insights into the evolution of plastomes at low taxonomic levels.

Contrasting patterns of taxonomic, phylogenetic and functional variation along a Costa Rican altitudinal gradient in the plant family Melastomataceae

The functional composition of plant communities in montane regions has been studied for decades, and most recent analyses find that environmentally favourable landscapes at lower altitudes tend to be dominated by species with resource-acquisitive traits, while more resource-conservative taxa dominate higher-altitude communities. However, it is unclear the extent to which this pattern is driven by co-gradient variation within clades or changes in clade representation across the gradient. To test for co-gradient variation, species composition, phylogenetic structure and functional traits were quantified for 97 species within the plant family Melastomataceae at five locations across a 2500-m altitudinal gradient along Volcán Barva in Costa Rica. Average melastome leaf force to punch, specific leaf area and leaf size vary with altitude, while four other functional traits do not. Taxonomic dissimilarity between communities was correlated with altitudinal difference, while phylogenetic dissimilarity was correlated with altitudinal dissimilarity only when measured with a metric that emphasizes shallow turnover of the tips of the phylogeny. These results highlight how species turnover may be more pronounced than functional or phylogenetic variation along altitudinal gradients. In addition, these results highlight the conservation value of lowland tropical forests, which here harbour a disproportionate amount of phylogenetic and functional diversity.

Comparative analyses of plastid genomes from fourteen Cornales species: inferences for phylogenetic relationships and genome evolution

Background

The Cornales is the basal lineage of the asterids, the largest angiosperm clade. Phylogenetic relationships within the order were previously not fully resolved. Fifteen plastid genomes representing 14 species, ten genera and seven families of Cornales were newly sequenced for comparative analyses of genome features, evolution, and phylogenomics based on different partitioning schemes and filtering strategies.

Results

All plastomes of the 14 Cornales species had the typical quadripartite structure with a genome size ranging from 156,567 bp to 158,715 bp, which included two inverted repeats (25,859–26,451 bp) separated by a large single-copy region (86,089–87,835 bp) and a small single-copy region (18,250–18,856 bp) region. These plastomes encoded the same set of 114 unique genes including 31 transfer RNA, 4 ribosomal RNA and 79 coding genes, with an identical gene order across all examined Cornales species. Two genes (rpl22 and ycf15) contained premature stop codons in seven and five species respectively. The phylogenetic relationships among all sampled species were fully resolved with maximum support. Different filtering strategies (none, light and strict) of sequence alignment did not have an effect on these relationships. The topology recovered from coding and noncoding data sets was the same as for the whole plastome, regardless of filtering strategy. Moreover, mutational hotspots and highly informative regions were identified.

Conclusions

Phylogenetic relationships among families and intergeneric relationships within family of Cornales were well resolved. Different filtering strategies and partitioning schemes do not influence the relationships. Plastid genomes have great potential to resolve deep phylogenetic relationships of plants.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4319-9) contains supplementary material, which is available to authorized users.

The complete chloroplast genome sequence of Barthea barthei (Melastomataceae), a shrub endemic to southern China

The plant genus Barthea is monotypic, comprising of only a single species (Barthea barthei), and is endemic to southern China. In this study, we report the complete chloroplast genome of B. barthei, assembled from whole-genome high-throughput sequencing data, as a resource for future studies on the taxonomy and evolution of Barthea. The chloroplast genome was 155,951 bp in length, with a large single-copy (LSC) region of 85,882 bp, a small single-copy (SSC) region of 16,445 bp, separated by two inverted repeat (IR) regions of 26,812 bp each. It was predicted to contain a total of 130 genes, with an overall GC content of 36.99%. Phylogenetic analysis placed B. barthei closest to Opisthocentra sp. in Melastomataceae.

The complete chloroplast genome sequence of Melastoma candidum (Melastomataceae)

The plant genus Melastoma is comprised of members estimated to have formed through recent species radiation. Natural hybridization among member species further complicates taxonomy within the genus. Herein, we report the complete chloroplast genome of M. candidum, assembled from partial data obtained from a parallel whole-genome Illumina paired-end sequencing effort on the species. The chloroplast genome was 156,682 bp in length, with a large single-copy (LSC) region of 86,084 bp, a small single-copy (SSC) region of 17,094 bp, separated by two inverted repeat (IR) regions of 26,752 bp each. It was predicted to contain a total of 129 genes, with an overall GC content of 37.17%. Phylogenetic analysis placed M. candidum in the same clade as species within the Melastomeae tribe of Melastomataceae.

Plastome Rearrangements in the "Adenocalymma-Neojobertia" Clade (Bignonieae, Bignoniaceae) and Its Phylogenetic Implications

The chloroplast is one of the most important organelles of plants. This organelle has a circular DNA with approximately 130 genes. The use of plastid genomic data in phylogenetic and evolutionary studies became possible with high-throughput sequencing methods, which allowed us to rapidly obtain complete genomes at a reasonable cost. Here, we use high-throughput sequencing to study the "Adenocalymma-Neojobertia" clade (Bignonieae, Bignoniaceae). More specifically, we use Hi-Seq Illumina technology to sequence 10 complete plastid genomes. Plastomes were assembled using selected plastid reads and de novo approach with SPAdes. The 10 assembled genomes were analyzed in a phylogenetic context using five different partition schemes: (1) 91 protein-coding genes ("coding"); (2) 76 introns and spacers with alignment manually edited ("non-coding edited"); (3) 76 non-coding regions with poorly aligned regions removed using T-Coffee ("non-coding filtered"); (4) 91 coding regions plus 76 non-coding regions edited ("coding + non-coding edited"); and, (5) 91 protein-coding regions plus the 76 filtered non-coding regions ("coding + non-coding filtered"). Fragmented regions were aligned using Mafft. Phylogenetic analyses were conducted using Maximum Likelihood (ML) and Bayesian Criteria (BC). The analyses of the individual plastomes consistently recovered an expansion of the Inverted Repeated (IRs) regions and a compression of the Small Single Copy (SSC) region. Major genomic translocations were observed at the Large Single Copy (LSC) and IRs. ML phylogenetic analyses of the individual datasets led to the same topology, with the exception of the analysis of the "non-coding filtered" dataset. Overall, relationships were strongly supported, with the highest support values obtained through the analysis of the "coding + non-coding edited" dataset. Four regions at the LSC, SSC, and IR were selected for primer development. The "Adenocalymma-Neojobertia" clade shows an unusual pattern of plastid structure variation, including four major genomic translocations. These rearrangements challenge the current view of conserved plastid genome architecture in terms of gene order. It also complicates both genomic assemblies using reference genomes and sequence alignments using whole plastomes. Therefore, strategies that employ de novo assemblies and manual evaluation of sequence alignments are required to prevent assembly and alignment errors.