Exceptionally high rates of positive selection on the rbcL gene in the genus Ilex (Aquifoliaceae)

Chloroplast genome structure analysis of Equisetum unveils phylogenetic relationships to ferns and mutational hotspot region

Equisetum is one of the oldest extant group vascular plants and is considered to be the key to understanding vascular plant evolution. Equisetum is distributed almost all over the world and has a high degree of adaptability to different environments. Despite the fossil record of horsetails (Equisetum, Equisetaceae) dating back to the Carboniferous, the phylogenetic relationship of this genus is not well, and the chloroplast evolution in Equisetum remains poorly understood. In order to fill this gap, we sequenced, assembled, and annotated the chloroplast genomes of 12 species of Equisetum, and compared them to 13 previously published vascular plants chloroplast genomes to deeply examine the plastome evolutionary dynamics of Equisetum. The chloroplast genomes have a highly conserved quadripartite structure across the genus, but these chloroplast genomes have a lower GC content than other ferns. The size of Equisetum plastomes ranges from 130,773 bp to 133,684 bp and they encode 130 genes. Contraction/expansion of IR regions and the number of simple sequences repeat regions underlie large genomic variations in size among them. Comparative analysis revealed we also identified 13 divergence hotspot regions. Additionally, the genes accD and ycf1 can be used as potential DNA barcodes for the identification and phylogeny of the genus Equisetum. Twelve photosynthesis-related genes were specifically selected in Equisetum. Comparative genomic analyses implied divergent evolutionary patterns between Equisetum and other ferns. Phylogenomic analyses and molecular dating revealed a relatively distant phylogenetic relationship between Equisetum and other ferns, supporting the division of pteridophyte into Lycophytes, Equisetaceae and ferns. The results show that the chloroplast genome can be used to solve phylogenetic problems within or between Equisetum species, and also provide genomic resources for the study of Equisetum systematics and evolution.

Plastome evolution in the East Asian lobelias (Lobelioideae) using phylogenomic and comparative analyses

Lobelia species, as rich source of the alkaloid lobeline which has been shown to have important biological activity, have been used in folk medicine throughout East Asia to treat various diseases. However, Lobelia is a complex and varied genus in East Asia and is thus difficult to identify. Genomic resources would aid identification, however the availability of such information is poor, preventing a clear understanding of their evolutionary history from being established. To close this gap in the available genomic data, in this study, 17 plastomes of East Asian lobelias were newly sequenced and assembled. Although the plastomes of Lobelia sect. Hypsela, L. sect. Speirema, and L. sect. Rhynchopetalum shared the gene structure, the inverted repeat (IR)/large single copy (LSC) boundaries, genome size, and the number of repeats were variable, indicating the non-conservative nature of plastome evolution within these sections. However, the genomes of the Lobelia sect. Delostemon and L. sect. Stenotium showed rearrangements, revealing that these two sections might have undergone different evolutionary histories. We assessed nine hotspot genes and 27-51 simple sequence repeat motifs, which will also serve as valuable DNA barcode regions in future population genetics studies and for the delineation of plant species. Our phylogenetic analysis resolved the evolutionary positions of the five sections in agreement with previous evolutionary trees based on morphological features. Although phylogenetic reconstruction of Lobelioideae based on the rpoC2 gene has rarely been performed, our results indicated that it contains a considerable amount of phylogenetic information and offers great promise for further phylogenetic analysis of Lobelioideae. Our site-specific model identified 173 sites under highly positive selections. The branch-site model exhibited 11 positive selection sites involving four genes in the East Asian branches. These four genes may play critical roles in the adaptation of East Asian taxa to diverse environments. Our study is the first to detect plastome organization, phylogenetic utility, and signatures of positive selection in the plastomes of East Asian lobelias, which will help to further advance taxonomic and evolutionary studies and the utilization of medicinal plant resources.

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.

Assessing Phenotypic Variability in Some Eastern European Insular Populations of the Climatic Relict Ilex aquifolium L

Through its natural or cultivated insular population distribution, Ilex aquifolium L. is a paramount species which is exceptionally suitable for studying phenotypic variability and plasticity through the assessment of morphological, physiological, biochemical and genomic features with respect to acclimation and/or adaptation efficiency. The current study is focused on four insular populations of Ilex aquifolium from Eastern Europe (i.e., in Romania, Hungary, Serbia and Bulgaria), and presents an initial evaluation of phenotypic variability in order to conclude our research on phylogenetic relationships and phytochemical profiles, including several descriptive and quantitative morphological traits. Taken together, the data from different methods in this paper indicate that the Bulgarian and Romanian populations can be distinguished from each other and from Serbian and Hungarian populations, while the latter show a higher level of resemblance with regards to their quantitative morphological traits. It is likely that these morphological traits are determined through some quantitative trait loci implicated in stress responses generated by light, temperature, soil water, soil fertility and salinity conditions that will need to be analysed in terms of their physiological, genomic and metabolomics traits in future studies.

Plastid genome evolution of a monophyletic group in the subtribe Lauriineae (Laureae, Lauraceae)

Litsea, a non-monophyletic group of the tribe Laureae (Lauraceae), plays important roles in the tropical and subtropical forests of Asia, Australia, Central and North America, and the islands of the Pacific. However, intergeneric relationships between Litsea and Laurus, Lindera, Parasassafras and Sinosassafras of the tribe Laureae remain unresolved. In this study, we present phylogenetic analyses of seven newly sequenced Litsea plastomes, together with 47 Laureae plastomes obtained from public databases, representing six genera of the Laureae. Our results highlight two highly supported monophyletic groups of Litsea taxa. One is composed of 16 Litsea taxa and two Lindera taxa. The 18 plastomes of these taxa were further compared for their gene structure, codon usage, contraction and expansion of inverted repeats, sequence repeats, divergence hotspots, and gene evolution. The complete plastome size of newly sequenced taxa varied between 152,377 bp (Litsea auriculata) and 154,117 bp (Litsea pierrei). Seven of the 16 Litsea plastomes have a pair of insertions in the IRa (trnL-trnH) and IRb (ycf2) regions. The 18 plastomes of Litsea and Lindera taxa exhibit similar gene features, codon usage, oligonucleotide repeats, and inverted repeat dynamics. The codons with the highest frequency among these taxa favored A/T endings and each of these plastomes had nine divergence hotspots, which are located in the same regions. We also identified six protein coding genes (accD, ndhJ, rbcL, rpoC2, ycf1 and ycf2) under positive selection in Litsea; these genes may play important roles in adaptation of Litsea species to various environments.

Molecular diversity of exotic durian (Durio spp.) germplasm: a case study of Kalimantan, Indonesia

Background

Durian of Indonesia, specifically Durio zibethinus, is a potential agricultural commodity for domestic and international markets. However, its quality is still less competitive or significantly lower to fulfill the export market, compared to a similar one from other countries. This study aimed to determine and analyze the genetic diversity and relationship of the exotic durian (Durio spp.) germplasm originally from Kalimantan, Indonesia, using the rbcL marker.

Results

Based on this marker, the durian germplasm has a low genetic diversity (π%=0.24). It may strongly correspond with the variability sites or mutation present in the region. In this case, the rbcL region of the durian germplasm has generated 23 variable sites with a transition/transversion (Ti/Tv) bias value of 1.00. However, following the phylogenetic and principal component analyses, this germplasm is separated into four main clades and six groups, respectively. In this case, D. zibethinus was very closely related to D. exleyanus. Meanwhile, D. lowianus and D. excelsus were the farthest. In further analysis, 29 durians were very closely related, and the farthest was shown by Durian Burung (D. acutifolius) and Kalih Haliyang (D. kutejensis) as well as Pampaken Burung Kecil (D. kutejensis) and Durian Burung (D. acutifolius) with a divergence coefficient of 0.011. The Pearson correlation analysis confirms that 20 pairs of individual durians have a strong relation, shown by, e.g., Maharawin Hamak and Durian Burung as well as Mantuala Batu Hayam and Durian Burung Besar.

Conclusion

While the durian has a low genetic diversity, the phylogenetic analyses revealed that this germplasm originally from Kalimantan, Indonesia, shows unique relationships. These findings may provide a beneficial task in supporting the durian genetic conservation and breeding practices in the future, locally and globally.

Evolution of Rosaceae Plastomes Highlights Unique Cerasus Diversification and Independent Origins of Fruiting Cherry

Rosaceae comprises numerous types of economically important fruits, ornamentals, and timber. The lack of plastome characteristics has blocked our understanding of the evolution of plastome and plastid genes of Rosaceae crops. Using comparative genomics and phylogenomics, we analyzed 121 Rosaceae plastomes of 54 taxa from 13 genera, predominantly including Cerasus (true cherry) and its relatives. To our knowledge, we generated the first comprehensive map of genomic variation across Rosaceae plastomes. Contraction/expansion of inverted repeat regions and sequence losses of the two single-copy regions underlie large genomic variations in size among Rosaceae plastomes. Plastid protein-coding genes were characterized with a high proportion (over 50%) of synonymous variants and insertion-deletions with multiple triplets. Five photosynthesis-related genes were specially selected in perennial woody trees. Comparative genomic analyses implied divergent evolutionary patterns between pomaceous and drupaceous trees. Across all examined plastomes, unique and divergent evolution was detected in Cerasus plastomes. Phylogenomic analyses and molecular dating highlighted the relatively distant phylogenetic relationship between Cerasus and relatives (Microcerasus, Amygdalus, Prunus, and Armeniaca), which strongly supported treating the monophyletic true cherry group as a separate genus excluding dwarf cherry. High genetic differentiation and distinct phylogenetic relationships implied independent origins and domestication between fruiting cherries, particularly between Prunus pseudocerasus (Cerasus pseudocerasus) and P. avium (C. avium). Well-resolved maternal phylogeny suggested that cultivated P. pseudocerasus originated from Longmenshan Fault zone, the eastern edge of Himalaya-Hengduan Mountains, where it was subjected to frequent genomic introgression between its presumed wild ancestors and relatives.

Gene Loss, Pseudogenization in Plastomes of Genus Allium (Amaryllidaceae), and Putative Selection for Adaptation to Environmental Conditions

Amaryllidaceae is a large family with more than 1,600 species, belonging to 75 genera. The largest genus—Allium—is vast, comprising about a thousand species. Allium species (as well as other members of the Amaryllidaceae) are widespread and diversified, they are adapted to a wide range of habitats from shady forests to open habitats like meadows, steppes, and deserts. The genes present in chloroplast genomes (plastomes) play fundamental roles for the photosynthetic plants. Plastome traits could thus be associated with geophysical abiotic characteristics of habitats. Most chloroplast genes are highly conserved and are used as phylogenetic markers for many families of vascular plants. Nevertheless, some studies revealed signatures of positive selection in chloroplast genes of many plant families including Amaryllidaceae. We have sequenced plastomes of the following nine Allium (tribe Allieae of Allioideae) species: A. zebdanense, A. moly, A. victorialis, A. macleanii, A. nutans, A. obliquum, A. schoenoprasum, A. pskemense, A. platyspathum, A. fistulosum, A. semenovii, and Nothoscordum bivalve (tribe Leucocoryneae of Allioideae). We compared our data with previously published plastomes and provided our interpretation of Allium plastome genes’ annotations because we found some noteworthy inconsistencies with annotations previously reported. For Allium species we estimated the integral evolutionary rate, counted SNPs and indels per nucleotide position as well as compared pseudogenization events in species of three main phylogenetic lines of genus Allium to estimate whether they are potentially important for plant physiology or just follow the phylogenetic pattern. During examination of the 38 species of Allium and the 11 of other Amaryllidaceae species we found that rps16, rps2, infA, ccsA genes have lost their functionality multiple times in different species (regularly evolutionary events), while the pseudogenization of other genes was stochastic events. We found that the “normal” or “pseudo” state of rps16, rps2, infA, ccsA genes correlates well with the evolutionary line of genus the species belongs to. The positive selection in various NADH dehydrogenase (ndh) genes as well as in matK, accD, and some others were found. Taking into account known mechanisms of coping with excessive light by cyclic electron transport, we can hypothesize that adaptive evolution in genes, coding subunits of NADH-plastoquinone oxidoreductase could be driven by abiotic factors of alpine habitats, especially by intensive light and UV radiation.

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.

Molecular Evolution of rbcL in Orthotrichales (Bryophyta): Site Variation, Adaptive Evolution, and Coevolutionary Patterns of Amino Acid Replacements

Molecular evolution of the large subunit of the RuBisCO enzyme is understudied in early diverging land plants. These groups show morphological and eco-physiological adaptations to the uneven and intermittent distribution of water in the terrestrial environment. This might have prompted a continuous fine-tuning of RuBisCO under a selective pressure modifying the species-specific optima for photosynthesis in contrasting microdistributions and environmental niches. To gain a better insight into the molecular evolution of RuBisCO large subunits, the aim of this study was to assess the pattern of evolutionary change in the amino acid residues in a monophyletic group of Bryophyta (Orthotrichaceae). Tests for positive, neutral, or purifying selection at the amino acid level were assessed by comparing rates (ω) of non-synonymous (d_N) and synonymous (d_S) nucleotide substitutions along a Maximum Likelihood phylogenetic tree. Molecular adaptation tests using likelihood ratio tests, reconstruction of ancestral amino acid sites, and intra-protein coevolution analyses were performed. Variable amino acid sites (39) were unevenly distributed across the LSU. The residues are located on rbcL sites that are highly variable in higher plants and close to key regions implying dimer-dimer (L₂L₂), RuBisCO-activase interactions, and conformational functions during catalysis. Ten rbcL sites (32, 33, 91, 230, 247, 251, 255, 424, 449 and 475) have been identified by the Bayesian Empirical Bayes inference to be under positive selection and under adaptive evolution under the M8 model. The pattern of amino acid variation suggests that it is not lineage specific, but rather representative of a case of convergent evolution, suggesting recurrent changes that potentially favor the same amino acid substitutions that are likely optimized the RuBisCO activity.

Comparative genomics of cetartiodactyla: energy metabolism underpins the transition to an aquatic lifestyle

Foraging disruption caused by human activities is emerging as a key issue in cetacean conservation because it can affect nutrient levels and the amount of energy available to individuals to invest into reproduction. Our ability to predict how anthropogenic stressors affect these ecological processes and ultimately population trajectory depends crucially on our understanding of the complex physiological mechanisms that detect nutrient availability and regulate energy metabolism, foraging behavior and life-history decisions. These physiological mechanisms are likely to differ considerably from terrestrial mammalian model systems. Here, we examine nucleotide substitution rates in cetacean and other artiodactyl genomes to identify signatures of selection in genes associated with nutrient sensing pathways. We also estimated the likely physiological consequences of adaptive amino acid substitutions for pathway functions. Our results highlight that genes involved in the insulin, mTOR and NF-ĸB pathways are subject to significant positive selection in cetaceans compared to terrestrial artiodactyla. These genes may have been positively selected to enable cetaceans to adapt to a glucose-poor diet, to overcome deleterious effects caused by hypoxia during diving (e.g. oxidative stress and inflammation) and to modify fat-depot signaling functions in a manner different to terrestrial mammals. We thus show that adaptation in cetaceans to an aquatic lifestyle significantly affected functions in nutrient sensing pathways. The use of fat stores as a condition index in cetaceans may be confounded by the multiple and critical roles fat has in regulating cetacean metabolism, foraging behavior and diving physiology.

Dissection for Floral Micromorphology and Plastid Genome of Valuable Medicinal Borages Arnebia and Lithospermum (Boraginaceae)

The genera Arnebia and Lithospermum (Lithospermeae-Boraginaceae) comprise 25-30 and 50-60 species, respectively. Some of them are economically valuable, as their roots frequently contain a purple-red dye used in the cosmetic industry. Furthermore, dried roots of Arnebia euchroma, A. guttata, and Lithospermum erythrorhizon, which have been designated Lithospermi Radix, are used as traditional Korean herbal medicine. This study is the first report on the floral micromorphology and complete chloroplast (cp) genome sequences of A. guttata (including A. tibetana), A. euchroma, and L. erythrorhizon. We reveal great diversity in floral epidermal cell patterns, gynoecium, and structure of trichomes. The cp genomes were 149,361-150,465 bp in length, with conserved quadripartite structures. In total, 112 genes were identified, including 78 protein-coding regions, 30 tRNA genes, and four rRNA genes. Gene order, content, and orientation were highly conserved and were consistent with the general structure of angiosperm cp genomes. Comparison of the four cp genomes revealed locally divergent regions, mainly within intergenic spacer regions (atpH-atpI, petN-psbM, rbcL-psaI, ycf4-cemA, ndhF-rpl32, and ndhC-trnV-UAC). To facilitate species identification, we developed molecular markers psaA- ycf3 (PSY), trnI-CAU- ycf2 (TCY), and ndhC-trnV-UAC (NCTV) based on divergence hotspots. High-resolution phylogenetic analysis revealed clear clustering and a close relationship of Arnebia to its Lithospermum sister group, which was supported by strong bootstrap values and posterior probabilities. Overall, gynoecium characteristics and genetic distance of cp genomes suggest that A. tibetana, might be recognized as an independent species rather than a synonym of A. guttata. The present morphological and cp genomic results provide useful information for future studies, such as taxonomic, phylogenetic, and evolutionary analysis of Boraginaceae.