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.

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.

Recent Advances in Vitamin D3 Intervention to Eradicate Helicobacter pylori Infection.

Helicobacter pylori (HP) infections affect approximately one-third of children worldwide. In China, the incidence of HP infection in children ranges from approximately 30% to 60%. In addition to damaging the gastrointestinal tract mucosa, HP infection in children can negatively affect their growth and development, hematology, respiratory and hepatobiliary system, skin, nutritional metabolism, and autoimmune system. However, the rate of HP eradication also fell considerably from the previous rate due to the presence of drug-resistant HP strains and the limited types of antibiotics that can be used in young patients. Vitamin D3 (VitD3) is a steroid hormone that can reduce inflammation in the stomach mucosa induced by HP and can alleviate and eradicate HP through a variety of pathways and mechanisms, including immune regulation and the stimulation of antimicrobial peptide (AMP) secretion and Ca2+ influx, to reestablish lysosomal acidification; thus, these results provide new strategies and ideas for the eradication of drug-resistant HP strains.

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.

25-hydroxyvitamin D levels in children of different ages and with varying degrees of Helicobacter pylori infection and immunological features.

Background: Helicobacter pylori (HP) is a major cause of upper digestive tract diseases. However, the relationship between HP infection and 25-hydroxyvitamin D [25(OH)D] levels in children has not been fully elucidated. This study investigated the levels of 25(OH)D in children of different ages and with varying degrees of HP infection and immunological features as well as the correlations between 25(OH)D levels in children infected with HP and their ages and degrees of infection. Materials and methods: Ninety-four children who underwent upper digestive endoscopy were divided into an HP-positive group without peptic ulcers (Group A), an HP-positive group with peptic ulcers (Group B) and an HP-negative control group (Group C). The serum levels of 25(OH)D and immunoglobulin and the percentages of lymphocyte subsets were determined. HP colonization, the degree of inflammation, and the degree of activity were further evaluated by HE staining and immunohistochemical staining in gastric mucosal biopsy. Results: The 25(OH)D level of the HP-positive groups (50.93 ± 16.51 nmol/L) was significantly lower than that of the HP-negative group (62.89 ± 19.18 nmol/L). The 25(OH)D level of Group B (47.79 ± 14.79 nmol/L) was lower than that of Group A (51.53 ± 17.05 nmol/L) and was significantly lower than that of Group C (62.89 ± 19.18 nmol/L). The 25(OH)D level decreased with increasing age, and there was a significant difference between Group C subjects who were ≤5 years old and those who were aged 6-9 years and ≥10 years. The 25(OH)D level was negatively correlated with HP colonization (r = -0.411, P < 0.01) and the degree of inflammation (r = -0.456, P < 0.01). The percentages of lymphocyte subsets and immunoglobulin levels among Groups A, B and C were not significantly different. Conclusions: The 25(OH)D level was negatively correlated with HP colonization and the degree of inflammation. As the age of the children increased, the level of 25(OH)D decreased, and the susceptibility to HP infection increased.

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.

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.

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

MAIN CONCLUSION: 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.

Notch-mediated lactate metabolism regulates MDSC development through the Hes1/MCT2/c-Jun axis.

Myeloid-derived suppressor cells (MDSCs) and tumor-associated macrophages (TAMs) play critical roles in tumorigenesis. However, the mechanisms underlying MDSC and TAM development and function remain unclear. In this study, we find that myeloid-specific activation of Notch/RBP-J signaling downregulates lactate transporter MCT2 transcription via its downstream molecule Hes1, leading to reduced intracellular lactate levels, blunted granulocytic MDSC (G-MDSC) differentiation, and enhanced TAM maturation. We identify c-Jun as a novel intracellular sensor of lactate in myeloid cells using liquid-chromatography-mass spectrometry (LC-MS) followed by CRISPR-Cas9-mediated gene disruption. Meanwhile, lactate interacts with c-Jun to protect from FBW7 ubiquitin-ligase-mediated degradation. Activation of Notch signaling and blockade of lactate import repress tumor progression by remodeling myeloid development. Consistently, the relationship between the Notch-MCT2/lactate-c-Jun axis in myeloid cells and tumorigenesis is also confirmed in clinical lung cancer biopsies. Taken together, our current study shows that lactate metabolism regulated by activated Notch signaling might participate in MDSC differentiation and TAM maturation.

Overexpression of SLC25A51 promotes hepatocellular carcinoma progression by driving aerobic glycolysis through activation of SIRT5.

Solute carrier family 25 member 20 (SLC25A51) is a newly identified mammalian mitochondrial NAD+ transporter. However, the clinicopathological and biological significance of SLC25A51 in human cancers, including hepatocellular carcinoma (HCC), remains unclear. The aim of this study was to define the role of SLC25A51 in HCC progression. Here we demonstrate that SLC25A51 is significantly overexpressed in human HCC specimens and cell lines, caused by, at least in partial, the decrease of miR-212-3p. SLC25A51 overexpression is positively correlated with the clinicopathological characteristics of vascular invasion and tumor diameter, as well as poor survival in patients with HCC. Knockdown of SLC25A51 attenuated, while overexpression of SLC25A51 enhanced the growth and metastasis of HCC cells both in vitro and in vivo. Mechanistically, glucose metabolism reprogramming from oxidative phosphorylation to glycolysis by activation of mitochondrial sirtuin 5 (SIRT5) was found to contribute to the promotion of growth and metastasis by SLC25A51 in HCC cells. Together, these findings reveal important roles of SLC25A51 in HCC tumorigenesis and suggest SLC25A51 as a promising prognostic marker and therapeutic target for treating HCC.

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.

New Genes Interacted With Recent Whole-Genome Duplicates in the Fast Stem Growth of Bamboos.

As drivers of evolutionary innovations, new genes allow organisms to explore new niches. However, clear examples of this process remain scarce. Bamboos, the unique grass lineage diversifying into the forest, have evolved with a key innovation of fast growth of woody stem, reaching up to 1 m/day. Here, we identify 1,622 bamboo-specific orphan genes that appeared in recent 46 million years, and 19 of them evolved from noncoding ancestral sequences with entire de novo origination process reconstructed. The new genes evolved gradually in exon-intron structure, protein length, expression specificity, and evolutionary constraint. These new genes, whether or not from de novo origination, are dominantly expressed in the rapidly developing shoots, and make transcriptomes of shoots the youngest among various bamboo tissues, rather than reproductive tissue in other plants. Additionally, the particularity of bamboo shoots has also been shaped by recent whole-genome duplicates (WGDs), which evolved divergent expression patterns from ancestral states. New genes and WGDs have been evolutionarily recruited into coexpression networks to underline fast-growing trait of bamboo shoot. Our study highlights the importance of interactions between new genes and genome duplicates in generating morphological innovation.

Comparative plastomic analysis and insights into the phylogeny of Salvia (Lamiaceae).

Salvia is the largest genus of Lamiaceae, with almost 1000 species, and has been divided into 11 subgenera. Salvia subg. Glutinaria, native to East Asia, is particularly important because of its potential medicinal value. However, the interspecific relationships of this subgenus have not been resolved and the plastomes of Salvia have rarely been studied. In the current study, we compared plastid genome structure and organization of 19 species of Salvia (14 newly sequenced and 5 previously published). Our comparative analysis showed that all Salvia plastomes examined have a quadripartite structure typical of most angiosperms and contain an identical set of 114 unique genes (80 protein-coding genes, 4 rRNA genes, and 30 tRNA genes). The plastome structure of all Salvia species is highly conserved like other Lamiaceae plastomes. Gene content, gene order, and GC content were highly similar in these plastomes. The inverted repeats/single copy region (IR/SC) boundaries of Salvia are highly conserved, and IR contraction only occurred in two species (Salvia mekongensis and S. rosmarinus). In Salvia, sequence divergence was higher in non-coding regions than in coding regions. We found that using large single copy (LSC) and small single copy regions (SSC) with exclusion of the rapidly evolving sites produced the highest resolution in phylogenetic analysis of Salvia, suggesting that using suitable informative sites to build trees is more conducive in phylogenetic research. This study assembled a powerful matrix data set for studying the phylogeny of Salvia, resolving the interspecific relationship of Salvia subg. Glutinaria. The newly sequenced plastid genomes will also enrich the plastome database of Salvia, providing the scientific basis for the development and utilization of germplasm resources of this large and important genus.

[Comparative analysis of binding multi-fold rib graft, iliac bone graft and titanium mesh graft during surgery of tuberculosis of thoracic vertebra].

OBJECTIVE: To compare the clinical effects of three different methods of binding multi-fold rib graft, iliac bone graft and titanium mesh graft in tuberculosis of thoracic vertebra by approach of transverse rib process. METHODS: A hundred and seven patients with tuberculosis of thoracic vertebra received surgical treatment from January 2010 to December 2016 were retrospectively analyzed. The patients were divided into three groups according different methods of bone graft. The surgical approach of the transverse rib process was used in all 107 patients, after thoroughly remove the necrotic tissue of tuberculosis, three different bone grafts were used respectively including iliac bone graft (36 cases, group A), binding multi-fold rib graft (35 cases, group B), titanium mesh bone graft (36 cases, group C). Perioperative indexes, the time required for bone graft during operation, intraoperation blood loss, the loss rate of the anterior edge of the lesion, Cobb angle, postoperative bone graft fusion time, spinal nerve recovery and Oswestry Disability Index were compared among three groups. RESULTS: All the patients were followed up for 13 to 24 months, and the operation time required for bone graft was (23.2±4.1) min in group A, (23.8± 4.4)min in group B, and (25.5±4.2) min in group C, with no statistically significant difference among three groups (P>0.05). Intraoperative blood loss was (541.6±35.3) ml in group A, (546.8±27.8) ml in group B, and (540.1±34.5) ml in group C, withno statistically significant difference among three groups(P>0.05). Preoperative anterior vertebral height loss rate was (46.0± 3.1)% in group A, (46.4±3.3)% in group B, and (45.3±3.6)% in group B;at the final follow up, the loss rate of anterior vertebral height among three groups was (8.6±5.0)%, (8.1±4.2)%, (9.4±4.3)%, respectively. There were no statistically significant differences before operation and final follow-up among three groups (P>0.05). Preoperative Cobb angle was (35.1±4.8)° in group A, (35.2±4.5)° in group B and (35.2±4.5)° in group C, with no statistically significant difference among three groups (P>0.05);postoperative at 3 days, Cobb angle in three groups was (15.1±3.6)°, (15.3±3.1)° and (15.2±3.4)°, respectively, there was no statistically significant difference among three groups (P>0.05);at the final follow-up, the Cobb angle among three groups was (17.7±3.3)°, (17.9±3.9)°, (18.6±3.6)°, respectively, with no statistically significant difference among three groups (P>0.05). The time of bone graft fusion was (5.6±0.5) months in group A, (5.6±0.6) months in group B and (5.8±0.6)months in group C, with no statistically significant difference among three groups (P>0.05). Frankel classification at the final follow up, 4 cases were grade B, 7 cases were grade C, 10 cases were grade D, and 86 cases were grade E. Spinal nerve function in all three groups recovered to a certain extent after treatment, with no statistically significant difference among three groups (P> 0.05). Oswestry Disability Index at the final follow-up showed no statistically significant difference among three groups (P> 0.05). CONCLUSION: The approach of transverse rib process for debridement of lesions can effectively treat tuberculosis of thoracic vertebra by binding multi-fold rib graft, iliac bone graft and titanium mesh graft, but binding multi-fold rib graft can effectively avoid iliac bone donor complications, and is an effective alternative to iliac bone graft, which is worth popularizing.

The Pharus latifolius genome bridges the gap of early grass evolution.

The grass family (Poaceae) includes all commercial cereal crops and is a major contributor to biomass in various terrestrial ecosystems. The ancestry of all grass genomes includes a shared whole-genome duplication (WGD), named rho (ρ) WGD, but the evolutionary significance of ρ-WGD remains elusive. We sequenced the genome of Pharus latifolius, a grass species (producing a true spikelet) in the subfamily Pharoideae, a sister lineage to the core Poaceae including the (Panicoideae, Arundinoideae, Chloridoideae, Micrairoideae, Aristidoideae, and Danthonioideae (PACMAD) and Bambusoideae, Oryzoideae, and Pooideae (BOP) clades. Our results indicate that the P. latifolius genome has evolved slowly relative to cereal grass genomes, as reflected by moderate rates of molecular evolution, limited chromosome rearrangements and a low rate of gene loss for duplicated genes. We show that the ρ-WGD event occurred approximately 98.2 million years ago (Ma) in a common ancestor of the Pharoideae and the PACMAD and BOP grasses. This was followed by contrasting patterns of diploidization in the Pharus and core Poaceae lineages. The presence of two FRIZZY PANICLE-like genes in P. latifolius, and duplicated MADS-box genes, support the hypothesis that the ρ-WGD may have played a role in the origin and functional diversification of the spikelet, an adaptation in grasses related directly to cereal yields. The P. latifolius genome sheds light on the origin and early evolution of grasses underpinning the biology and breeding of cereals.

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.

Parallel ddRAD and Genome Skimming Analyses Reveal a Radiative and Reticulate Evolutionary History of the Temperate Bamboos.

Rapid evolutionary radiations are among the most challenging phylogenetic problems, wherein different types of data (e.g., morphology and molecular) or genetic markers (e.g., nuclear and organelle) often yield inconsistent results. The tribe Arundinarieae, that is, the temperate bamboos, is a clade of tetraploid originated 22 Ma and subsequently radiated in East Asia. Previous studies of Arundinarieae have found conflicting relationships and/or low support. Here, we obtain nuclear markers from ddRAD data for 213 Arundinarieae taxa and parallel sampling of chloroplast genomes from genome skimming for 147 taxa. We first assess the feasibility of using ddRAD-seq data for phylogenetic estimates of paleopolyploid and rapidly radiated lineages, optimize clustering thresholds, and analysis workflow for orthology identification. Reference-based ddRAD data assembly approaches perform well and yield strongly supported relationships that are generally concordant with morphology-based taxonomy. We recover five major lineages, two of which are notable (the pachymorph and leptomorph lineages), in that they correspond with distinct rhizome morphologies. By contrast, the phylogeny from chloroplast genomes differed significantly. Based on multiple lines of evidence, the ddRAD tree is favored as the best species tree estimation for temperate bamboos. Using a time-calibrated ddRAD tree, we find that Arundinarieae diversified rapidly around the mid-Miocene corresponding with intensification of the East Asian monsoon and the evolution of key innovations including the leptomorph rhizomes. Our results provide a highly resolved phylogeny of Arundinarieae, shed new light on the radiation and reticulate evolutionary history of this tribe, and provide an empirical example for the study of recalcitrant plant radiations. [Arundinarieae; ddRAD; paleopolyploid; genome skimming; rapid diversification; incongruence.].

ddRAD analyses reveal a credible phylogenetic relationship of the four main genera of Bambusa-Dendrocalamus-Gigantochloa complex (Poaceae: Bambusoideae).

The Bambusa-Dendrocalamus-Gigantochloa complex (BDG complex) is the most diversified and phylogenetically recalcitrant group of the paleotropical woody bamboos. Species of this complex occur in tropical and subtropical Asia and most of them are of great economic, cultural and ecological value. The lack of resolution achieved through the analyses of previous molecular datasets has long confounded its phylogenetic estimation and generic delimitation. Here, we adopted a ddRAD-seq strategy to investigate phylogenetic relationships of the four main genera (Bambusa, Dendrocalamus, Gigantochloa, and Melocalamus) in the BDG complex. A total of 102 species were sampled, and SNP data were generated. Both MP and ML analyses of the ddRAD-seq data resulted in a well-resolved topology with Gigantochloa and Melocalamus confirmed as monophyletic, and Melocalamus resolved as sister to the rest of the complex. Bambusa and Dendrocalamus were both resolved as paraphyletic. The phylogenetic relationships were mostly supported by morphological evidence including characters of the branch complement, rachilla, lodicules, filaments and stigma. We also generated and assembled complete plastid genomes of 48 representative species. There were conflicts between the plastome and the ddRAD topologies. Our study demonstrated that RAD-seq can be used to reconstruct evolutionary history of lineages such as the bamboos where ancient hybridization and polyploidy play a significant role. The four genera of the BDG complex have a complex evolutionary history which is likely a product of ancient introgression events.

Using nuclear loci and allelic variation to disentangle the phylogeny of Phyllostachys (Poaceae, Bambusoideae).

With the development of sequencing technologies, the use of multiple nuclear genes has become conventional for resolving difficult phylogenies. However, this technique also presents challenges due to gene-tree discordance, as a result of incomplete lineage sorting (ILS) and reticulate evolution. Although alleles can show sequence variation within individuals, which contain information regarding the evolution of organisms, they continue to be ignored in almost all phylogenetic analyses using randomly phased genome sequences. Here, we tried to incorporate alleles from multiple nuclear loci to study the phylogeny of the economically important bamboo genus Phyllostachys (Poaceae, Bambusoideae). Obtaining a total of 3926 sequences, we documented extensive allelic variation for 61 genes from 39 sampled species. Using datasets consisting of selected alleles, we demonstrated substantial discordance among phylogenetic relationships inferred from different alleles, as well as between concatenation and coalescent methods. Furthermore, ILS and hybridization were suggested to be underlying causes of the discordant phylogenetic signals. Taking these possible causes for conflicting phylogenetic results into consideration, we recovered the monophyly of Phyllostachys and its two morphology-defined sections. Our study also suggests that alleles deserve more attention in phylogenetic studies, since ignoring them can yield highly supported but spurious phylogenies. Meanwhile, alleles are helpful for unraveling complex evolutionary processes, particularly hybridization.