Diversification, biogeographic pattern, and demographic history of Taiwanese Scutellaria species inferred from nuclear and chloroplast DNA

A leading-edge scenario in the phylogeography and evolutionary history of East Asian insular Taxus in Taiwan and the Philippines

The cool temperate origin of gymnosperm Taxus species in East Asia is specifically diverse and widespread. Certain lineages have managed to extend their distribution further south to subtropical and tropical islands such as Taiwan and the Philippines. To address questions including whether these insular lineages, recently identified as T. phytonii, have become genetically distinct from each other and from their continental relatives, and when and how they colonized their residing islands, we sampled over 11 populations, covering 179 Taxus individuals from Taiwan and the Philippines. Using four cpDNA and one nuclear marker, we showed in population genetic and genealogical analyses that the two insular lineages were genetically distinct from each other and also from other continental Taxus and that they represented each other's closest relative. Estimated with the coalescent-based multi-type tree (MTT) analyses, we inferred an origin of Taiwanese T. phytonii more ancient than 2.49 Mya and that of Philippine T. phytonii more ancient than 1.08 Mya. In addition, the divergence demographic history revealed by both MTT and isolation with migration (IM) analyses indicated the presence of recent post-split migrations from a continental taxon, T. mairei, to Taiwanese T. phytonii, as well as from Taiwanese T. phytonii to Philippine T. phytonii. Overall, this study suggests Taiwan as a stepping stone through which the temperate-origin yew trees can extend their distributions to tropical regions such as the Philippines.

Comparative analysis of medicinal plants Scutellaria baicalensis and common adulterants based on chloroplast genome sequencing

BACKGROUND

Scutellaria baicalensis Georgi has been extensively used as a medicinal herb in China for over 2000 years. They may be intentionally or inadvertently substituted or blended with comparable species in the local market, threatening clinical medication safety. Molecular markers are effective tools to prevent misidentification and eliminate doping and falsification among Scutellaria plants. This study screened four highly variable regions to identify Scutellaria and its adulterants. In addition, a phylogenetic analysis was performed using the complete cp genome combined with published Scutellaria species samples. Moreover, a comparative analysis of the cp genomes was conducted to investigate the cp genome evolution of S. baicalensis.

RESULTS

The complete cp genome of five species of Scutellaria was sequenced for the first time, and four previously published Scutellaria species were re-sequenced. They all exhibited a conserved quadripartite structure in their cp genomes, including two distinct regions, namely a small and large single copy region, respectively, and two inverted repeats encompassing the majority of ribosomal RNA genes. Furthermore, the nine species exhibited high conservation from aspects of the genome structure, codon usage, repeat sequences, and gene content. Four highly variable regions (matK-rps16, ndhC-trnV-UAC, psbE-petL, and rps16-trnQ-UUG) may function as potential molecular markers for differentiating S. baicalensis from its adulterants. Additionally, the monophyly of Scutellaria was ascertained and could be reclassified into two subgenera, subgenus Anaspis and subgenus Scutellaria, as evidenced by the phylogenetic analyses on sequences of cp genome and shared protein-coding sequences. According to the molecular clock analysis, it has been inferred that the divergence of Scutellaria occurred at approximately 4.0 Mya during the Pliocene Epoch.

CONCLUSION

Our study provides an invaluable theoretical basis for further Scutellaria species identification, phylogenetics, and evolution analysis.

Ghost introgression facilitates genomic divergence of a sympatric cryptic lineage in Cycas revoluta

A cryptic lineage is a genetically diverged but morphologically unrecognized variant of a known species. Clarifying cryptic lineage evolution is essential for quantifying species diversity. In sympatric cryptic lineage divergence compared with allopatric divergence, the forces of divergent selection and mating patterns override geographical isolation. Introgression, by supplying preadapted or neutral standing genetic variations, can promote sympatric cryptic lineage divergence via selection. However, most studies concentrated on extant species introgression, ignoring the genetic legacy of introgression from extinct or unsampled lineages ("ghost introgression"). Cycads are an ideal plant for studying the influence of ghost introgression because of their common interspecific gene flow and past high extinction rate. Here, we utilized reference-based ddRADseq to clarify the role of ghost introgression in the evolution of a previously identified sympatric cryptic lineage in Cycas revoluta. After re-evaluating the evolutionary independency of cryptic lineages, the group-wise diverged single-nucleotide polymorphisms among sympatric and allopatric lineages were compared and functionally annotated. Next, we employed an approximate Bayesian computation method for hypothesis testing to clarify the cryptic lineage evolution and ghost introgression effect. SNPs with the genomic signatures of ghost introgression were further annotated. Our results reconfirmed the evolutionary independency of cryptic lineage among C. revoluta and demonstrated that ghost introgression to the noncryptic lineage facilitated their divergence. Gene function related to heat stress and disease resistance implied ecological adaptation of the main extant populations of C. revoluta.

Identification of three cultivated varieties of Scutellaria baicalensis using the complete chloroplast genome as a super-barcode

Scutellaria baicalensis has been one of the most commonly used traditional Chinese medicinal plants in China for more than 2000 years. The three new varieties cultivated could not be distinguished by morphology before flowering. It will hinder the promotion of later varieties. Chloroplast DNA has been widely used in species identification. Moreover, previous studies have shown that complete chloroplast genome sequences have been suggested as super barcodes for identifying plants. Therefore, we sequenced and annotated the complete chloroplast genomes of three cultivated varieties. The chloroplast genomes of SBW, SBR, and SBP were 151,702 bp, 151,799 bp, and 151,876 bp, which contained 85 protein-coding genes, 36 tRNA genes, and 8 rRNA genes. The analysis of the repeat sequences, codon usage, and comparison of chloroplast genomes shared a high degree of conservation. However, the sliding window results show significant differences among the three cultivated varieties in matK-rps16 and petA-psbJ. And we found that the matK-rps16 sequence can be used as a barcode for the identification of three varieties. In addition, the complete chloroplast genome contains more variations and can be used as a super-barcode to identify these three cultivated varieties. Based on the protein-coding genes, the phylogenetic tree demonstrated that SBP was more closely related to SBW, in the three cultivated varieties. Interestingly, we found that S. baicalensis and S. rehderiana are closely related, which provides new ideas for the development of S. baicalensis. The divergence time analysis showed that the three cultivated varieties diverged at about 0.10 Mya. Overall, this study showed that the complete chloroplast genome could be used as a super-barcode to identify three cultivated varieties of S. baicalensis and provide biological information, and it also contributes to bioprospecting.

Morphological and Genetic Diversity of Scutellaria tuvensis Juz., an Endemic of Desert Steppes

For the first time, an assessment of phenotypic variability and genetic polymorphism was performed on endemic plants Scutellaria tuvensis Juz. growing in Tuva (five populations; Russia). Based on morphological traits of individuals, principal component analysis clustered the individuals into three groups depending on characteristics of their habitats: group 1 turned out to be sampled from beach gravel, group 2 from a detrital cone, and group 3 from coarse rock fragments; this finding was confirmed by specific features of the development of the individuals in these habitats. Using inter-simple sequence repeat markers, high genetic polymorphism was identified at the population level: the proportion of polymorphic loci was 95%, expected heterozygosity 0.221, the absolute number of alleles 1.533, and the effective number of alleles 1.376. Population 3 (P 3) was the most genetically homogeneous; P 5 was characterized by the highest genetic diversity. In an unweighted pair group method with arithmetic mean dendrogram, the studied populations formed two major groups: the first cluster included P 4 and P 5, and the second cluster contained plants collected in P 1, P 2, and P 3. An analysis of the population structure using the STRUCTURE software showed the same result, dividing the sample under study into two subpopulations. The genetic differentiation index among populations was 0.232, and gene flow 1.655. According to analysis of molecular variance, intrapopulation differences accounted for 73% of total genetic diversity.

Traditional uses, ten-years research progress on phytochemistry and pharmacology, and clinical studies of the genus Scutellaria

ETHNOPHARMACOLOGICAL RELEVANCE

Scutellaria (Lamiaceae), which includes approximately 360-469 accepted species, is widespread in Europe, North America, East Asia, and South America. Several species have a long history being used as traditional medicines to treat respiratory, peptic, neurological, and hepatic and gall diseases. The phytochemistry and pharmacology of the genus Scutellaria have been developed dramatically in the past ten years, and the traditional uses and clinical studies of the genus have not been systematically summarized. Therefore, it is especially valuable to review the current state of knowledge to provide a basis for further exploration of its medicinal potential.

AIM OF THE REVIEW

The review aims to provide updated information on the ethnopharmacology, the ten-year research progress of phytochemistry and pharmacology, and clinical studies of Scutellaria and to explore the potential medicinal values and further studies of Scutellaria.

MATERIALS AND METHODS

This review is based on published studies and books from the library and electronic sources, including SciFinder, Scopus, PubMed, Web of Science, Baidu Scholar, CNKI, the online ethnobotanical database, and ethnobotanical monographs. This literature is related to ethnopharmacology, the ten-year research progress on the phytochemistry and pharmacology, and clinical studies of Scutellaria.

RESULTS

A total of 50 species, 5 subspecies and 17 varieties of the genus Scutellaria are used as traditional medicine with various biological activities. In the past ten years, 208 chemical constituents have been identified from 16 species and 1 variety of the genus Scutellaria, such as neo-clerodane diterpenoids, sesterterpenoids, terpenoids, flavonoids. Pharmacological research has demonstrated that the extracts and compounds identified from this genus exhibit extensive biological activities, including anticancer, antioxidant, anti-inflammatory, antiviral and antibacterial activities, effects on cardiovascular, cerebrovascular diseases as well as hepatoprotective and neuroprotective effects. The species S. baicalensis, S. barbata, and S. lateriflora and the main compounds baicalein, baicalin and wogonin are involved in clinical trials, which point the way for us to conduct further studies, such as study on the anticancer, antihypertensive, anti-infective, anti-inflammatory, neuroprotective and other effects of Scutellaria.

CONCLUSIONS

The species included in the genus Scutellaria can be used to treat cancer, infection, hepatic disorders, cardiovascular and cerebrovascular diseases, neurodegenerative diseases, and other diseases. Some indications in traditional medicines have been confirmed by modern pharmacological studies, such as anticancer, anti-inflammatory, anti-infective activity, and hepatoprotective and neuroprotective effects. The available literature indicated that most of the bioactivities could be attributed to flavonoids and neo-clerodane diterpenoids. Although there are some uses of Scutellaria in clinical practice, the existing research on this genus is still limited. In order to expand the development of medicinal resources of Scutellaria, the already studied species in this genus are recommended for more comprehensive investigation on their active substances, pharmacological mechanisms, quality control, clinical use and new drug research. Additionally, it is necessary to study species that their chemical composition or pharmacological activity have not yet been investigated, especially those used in folk medicine.

Characterization and comparative analysis among plastome sequences of eight endemic Rubus (Rosaceae) species in Taiwan

Genus Rubus represents the second largest genus of the family Rosaceae in Taiwan, with 41 currently recognized species across three subgenera (Chamaebatus, Idaoeobatus, and Malochobatus). Despite previous morphological and cytological studies, little is known regarding the overall phylogenetic relationships among the Rubus species in Taiwan, and their relationships to congeneric species in continental China. We characterized eight complete plastomes of Taiwan endemic Rubus species: subg. Idaeobatus (R. glandulosopunctatus, R. incanus, R. parviaraliifolius, R rubroangustifolius, R. taitoensis, and R. taiwanicolus) and subg. Malachobatus (R. kawakamii and R. laciniastostipulatus) to determine their phylogenetic relationships. The plastomes were highly conserved and the size of the complete plastome sequences ranged from 155,566 to 156,236 bp. The overall GC content ranged from 37.0 to 37.3%. The frequency of codon usage showed similar patterns among species, and 29 of the 73 common protein-coding genes were positively selected. The comparative phylogenomic analysis identified four highly variable intergenic regions (rps16/trnQ, petA/psbJ, rpl32/trnL-UAG, and trnT-UGU/trnL-UAA). Phylogenetic analysis of 31 representative complete plastomes within the family Rosaceae revealed three major lineages within Rubus in Taiwan. However, overall phylogenetic relationships among endemic species require broader taxon sampling to gain new insights into infrageneric relationships and their plastome evolution.

An updated tribal classification of Lamiaceae based on plastome phylogenomics

BACKGROUND

A robust molecular phylogeny is fundamental for developing a stable classification and providing a solid framework to understand patterns of diversification, historical biogeography, and character evolution. As the sixth largest angiosperm family, Lamiaceae, or the mint family, consitutes a major source of aromatic oil, wood, ornamentals, and culinary and medicinal herbs, making it an exceptionally important group ecologically, ethnobotanically, and floristically. The lack of a reliable phylogenetic framework for this family has thus far hindered broad-scale biogeographic studies and our comprehension of diversification. Although significant progress has been made towards clarifying Lamiaceae relationships during the past three decades, the resolution of a phylogenetic backbone at the tribal level has remained one of the greatest challenges due to limited availability of genetic data.

RESULTS

We performed phylogenetic analyses of Lamiaceae to infer relationships at the tribal level using 79 protein-coding plastid genes from 175 accessions representing 170 taxa, 79 genera, and all 12 subfamilies. Both maximum likelihood and Bayesian analyses yielded a more robust phylogenetic hypothesis relative to previous studies and supported the monophyly of all 12 subfamilies, and a classification for 22 tribes, three of which are newly recognized in this study. As a consequence, we propose an updated phylogenetically informed tribal classification for Lamiaceae that is supplemented with a detailed summary of taxonomic history, generic and species diversity, morphology, synapomorphies, and distribution for each subfamily and tribe.

CONCLUSIONS

Increased taxon sampling conjoined with phylogenetic analyses based on plastome sequences has provided robust support at both deep and shallow nodes and offers new insights into the phylogenetic relationships among tribes and subfamilies of Lamiaceae. This robust phylogenetic backbone of Lamiaceae will serve as a framework for future studies on mint classification, biogeography, character evolution, and diversification.

Comparative Genome Analysis of Scutellaria baicalensis and Scutellaria barbata Reveals the Evolution of Active Flavonoid Biosynthesis

Scutellaria baicalensis (S. baicalensis) and Scutellaria barbata (S. barbata) are common medicinal plants of the Lamiaceae family. Both produce specific flavonoid compounds, including baicalein, scutellarein, norwogonin, and wogonin, as well as their glycosides, which exhibit antioxidant and antitumor activities. Here, we report chromosome-level genome assemblies of S. baicalensis and S. barbata with quantitative chromosomal variation (2n = 18 and 2n = 26, respectively). The divergence of S. baicalensis and S. barbata occurred far earlier than previously reported, and a whole-genome duplication (WGD) event was identified. The insertion of long terminal repeat elements after speciation might be responsible for the observed chromosomal expansion and rearrangement. Comparative genome analysis of the congeneric species revealed the species-specific evolution of chrysin and apigenin biosynthetic genes, such as the S. baicalensis-specific tandem duplication of genes encoding phenylalanine ammonia lyase and chalcone synthase, and the S. barbata-specific duplication of genes encoding 4-CoA ligase. In addition, the paralogous duplication, colinearity, and expression diversity of CYP82D subfamily members revealed the functional divergence of genes encoding flavone hydroxylase between S. baicalensis and S. barbata. Analyzing these Scutellaria genomes reveals the common and species-specific evolution of flavone biosynthetic genes. Thus, these findings would facilitate the development of molecular breeding and studies of biosynthesis and regulation of bioactive compounds.

Leveraging plastomes for comparative analysis and phylogenomic inference within Scutellarioideae (Lamiaceae)

Scutellaria, or skullcaps, are medicinally important herbs in China, India, Japan, and elsewhere. Though Scutellaria is the second largest and one of the more taxonomically challenging genera within Lamiaceae, few molecular systematic studies have been undertaken within the genus; in part due to a paucity of available informative markers. The lack of informative molecular markers for Scutellaria hinders our ability to accurately and robustly reconstruct phylogenetic relationships, which hampers our understanding of the diversity, phylogeny, and evolutionary history of this cosmopolitan genus. Comparative analyses of 15 plastomes, representing 14 species of subfamily Scutellarioideae, indicate that plastomes within Scutellarioideae contain about 151,000 nucleotides, and possess a typical quadripartite structure. In total, 590 simple sequence repeats, 489 longer repeats, and 16 hyper-variable regions were identified from the 15 plastomes. Phylogenetic relationships among the 14 species representing four of the five genera of Scutellarioideae were resolved with high support values, but the current infrageneric classification of Scutellaria was not supported in all analyses. Complete plastome sequences provide better resolution at an interspecific level than using few to several plastid markers in phylogenetic reconstruction. The data presented here will serve as a foundation to facilitate DNA barcoding, species identification, and systematic research within Scutellaria, which is an important medicinal plant resource worldwide.

Taxonomic notes on Scutellaria taipeiensis (Lamiaceae) from morphological and molecular data

The genus Scutellaria comprises eight species distributed from 50 to 2000 m in Taiwan. Amongst them, S. barbata and S. taipeiensis are very similar on the basis of morphological and plastid DNA sequence information. Therefore, a comprehensive study of the taxonomic status of S. taipeiensis is necessary. We reviewed the herbarium sheets, related literature and protologues and compared morphologies of these two species, as well as their phylogenetic relationships. All evidence, including the diagnostic characters between S. taipeiensis and S. barbata, suggest that they belonged to a single species rather than two. As a result, S. taipeiensis is treated as a synonym of S. barbata.

Nuclear and chloroplast DNA phylogeography suggests an Early Miocene southward expansion of Lithocarpus (Fagaceae) on the Asian continent and islands

BACKGROUND

Most genera of Fagaceae are thought to have originated in the temperate regions except for the genus Lithocarpus, the stone oaks. Lithocarpus is distributed in subtropical and tropical Asia, and its ancestral population is hypothesized to be distributed in tropical regions in Borneo and Indochina. Borneo and the nearby islands (the Greater Sunda Islands) were connected to the Malay Peninsula and Indochina prior to the Pliocene epoch and formed the former Sundaland continent. The Southeast Asian Lithocarpus, is thought to have dispersed between continental Asia and the present Sundaland. The drastic climate changes during the Pliocene and Pleistocene epochs which caused periodic sea-level changes is often used to explain the cause of its diversity. The aim of this study was to establish phylogenetic relationships by analyzing nuclear (nrDNA) and chloroplast (cpDNA) DNA in order to describe and analyze the origin, causes of diversification and historical biogeography of Lithocarpus.

RESULTS

Phylogeny reconstructed through the multiple-species coalescent method with nrDNA and cpDNA revealed that the continental-Asian taxa were clustered at the basal lineages. The derived lineages of tropical Lithocarpus, with the inference of a subtropical ancestral state, imply a southward migration in the Early Miocene period with subsequent in situ diversification in the Greater Sunda Islands. The gradual decrease in temperature since the Middle Miocene period is proposed as a cause of the increase in the net diversification rate.

CONCLUSIONS

The historical ancestral origin of Lithocarpus has been suggested to be mainland Asia. Southward migration in the Early Miocene period with subsequent in situ diversification could explain the current diversity of stone oaks in Southeast Asia. This study also considered the multiple origins of stone oaks currently indigenous to the subtropical islands offshore and near mainland China. Our results provide phylogenetic evidence for a subtropical origin of Asian stone oaks and reveal the process of diversification and how it fits into the timeline of major geologic and climatic events rather than local, episodic, rate-shifting events.

Interspecific Divergence of Two Sinalliaria (Brassicaceae) Species in Eastern China

How endemic species originated in eastern Asia has interested botanists for a long time. In this study, we combined experimental and computational modeling approaches to examine the morphological and genetic divergence and reproductive isolation of two tentative species of Sinalliaria (Brassicaceae) endemic to eastern China, S. limprichtiana and S. grandifolia. Most of the examined morphological characters (including hairs of leaf blades and stems, corolla length and width, and flower stalk length) were well-delineated between two species at the same ploidy level, and there was clear evidence of reproductive isolation between them (mainly due to post-pollination barriers) in the common garden environment. There were also strong and consistent divergences in the population genetic data. Coalescent simulations based on sequence variation of the nuclear genes suggest that interspecific divergence began during the Pleistocene when the climate oscillated in eastern Asia. Gene flow between two species appears to have been very limited and asymmetrical. Our results suggested that both species are well-differentiated and that the fast divergence between them might have been together shaped by both stochastic processes and habitat selection pressures.

Continuation of the genetic divergence of ecological speciation by spatial environmental heterogeneity in island endemic plants

Divergent selection plays a critical role not only as a speciation driver but also in maintaining post-speciation divergence. In the absence of direct evidence, ancestral interspecific gene flow between incipient species can reflect ancient selective pressure for ecological speciation. In the present study, two late-Pleistocene diverged species endemic to Taiwan, Scutellaria playfairii and S. tashiroi, were spatially and ecologically partitioned with partial overlap. Multilocus genome-scan analyses and in silico evaluation revealed ancestral interspecific gene flow but distinct genetic compositions, implying that adaptive divergence contributed to their speciation. Ecological niche modeling and principal component analysis suggested incomplete divergent niches between the two species; the species distribution is therefore consistent with Hutchinson's metaphor of multidimensional hypervolume niches rather than attributable to a single factor. Constraint ordination analysis supported this inference of a combination of variables explaining the genetic structure. The rare occurrence of hybrids in the sympatric population suggested hybrid breakdown, providing further evidence of divergent selection blocking gene flow. The correlation of environmental variables with integrated genetic components demonstrated that environmental heterogeneity maintains the species and population differentiation. This study highlights the importance of environmental heterogeneity and divergent selection for the rapid speciation and recent diversification of island plants.

Phylogeography and genetic structure of the oriental river prawn Macrobrachium nipponense (Crustacea: Decapoda: Palaemonidae) in East Asia

The oriental river prawn (Macrobrachium nipponense) is mainly distributed in East Asia. The phylogeography, population genetic structure and historical demography of this species in the East Asia were examined by using partial sequences of the cytochrome oxidase subunit I (COI) and 16S rRNA in mitochondrial DNA. Ten populations that included 239 individuals were collected from Taiwan (Shihmen Reservoir, SMR, Mingte Reservoir, MTR and Chengching Lake Reservoir, CLR), mainland China (Taihu Lake, TLC, Min River, MRC, Jiulong River, JRC and Shenzhen Reservoir, SRC), Japan (Biwa Lake, BLJ and Kasumigaura Lake, KLJ) and Korea (Han River, HRK). The nucleotide diversity (π) of all individuals was 0.01134, with values ranging from 0.0089 (BLJ, Japan) to 0.01425 (MTR, Taiwan). A total of 83 haplotypes were obtained, and the haplotypes were divided into 2 main lineages: lineage A included the specimens from BLJ, KLJ, CLR, MTR, TLC, MRC and JRC, and lineage B comprised the ones from HRK, SRC, SMR, MTR, TLC, MRC and JRC. Lineage A could be further divided two sub-lineages (A₁ and A₂). Individuals of lineage A₂ were only from TLC. Demographic expansion was observed in each lineage, starting within the second-to-latest interglacial period for lineage A and within the last glacial period for lineage B. All F_ST values among the ten populations were significantly different, except for the values between MRC and JRC, and SMR and SRC. The phylogeography and genetic structure of M. nipponense in East Asia might be influenced by Pleistocene glacial cycles, lake isolation and human introduction. The possible dispersal routes of M. nipponense in the East Asia were also discussed.

Local Climate Heterogeneity Shapes Population Genetic Structure of Two Undifferentiated Insular Scutellaria Species

Spatial climate heterogeneity may not only affect adaptive gene frequencies but could also indirectly shape the genetic structure of neutral loci by impacting demographic dynamics. In this study, the effect of local climate on population genetic variation was tested in two phylogenetically close Scutellaria species in Taiwan. Scutellaria taipeiensis, which was originally assumed to be an endemic species of Taiwan Island, is shown to be part of the widespread species S. barbata based on the overlapping ranges of genetic variation and climatic niches as well as their morphological similarity. Rejection of the scenario of "early divergence with secondary contact" and the support for multiple origins of populations of S. taipeiensis from S. barbata provide strong evolutionary evidence for a taxonomic revision of the species combination. Further tests of a climatic effect on genetic variation were conducted. Regression analyses show nonlinear correlations among any pair of geographic, climatic, and genetic distances. However, significantly, the bioclimatic variables that represent the precipitation from late summer to early autumn explain roughly 13% of the genetic variation of our sampled populations. These results indicate that spatial differences of precipitation in the typhoon season may influence the regeneration rate and colonization rate of local populations. The periodic typhoon episodes explain the significant but nonlinear influence of climatic variables on population genetic differentiation. Although, the climatic difference does not lead to species divergence, the local climate variability indeed impacts the spatial genetic distribution at the population level.

Imbalanced positive selection maintains the functional divergence of duplicated DIHYDROKAEMPFEROL 4-REDUCTASE genes

Gene duplication could be beneficial by functional division but might increase the risk of genetic load. The dynamics of duplicated paralogs number could involve recombination, positive selection, and functional divergence. Duplication of DIHYDROFLAVONOL 4-REDUCTASE (DFR) has been reported in several organisms and may have been retained by escape from adaptive conflict (EAC). In this study, we screened the angiosperm DFR gene focusing on a diversified genus Scutellaria to investigate how these duplicated genes are retained. We deduced that gene duplication involved multiple independent events in angiosperms, but the duplication of DFR was before the divergence of Scutellaria. Asymmetric positive selective pressures resulted in different evolutionary rates between the duplicates. Different numbers of regulatory elements, differential codon usages, radical amino acid changes, and differential gene expressions provide evidences of functional divergence between the two DFR duplicates in Scutellaria, implying adaptive subfunctionalization between duplicates. The discovery of pseudogenes accompanying a reduced replacement rate in one DFR paralogous gene suggested possibly leading to “loss of function” due to dosage imbalance after the transient adaptive subfunctionalization in the early stage of duplication. Notwithstanding, episodic gene duplication and functional divergence may be relevant to the diversification of ecological function of DFR gene in Scutellaria.

Diversifying selection of the anthocyanin biosynthetic downstream gene UFGT accelerates floral diversity of island Scutellaria species

BACKGROUND

Adaptive divergence, which usually explains rapid diversification within island species, might involve the positive selection of genes. Anthocyanin biosynthetic pathway (ABP) genes are important for floral diversity, and are related to stress resistance and pollination, which could be responsible for species diversification. Previous studies have shown that upstream genes of ABP are subject to selective constraints and have a slow evolutionary rate, while the constraints on downstream genes are lower.

RESULTS

In this study, we confirmed these earlier observations of heterogeneous evolutionary rate in upstream gene CHS and the downstream gene UFGT, both of which were expressed in Scutellaria sp. inflorescence buds. We found a higher evolutionary rate and positive selection for UFGT. The codons under positive selection corresponded to the diversified regions, and the presence or absence of an α-helix might produce conformation changes in the Rossmann-like fold structure. The significantly high evolutionary rates for UFGT genes in Taiwanese lineages suggested rapid accumulation of amino acid mutations in island species. The results showed positive selection in closely related species and explained the high diversification of floral patterns in these recently diverged species. In contrast, non-synonymous mutation rate decreases in long-term divergent species could reduce mutational load and prevent the accumulation of deleterious mutations.

CONCLUSIONS

Together with the positive selection of transcription factors for ABP genes described in a previous study, these results confirmed that positive selection takes place at a translational level and suggested that the high divergence of ABP genes is related to the floral diversity in island Scutellaria species.

Biogeography of the Phalaenopsis amabilis species complex inferred from nuclear and plastid DNAs

BACKGROUND

Phalaenopsis is one of the important commercial orchids in the world. Members of the P. amabilis species complex represent invaluable germplasm for the breeding program. However, the phylogeny of the P. amabilis species complex is still uncertain. The Phalaenopsis amabilis species complex (Orchidaceae) consists of subspecies amabilis, moluccana, and rosenstromii of P. amabilis, as well as P. aphrodite ssp. aphrodite, P. ap. ssp. formosana, and P. sanderiana. The aims of this study were to reconstruct the phylogeny and biogeographcial patterns of the species complex using Neighbor Joining (NJ), Maxinum Parsimony (MP), Bayesian Evolutionary Analysis Sampling Trees (BEAST) and Reconstruct Ancestral State in Phylogenies (RASP) analyses based on sequences of internal transcribed spacers 1 and 2 from the nuclear ribosomal DNA and the trnH-psbA spacer from the plastid DNA.

RESULTS

A pattern of vicariance, dispersal, and vicariance + dispersal among disjunctly distributed taxa was uncovered based on RASP analysis. Although two subspecies of P. aphrodite could not be differentiated from each other in dispersal state, they were distinct from P. amabilis and P. sanderiana. Within P. amabilis, three subspecies were separated phylogenetically, in agreement with the vicariance or vicariance + dispersal scenario, with geographic subdivision along Huxley's, Wallace's and Lydekker's Lines. Molecular dating revealed such subdivisions among taxa of P. amabilis complex dating back to the late Pleistocene. Population-dynamic analyses using a Bayesian skyline plot suggested that the species complex experienced an in situ range expansion and population concentration during the late Last Glacial Maximum (LGM).

CONCLUSIONS

Taxa of the P. amabilis complex with disjunct distributions were differentiated due to vicariance or vicariance + dispersal, with events likely occurring in the late Pleistocene. Demographic growth associated with the climatic oscillations in the Würm glacial period followed the species splits. Nevertheless, a subsequent population slowdown occurred in the late LGM due to extinction of regional populations. The reduction of suitable habitats resulted in geographic fragmenttation of the remaining taxa.

Positive selection and functional divergence of R2R3-MYB paralogous genes expressed in inflorescence buds of Scutellaria species (Labiatae)

Anthocyanin is the main pigment forming floral diversity. Several transcription factors that regulate the expression of anthocyanin biosynthetic genes belong to the R2R3-MYB family. Here we examined the transcriptomes of inflorescence buds of Scutellaria species (skullcaps), identified the expression R2R3-MYBs, and detected the genetic signatures of positive selection for adaptive divergence across the rapidly evolving skullcaps. In the inflorescence buds, seven R2R3-MYBs were identified. MYB11 and MYB16 were detected to be positively selected. The signature of positive selection on MYB genes indicated that species diversification could be affected by transcriptional regulation, rather than at the translational level. When comparing among the background lineages of Arabidopsis, tomato, rice, and Amborella, heterogeneous evolutionary rates were detected among MYB paralogs, especially between MYB13 and MYB19. Significantly different evolutionary rates were also evidenced by type-I functional divergence between MYB13 and MYB19, and the accelerated evolutionary rates in MYB19, implied the acquisition of novel functions. Another paralogous pair, MYB2/7 and MYB11, revealed significant radical amino acid changes, indicating divergence in the regulation of different anthocyanin-biosynthetic enzymes. Our findings not only showed that Scutellaria R2R3-MYBs are functionally divergent and positively selected, but also indicated the adaptive relevance of regulatory genes in floral diversification.