Young evolutionary origins of dioecy in the genus Asparagus.

PREMISE: Dioecy (separate sexes) has independently evolved numerous times across the angiosperm phylogeny and is recently derived in many lineages. However, our understanding is limited regarding the evolutionary mechanisms that drive the origins of dioecy in plants. The recent and repeated evolution of dioecy across angiosperms offers an opportunity to make strong inferences about the ecological, developmental, and molecular factors influencing the evolution of dioecy, and thus sex chromosomes. The genus Asparagus (Asparagaceae) is an emerging model taxon for studying dioecy and sex chromosome evolution, yet estimates for the age and origin of dioecy in the genus are lacking. METHODS: We use plastome sequences and fossil time calibrations in phylogenetic analyses to investigate the age and origin of dioecy in the genus Asparagus. We also review the diversity of sexual systems present across the genus to address contradicting reports in the literature. RESULTS: We estimate that dioecy evolved once or twice approximately 2.78-3.78 million years ago in Asparagus, of which roughly 27% of the species are dioecious and the remaining are hermaphroditic with monoclinous flowers. CONCLUSIONS: Our findings support previous work implicating a young age and the possibility of two origins of dioecy in Asparagus, which appear to be associated with rapid radiations and range expansion out of Africa. Lastly, we speculate that paleoclimatic oscillations throughout northern Africa may have helped set the stage for the origin(s) of dioecy in Asparagus approximately 2.78-3.78 million years ago.

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.

High-quality Cymbidium mannii genome and multifaceted regulation of crassulacean acid metabolism in epiphytes.

Epiphytes with crassulacean acid metabolism (CAM) photosynthesis are widespread among vascular plants, and repeated evolution of CAM photosynthesis is a key innovation for micro-ecosystem adaptation. However, we lack a complete understanding of the molecular regulation of CAM photosynthesis in epiphytes. Here, we report a high-quality chromosome-level genome assembly of a CAM epiphyte, Cymbidium mannii (Orchidaceae). The 2.88-Gb orchid genome with a contig N50 of 22.7 Mb and 27 192 annotated genes was organized into 20 pseudochromosomes, 82.8% of which consisted of repetitive elements. Recent expansions of long terminal repeat retrotransposon families have made a major contribution to the evolution of genome size in Cymbidium orchids. We reveal a holistic scenario of molecular regulation of metabolic physiology using high-resolution transcriptomics, proteomics, and metabolomics data collected across a CAM diel cycle. Patterns of rhythmically oscillating metabolites, especially CAM-related products, reveal circadian rhythmicity in metabolite accumulation in epiphytes. Genome-wide analysis of transcript and protein level regulation revealed phase shifts during the multifaceted regulation of circadian metabolism. Notably, we observed diurnal expression of several core CAM genes (especially ßCA and PPC) that may be involved in temporal fixation of carbon sources. Our study provides a valuable resource for investigating post-transcription and translation scenarios in C. mannii, an Orchidaceae model for understanding the evolution of innovative traits in epiphytes.

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.

DNA barcoding of Cymbidium by genome skimming: Call for next-generation nuclear barcodes.

Cymbidium is an orchid genus that has undergone rapid radiation and has high ornamental, economic, ecological and cultural importance, but its classification based on morphology is controversial. The plastid genome (plastome), as an extension of plant standard DNA barcodes, has been widely used as a potential molecular marker for identifying recently diverged species or complicated plant groups. In this study, we newly generated 237 plastomes of 50 species (at least two individuals per species) by genome skimming, covering 71.4% of members of the genus Cymbidium. Sequence-based analyses (barcoding gaps and automatic barcode gap discovery) and tree-based analyses (maximum likelihood, Bayesian inference and multirate Poisson tree processes model) were conducted for species identification of Cymbidium. Our work provides a comprehensive DNA barcode reference library for Cymbidium species identification. The results show that compared with standard DNA barcodes (rbcL + matK) as well as the plastid trnH-psbA, the species identification rate of the plastome increased moderately from 58% to 68%. At the same time, we propose an optimized identification strategy for Cymbidium species. The plastome cannot completely resolve the species identification of Cymbidium, the main reasons being incomplete lineage sorting, artificial cultivation, natural hybridization and chloroplast capture. To further explore the potential use of nuclear data in identifying species, the Skmer method was adopted and the identification rate increased to 72%. It appears that nuclear genome data have a vital role in species identification and are expected to be used as next-generation nuclear barcodes.

[Advances and Prospect of Sampling Techniques and Analytical Methods for Trace Elements in the Ocean: Progress of Trace Element Platform Construction in Xiamen University].

Trace elements, which are important chemical components in the ocean, generally refer to those chemical elements with concentrations below 10 µmol·kg-1in seawater. Some trace elements, such as Fe and Zn, serve as essential micronutrients for marine organisms, which regulate marine primary productivity and are closely related to the biogeochemical cycle of carbon and nitrogen and therefore affect the global environment and climate change. In contrast, some elements, such as Pb, are anthropogenic pollutants largely released by human activities. In addition, some trace elements and their isotopes can be used as tracers for oceanographic processes and proxies for paleoceanography. However, the high saline matrix and extremely low trace element concentrations in seawater, as well as the contamination from research vessels, sampling equipment, and the surrounding environment during the process of sample collection, pretreatment, and analysis, have restricted researchers from obtaining reliable trace element data in the ocean for a long period of time. Nevertheless, high quality samples and accurate data are prerequisites for investigating the biogeochemical and environmental behavior of marine trace elements. This paper reviews the development of sampling techniques and analytical methods for trace elements in seawater, introduces the research history and platform construction activities in Xiamen University in this field, summarizes the advantages and disadvantages of various sampling and analytical techniques and methods, and presents the perspectives on future developments in the research on trace elements in the ocean.

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.

Blood Lead and High-Density Lipoprotein Concentrations in Relation to Human Blood Pressure: A Cross Sectional Study.

Background: While the relationship between blood pressure and blood lead has been studied more extensively, the effect of high-density lipoprotein (HDL) concentration on this relationship remains uncertain. Therefore, this study aimed to determine the effect of HDL concentration on the relationship between blood lead and blood pressure. Methods: The research used cross-sectional data from the 2005 to 2014 National Health and Nutrition Examination Survey (NHANES), which included 16,451 participants aged 20-60 years. Multivariable linear regression was used to evaluate the correlation among blood lead, systolic blood pressure (SBP), and diastolic blood pressure (DBP). HDL concentration was determined by low HDL concentration (≤ 49 mg/dl) and high HDL concentration (>49 mg/dl) stratified. The effect of HDL concentration was assessed by an interaction test between blood lead and SBP in multivariable linear regression. Results: In this cross-sectional research, we identified a positive correlation between blood lead and SBP, but not DBP. The relationship between blood lead and SBP was different in the group with low and high HDL concentrations (ß: 0.21 95% Cl:-0.05-0.46 vs. ß:0.47 95% Cl: 0.15-0.79). In addition, high HDL significantly altered the positive correlation between blood lead and SBP (P-value of interaction < 0.001). Conclusion: The study suggests an interaction between HDL and blood lead in elevating SBP, which may have important clinical implications.

Species discrimination in Schima (Theaceae): Next-generation super-barcodes meet evolutionary complexity.

Plastid genome and nuclear ribosomal DNA (nrDNA) arrays, proposed recently as "super-barcodes," might provide additional discriminatory power and overcome the limitations of traditional barcoding loci, yet super-barcodes need to be tested for their effectiveness in more plant groups. Morphological homoplasy among Schima species makes the genus a model for testing the efficacy of super-barcodes. In this study, we generated multiple data sets comprising standard DNA barcodes (matK, rbcL, trnH-psbA, nrITS) and super-barcodes (plastid genome, nrDNA arrays) across 58 individuals from 12 out of 13 species of Schima from China. No samples were correctly assigned to species using standard DNA barcodes and nrDNA arrays, while only 27.27% of species with multiple accessions were distinguished using the plastid genome and its partitioned data sets-the lowest estimated rate of super-barcode success in the literature so far. For Schima and other taxa with similarly recently divergence and low levels of genetic variation, incomplete lineage sorting, hybridization or taxonomic oversplitting are all possible causes of the failure. Taken together, our study suggests that by no means are super-barcodes immune to the challenges imposed by evolutionary complexity. We therefore call for developing multilocus nuclear markers for species discrimination in plant groups.

Testing the Complete Plastome for Species Discrimination, Cryptic Species Discovery and Phylogenetic Resolution in Cephalotaxus (Cephalotaxaceae).

Species of Cephalotaxus have great economic and ecological values. However, the taxonomy and interspecific phylogenetic relationships within the genus have been controversial and remained not fully resolved until now. To date, no study examined the efficiency of the complete plastome as super-barcode across Cephalotaxus species with multiple samples per taxon. In this study, we have evaluated the complete plastome in species discrimination and phylogenetic resolution in Cephalotaxus by including 32 individuals of all eight recognized species and five varieties following Farjon's classification (2010) with multiple samples per taxon. Our results indicated that not all species recognized in recent taxonomic revisions of Cephalotaxus could be distinguished and not all were monophyletic. Based on the plastome phylogeny, a new taxonomic classification for the genus comprising nine species and two varieties, including a cryptic species, was proposed. The phylogeny also resolved all interspecific relationships. Compared to the plastome based classification, standard DNA barcodes, alone or in combination, only recognized a maximum of seven out of the nine species. Moreover, two highly variable single loci, ycf1 and rps16, each alone achieved full species discrimination. With the moderate length of 1079 bp, rps16 is proposed as a specific barcode to discriminate Cephalotaxus species. The super-barcodes and specific barcode candidates will aid in the identification of endangered Cephalotaxus species, and to help focus conservation measures.

Genome Skimming Contributes to Clarifying Species Limits in Paris Section Axiparis (Melanthiaceae).

Paris L. section Axiparis H. Li (Melanthiaceae) is a taxonomically perplexing taxon with considerable confusion regarding species delimitation. Based on the analyses of morphology and geographic distribution of each species currently recognized in the taxon, we propose a revision scheme that reduces the number of species in P. sect. Axiparis from nine to two. To verify this taxonomic proposal, we employed a genome skimming approach to recover the plastid genomes (plastomes) and nuclear ribosomal DNA (nrDNA) regions of 51 individual plants across the nine described species of P. sect. Axiparis by sampling multiple accessions per species. The species boundaries within P. sect. Axiparis were explored using phylogenetic inference and three different sequence-based species delimitation methods (ABGD, mPTP, and SDP). The mutually reinforcing results indicate that there are two species-level taxonomic units in P. sect. Axiparis (Paris forrestii s.l. and P. vaniotii s.l.) that exhibit morphological uniqueness, non-overlapping distribution, genetic distinctiveness, and potential reproductive isolation, providing strong support to the proposed species delimitation scheme. This study confirms that previous morphology-based taxonomy overemphasized intraspecific and minor morphological differences to delineate species boundaries, therefore resulting in an overestimation of the true species diversity of P. sect. Axiparis. The findings clarify species limits and will facilitate robust taxonomic revision in P. sect. Axiparis.

Cryptic Species Diversification of the Pedicularis siphonantha Complex (Orobanchaceae) in the Mountains of Southwest China Since the Pliocene.

Morphological approaches often fail to delimit species in recently derived species complexes. This can be exacerbated in historical collections which may have lost key features in specimen preparation and preservation. Here, we examine the Pedicularis siphonantha complex, endemic to the Mountains of Southwest China. This complex is characterized by its red/purple/pink and long-tubular corolla, and twisted, beaked galea. However, herbarium specimens are often difficult to identify to species. Molecular approaches using nrITS or nuclear ribosomal internal transcribed spacer (nrITS) + plastid DNA (ptDNA) have been successfully used for species identification in Pedicularis. To resolve taxonomic confusion in the Pedicularis siphonantha complex, we reconstructed phylogenies of the complex using nrITS and four plastid DNA loci (matK, rbcL, trnH-psbA, and trnL-F). To recover as much of the phylogenetic history as possible, we sampled individuals at the population level. Topological incongruence between the nrITS and ptDNA datasets was recovered in clades including two widely distributed species, Pedicularis milliana and Pedicularis tenuituba. Based on morphological, geographical, and genetic evidence, we suggest that hybridization/introgression has occurred between P. milliana and Pedicularis sigmoidea/Pedicularis sp. 1 in the Yulong Snow Mountain of Lijiang, northwest Yunnan, and between P. tenuituba and Pedicularis leptosiphon in Ninglang, northwest Yunnan. After removing conflicting DNA regions in Pedicularis dolichosiphon (nrITS) and P. milliana (ptDNA), the concatenated nrITS and ptDNA phylogenies distinguish 11 species in the P. siphonantha complex, including two undescribed species, from the Jiaozi and Yulong Snow Mountains, respectively. Phylogeographical analyses indicate that the P. siponantha complex originated from south of the Hengduan Mountains, expanding north to the Himalayas and the Yunnan-Guizhou Plateau. Moreover, the uplift of the Qinghai-Tibet Plateau and climate oscillations may have driven further diversification in the complex.

Testing genome skimming for species discrimination in the large and taxonomically difficult genus Rhododendron.

Standard plant DNA barcodes based on 2-3 plastid regions, and nrDNA ITS show variable levels of resolution, and fail to discriminate among species in many plant groups. Genome skimming to recover complete plastid genome sequences and nrDNA arrays has been proposed as a solution to address these resolution limitations. However, few studies have empirically tested what gains are achieved in practice. Of particular interest is whether adding substantially more plastid and nrDNA characters will lead to an increase in discriminatory power, or whether the resolution limitations of standard plant barcodes are fundamentally due to plastid genomes and nrDNA not tracking species boundaries. To address this, we used genome skimming to recover near-complete plastid genomes and nuclear ribosomal DNA from Rhododendron species and compared discrimination success with standard plant barcodes. We sampled 218 individuals representing 145 species of this species-rich and taxonomically difficult genus, focusing on the global biodiversity hotspots of the Himalaya-Hengduan Mountains. Only 33% of species were distinguished using ITS+matK+rbcL+trnH-psbA. In contrast, 55% of species were distinguished using plastid genome and nrDNA sequences. The vast majority of this increase is due to the additional plastid characters. Thus, despite previous studies showing an asymptote in discrimination success beyond 3-4 plastid regions, these results show that a demonstrable increase in discriminatory power is possible with extensive plastid genome data. However, despite these gains, many species remain unresolved, and these results also reinforce the need to access multiple unlinked nuclear loci to obtain transformative gains in species discrimination in plants.

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.

Organelle Genomes and Transcriptomes of Nymphaea Reveal the Interplay between Intron Splicing and RNA Editing.

Posttranscriptional modifications, including intron splicing and RNA editing, are common processes during regulation of gene expression in plant organelle genomes. However, the intermediate products of intron-splicing, and the interplay between intron-splicing and RNA-editing were not well studied. Most organelle transcriptome analyses were based on the Illumina short reads which were unable to capture the full spectrum of transcript intermediates within an organelle. To fully investigate the intermediates during intron splicing and the underlying relationships with RNA editing, we used PacBio DNA-seq and Iso-seq, together with Illumina short reads genome and transcriptome sequencing data to assemble the chloroplast and mitochondrial genomes of Nymphaea 'Joey Tomocik' and analyze their posttranscriptional features. With the direct evidence from Iso-seq, multiple intermediates partially or fully intron-spliced were observed, and we also found that both cis- and trans-splicing introns were spliced randomly. Moreover, by using rRNA-depleted and non-Oligo(dT)-enrichment strand-specific RNA-seq data and combining direct SNP-calling and transcript-mapping methods, we identified 98 and 865 RNA-editing sites in the plastome and mitogenome of N. 'Joey Tomocik', respectively. The target codon preference, the tendency of increasing protein hydrophobicity, and the bias distribution of editing sites are similar in both organelles, suggesting their common evolutionary origin and shared editing machinery. The distribution of RNA editing sites also implies that the RNA editing sites in the intron and exon regions may splice synchronously, except those exonic sites adjacent to intron which could only be edited after being intron-spliced. Our study provides solid evidence for the multiple intermediates co-existing during intron-splicing and their interplay with RNA editing in organelle genomes of a basal angiosperm.

Plastome Evolution in the Hyperdiverse Genus Euphorbia (Euphorbiaceae) Using Phylogenomic and Comparative Analyses: Large-Scale Expansion and Contraction of the Inverted Repeat Region.

With c. 2,000 species, Euphorbia is one of the largest angiosperm genera, yet a lack of chloroplast genome (plastome) resources impedes a better understanding of its evolution. In this study, we assembled and annotated 28 plastomes from Euphorbiaceae, of which 15 were newly sequenced. Phylogenomic and comparative analyses of 22 plastome sequences from all four recognized subgenera within Euphorbia revealed that plastome length in Euphorbia is labile, presenting a range of variation c. 42 kb. Large-scale expansions of the inverted repeat (IR) region were identified, and at the extreme opposite, the near-complete loss of the IR region (with only 355 bp left) was detected for the first time in Euphorbiaceae. Other structural variations, including gene inversion and duplication, and gene loss/pseudogenization, were also observed. We screened the most promising molecular markers from both intergenic and coding regions for phylogeny-based utilities, and estimated maximum likelihood and Bayesian phylogenies from four datasets including whole plastome sequences. The monophyly of Euphorbia is supported, and its four subgenera are recovered in a successive sister relationship. Our study constitutes the first comprehensive investigation on the plastome structural variation in Euphorbia and it provides resources for phylogenetic research in the genus, facilitating further studies on its taxonomy, evolution, and conservation.

Gene duplications and phylogenomic conflict underlie major pulses of phenotypic evolution in gymnosperms.

Inferring the intrinsic and extrinsic drivers of species diversification and phenotypic disparity across the tree of life is a major challenge in evolutionary biology. In green plants, polyploidy (or whole-genome duplication, WGD) is known to play a major role in microevolution and speciation, but the extent to which WGD has shaped macroevolutionary patterns of diversification and phenotypic innovation across plant phylogeny remains an open question. Here, we examine the relationship of various facets of genomic evolution-including gene and genome duplication, genome size, and chromosome number-with macroevolutionary patterns of phenotypic innovation, species diversification, and climatic occupancy in gymnosperms. We show that genomic changes, such as WGD and genome-size shifts, underlie the origins of most major extant gymnosperm clades, and notably, our results support an ancestral WGD in the gymnosperm lineage. Spikes of gene duplication typically coincide with major spikes of phenotypic innovation, while increased rates of phenotypic evolution are typically found at nodes with high gene-tree conflict, representing historic population-level dynamics during speciation. Most shifts in gymnosperm diversification since the rise of angiosperms are decoupled from putative WGDs and instead are associated with increased rates of climatic occupancy evolution, particularly in cooler and/or more arid climatic conditions, suggesting that ecological opportunity, especially in the later Cenozoic, and environmental heterogeneity have driven a resurgence of gymnosperm diversification. Our study provides critical insight on the processes underlying diversification and phenotypic evolution in gymnosperms, with important broader implications for the major drivers of both micro- and macroevolution in plants.

Plastid NDH Pseudogenization and Gene Loss in a Recently Derived Lineage from the Largest Hemiparasitic Plant Genus Pedicularis (Orobanchaceae).

The plastid genome (plastome) is highly conserved in both gene order and content and has a lower mutation rate than the nuclear genome. However, the plastome is more variable in heterotrophic plants. To date, most such studies have investigated just a few species or only holoheterotrophic groups, and few have examined plastome evolution in recently derived lineages at an early stage of transition from autotrophy to heterotrophy. In this study, we investigated the evolutionary dynamics of plastomes in the monophyletic and recently derived Pedicularis sect. Cyathophora (Orobanchaceae). We obtained 22 new plastomes, 13 from the six recognized species of section Cyathophora, six from hemiparasitic relatives and three from autotrophic relatives. Comparative analyses of gene content, plastome structure and selection pressure showed dramatic differences among species in section Cyathophora and in Pedicularis as a whole. In comparison with autotrophic relatives and other Pedicularis spp., we found that the inverted repeat (IR) region in section Cyathophora had expansions to the small single-copy region, with a large expansion event and two independent contraction events. Moreover, NA(D)H dehydrogenase, accD and ccsA have lost function multiple times, with the function of accD being replaced by nuclear copies of an accD-like gene in Pedicularis spp. The ccsA and ndhG genes may have evolved under selection in association with IR expansion/contraction events. This study is the first to report high plastome variation in a recently derived lineage of hemiparasitic plants and therefore provides evidence for plastome evolution in the transition from autotrophy to heterotrophy.

Organelle Phylogenomics and Extensive Conflicting Phylogenetic Signals in the Monocot Order Poales.

The Poales is one of the largest orders of flowering plants with significant economic and ecological values. Reconstructing the phylogeny of the Poales is important for understanding its evolutionary history that forms the basis for biological studies. However, due to sparse taxon sampling and limited molecular data, previous studies have resulted in a variety of contradictory topologies. In particular, there are three nodes surrounded by incongruence: the phylogenetic ambiguity near the root of the Poales tree, the sister family of Poaceae, and the delimitation of the xyrid clade. We conducted a comprehensive sampling and reconstructed the phylogenetic tree using plastid and mitochondrial genomic data from 91 to 66 taxa, respectively, representing all the 16 families of Poales. Our analyses support the finding of Bromeliaceae and Typhaceae as the earliest diverging groups within the Poales while having phylogenetic relationships with the polytomy. The clade of Ecdeiocoleaceae and Joinvilleaceae is recovered as the sister group of Poaceae. The three families, Mayacaceae, Eriocaulaceae, and Xyridaceae, of the xyrid assembly diverged successively along the backbone of the Poales phylogeny, and thus this assembly is paraphyletic. Surprisingly, we find substantial phylogenetic conflicts within the plastid genomes of the Poales, as well as among the plastid, mitochondrial, and nuclear data. These conflicts suggest that the Poales could have a complicated evolutionary history, such as rapid radiation and polyploidy, particularly allopolyploidy through hybridization. In sum, our study presents a new perspicacity into the complex phylogenetic relationships and the underlying phylogenetic conflicts within the Poales.

De novo assembling a complete mitochondrial genome of Pedicularis rex (Orobanchaceae) using GetOrganelle toolkit.

We reported the first mitogenome of Pedicularis from P. rex (Orobanchaceae), which is endemic to SW China. The complete mitochondrial genome (mitogenome or chondriome) was a single circular chromosome that was 219,859 bp in length. It contains 56 genes, including 34 protein-coding (cox2 and atp9 with two copies), 19 transfer RNA (tRNA), and three ribosomal RNA (rRNA) genes. Phylogenetic analysis showed that Pedicularis rex was closely related to Castilleja paramensis.