Distinctive plastome evolution in carnivorous angiosperms.

BACKGROUND: Independent origins of carnivory in multiple angiosperm families are fabulous examples of convergent evolution using a diverse array of life forms and habitats. Previous studies have indicated that carnivorous plants have distinct evolutionary trajectories of plastid genome (plastome) compared to their non-carnivorous relatives, yet the extent and general characteristics remain elusive. RESULTS: We compared plastomes from 9 out of 13 carnivorous families and their non-carnivorous relatives to assess carnivory-associated evolutionary patterns. We identified inversions in all sampled Droseraceae species and four species of Utricularia, Pinguicula, Darlingtonia and Triphyophyllum. A few carnivores showed distinct shifts in inverted repeat boundaries and the overall repeat contents. Many ndh genes, along with some other genes, were independently lost in several carnivorous lineages. We detected significant substitution rate variations in most sampled carnivorous lineages. A significant overall substitution rate acceleration characterizes the two largest carnivorous lineages of Droseraceae and Lentibulariaceae. We also observe moderate substitution rates acceleration in many genes of Cephalotus follicularis, Roridula gorgonias, and Drosophyllum lusitanicum. However, only a few genes exhibit significant relaxed selection. CONCLUSION: Our results indicate that the carnivory of plants have different effects on plastome evolution across carnivorous lineages. The complex mechanism under carnivorous habitats may have resulted in distinctive plastome evolution with conserved plastome in the Brocchinia hechtioides to strongly reconfigured plastomes structures in Droseraceae. Organic carbon obtained from prey and the efficiency of utilizing prey-derived nutrients might constitute possible explanation.

Genetic imprints of grafting in wild iron walnut populations in southwestern China.

BACKGROUND: Anthropogenic activities are causing unprecedented loss of genetic diversity in many species. However, the effects on genetic diversity from large-scale grafting onto wild plants of crop species are largely undetermined. Iron walnut (Juglans sigillata Dode) is a deciduous nut tree crop endemic to southwestern China with a long history of cultivation. Due to the rapid expansion of the walnut industry, many natural populations are now being replaced by cultivars grafted onto wild rootstocks. However, little is known about the potential genetic consequences of such action on natural populations. RESULTS: We sampled the scion and the rootstock from each of 149 grafted individuals within nine wild populations of J. sigillata from Yunnan Province which is the center of walnut diversity and cultivation in China, and examined their genetic diversity and population structure using 31 microsatellite loci. Scions had lower genetic diversity than rootstocks, and this pattern was repeated in seven of the nine examined populations. Among those seven populations, AMOVA and clustering analyses showed a clear genetic separation between all rootstocks and all scions. However, the two remaining populations, both from northern Yunnan, showed genetic similarity between scions and rootstocks, possibly indicating that wild populations here are derived from feralized local cultivars. Moreover, our data indicated probable crop-to-wild gene flow between scions and rootstocks, across all populations. CONCLUSIONS: Our results indicate that large-scale grafting has been causing genetic diversity erosion and genetic structure breakdown in the wild material of J. sigillata within Yunnan. To mitigate these effects, we caution against the overuse of grafting in wild populations of iron walnut and other crop species and recommend the preservation of natural genotypes through in situ  and ex situ conservation.

Multiple paternally inherited chloroplast capture events associated with Taxus speciation in the Hengduan Mountains.

Mountainous regions provide a multitude of habitats and opportunities for complex speciation scenarios. Hybridization leading to chloroplast capture, which can be revealed by incongruent phylogenetic trees, is one possible outcome. Four allopatric Taxus lineages (three species and an undescribed lineage) from the Hengduan Mountains, southwest China, exhibit conflicting phylogenetic relationships between nuclear and chloroplast phylogenies. Here, we use multi-omic data at the population level to investigate their historical speciation processes. Population genomic analysis based on ddRAD-seq data revealed limited contemporary inter-specific gene flow involving only populations located close to another species. In a historical context, chloroplast and nuclear data (transcriptome) consistently showed conflicting phylogenetic relationships for T. florinii and the Emei type lineage. ILS and chloroplast recombination were excluded as possible causes, and transcriptome and ddRAD-seq data revealed an absence of the mosaic nuclear genomes that characterize hybrid origin scenarios. Therefore, T. florinii appears to have originated when a lineage of T. florinii captured the T. chinensis plastid type, whereas plastid introgression in the opposite direction generated the Emei Type. All four species have distinct ecological niche based on community investigations and ecological niche analyses. We propose that the origins of both species represent very rare examples of chloroplast capture events despite the paternal cpDNA inheritance of gymnosperms. Specifically, allopatrically and/or ecologically diverged parental species experienced a rare secondary contact, subsequent hybridization and reciprocal chloroplast capture, generating two new lineages, each of which acquired a unique ecological niche. These events might have been triggered by orogenic activities of the Hengduan Mountains and an intensification of the Asian monsoon in the late Miocene, and may represent a scenario more common in these mountains than presently known.

Phylogenomics and morphological evolution of the mega-diverse genus Artemisia (Asteraceae: Anthemideae): implications for its circumscription and infrageneric taxonomy.

BACKGROUND AND AIMS: Artemisia is a mega-diverse genus consisting of ~400 species. Despite its medicinal importance and ecological significance, a well-resolved phylogeny for global Artemisia, a natural generic delimitation and infrageneric taxonomy remain missing, owing to the obstructions from limited taxon sampling and insufficient information on DNA markers. Its morphological characters, such as capitulum, life form and leaf, show marked variations and are widely used in its infrageneric taxonomy. However, their evolution within Artemisia is poorly understood. Here, we aimed to reconstruct a well-resolved phylogeny for global Artemisia via a phylogenomic approach, to infer the evolutionary patterns of its key morphological characters and to update its circumscription and infrageneric taxonomy. METHODS: We sampled 228 species (258 samples) of Artemisia and its allies from both fresh and herbarium collections, covering all the subgenera and its main geographical areas, and conducted a phylogenomic analysis based on nuclear single nucleotide polymorphisms (SNPs) obtained from genome skimming data. Based on the phylogenetic framework, we inferred the possible evolutionary patterns of six key morphological characters widely used in its previous taxonomy. KEY RESULTS: The genus Kaschgaria was revealed to be nested in Artemisia with strong support. A well-resolved phylogeny of Artemisia consisting of eight highly supported clades was recovered, two of which were identified for the first time. Most of the previously recognized subgenera were not supported as monophyletic. Evolutionary inferences based on the six morphological characters showed that different states of these characters originated independently more than once. CONCLUSIONS: The circumscription of Artemisia is enlarged to include the genus Kaschgaria. The morphological characters traditionally used for the infrageneric taxonomy of Artemisia do not match the new phylogenetic tree. They experienced a more complex evolutionary history than previously thought. We propose a revised infrageneric taxonomy of the newly circumscribed Artemisia, with eight recognized subgenera to accommodate the new results.

Phylogenetic Relationships and Next-Generation Barcodes in the Genus Torreya Reveal a High Proportion of Misidentified Cultivated Plants.

Accurate species identification is key to conservation and phylogenetic inference. Living plant collections from botanical gardens/arboretum are important resources for the purpose of scientific research, but the proportion of cultivated plant misidentification are un-tested using DNA barcodes. Here, we assembled the next-generation barcode (complete plastid genome and complete nrDNA cistron) and mitochondrial genes from genome skimming data of Torreya species with multiple accessions for each species to test the species discrimination and the misidentification proportion of cultivated plants used in Torreya studies. A total of 38 accessions were included for analyses, representing all nine recognized species of genus Torreya. The plastid phylogeny showed that all 21 wild samples formed species-specific clades, except T. jiulongshanensis. Disregarding this putative hybrid, seven recognized species sampled here were successfully discriminated by the plastid genome. Only the T. nucifera accessions grouped into two grades. The species identification rate of the nrDNA cistron was 62.5%. The Skmer analysis based on nuclear reads from genome skims showed promise for species identification with seven species discriminated. The proportion of misidentified cultivated plants from arboreta/botanical gardens was relatively high with four accessions (23.5%) representing three species. Interspecific relationships within Torreya were fully resolved with maximum support by plastomes, where Torreya jackii was on the earliest diverging branch, though sister to T. grandis in the nrDNA cistron tree, suggesting that this is likely a hybrid species between T. grandis and an extinct Torreya ancestor lineage. The findings here provide quantitative insights into the usage of cultivated samples for phylogenetic study.

Phylogenomics and the flowering plant tree of life.

The advances accelerated by next-generation sequencing and long-read sequencing technologies continue to provide an impetus for plant phylogenetic study. In the past decade, a large number of phylogenetic studies adopting hundreds to thousands of genes across a wealth of clades have emerged and ushered plant phylogenetics and evolution into a new era. In the meantime, a roadmap for researchers when making decisions across different approaches for their phylogenomic research design is imminent. This review focuses on the utility of genomic data (from organelle genomes, to both reduced representation sequencing and whole-genome sequencing) in phylogenetic and evolutionary investigations, describes the baseline methodology of experimental and analytical procedures, and summarizes recent progress in flowering plant phylogenomics at the ordinal, familial, tribal, and lower levels. We also discuss the challenges, such as the adverse impact on orthology inference and phylogenetic reconstruction raised from systematic errors, and underlying biological factors, such as whole-genome duplication, hybridization/introgression, and incomplete lineage sorting, together suggesting that a bifurcating tree may not be the best model for the tree of life. Finally, we discuss promising avenues for future plant phylogenomic studies.

Phylotranscriptomics of Theaceae: generic-level relationships, reticulation and whole-genome duplication.

BACKGROUND AND AIMS: Theaceae, with three tribes, nine genera and more than 200 species, are of great economic and ecological importance. Recent phylogenetic analyses based on plastomic data resolved the relationships among the three tribes and the intergeneric relationships within two of those tribes. However, generic-level relationships within the largest tribe, Theeae, were not fully resolved. The role of putative whole-genome duplication (WGD) events in the family and possible hybridization events among genera within Theeae also remain to be tested further. METHODS: Transcriptomes or low-depth whole-genome sequencing of 57 species of Theaceae, as well as additional plastome sequence data, were generated. Using a dataset of low-copy nuclear genes, we reconstructed phylogenetic relationships using concatenated, species tree and phylogenetic network approaches. We further conducted molecular dating analyses and inferred possible WGD events by examining the distribution of the number of synonymous substitutions per synonymous site (Ks) for paralogues in each species. For plastid protein-coding sequences , phylogenies were reconstructed for comparison with the results obtained from analysis of the nuclear dataset. RESULTS: Based on the 610 low-copy nuclear genes (858 606 bp in length) investigated, Stewartieae was resolved as sister to the other two tribes. Within Theeae, the Apterosperma-Laplacea clade grouped with Pyrenaria, leaving Camellia and Polyspora as sister. The estimated ages within Theaceae were largely consistent with previous studies based mainly on plastome data. Two reticulation events within Camellia and one between the common ancestor of Gordonia and Schima were found. All members of the tea family shared two WGD events, an older At-γ and a recent Ad-ß; both events were also shared with the outgroups (Diapensiaceae, Pentaphylacaceae, Styracaceae and Symplocaceae). CONCLUSIONS: Our analyses using low-copy nuclear genes improved understanding of phylogenetic relationships at the tribal and generic levels previously proposed based on plastome data, but the phylogenetic position of the Apterosperma-Laplacea clade needs more attention. There is no evidence for extensive intergeneric hybridization within Theeae or for a Theaceae-specific WGD event. Land bridges (e.g. the Bering land bridge) during the Late Oligocene may have permitted the intercontinental plant movements that facilitated the putative ancient introgression between the common ancestor of Gordonia and Schima.

Resolution, conflict and rate shifts: insights from a densely sampled plastome phylogeny for Rhododendron (Ericaceae).

BACKGROUND AND AIMS: Rhododendron is a species-rich and taxonomically challenging genus due to recent adaptive radiation and frequent hybridization. A well-resolved phylogenetic tree would help to understand the diverse history of Rhododendron in the Himalaya-Hengduan Mountains where the genus is most diverse. METHODS: We reconstructed the phylogeny based on plastid genomes with broad taxon sampling, covering 161 species representing all eight subgenera and all 12 sections, including ~45 % of the Rhododendron species native to the Himalaya-Hengduan Mountains. We compared this phylogeny with nuclear phylogenies to elucidate reticulate evolutionary events and clarify relationships at all levels within the genus. We also estimated the timing and diversification history of Rhododendron, especially the two species-rich subgenera Rhododendron and Hymenanthes that comprise >90 % of Rhododendron species in the Himalaya-Hengduan Mountains. KEY RESULTS: The full plastid dataset produced a well-resolved and supported phylogeny of Rhododendron. We identified 13 clades that were almost always monophyletic across all published phylogenies. The conflicts between nuclear and plastid phylogenies suggested strongly that reticulation events may have occurred in the deep lineage history of the genus. Within Rhododendron, subgenus Therorhodion diverged first at 56 Mya, then a burst of diversification occurred from 23.8 to 17.6 Mya, generating ten lineages among the component 12 clades of core Rhododendron. Diversification in subgenus Rhododendron accelerated c. 16.6 Mya and then became fairly continuous. Conversely, Hymenanthes diversification was slow at first, then accelerated very rapidly around 5 Mya. In the Himalaya-Hengduan Mountains, subgenus Rhododendron contained one major clade adapted to high altitudes and another to low altitudes, whereas most clades in Hymenanthes contained both low- and high-altitude species, indicating greater ecological plasticity during its diversification. CONCLUSIONS: The 13 clades proposed here may help to identify specific ancient hybridization events. This study will help to establish a stable and reliable taxonomic framework for Rhododendron, and provides insight into what drove its diversification and ecological adaption. Denser sampling of taxa, examining both organelle and nuclear genomes, is needed to better understand the divergence and diversification history of Rhododendron.

Plastid phylogenomic insights into relationships of all flowering plant families.

BACKGROUND: Flowering plants (angiosperms) are dominant components of global terrestrial ecosystems, but phylogenetic relationships at the familial level and above remain only partially resolved, greatly impeding our full understanding of their evolution and early diversification. The plastome, typically mapped as a circular genome, has been the most important molecular data source for plant phylogeny reconstruction for decades. RESULTS: Here, we assembled by far the largest plastid dataset of angiosperms, composed of 80 genes from 4792 plastomes of 4660 species in 2024 genera representing all currently recognized families. Our phylogenetic tree (PPA II) is essentially congruent with those of previous plastid phylogenomic analyses but generally provides greater clade support. In the PPA II tree, 75% of nodes at or above the ordinal level and 78% at or above the familial level were resolved with high bootstrap support (BP ≥ 90). We obtained strong support for many interordinal and interfamilial relationships that were poorly resolved previously within the core eudicots, such as Dilleniales, Saxifragales, and Vitales being resolved as successive sisters to the remaining rosids, and Santalales, Berberidopsidales, and Caryophyllales as successive sisters to the asterids. However, the placement of magnoliids, although resolved as sister to all other Mesangiospermae, is not well supported and disagrees with topologies inferred from nuclear data. Relationships among the five major clades of Mesangiospermae remain intractable despite increased sampling, probably due to an ancient rapid radiation. CONCLUSIONS: We provide the most comprehensive dataset of plastomes to date and a well-resolved phylogenetic tree, which together provide a strong foundation for future evolutionary studies of flowering plants.

Natural hybridization among three Rhododendron species (Ericaceae) revealed by morphological and genomic evidence.

BACKGROUND: Natural hybridization can influence the adaptive response to selection and accelerate species diversification. Understanding the composition and structure of hybrid zones may elucidate patterns of hybridization processes that are important to the formation and maintenance of species, especially for taxa that have experienced rapidly adaptive radiation. Here, we used morphological traits, ddRAD-seq and plastid DNA sequence data to investigate the structure of a Rhododendron hybrid zone and uncover the hybridization patterns among three sympatric and closely related species. RESULTS: Our results show that the hybrid zone is complex, where bi-directional hybridization takes place among the three sympatric parental species: R. spinuliferum, R. scabrifolium, and R. spiciferum. Hybrids between R. spinuliferum and R. spiciferum (R. ×duclouxii) comprise multiple hybrid classes and a high proportion of F1 generation hybrids, while a novel hybrid taxon between R. spinuliferum and R. scabrifolium dominated the F2 generation, but no backcross individuals were detected. The hybrid zone showed basically coincident patterns of population structure between genomic and morphological data. CONCLUSIONS: Natural hybridization exists among the three Rhododendron species in the hybrid zone, although patterns of hybrid formation vary between hybrid taxa, which may result in different evolutionary outcomes. This study represents a unique opportunity to dissect the ecological and evolutionary mechanisms associated with adaptive radiation of Rhododendron species in a biodiversity hotspot.

Differential expressions of anthocyanin synthesis genes underlie flower color divergence in a sympatric Rhododendron sanguineum complex.

BACKGROUND: The Rhododendron sanguineum complex is endemic to alpine mountains of northwest Yunnan and southeast Tibet of China. Varieties in this complex exhibit distinct flower colors even at the bud stage. However, the underlying molecular regulations for the flower color variation have not been well characterized. Here, we investigated this via measuring flower reflectance profiles and comparative transcriptome analyses on three coexisting varieties of the R. sanguineum complex, with yellow flush pink, bright crimson, and deep blackish crimson flowers respectively. We compared the expression levels of differentially-expressed-genes (DEGs) of the anthocyanin / flavonoid biosynthesis pathway using RNA-seq and qRT-PCR data. We performed clustering analysis based on transcriptome-derived Single Nucleotide Polymorphisms (SNPs) data, and finally analyzed the promoter architecture of DEGs. RESULTS: Reflectance spectra of the three color morphs varied distinctively in the range between 400 and 700 nm, with distinct differences in saturation, brightness, hue, and saturation/hue ratio, an indirect measurement of anthocyanin content. We identified 15,164 orthogroups that were shared among the three varieties. The SNP clustering analysis indicated that the varieties were not monophyletic. A total of 40 paralogous genes encoding 12 enzymes contributed to the flower color polymorphism. These anthocyanin biosynthesis-related genes were associated with synthesis, modification and transportation properties (RsCHS, RsCHI, RsF3H, RsF3'H, RsFLS, RsANS, RsAT, RsOMT, RsGST), as well as genes involved in catabolism and degradation (RsBGLU, RsPER, RsCAD). Variations in sequence and cis-acting elements of these genes might correlate with the anthocyanin accumulation, thus may contribute to the divergence of flower color in the R. sanguineum complex. CONCLUSIONS: Our results suggested that the varieties are very closely related and flower color variations in the R. sanguineum complex correlate tightly with the differential expression levels of genes involved in the anabolic and catabolic synthesis network of anthocyanin. Our study provides a scenario involving intricate relationships between genetic mechanisms for floral coloration accompanied by gene flow among the varieties that may represent an early case of pollinator-mediated incipient sympatric speciation.

Repeated intercontinental migrations and recurring hybridizations characterise the evolutionary history of yew (Taxus L.).

The genus Taxus (Taxaceae) consists of 16 genetically well-defined lineages that are predominantly distributed across the Northern hemisphere. We investigated its biogeographic origin and evolutionary history by sampling 13 chloroplast gene sequences, the nuclear internal transcribed spacers (ITS) and NEEDLY sequences for all 16 lineages. We applied Maximum Parsimony and Bayesian Inference analyses to infer their phylogenetic relationships, time-calibrated phylogenies using BEAST and inferred the ancestral area of occupancy with BioGeoBEARS. We found strong evidence for the hybrid origin of three lineages and dated these events to a rather narrow time window of 6.8-4.9 million years ago (Mya). The dated phylogenies inferred an Upper Cretaceous origin of the genus, with the extant lineages diversifying in North America much later during the Oligocene/early Miocene. Repeated migrations via the Bering land bridge to Eurasia and back were further inferred, with the return to North America as a possible result of vicariance. The diversification in Eurasia (from ~8 Mya onwards) coincided with the orogeny of the Hengduan Mountains, the intensification of the East Asian summer monsoon and the occupancy of ecological niches by lineages that experienced secondary contacts and hybridizations in the Hengduan Mountains and Qinling Mountain, especially around the Sichuan basin. We provide a hypothesis for the evolution of extant lineages of Taxus, a genus with an old and complex evolutionary history. The study highlights that the history of complex species can be unravelled with a careful dissection of phylogenetic signals.

Greater than the sum of the parts: how the species composition in different forest strata influence ecosystem function.

The mechanisms underpinning forest biodiversity-ecosystem function relationships remain unresolved. Yet, in heterogeneous forests, ecosystem function of different strata could be associated with traits or evolutionary relationships differently. Here, we integrate phylogenies and traits to evaluate the effects of elevational diversity on above-ground biomass across forest strata and spatial scales. Community-weighted means of height and leaf phosphorous concentration and functional diversity in specific leaf area exhibited positive correlations with tree biomass, suggesting that both positive selection effects and complementarity occur. However, high shrub biomass is associated with greater dissimilarity in seed mass and multidimensional trait space, while species richness or phylogenetic diversity is the most important predictor for herbaceous biomass, indicating that species complementarity is especially important for understory function. The strength of diversity-biomass relationships increases at larger spatial scales. We conclude that strata- and scale- dependent assessments of community structure and function are needed to fully understand how biodiversity influences ecosystem function.

Direct and indirect effects of climate on richness drive the latitudinal diversity gradient in forest trees.

Climate is widely recognised as an important determinant of the latitudinal diversity gradient. However, most existing studies make no distinction between direct and indirect effects of climate, which substantially hinders our understanding of how climate constrains biodiversity globally. Using data from 35 large forest plots, we test hypothesised relationships amongst climate, topography, forest structural attributes (stem abundance, tree size variation and stand basal area) and tree species richness to better understand drivers of latitudinal tree diversity patterns. Climate influences tree richness both directly, with more species in warm, moist, aseasonal climates and indirectly, with more species at higher stem abundance. These results imply direct limitation of species diversity by climatic stress and more rapid (co-)evolution and narrower niche partitioning in warm climates. They also support the idea that increased numbers of individuals associated with high primary productivity are partitioned to support a greater number of species.

Plastid phylogenomics and biogeographic analysis support a trans-Tethyan origin and rapid early radiation of Cornales in the Mid-Cretaceous.

The Cornales is a relatively small but morphologically diverse order in the basal position of the Asterids clade. Previous study hypothesized that the order might have undergone ancient rapid radiation during the Cretaceous when major angiosperm lineages were established. We conducted the phylogenomic analysis of Cornales using 81 plastid genome sequences with 67 newly generated in this study to test the hypothesis. This sampling represents all the families and 31 out of 48 genera in the order. Phylogenetic analyses were conducted using different datasets to examine the effects of different coding positions and character coding methods. We further conducted divergence time, diversification rate, and biogeographic analyses to understand the early evolutionary history of Cornales in space and time. Our phylogenetic analyses of four datasets (the amino acid characters, the 1st and 2nd codon positions of protein coding genes, nucleotide characters with degenerated coding method, and noncoding regions) resulted in a robust phylogeny congruent with results of previous studies, showing (((Cornaceae-Alangiaceae)-(Curtisiaceae-Grubbiaceae))-(((Nyssaceae-Davidiaceae)-Mastixiaceae)-((Hydrostachyaceae-(Hydrangeaceae-Loasaceae)))). Phylogenetic relationships within families were also well resolved. Conflicts in the placement of Hydrostachyaceae were found from analyses of two datasets, the nucleotide characters of all codon position and the 3rd codon positions, where the family was united with Loasaceae, but not strongly supported. Results from divergence time analyses suggested a mid-Cretaceous origin of Cornales followed by rapid early diversification into major clades/families within 10 million years. The early diversification of Cornales may have been facilitated by divergence in habitat and morphology following geographic dispersals. The ancestral distribution of the order was inferred as a widespread range covering Asia, Europe, North America, and Africa when including fossils in the analyses, suggesting an origin of the order likely along the Tethys Seaway where the areas were connected in the mid-Cretaceous. Inferred geographic origins of each family differed to some extent between analyses including fossils vs excluding fossils. In the analysis with extant and fossil species, the origins of the African Hydrostachyaceae and Grubbiaceae-Curtisiaceae clade were inferred to have involved two independent events, an intercontinental dispersal from the northern hemisphere to Africa and an intercontinental vicariance between the northern hemisphere and Africa, respectively. Other families were inferred to have evolved in the northern hemisphere with subsequent intercontinental dispersal(s) to other areas including to Central and South America, during their subsequent diversification. Net diversification rate analysis based on treePL dated phylogeny using MEDUSA detected a nearly 5-fold decrease in the African endemic Curtisiaceae-Grubbiaceae (CuG) clade and an increase of rate in the Hydrangeaceae-Loasaceae (HL) clade. Within HL, a decrease in the Fendlera-Jamesia clade and an increase in the Philadelphus clade were also detected. The findings are also consistent with the level of present species diversity in these lineages. Our study demonstrated the value of plastid genome in phylogenomic study, but posed an old challenge of biogeographic study with fossil data and raised caution for the synonymous substitution sites of plastid genome in phylogenomics studies.

Development of polymorphic microsatellite markers for tree peony Paeonia delavayi (Paeoniaceae) using ddRAD-seq data.

Microsatellite markers were developed for the tree peony Paeonia delavayi to investigate fine scale population genetics of this species. Using ddRAD-seq data from twenty individuals of P. delavayi, we identified 529 polymorphic microsatellite loci, of which 195 were suitable for designing microsatellite primers. Of the 120 microsatellite loci selected for validation, 20 were successfully amplified with clear peaks and displayed polymorphism. Three populations were genotyped using the 20 polymorphic microsatellites. The number of alleles per locus ranged from two to thirteen. Observed and expected heterozygosity ranged from 0 to 0.941 and 0 to 0.834 respectively. The cross-species amplification test using five individuals from a population of P. ludlowii showed that 15 of the 20 polymorphic loci were successfully amplified, and four loci showed polymorphism. Among the 22 alleles occurring in P. ludlowii across fifteen loci, eight alleles across five loci were exclusive to P. ludlowii. The results demonstrate that ddRAD-seq is an efficient method for the development of microsatellite markers for non-model organisms with large genomes. The newly developed markers will be valuable tools to investigate the genetic diversity, genetic structure, and gene flow of P. delavayi from local to regional spatial scales.

Incomplete reproductive isolation between Rhododendron taxa enables hybrid formation and persistence.

The evolutionary consequences of hybridization ultimately depend on the magnitude of reproductive isolation between hybrids and their parents. We evaluated the relative contributions of pre- and post-zygotic barriers to reproduction for hybrid formation, hybrid persistence and potential for reproductive isolation of hybrids formed between two Rhododendron species, R. spiciferum and R. spinuliferum. Our study established that incomplete reproductive isolation promotes hybrid formation and persistence and delays hybrid speciation. All pre-zygotic barriers to reproduction leading to hybrid formation are incomplete: parental species have overlapping flowering; they share the same pollinators; reciprocal assessments of pollen tube germination and growth do not differ among parents. The absence of post-zygotic barriers between parental taxa indicates that the persistence of hybrids is likely. Reproductive isolation was incomplete between hybrids and parents in all cases studied, although asymmetric differences in reproductive fitness were prevalent and possibly explain the genetic structure of natural hybrid swarms where hybridization is known to be bidirectional but asymmetric. Introgression, rather than speciation, is a probable evolutionary outcome of hybridization between the two Rhododendron taxa. Our study provides insights into understanding the evolutionary implications of natural hybridization in woody plants.

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.

Insights into the historical assembly of East Asian subtropical evergreen broadleaved forests revealed by the temporal history of the tea family.

Subtropical evergreen broadleaved forests (EBLFs) inhabit large areas of East Asia. Although paleovegetation reconstructions have revealed that the subtropical EBLFs existed in Southwest China during the Miocene, the historical construction of these forests remains poorly known. Here, we used the tea family (Theaceae), a characteristic component of the subtropical EBLFs, to gain new insights into the assembly of this important biome. Using a robust phylogenetic framework of Theaceae based on plastome and nuclear ribosomal DNA sequence data, the temporal history of the family was reconstructed. Data from other characteristic components of subtropical EBLFs, including Fagaceae, Lauraceae and Magnoliaceae, were also integrated. Most of the essential elements of the subtropical EBLFs appear to have originated around the Oligocene-Miocene (O-M) boundary. However, small woody lineages (e.g. Camellia, Hartia) from Theaceae were dated to the late Miocene. Accelerated net diversification rates within Theaceae were also detected near the O-M transition period and the late Miocene. Our results suggest that two independent intensifications of the East Asian summer monsoon (EASM) around the O-M boundary and the late Miocene may have facilitated the historical assembly of the subtropical EBLFs in East Asia.

Asymmetrical natural hybridization varies among hybrid swarms between two diploid Rhododendron species.

Background and Aims: The extent to which hybridization leads to gene flow between plant species depends on the structure of hybrid populations. However, if this varies between locations, species barriers might prove permeable in some locations but not in others. To assess possible variation in hybrid population structure, the magnitude and direction of natural hybridization between two Chinese endemic species, Rhododendron spiciferum and Rhododendron spinuliferum , were evaluated. Methods: Thirteen nuclear microsatellite markers were employed to characterize 566 individuals collected from 15 non-allopatric populations and nine allopatric parental populations. Chloroplast DNA (cpDNA) sequences were obtained from a subset of samples. Genetic structure and direction of gene flow was determined using a combination of STRUCTURE and NEWHYBRIDS analysis. Key Results: Nuclear analysis revealed that parental taxa formed two genetically distinct clusters and hybrids shared the genetic background of both parents and did not form a separate genetic lineage. Overall, hybrid swarms were dominated by early- and later-generation hybrids, with a significantly higher proportion of hybrids (59·6 %) possessing >50 % R. spiciferum-like nuclear germplasm. The cpDNA analysis further indicated that a significantly greater proportion of hybrids (61·1 %) possessed the R. spiciferum cpDNA haplotype. Conclusions: Gene flow between R. spiciferum and R. spinuliferum was found to be bidirectional in 14 of the 15 hybrid swarms and asymmetrical in six hybrid swarms. Asymmetrical gene flow was evident for only nuclear DNA (nDNA) in two populations, for only cpDNA in three populations, and for both nDNA and cpDNA in one population. Collectively, the variation in genetic structure found among the 15 hybrid swarms suggests that introgression rather than hybrid speciation is a more likely outcome of hybridization between these hybridizing taxa.