Unprecedented variation pattern of plastid genomes and the potential role in adaptive evolution in Poales

BACKGROUND

The plastid is the photosynthetic organelle in plant cell, and the plastid genomes (plastomes) are generally conserved in evolution. As one of the most economically and ecologically important order of angiosperms, Poales was previously documented to exhibit great plastomic variation as an order of photoautotrophic plants.

RESULTS

We acquired 93 plastomes, representing all the 16 families and 5 major clades of Poales to reveal the extent of their variation and evolutionary pattern. Extensive variation including the largest one in monocots with 225,293 bp in size, heterogeneous GC content, and a wide variety of gene duplication and loss were revealed. Moreover, rare occurrences of three inverted repeat (IR) copies in angiosperms and one IR loss were observed, accompanied by short IR (sIR) and small direct repeat (DR). Widespread structural heteroplasmy, diversified inversions, and unusual genomic rearrangements all appeared in Poales, occasionally within a single species. Extensive repeats in the plastomes were found to be positively correlated with the observed inversions and rearrangements. The variation all showed a "small-large-moderate" trend along the evolution of Poales, as well as for the sequence substitution rate. Finally, we found some positively selected genes, mainly in C4 lineages, while the closely related lineages of those experiencing gene loss tended to have undergone more relaxed purifying selection.

CONCLUSIONS

The variation of plastomes in Poales may be related to its successful diversification into diverse habitats and multiple photosynthetic pathway transitions. Our order-scale analyses revealed unusual evolutionary scenarios for plastomes in the photoautotrophic order of Poales and provided new insights into the plastome evolution in angiosperms as a whole.

Phylogenomics and the rise of the angiosperms

Angiosperms are the cornerstone of most terrestrial ecosystems and human livelihoods1,2. A robust understanding of angiosperm evolution is required to explain their rise to ecological dominance. So far, the angiosperm tree of life has been determined primarily by means of analyses of the plastid genome3,4. Many studies have drawn on this foundational work, such as classification and first insights into angiosperm diversification since their Mesozoic origins5-7. However, the limited and biased sampling of both taxa and genomes undermines confidence in the tree and its implications. Here, we build the tree of life for almost 8,000 (about 60%) angiosperm genera using a standardized set of 353 nuclear genes8. This 15-fold increase in genus-level sampling relative to comparable nuclear studies9 provides a critical test of earlier results and brings notable change to key groups, especially in rosids, while substantiating many previously predicted relationships. Scaling this tree to time using 200 fossils, we discovered that early angiosperm evolution was characterized by high gene tree conflict and explosive diversification, giving rise to more than 80% of extant angiosperm orders. Steady diversification ensued through the remaining Mesozoic Era until rates resurged in the Cenozoic Era, concurrent with decreasing global temperatures and tightly linked with gene tree conflict. Taken together, our extensive sampling combined with advanced phylogenomic methods shows the deep history and full complexity in the evolution of a megadiverse clade.

Evolution of the FLOWERING LOCUS T-like genes in angiosperms: a core-Lamiales-specific diversification

Plant life-history is determined by two transitions, the germination and the flowering times, in which the phosphatidylethanolamine-binding proteins (PEBP) FLOWERING LOCUS T (FT) and TERMINAL FLOWER1 (TFL1) play key regulatory roles. Compared to the highly conserved TFL1-likes, FT-like genes vary in copy numbers significantly in gymnosperms and monocots of the angiosperms, while sporadic duplications can be observed in eudicots. Here, via a systematic analysis of the PEBPs in angiosperms with a special focus on twelve representative species featuring high-quality genomes in the Lamiales order, we identified a successive lineage-specific but systematic expansion of FT-like genes in the families of core Lamiales. The first expansion event generated FT1-likes mainly via a core-Lamiales-specific whole-genome-duplication (cL-WGD), while on the other hand, a likely random duplication produced the FT2-likes in the lineages containing Scrophulariaceae and rest of the core Lamiales. Both FT1- and FT2-like genes were further amplified tandemly in some families. These expanded FT-likes featured highly diverged expression patterns and structural variation, indicating functional diversification. Intriguingly, some core Lamiales contained the relict MOTHER OF FT AND TFL1 like 2 (MFT2) that likely expanded in the common ancestor of angiosperms. Our data showcase the highly dynamic lineage-specific expansion of the FT-like genes, thus provide important and fresh evolutionary insights into the gene-regulatory-network underpinning flowering time diversity in Lamiales, and more generally, in angiosperms.

Phylogenomic profiles of whole-genome duplications in Poaceae and landscape of differential duplicate retention and losses among major Poaceae lineages

Poaceae members shared a whole-genome duplication called rho. However, little is known about the evolutionary pattern of the rho-derived duplicates among Poaceae lineages and implications in adaptive evolution. Here we present phylogenomic/phylotranscriptomic analyses of 363 grasses covering all 12 subfamilies and report nine previously unknown whole-genome duplications. Furthermore, duplications from a single whole-genome duplication were mapped to multiple nodes on the species phylogeny; a whole-genome duplication was likely shared by woody bamboos with possible gene flow from herbaceous bamboos; and recent paralogues of a tetraploid Oryza are implicated in tolerance of seawater submergence. Moreover, rho duplicates showing differential retention among subfamilies include those with functions in environmental adaptations or morphogenesis, including ACOT for aquatic environments (Oryzoideae), CK2β for cold responses (Pooideae), SPIRAL1 for rapid cell elongation (Bambusoideae), and PAI1 for drought/cold responses (Panicoideae). This study presents a Poaceae whole-genome duplication profile with evidence for multiple evolutionary mechanisms that contribute to gene retention and losses.

Variations and reduction of plastome are associated with the evolution of parasitism in Convolvulaceae

Parasitic lifestyle can often relax the constraint on the plastome, leading to gene pseudogenization and loss, and resulting in diverse genomic structures and rampant genome degradation. Although several plastomes of parasitic Cuscuta have  been reported, the evolution of parasitism in the family Convolvulaceae which is linked to structural variations and reduction of plastome has not been well investigated. In this study, we assembled and collected 40 plastid genomes belonging to 23 species representing four subgenera of Cuscuta and ten species of autotrophic Convolvulaceae. Our findings revealed nine types of structural variations and six types of inverted repeat (IR) boundary variations in the plastome of Convolvulaceae spp. These structural variations were associated with the shift of parasitic lifestyle, and IR boundary shift, as well as the abundance of long repeats. Overall, the degradation of Cuscuta plastome proceeded gradually, with one clade exhibiting an accelerated degradation rate. We observed five stages of gene loss in Cuscuta, including NAD(P)H complex → PEP complex → Photosynthesis-related → Ribosomal protein subunits → ATP synthase complex. Based on our results, we speculated that the shift of parasitic lifestyle in early divergent time promoted relaxed selection on plastomes, leading to the accumulation of microvariations, which ultimately resulted in the plastome reduction. This study provides new evidence towards a better understanding of plastomic evolution, variation, and reduction in the genus Cuscuta.

[The efficacy and safety of immunotherapy combined with chemotherapy neoadjuvant in locally advanced resectable hypopharyngeal squamous cell carcinoma]

Objective: To explore the efficacy and safety of immunoneoadjuvant therapy with pembrolizumab combined with chemotherapy in locally advanced resectable hypopharyngeal squamous cell carcinoma patients. Methods: This study was a prospective, single arm, single center clinical study that was opened for enrollment in April 2021. Patients who met the inclusion criteria at the Cancer Hospital of the Chinese Academy of Medical Sciences were treated with neoadjuvant therapy of pembrolizumab combined with cisplatin and paclitaxel, and after treatments, received surgery and postoperative adjuvant therapy. The main endpoint of this study was postoperative pathological complete response (pCR), and other observations included adverse reactions and long-term prognoses of patients after neoadjuvant therapy. Results: By September 2023, a total of 23 patients who underwent neoadjuvant therapy and surgery were enrolled in the study and all patients were males aged 49-74 years. All patients were locally advanced stage, including 3 patients in stage Ⅲ and 20 patients in stage Ⅳ. There were 12 cases of primary lesions with posterior ring involvement accompanied by fixation of one vocal cord and 20 cases of regional lymph node metastases classified as N2. Eighteen cases received a two cycle regimen and 5 cases received a three cycle regimen for neoadjuvant therapy. The postoperative pCR rate was 26.1% (6/23), with no surgical delay caused by adverse drug reactions. The laryngeal preservation rate was 87.0% (20/23). Pharyngeal fistula was the main surgical complication, with an incidence of 21.7% (5/23). The median follow-up time was 15 months, and 3 patients experienced local recurrence. Conclusions: The immunoneoadjuvant therapy of pembrolizumab combined with chemotherapy has a high pCR rate in locally advanced resectable hypopharyngeal squamous cell carcinoma, with increased laryngeal preservation rate and no significant impact on surgical safety.

Parallel evolution of morphological and genomic selfing syndromes accompany the breakdown of heterostyly

Evolutionary transitions from outcrossing to selfing in flowering plants have convergent morphological and genomic signatures and can involve parallel evolution within related lineages. Adaptive evolution of morphological traits is often assumed to evolve faster than nonadaptive features of the genomic selfing syndrome. We investigated phenotypic and genomic changes associated with transitions from distyly to homostyly in the Primula oreodoxa complex. We determined whether the transition to selfing occurred more than once and investigated stages in the evolution of morphological and genomic selfing syndromes using 22 floral traits and both nuclear and plastid genomic data from 25 populations. Two independent transitions were detected representing an earlier and a more recently derived selfing lineage. The older lineage exhibited classic features of the morphological and genomic selfing syndrome. Although features of both selfing syndromes were less developed in the younger selfing lineage, they exhibited parallel development with the older selfing lineage. This finding contrasts with the prediction that some genomic changes should lag behind adaptive changes to morphological traits. Our findings highlight the value of comparative studies on the timing and extent of transitions from outcrossing to selfing between related lineages for investigating the tempo of morphological and molecular evolution.

Genome assemblies of 11 bamboo species highlight diversification induced by dynamic subgenome dominance

Polyploidy (genome duplication) is a pivotal force in evolution. However, the interactions between parental genomes in a polyploid nucleus, frequently involving subgenome dominance, are poorly understood. Here we showcase analyses of a bamboo system (Poaceae: Bambusoideae) comprising a series of lineages from diploid (herbaceous) to tetraploid and hexaploid (woody), with 11 chromosome-level de novo genome assemblies and 476 transcriptome samples. We find that woody bamboo subgenomes exhibit stunning karyotype stability, with parallel subgenome dominance in the two tetraploid clades and a gradual shift of dominance in the hexaploid clade. Allopolyploidization and subgenome dominance have shaped the evolution of tree-like lignified culms, rapid growth and synchronous flowering characteristic of woody bamboos as large grasses. Our work provides insights into genome dominance in a remarkable polyploid system, including its dependence on genomic context and its ability to switch which subgenomes are dominant over evolutionary time.

Genomes of Meniocus linifolius and Tetracme quadricornis reveal the ancestral karyotype and genomic features of core Brassicaceae

Brassicaceae represents an important plant family from both a scientific and economic perspective. However, genomic features related to the early diversification of this family have not been fully characterized, especially upon the uplift of the Tibetan Plateau, which was followed by increasing aridity in the Asian interior, intensifying monsoons in Eastern Asia, and significantly fluctuating daily temperatures. Here, we reveal the genomic architecture that accompanied early Brassicaceae diversification by analyzing two high-quality chromosome-level genomes for Meniocus linifolius (Arabodae; clade D) and Tetracme quadricornis (Hesperodae; clade E), together with genomes representing all major Brassicaceae clades and the basal Aethionemeae. We reconstructed an ancestral core Brassicaceae karyotype (CBK) containing 9 pseudochromosomes with 65 conserved syntenic genomic blocks and identified 9702 conserved genes in Brassicaceae. We detected pervasive conflicting phylogenomic signals accompanied by widespread ancient hybridization events, which correlate well with the early divergence of core Brassicaceae. We identified a successive Brassicaceae-specific expansion of the class I TREHALOSE-6-PHOSPHATE SYNTHASE 1 (TPS1) gene family, which encodes enzymes with essential regulatory roles in flowering time and embryo development. The TPS1s were mainly randomly amplified, followed by expression divergence. Our results provide fresh insights into historical genomic features coupled with Brassicaceae evolution and offer a potential model for broad-scale studies of adaptive radiation under an ever-changing environment.

Positive selection and relaxed purifying selection contribute to rapid evolution of male-biased genes in a dioecious flowering plant

Sex-biased genes offer insights into the evolution of sexual dimorphism. Sex-biased genes, especially those with male bias, show elevated evolutionary rates of protein sequences driven by positive selection and relaxed purifying selection in animals. Although rapid sequence evolution of sex-biased genes and evolutionary forces have been investigated in animals and brown algae, less is known about evolutionary forces in dioecious angiosperms. In this study, we separately compared the expression of sex-biased genes between female and male floral buds and between female and male flowers at anthesis in dioecious Trichosanthes pilosa (Cucurbitaceae). In floral buds, sex-biased gene expression was pervasive, and had significantly different roles in sexual dimorphism such as physiology. We observed higher rates of sequence evolution for male-biased genes in floral buds compared to female-biased and unbiased genes. Male-biased genes under positive selection were mainly associated with functions to abiotic stress and immune responses, suggesting that high evolutionary rates are driven by adaptive evolution. Additionally, relaxed purifying selection may contribute to accelerated evolution in male-biased genes generated by gene duplication. Our findings, for the first time in angiosperms, suggest evident rapid evolution of male-biased genes, advance our understanding of the patterns and forces driving the evolution of sexual dimorphism in dioecious plants.

Anthropogenic disturbance driving population decline of a dominant tree in East Asia evergreen broadleaved forests over the last 11,000 years

Current biodiversity loss is generally considered to have been caused by anthropogenic disturbance, but it is unclear when anthropogenic activities began to affect biodiversity loss. One hypothesis suggests it began with the Industrial Revolution, whereas others propose that anthropogenic disturbance has been associated with biodiversity decline since the early Holocene. To test these hypotheses, we examined the unique vegetation of evergreen broadleaved forests (EBLFs) in East Asia, where humans have affected landscapes since the early Holocene. We adopted a genomic approach to infer the demographic history of a dominant tree (Litsea elongata) of EBLFs. We used Holocene temperature and anthropogenic disturbance factors to calculate the correlation between these variables and the historical effective population size of L. elongata with Spearman statistics and integrated the maximum-entropy niche model to determine the impact of climate change and anthropogenic disturbance on fluctuation in its effective population size. We identified 9 well-defined geographic clades for the populations of L. elongata. Based on the estimated historical population sizes of these clades, all the populations contracted, indicating persistent population decline over the last 11,000 years. Demographic history of L. elongata and human population change, change in cropland use, and change in irrigated rice area were significantly negatively correlated, whereas climate change in the Holocene was not correlated with demographic history. Our results support the early human impact hypothesis and provide comprehensive evidence that early anthropogenic disturbance may contribute to the current biodiversity crisis in East Asia.

Remarkable mitochondrial genome heterogeneity in Meniocus linifolius (Brassicaceae)

KEY MESSAGE

Detailed analyses of 16 genomes identified a remarkable acceleration of mutation rate, hence mitochondrial sequence and structural heterogeneity, in Meniocus linifolius (Brassicaceae). The powerhouse, mitochondria, in plants feature high levels of structural variation, while the encoded genes are normally conserved. However, the substitution rates and spectra of mitochondria DNA within the Brassicaceae, a family with substantial scientific and economic importance, have not been adequately deciphered. Here, by analyzing three newly assembled and 13 known mitochondrial genomes (mitogenomes), we report the highly variable genome structure and mutation rates in Brassicaceae. The genome sizes and GC contents are 196,604 bp and 46.83%, 288,122 bp and 44.79%, and 287,054 bp and 44.93%, for Meniocus linifolius (Mli), Crucihimalaya lasiocarpa (Cla), and Lepidium sativum (Lsa), respectively. In total, 29, 33, and 34 protein-coding genes (PCGs) and 14, 18, and 18 tRNAs are annotated for Mli, Cla, and Lsa, respectively, while all mitogenomes contain one complete circular molecule with three rRNAs and abundant RNA editing sites. The Mli mitogenome features four conformations likely mediated by the two pairs of long repeats, while at the same time seems to have an unusual evolutionary history due to higher GC content, loss of more genes and sequences, but having more repeats and plastid DNA insertions. Corroborating with these, an ambiguous phylogenetic position with long branch length and elevated synonymous substitution rate in nearly all PCGs are observed for Mli. Taken together, our results reveal a high level of mitogenome heterogeneity at the family level and provide valuable resources for further understanding the evolutionary pattern of organelle genomes in Brassicaceae.

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.

A comparative full-length transcriptomic resource provides insight into the perennial monocarpic mass flowering

Perennial monocarpic mass flowering represents as a key developmental innovation in flowering time diversity in several biological and economical essential families, such as the woody bamboos and the shrubby Strobilanthes. However, molecular and genetic mechanisms underlying this important biodiversity remain poorly investigated. Here, we generated a full-length transcriptome resource incorporated into the BlueOmics database (http://blueomics.iflora.cn) for two Strobilanthes species, which feature contrasting flowering time behaviors. Using about 112 and 104 Gb Iso-seq reads together with ~185 and ~75 Gb strand-specific RNA seq data, we annotated 80 971 and 79 985 non-redundant full-length transcripts for the perennial polycarpic Strobilanthes tetrasperma and the perennial monocarpic Strobilanthes biocullata, respectively. In S. tetrasperma, we identified 8794 transcripts showing spatiotemporal expression in nine tissues. In leaves and shoot apical meristems at two developmental stages, 977 and 1121 transcripts were differentially accumulated in S. tetrasperma and S. biocullata, respectively. Interestingly, among the 33 transcription factors showing differential expression in S. tetrasperma but without differential expression in S. biocullata, three were involved potentially in the photoperiod and circadian-clock pathway of flowering time regulation (FAR1 RELATED SEQUENCE 12, FRS12; NUCLEAR FACTOR Y A1, NFYA1; PSEUDO-RESPONSE REGULATOR 5, PRR5), hence provides an important clue in deciphering the flowering diversity mechanisms. Our data serve as a key resource for further dissection of molecular and genetic mechanisms underpinning key biological innovations, here, the perennial monocarpic mass flowering.

Dynamics of imprinted genes and their epigenetic mechanisms in castor bean seed with persistent endosperm

Genomic imprinting refers to parent-of-origin-dependent gene expression and primarily occurs in the endosperm of flowering plants, but its functions and epigenetic mechanisms remain to be elucidated in eudicots. Castor bean, a eudicot with large and persistent endosperm, provides an excellent system for studying the imprinting. Here, we identified 131 imprinted genes in developing endosperms and endosperm at seed germination phase of castor bean, involving into the endosperm development, accumulation of storage compounds and specially seed germination. Our results showed that the transcriptional repression of maternal allele of DNA METHYLTRANSFERASE 1 (MET1) may be required for maternal genome demethylation in the endosperm. DNA methylation analysis showed that only a small fraction of imprinted genes was associated with allele-specific DNA methylation, and most of them were closely associated with constitutively unmethylated regions (UMRs), suggesting a limited role for DNA methylation in controlling genomic imprinting. Instead, histone modifications can be asymmetrically deposited in maternal and paternal genomes in a DNA methylation-independent manner to control expression of most imprinted genes. These results expanded our understanding of the occurrence and biological functions of imprinted genes and showed the evolutionary flexibility of the imprinting machinery and mechanisms in plants.

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.

Five new and noteworthy species of Epidendroideae (Orchidaceae) from southwestern China based on morphological and phylogenetic evidence

Five new orchid species from southwestern China's Yunnan Province and the Tibetan Autonomous Region, Neottialihengiae, Neottiachawalongensis, Papilionanthemotuoensis, Gastrochiluslihengiae, and Gastrochilusbernhardtianus, are described and illustrated. To confirm their identities, and to resolve phylogenetic relationships, we sequenced the complete plastomes of these taxa with their congeneric species, adding new plastomes of three Neottia species, two Papilionanthe species and nine Gastrochilus species. Combined with published plastid sequences, our well-resolved phylogeny supported the alliance of N.lihengiae with the the N.grandiflora + N.pinetorum clade. Neottiachawalongensis is now sister to N.alternifolia, while P.motuoensis is closely related to P.subulata + P.teres. Conversely, phylogenetic analyses based on complete plastomes and plastid sequences showed inconsistent relationships among taxa in the genus Gastrochilus, but the two new species, G.lihengiae and G.bernhardtianus were supported by all datasets.

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.

[Exploration of the clinical application of combined endoscopic and laparoscopic surgery in early gastric cancer: 15 cases]

Objective: To investigate the application of combined gastroscopy and laparoscopy (dual scope) in the treatment of early gastric cancer. Methods: In this descriptive case series study, we retrospectively collected data on 15 patients with cT1b stage gastric cancer who had undergone combined laparoscopic and endoscopic surgery in the 900th Hospital of the People's Liberation Army of China from May 2020 to October 2022. The study cohort comprised nine men and six women of median age 59 (range: 47-76) years and median body mass index 20.9 (range: 18.3-26.2) kg/m2. Seven of the lesions were located on the lesser curvature of the gastric antrum and eight in the gastric angle. All lesions were biopsied for pathological examination and evaluated by endoscopic ultrasonography, followed by endoscopic submucosal dissection (ESD) and laparoscopic regional lymph node dissection. Studied variables included surgical and pathological features, postoperative factors, and outcomes. Results: In this group of patients, the median (range) operative time for ESD was 45 (30-82) minutes, the duration of laparoscopic lymph node dissection (45.1±8.6) minutes, and the median (range) intraoperative blood loss during lymph node dissection 30 (10-80) mL. Of the 13 patients with negative postoperative horizontal margins, four were stage SM1 and had no lymph node metastases (Stage SM1) and nine were Stage SM2, of which had one positive regional lymph node and two received additional standard distal gastrectomy with D2 lymphadenectomy concurrently because of positive ESD specimens (lymph node negative). No lymph node metastases were found in the surgical specimens of these patients. The remaining two patients had positive vertical margins; both had undergone concurrent standard distal gastrectomy with D2 lymphadenectomy. One of them was found to be lymph node positive (No. 3, one node). Four patients had impaired gastric emptying after dual-scope treatment, all of whom recovered well with symptomatic management; one patient with a suspected lymphatic leak was also managed conservatively. There were no cases of postoperative bleeding, abdominal infection, or incisional infection. At a median follow-up of 14 (6-26) months, no tumor recurrence or metastasis had been identified in any of the patients. Three patients had a grade B nutrition score 3 to 6 months after surgery, all of whom had undergone major gastrectomy, and two patients who had undergone dual-scope surgery reported an increase in acid reflux and belching after surgery compared with the preoperative period. Conclusion: A combined technique is safe and feasible for the treatment of early gastric cancer and is worthy of further exploration.

[Molecular epidemiological characterization of hypervirulent carbapenem-resistant Klebsiella pneumoniae in a hospital in Henan Province from 2020 to 2022]

Objective: The study investigated the clinical distribution, antimicrobial resistance and epidemiologic characteristics of hypervirulent Carbapenem-resistant Klebsiella pneumoniae (hv-CRKP) in a hospital in Henan Province to provide a scientific basis for antibiotic use and nosocomial infection prevention and control. Methods: A retrospective analysis of the clinical data from the cases was carried out in this study. Clinical data of patients infected with the CRKP strain isolated from the clinical microbiology laboratory of Henan Provincial Hospital of Traditional Chinese Medicine from January 2020 to December 2022 were retrospectively analyzed. A string test, virulence gene screening, serum killing, and a G. mellonella infection model were used to screen hv-CRKP isolates. The clinical characteristics of hv-CRKP and the drug resistance rate of hv-CRKP to twenty-five antibiotics were analyzed using WHONET 5.6. Carbapenemase phenotypic characterization of the hv-CRKP was performed by colloidal gold immunochromatographic assay, and Carbapenemase genotyping, multi-locus sequence typing (MLST) and capsular serotyping of hv-CRKP isolates were performed by PCR and Sanger sequencing. Results: A total of non-duplicate 264 CRKP clinical isolates were detected in the hospital from 2020 to 2022, and 23 hv-CRKP isolates were detected, so the corresponding detection rate of hv-CRKP was 8.71% (23/264). The hv-CRKP isolates in this study were mainly from the intensive care unit (10/23) and neurosurgery department (8/23), and the main sources of hv-CRKP isolates were sputum (10/23) and bronchoalveolar lavage fluid (6/23). The hv-CRKP isolates in this study were highly resistant to β-lactam antibiotics, fluoroquinolones and aminoglycosides, and were only susceptible to colistin, tigecycline and ceftazidime/avibactam. The detection rate of the blaKPC-2 among 23 hv-CRKP isolates was 91.30% (21/23) and none of the class B and class D carbapenemases were detected. Results of MLST and capsular serotypes showed that ST11 type hv-CRKP was the dominant strain in the hospital, accounting for 56.52% (13/23), and K64 (9/13) and KL47 (4/13) were the major capsular serotypes. Conclusion: The hv-CRKP isolates from the hospital are mainly from lower respiratory tract specimens from patients admitted to the intensive care department and the drug resistance is relatively severe. The predominant strains with certain polymorphisms are mainly composed of the KPC-2-producing ST11-K64 and ST11-KL47 hv-CRKP isolates in the hospital.

[A medium- and long-term comparative observation on volumetric changes of cervical disc herniation after symmetrically or asymmetrically decompression and conservative treatment for cervical spondylotic myelopathy]

Objective: To compare the volumetric changes of cervical disc herniation (CDH) after cervical microendoscopic laminoplasty(CMEL),expansive open-door laminoplasty (EOLP) and conservative treatment. Methods: A retrospective study was conducted involving 101 patients with cervical spondylotic myelopathy(CSM),at the Department of Orthopaedic Surgery,the First Affiliated Hospital of Zhengzhou University from April 2012 to April 2021. The patients included 52 males and 49 females with an age of (54.7±11.8) years(range:25 to 86 years). Among them, 35 patients accepted CMEL treatment,33 patients accepted EOLP treatment,while 33 patients accepted conservative treatment. Volume data of CDH were measured by three-dimensional analysis of the initial and follow-up MRI images. The absorption rate and reprotrusion rate of CDH were calculated. The happening of resorption or reprotrusion was defined when the ratio was greater than 5%. The clinical outcomes and quality of life were evaluated by the Japanese Orthopaedic Association (JOA) score and the neck disability index (NDI).Quantitative data was analyzed by one-way ANOVA with post LSD-t test (multiple comparison) or Kruskal-Wallis test. Categorical data was analyzed by χ2 test. Results: The follow-up time of the CMEL group,EOLP group and the conservative treatment group were (27.6±18.8)months,(21.6±6.9)months and(24.9±16.3)months respectively with no significant difference(P>0.05). Changes of CDH volume in patients:(1) There were 96 CDH of 35 patients in the CMEL group,among which 78 showed absorption. The absorption frequency was 81.3%(78/96) and the absorption rate was ranged 5.9% to 90.9%;9 CDH showed reprotrusion,the reprotrusion frequency was 9.4% (9/96) and the reprotrusion rate was 5.9% to 13.3%;(2) There were 94 CDH of 33 patients in the EOLP group,of which 45 showed absorption. The absorption prevalence was 47.9% (45/94) and the absorption rate was 5.0% to 26.7%;20 CDH showed reprotruded,with the reprotrusion frequency of 21.3% (20/94) and the reprotrusion rate was 5.8% to 28.3%;(3) There were 102 CDH in 33 patients of the conservative group. Among them, 5 showed absorption. The absorption frequency was 4.9% (5/102),and the absorption rate was 7.2% to 14.3%;58 CDH showed reprotruded with the re-protrusion ratio of 56.9% (58/102) and the re-protrusion rate was 5.4% to 174.1%. The absorption ratio and reprotrusion ratio of the CMEL group were statistically different from EOLP group or the conservative group (P<0.01).The absorption ratio and reprotrusion ratio of the EOLP group was different from conservative group (all P<0.01). In terms of clinical outcomes, the excellent/good rate of the JOA score and NDI scores in the CMEL group were different from that of conservative group (all P<0.01) but not from that of the EOLP group(P>0.05). Conclusions: CMEL is an effective method for the treatment of CSM,making CDH easier to resorption compared to the EOLP or conservative treatment,thus making a better decompression effect on the nerves. This study enlightened on a new strategy for the clinical treatment of CSM.

Burning lignin: overlooked cues for post-fire seed germination

Information about smoke cues for seed germination is fundamental to understanding fire adaptation. Recently, lignin-derived syringaldehyde (SAL) was identified as a new smoke cue for seed germination, which challenges the assumption that cellulose-derived karrikins are the primary smoke cues. We highlight the overlooked association between lignin and the fire adaptation of plants.

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

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

Variation in gene expression along an elevation gradient of Rhododendron sanguineum var. haemaleum assessed in a comparative transcriptomic analysis

Selection along environmental gradients may play a vital role in driving adaptive evolution. Nevertheless, genomic variation and genetic adaptation along environmental clines remains largely unknown in plants in alpine ecosystems. To close this knowledge gap, we assayed transcriptomic profiles of late flower bud and early leaf bud of Rhododendron sanguineum var. haemaleum from four different elevational belts between 3,000 m and 3,800 m in the Gaoligong Mountains. By comparing differences in gene expression of these samples, a gene co-expression network (WGCNA) was constructed to identify candidate genes related to elevation. We found that the overall gene expression patterns are organ-specific for the flower and leaf. Differentially expressed unigenes were identified in these organs. In flowers, these were mainly related to terpenoid metabolism (RsHMGR, RsTPS), while in leaves mainly related to anthocyanin biosynthesis (RsCHS, RsF3'5'H). Terpenoids are the main components of flower scent (fragrance) likely attracting insects for pollination. In response to fewer pollinators at higher elevation zone, it seems relatively less scent is produced in flower organs to reduce energy consumption. Secondary metabolites in leaves such as anthocyanins determine the plants' alternative adaptive strategy to extreme environments, such as selective pressures of insect herbivory from environmental changes and substrate competition in biosynthesis pathways at high elevations. Our findings indicated that the gene expression profiles generated from flower and leaf organs showed parallel expression shifts but with different functionality, suggesting the existence of flexibility in response strategies of plants exposed to heterogeneous environments across elevational gradients. The genes identified here are likely to be involved in the adaptation of the plants to these varying mountainous environments. This study thus contributes to our understanding of the molecular mechanisms of adaptation in response to environmental change.

The identity of Dinochloa species and enumeration of Melocalamus (Poaceae: Bambusoideae) in China

Three woody bamboo species collected in Hainan, China in 1940 have been described as Dinochloa based on vegetative specimens. However, the identity of these species has long been in doubt, largely because the vegetative phase in species of Dinochloa is morphologically similar to that in species of Melocalamus, a climbing or scrambling bamboo genus of the paleotropical woody bamboos (Poaceae: Bambusoideae) that consists of about 15 species and one variety. To determine the phylogenetic affinity of the three Dinochloa species from Hainan, we sampled almost all recognized Chinese species of Melocalamus and representative species of Dinochloa as well as other closely related genera, performed molecular phylogenetic analysis, and compared their morphology based on herbarium and fieldwork investigation. Our ddRAD data indicate that the three species from Hainan are closely related to Melocalamus, not Dinochloa. Morphological analysis showed that these three species have a climbing habit but do not grow spirally, their culm leaves have smooth bases, and there is a ring of powder and/or tomenta above and below the nodes. Taken together our findings indicate that the three species from Hainan originally published in Dinochloa should be transferred to Melocalamus, i.e., Melocalamus orenudus (McClure) D.Z. Li & J.X. Liu, Melocalamus puberulus (McClure) D.Z. Li & J.X. Liu, and Melocalamus utilis (McClure) D.Z. Li & J.X. Liu, respectively. This study concludes with an enumeration of Chinese species of Melocalamus, with a key to nine recognized species and one variety, and a lectotypification for M. compatiflorus.

The establishment of plants following long-distance dispersal

Long-distance dispersal (LDD) beyond the range of a species is an important driver of ecological and evolutionary patterns, but insufficient attention has been given to postdispersal establishment. In this review, we summarize current knowledge of the post-LDD establishment phase in plant colonization, identify six key determinants of establishment success, develop a general quantitative framework for post-LDD establishment, and address the major challenges and opportunities in future research. These include improving detection and understanding of LDD using novel approaches, investigating mechanisms determining post-LDD establishment success using mechanistic modeling and inference, and comparison of establishment between past and present. By addressing current knowledge gaps, we aim to further our understanding of how LDD affects plant distributions, and the long-term consequences of LDD events.

Testing complete plastomes and nuclear ribosomal DNA sequences for species identification in a taxonomically difficult bamboo genus Fargesia

Fargesia, the largest genus within the temperate bamboo tribe Arundinarieae, has more than 90 species mainly distributed in the mountains of Southwest China. The Fargesia bamboos are important components of the subalpine forest ecosystems that provide food and habitat for many endangered animals, including the giant panda. However, species-level identification of Fargesia is difficult. Moreover, the rapid radiation and slow molecular evolutionary rate of Fargesia pose a significant challenge to using DNA barcoding with standard plant barcodes (rbcL, matK, and ITS) in bamboos. With progress in the sequencing technologies, complete plastid genomes (plastomes) and nuclear ribosomal DNA (nrDNA) sequences have been proposed as organelle barcodes for species identification; however, these have not been tested in bamboos. We collected 196 individuals representing 62 species of Fargesia to comprehensively evaluate the discriminatory power of plastomes and nrDNA sequences compared to standard barcodes. Our analysis indicates that complete plastomes have substantially higher discriminatory power (28.6%) than standard barcodes (5.7%), whereas nrDNA sequences show a moderate improvement (65.4%) compared to ITS (47.2%). We also found that nuclear markers performed better than plastid markers, and ITS alone had higher discriminatory power than complete plastomes. The study also demonstrated that plastomes and nrDNA sequences can contribute to intrageneric phylogenetic resolution in Fargesia. However, neither of these sequences were able to discriminate all the sampled species, and therefore, more nuclear markers need to be identified.

New insights into intergeneric relationships of Hickeliinae (Poaceae: Bambusoideae) revealed by complete plastid genomes

The Hickeliinae (Poaceae: Bambusoideae) is an ecologically and economically significant subtribe of tropical bamboos restricted to Madagascar, Comoros, Reunion Island, and a small part of continental Africa (Tanzania). Because these bamboos rarely flower, field identification is challenging, and inferring the evolutionary history of Hickeliinae from herbarium specimens is even more so. Molecular phylogenetic work is critical to understanding this group of bamboos. Here, comparative analysis of 22 newly sequenced plastid genomes showed that members of all genera of Hickeliinae share evolutionarily conserved plastome structures. We also determined that Hickeliinae plastome sequences are informative for phylogenetic reconstructions. Phylogenetic analysis showed that all genera of Hickeliinae are monophyletic, except for Nastus, which is paraphyletic and forms two distant clades. The type species of Nastus (Clade II) is endemic to Reunion Island and is not closely related to other sampled species of Nastus endemic to Madagascar (Clade VI). Clade VI (Malagasy Nastus) is sister to the Sokinochloa + Hitchcockella clade (Clade V), and both clades have a clumping habit with short-necked pachymorph rhizomes. The monotypic Decaryochloa is remarkable in having the longest floret in Bambuseae and forms a distinct Clade IV. Clade III, which has the highest generic diversity, consists of Cathariostachys, Perrierbambus, Sirochloa, and Valiha, which are also morphologically diverse. This work provides significant resources for further genetic and phylogenomic studies of Hickeliinae, an understudied subtribe of bamboo.

Phylogenomic insights into the origin and evolutionary history of evergreen broadleaved forests in East Asia under Cenozoic climate change

The evergreen versus deciduous leaf habit is an important functional trait for adaptation of forest trees and has been hypothesized to be related to the evolutionary processes of the component species under paleoclimatic change, and potentially reflected in the dynamic history of evergreen broadleaved forests (EBLFs) in East Asia. However, knowledge about the shift of evergreen versus deciduous leaf with the impact of paleoclimatic change using genomic data remains rare. Here, we focus on the Litsea complex (Lauraceae), a key lineage with dominant species of EBLFs, to gain insights into how evergreen versus deciduous trait shifted, providing insights into the origin and historical dynamics of EBLFs in East Asia under Cenozoic climate change. We reconstructed a robust phylogeny of the Litsea complex using genome-wide single-nucleotide variants (SNVs) with eight clades resolved. Fossil-calibrated analyses, diversification rate shifts, ancestral habit, ecological niche modelling and climate niche reconstruction were employed to estimate its origin and diversification pattern. Taking into account studies on other plant lineages dominating EBLFs of East Asia, it was revealed that the prototype of EBLFs in East Asia probably emerged in the Early Eocene (55-50 million years ago [Ma]), facilitated by the greenhouse warming. As a response to the cooling and drying climate in the Middle to Late Eocene (48-38 Ma), deciduous habits were evolved in the dominant lineages of the EBLFs in East Asia. Up to the Early Miocene (23 Ma), the prevailing of East Asian monsoon increased the extreme seasonal precipitation and accelerated the emergence of evergreen habits of the dominant lineages, and ultimately shaped the vegetation resembling that of today.

Phylotranscriptomic analyses reveal multiple whole-genome duplication events, the history of diversification and adaptations in the Araceae

BACKGROUND AND AIMS

The Araceae are one of the most diverse monocot families with numerous morphological and ecological novelties. Plastid and mitochondrial genes have been used to investigate the phylogeny and to interpret shifts in the pollination biology and biogeography of the Araceae. In contrast, the role of whole-genome duplication (WGD) in the evolution of eight subfamilies remains unclear.

METHODS

New transcriptomes or low-depth whole-genome sequences of 65 species were generated through Illumina sequencing. We reconstructed the phylogenetic relationships of Araceae using concatenated and species tree methods, and then estimated the age of major clades using TreePL. We inferred the WGD events by Ks and gene tree methods. We investigated the diversification patterns applying time-dependent and trait-dependent models. The expansions of gene families and functional enrichments were analysed using CAFE and InterProScan.

KEY RESULTS

Gymnostachydoideae was the earliest diverging lineage followed successively by Orontioideae, Lemnoideae and Lasioideae. In turn, they were followed by the clade of 'bisexual climbers' comprised of Pothoideae and Monsteroideae, which was resolved as the sister to the unisexual flowers clade of Zamioculcadoideae and Aroideae. A special WGD event ψ (psi) shared by the True-Araceae clade occurred in the Early Cretaceous. Net diversification rates first declined and then increased through time in the Araceae. The best diversification rate shift along the stem lineage of the True-Araceae clade was detected, and net diversification rates were enhanced following the ψ-WGD. Functional enrichment analyses revealed that some genes, such as those encoding heat shock proteins, glycosyl hydrolase and cytochrome P450, expanded within the True-Araceae clade.

CONCLUSIONS

Our results improve our understanding of aroid phylogeny using the large number of single-/low-copy nuclear genes. In contrast to the Proto-Araceae group and the lemnoid clade adaption to aquatic environments, our analyses of WGD, diversification and functional enrichment indicated that WGD may play a more important role in the evolution of adaptations to tropical, terrestrial environments in the True-Araceae clade. These insights provide us with new resources to interpret the evolution of the Araceae.

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.

Inflorescences of Fargesiaangustissima T.P. Yi and Yushaniapauciramificans T.P. Yi (Poaceae, Bambusoideae) shed light on the taxonomy of the Sino-Himalayan alpine bamboos

The taxonomy of the Sino-Himalayan alpine bamboos is controversial due to their complex evolutionary history and further complicated by the scarcity of inflorescence. Here, we supplement the description of the inflorescence of Fargesiaangustissima T.P. Yi and Yushaniapauciramificans T.P. Yi, which shed light on the taxonomy of Fargesia Franchet, Borinda Stapleton and Yushania Keng. F.angustissima has compressed inflorescence unilateral stretching out from reduced spathe, showing a transitional state between species with condensed inflorescence embraced by spathe-like bracts and species with open inflorescence without bracts. Considering that extensive gene flow existed between several clades of Fargesia found in recent studies, a broadly-defined Fargesia s. l. should be adopted. Meanwhile, the inflorescence of Y.pauciramificans has typical characteristics of Yushania, such as axilla with tuberculate glands, rachilla internodes ciliate and cylindrical florets, supporting the delimitation of Yushania.

Monsoon intensification in East Asia triggered the evolution of its flora

INTRODUCTION

East Asia (EA), which falls within the region of the Asian monsoon that is composed of the East Asia monsoon (EAM) and the Indian monsoon (IM), is known for its high species diversity and endemism. This has been attributed to extreme physiographical heterogeneity in conjunction with climate and sea-level changes during the Pleistocene, this hypothesis has been widely proven by phylogeographic studies. Recently, dated phylogenies have indicated that the origins (stem age) of the flora occurred after the Oligocene-Miocene boundary and are related to the establishment of the EAM.

METHODS

Hence, this study further examined whether the strengthening of the monsoons triggered floral evolution via a meta-analysis of the tempo-spatial pattern of evolutionary radiation dates (crown ages) of 101 endemic seed plant genera.

RESULTS

Taxonomic diversification began during the late Eocene, whereas the accumulated number of diversifications did not significantly accelerate until the late Miocene. The distribution of the weighted mean and the average divergence times in the EAM, IM, or transitional regions all fall within the mid-late Miocene. Fossils of the Tertiary relict genera are mostly and widely distributed outside EA and only half of the earliest fossils in the EA region are not older than Miocene, while their divergence times are mostly after the late Miocene. The pattern of divergence time of monotypic and polytypic taxa suggest the climatic changes after the late Pliocene exert more influence on monotypic taxa.

DISCUSSION

The two key stages of floral evolution coincide with the intensifications of the EAM and IM, especially the summer monsoon which brings a humid climate. An integrated review of previous studies concerning flora, genus, and species levels further supports our suggestion that monsoon intensification in EA triggered the evolution of its flora.

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.

Variations in genetic diversity in cultivated Pistacia chinensis

Identification of the evolution history and genetic diversity of a species is important in the utilization of novel genetic variation in this species, as well as for its conservation. Pistacia chinensis is an important biodiesel tree crop in China, due to the high oil content of its fruit. The aim of this study was to uncover the genetic structure of P. chinensis and to investigate the influence of intraspecific gene flow on the process of domestication and the diversification of varieties. We investigated the genetic structure of P. chinensis, as well as evolution and introgression in the subpopulations, through analysis of the plastid and nuclear genomes of 39 P. chinensis individuals from across China. High levels of variation were detected in the P. chinensis plastome, and 460 intraspecific polymorphic sites, 104 indels and three small inversions were identified. Phylogenetic analysis and population structure using the plastome dataset supported five clades of P. chinensis. Population structure analysis based on the nuclear SNPs showed two groups, clearly clustered together, and more than a third of the total individuals were classified as hybrids. Discordance between the plastid and nuclear genomes suggested that hybridization events may have occurred between highly divergent samples in the P. chinensis subclades. Most of the species in the P. chinensis subclade diverged between the late Miocene and the mid-Pliocene. The processes of domestication and cultivation have decreased the genetic diversity of P. chinensis. The extensive variability and structuring of the P. chinensis plastid together with the nuclear genomic variation detected in this study suggests that much unexploited genetic diversity is available for improvement in this recently domesticated species.

Plastome phylogenomic analysis reveals evolutionary divergences of Polypodiales suborder Dennstaedtiineae

BACKGROUND

Polypodiales suborder Dennstaedtiineae contain a single family Dennstaedtiaceae, eleven genera, and about 270 species, and include some groups that were previously placed in Dennstaedtiaceae, Hypolepidaceae, Monachosoraceae, and Pteridaceae. The classification and phylogenetic relationships among these eleven genera have been poorly understood. To explore the deep relationships within suborder Dennstaedtiineae and estimate the early diversification of this morphologically heterogeneous group, we analyzed complete plastomes of 57 samples representing all eleven genera of suborder Dennstaedtiineae using maximum likelihood and Bayesian inference.

RESULTS

The phylogenetic relationships of all the lineages in the bracken fern family Dennstaedtiaceae were well resolved with strong support values. All six genera of Hypolepidoideae were recovered as forming a monophyletic group with full support, and Pteridium was fully supported as sister to all the other genera in Hypolepidoideae. Dennstaedtioideae (Dennstaedtia s.l.) fell into four clades with full support: the Microlepia clade, the northern Dennstaedtia clade, the Dennstaedtia globulifera clade, and the Dennstaedtia s.s. clade. Monachosorum was strongly resolved as sister to all the remaining genera of suborder Dennstaedtiineae. Based on the well resolved relationships among genera, the divergence between Monachosorum and other groups of suborder Dennstaedtiineae was estimated to have occurred in the Early Cretaceous, and all extant genera (and clades) in Dennstaedtiineae, were inferred to have diversified since the Late Oligocene.

CONCLUSION

This study supports reinstating a previously published family Monachosoraceae as a segregate from Dennstaedtiaceae, based on unique morphological evidence, the shady habitat, and the deep evolutionary divergence from its closest relatives.

Genetic analysis of walnut cultivars from southwest China: Implications for germplasm improvement

Walnuts are highly valued for their rich nutritional profile and wide medicinal applications. This demand has led to the intensification of breeding activities in major walnut production areas such as southwest China, in order to develop more superior cultivars. With the increasing number of cultivars, accurate identification becomes fundamental to selecting the right cultivar for grafting, industrial processing or development of new cultivars. To ensure proper identification of cultivars and understand the genetic structure of wild and cultivated material, we genotyped 362 cultivated and wild individuals of walnut trees from southwest China (with two additional populations from Xinjiang, plus three cultivars from Canada, France and Belgium) using 36 polymorphic microsatellite loci. We found relatively low indices of genetic diversity (H _O = 0.570, H _E = 0.404, N _A = 2.345) as well as a high level of clonality (>85% of cultivars), indicating reliance on genetically narrow sources of parental material for breeding. Our STRUCTURE and PCoA analyses generally delineated the two species, though considerable levels of introgression were also evident. More significantly, we detected a distinct genetic group of cultivated Juglans sigillata, which mainly comprised individuals of the popular 'Yangbidapao' landrace. Finally, a core set of 18 SSR loci was selected, which was capable of identifying 32 cultivars. In a nutshell, our results call for more utilization of genetically disparate material, including wild walnut trees, as parental sources to breed for more cultivars. The data reported herein will significantly contribute towards the genetic improvement and conservation of the walnut germplasm in southwest China.

Fragaria mitogenomes evolve rapidly in structure but slowly in sequence and incur frequent multinucleotide mutations mediated by microinversions

Plant mitochondrial DNA has been described as evolving rapidly in structure but slowly in sequence. However, many of the noncoding portions of plant mitogenomes are not homologous among species, raising questions about the rate and spectrum of mutations in noncoding regions. Recent studies have suggested that the lack of homology in noncoding regions could be due to increased sequence divergence. We compared 30 kb of coding and 200 kb of noncoding DNA from 13 sequenced Fragaria mitogenomes, followed by analysis of the rate of sequence divergence, microinversion events and structural variations. Substitution rates in synonymous sites and nongenic sites are nearly identical, suggesting that the genome-wide point mutation rate is generally consistent. A surprisingly high number of large multinucleotide substitutions were detected in Fragaria mitogenomes, which may have resulted from microinversion events and could affect phylogenetic signal and local rate estimates. Fragaria mitogenomes preferentially accumulate deletions relative to insertions and substantial genomic arrangements, whereas mutation rates could positively associate with these sequence and structural changes among species. Together, these observations suggest that plant mitogenomes exhibit low point mutations genome-wide but exceptionally high structural variations, and our results favour a gain-and-loss model for the rapid loss of homology among plant mitogenomes.

[ABC prognostic classification and MELD 3.0 and COSSH-ACLF Ⅱ prognostic evaluation in acute-on-chronic liver failure]

Objective: To investigate the ABC prognostic classification and the updated version of Model for End-stage Liver Disease (MELD) score 3.0 and Chinese Group on the Study of Severe Hepatitis B ACLF Ⅱ score (COSSH-ACLF Ⅱ score) to evaluate the prognostic value in acute-on-chronic liver failure (ACLF). Methods: ABC classification was performed on a 1 409 follow-up cohorts. The area under the receiver operating characteristic curve (AUROC) was used to analyze MELD, MELD 3.0, COSSH-Ⅱ and COSSH-Ⅱ score after 3 days of hospitalization (COSSH-Ⅱ-3d). The prognostic predictive ability of patients were evaluated for 360 days, and the prediction differences of different classifications and different etiologies on the prognosis of ACLF were compared. Results: The survival curve of 1 409 cases with ACLF showed that the difference between class A, B, and C was statistically significant, Log Rank (Mantel-Cox) χ²=80.133, P<0.01. Compared with class A and C, χ²=76.198, P<0.01, the difference between class B and C, was not statistically significant χ²=3.717, P>0.05. AUROC [95% confidence interval (CI)] analyzed MELD, MELD 3.0, COSSH-Ⅱ and COSSH-Ⅱ-3d were 0.644, 0.655, 0.817 and 0.839, respectively (P<0.01). COSSH-Ⅱ had better prognostic predictive ability with class A ACLF and HBV-related ACLF (HBV-ACLF) for 360-days, and AUROC (95% CI) were 0.877 and 0.881, respectively (P<0.01), while MELD 3.0 prognostic predictive value was not better than MELD. Conclusion: ACLF prognosis is closely related to ABC classification. COSSH-Ⅱ score has a high predictive value for the prognostic evaluation of class A ACLF and HBV-ACLF. COSSH-Ⅱ score has a better prognostic evaluation value after 3 days of hospitalization, suggesting that attention should be paid to the treatment of ACLF in the early stage of admission.

Single-base methylome analysis reveals dynamic changes of genome-wide DNA methylation associated with rapid stem growth of woody bamboos

CG and CHG methylation levels in the rapid shoot growth stages (ST2-ST4) of woody bamboos were obviously decreased, which might regulate the internode elongation during rapid shoot growth, while CHH methylation was strongly associated with shoot developmental time or age. DNA methylation plays a critical role in the regulation of plant growth and development. Woody bamboos have a unique trait of rapid stem growth resulted from internode elongation at the shooting period. However, it is still unclear whether DNA methylation significantly controls the bamboo rapid stem growth. Here we present whole-genome DNA methylation profiles of the paleotropical woody bamboo Bonia amplexicaulis at five newly defined stages of shoot growth, named ST1-ST5. We found that CG and CHG methylation levels in the rapid shoot growth stages (ST2-ST4) were significantly lower than in the incubation (ST1) and plateau stages (ST5). The changes in methylation levels mainly occurred in flanking regions of genes and gene body regions, and 23647 differentially methylated regions (DMRs) were identified between ST1 and rapid shoot growth stages (ST2-ST4). Combined with transcriptome analysis, we found that DMR-related genes enriched in the auxin and jasmonic acid (JA) signal transduction, and other pathways closely related to plant growth. Intriguingly, CHH methylation was not involved in the rapid shoot growth, but strongly associated with shoot developmental time by gradually accumulating in transposable elements (TEs) regions. Overall, our results reveal the importance of DNA methylation in regulating the bamboo rapid shoot growth and suggest a role of DNA methylation associated with development time or age in woody bamboos.

Phylogenomics and evolutionary history of Oreocnide (Urticaceae) shed light on recent geological and climatic events in SE Asia

Climate change and geological events have long been known to shape biodiversity, implying that these can likewise be viewed from a biological perspective. To study whether plants can shed light on this, and how they responded to climate change there, we examined Oreocnide, a genus widely distributed in SE Asia. Based on broad geographic sampling with genomic data, we employed an integrative approach of phylogenomics, molecular dating, historical biogeography, and ecological analyses. We found that Oreocnide originated in mainland East Asia and began to diversify ∼6.06 Ma, probably in response to a distinct geographic and climatic transition in East Asia at around that time, implying that the last important geological change in mainland SE Asia might be 1 Ma older than previously suggested. Around four immigration events to the islands of Malesia followed, indicating that immigration from the mainland could be an underestimated factor in the assembly of biotic communities in the region. Two detected increases of diversification rate occurred 3.13 and 1.19 Ma, which strongly implicated climatic rather than geological changes as likely drivers of diversification, with candidates being the Pliocene intensification of the East Asian monsoons, and Pleistocene climate and sea level fluctuations. Distribution modelling indicated that Pleistocene sea level and climate fluctuations were inferred to enable inter-island dispersal followed by allopatric separation, underpinning radiation in the genus. Overall, our study, based on multiple lines of evidence, linked plant diversification to the most recent climatic and geological events in SE Asia. We highlight the importance of immigration in the assembly and diversification of the SE Asian flora, and underscore the utility of plant clades, as independent lines of evidence, for reconstructing recent climatic and geological events in the SE Asian region.

Multitrophic diversity and biotic associations influence subalpine forest ecosystem multifunctionality

Biodiversity across multiple trophic levels is required to maintain multiple ecosystem functions. Yet, it remains unclear how multitrophic diversity and species interactions regulate ecosystem multifunctionality. Here, combining data from nine different trophic groups (including trees, shrubs, herbs, leaf mites, small mammals, bacteria, pathogenic fungi, saprophytic fungi and symbiotic fungi) and 13 ecosystem functions related to supporting, provisioning and regulating services, we used a multitrophic perspective to evaluate the effects of elevation, diversity and network complexity on scale-dependent subalpine forest multifunctionality. Our results demonstrate that elevation and soil pH significantly modified species composition and richness across multitrophic groups and influenced multiple functions simultaneously. We provide evidence that species richness across multiple trophic groups had stronger effects on multifunctionality than species richness at any single trophic level. Moreover, biotic associations, indicating the complexity of trophic networks, were positively associated with multifunctionality. The relative effects of diversity on multifunctionality increased at the scale of the larger community compared to a scale accounting for neighbouring interactions. Our results highlight the paramount importance of scale- and context- dependent multitrophic diversity and interactions for a better understanding of mountain ecosystem multifunctionality in a changing world. This article is protected by copyright. All rights reserved.

Structural Variation of Plastomes Provides Key Insight Into the Deep Phylogeny of Ferns

Structural variation of plastid genomes (plastomes), particularly large inversions and gene losses, can provide key evidence for the deep phylogeny of plants. In this study, we investigated the structural variation of fern plastomes in a phylogenetic context. A total of 127 plastomes representing all 50 recognized families and 11 orders of ferns were sampled, making it the most comprehensive plastomic analysis of fern lineages to date. The samples included 42 novel plastomes of 15 families with a focus on Hymenophyllales and Gleicheniales. We reconstructed a well-supported phylogeny of all extant fern families, detected significant structural synapomorphies, including 9 large inversions, 7 invert repeat region (IR) boundary shifts, 10 protein-coding gene losses, 7 tRNA gene losses or anticodon changes, and 19 codon indels (insertions or deletions) across the deep phylogeny of ferns, particularly on the backbone nodes. The newly identified inversion V5, together with the newly inferred expansion of the IR boundary R5, can be identified as a synapomorphy of a clade composed of Dipteridaceae, Matoniaceae, Schizaeales, and the core leptosporangiates, while a unique inversion V4, together with an expansion of the IR boundary R4, was verified as a synapomorphy of Gleicheniaceae. This structural evidence is in support of our phylogenetic inference, thus providing key insight into the paraphyly of Gleicheniales. The inversions of V5 and V7 together filled the crucial gap regarding how the "reversed" gene orientation in the IR region characterized by most extant ferns (Schizaeales and the core leptosporangiates) evolved from the inferred ancestral type as retained in Equisetales and Osmundales. The tRNA genes trnR-ACG and trnM-CAU were assumed to be relicts of the early-divergent fern lineages but intact in most Polypodiales, particularly in eupolypods; and the loss of the tRNA genes trnR-CCG, trnV-UAC, and trnR-UCU in fern plastomes was much more prevalent than previously thought. We also identified several codon indels in protein-coding genes within the core leptosporangiates, which may be identified as synapomorphies of specific families or higher ranks. This study provides an empirical case of integrating structural and sequence information of plastomes to resolve deep phylogeny of plants.

Determinants of Genetic Structure in a Highly Heterogeneous Landscape in Southwest China

Intra-specific genetic diversity is a fundamental component of biodiversity, and is key to species adaptation and persistence. However, significant knowledge gaps still exist in our understanding of the patterns of genetic diversity and their key determinants. Most previous investigations mainly utilized single-species and/or a limited number of explanatory variables; so here we mapped the patterns of plastid genetic diversity within 15 plant species, and explored the key determinants shaping these patterns using a wide range of variables. Population-level cpDNA sequence data for 15 plant species from the Longitudinal Range Gorge Region (LRGR), southwest China, were retrieved from literature and used to estimate haplotype diversity (H _D) and population pairwise genetic differentiation (F _ST) indices. Genetic diversity and divergence landscape surfaces were then generated based on the H _D and F _ST, respectively, to clarify the patterns of genetic structure in the region. Subsequently, we analyzed the relationships between plastid genetic diversity and 16 explanatory variables (classified as anthropogenic, climatic, and topographic). We found that the highest genetic diversity occurred in the Yulong Mountain region, with a significant proportion (~74.81%) of the high diversity land area being located outside of protected areas. The highest genetic divergence was observed approximately along the 25°N latitudinal line, with notable peaks in the western and eastern edges of the LRGR. Genetic diversity (H _D) was weakly but significantly positively correlated with both Latitude (lat) and Annual Mean Wet Day Frequency (wet), yet significantly negatively correlated with all of Longitude (long), Annual Mean Cloud Cover Percent (cld), Annual Mean Anthropogenic Flux (ahf), and Human Footprint Index (hfp). A combination of climatic, topographic, and anthropogenic factors explained a significant proportion (78%) of genetic variation, with topographic factors (lat and long) being the best predictors. Our analysis identified areas of high genetic diversity (genetic diversity "hotspots") and divergence in the region, and these should be prioritized for conservation. This study contributes to a better understanding of the features that shape the distribution of plastid genetic diversity in the LRGR and thus would inform conservation management efforts in this species-rich, but vulnerable region.

A well-supported nuclear phylogeny of Poaceae and implications for the evolution of C4 photosynthesis

Poaceae (the grasses) includes rice, maize, wheat, and other crops, and is the most economically important angiosperm family. Poaceae is also one of the largest plant families, consisting of over 11 000 species with a global distribution that contributes to diverse ecosystems. Poaceae species are classified into 12 subfamilies, with generally strong phylogenetic support for their monophyly. However, many relationships within subfamilies, among tribes and/or subtribes, remain uncertain. To better resolve the Poaceae phylogeny, we generated 342 transcriptomic and seven genomic datasets; these were combined with other genomic and transcriptomic datasets to provide sequences for 357 Poaceae species in 231 genera, representing 45 tribes and all 12 subfamilies. Over 1200 low-copy nuclear genes were retrieved from these datasets, with several subsets obtained using additional criteria, and used for coalescent analyses to reconstruct a Poaceae phylogeny. Our results strongly support the monophyly of 11 subfamilies; however, the subfamily Puelioideae was separated into two non-sister clades, one for each of the two previously defined tribes, supporting a hypothesis that places each tribe in a separate subfamily. Molecular clock analyses estimated the crown age of Poaceae to be ∼101 million years old. Ancestral character reconstruction of C₃/C₄ photosynthesis supports the hypothesis of multiple independent origins of C₄ photosynthesis. These origins are further supported by phylogenetic analysis of the ppc gene family that encodes the phosphoenolpyruvate carboxylase, which suggests that members of three paralogous subclades (ppc-aL1a, ppc-aL1b, and ppc-B2) were recruited as functional C₄ppc genes. This study provides valuable resources and a robust phylogenetic framework for evolutionary analyses of the grass family.