Characterizing genetic intra-tumor heterogeneity across 2,658 human cancer genomes.

Intra-tumor heterogeneity (ITH) is a mechanism of therapeutic resistance and therefore an important clinical challenge. However, the extent, origin, and drivers of ITH across cancer types are poorly understood. To address this, we extensively characterize ITH across whole-genome sequences of 2,658 cancer samples spanning 38 cancer types. Nearly all informative samples (95.1%) contain evidence of distinct subclonal expansions with frequent branching relationships between subclones. We observe positive selection of subclonal driver mutations across most cancer types and identify cancer type-specific subclonal patterns of driver gene mutations, fusions, structural variants, and copy number alterations as well as dynamic changes in mutational processes between subclonal expansions. Our results underline the importance of ITH and its drivers in tumor evolution and provide a pan-cancer resource of comprehensively annotated subclonal events from whole-genome sequencing data.

mEnrich-seq: methylation-guided enrichment sequencing of bacterial taxa of interest from microbiome.

Metagenomics has enabled the comprehensive study of microbiomes. However, many applications would benefit from a method that sequences specific bacterial taxa of interest, but not most background taxa. We developed mEnrich-seq (in which 'm' stands for methylation and seq for sequencing) for enriching taxa of interest from metagenomic DNA before sequencing. The core idea is to exploit the self versus nonself differentiation by natural bacterial DNA methylation and rationally choose methylation-sensitive restriction enzymes, individually or in combination, to deplete host and background taxa while enriching targeted taxa. This idea is integrated with library preparation procedures and applied in several applications to enrich (up to 117-fold) pathogenic or beneficial bacteria from human urine and fecal samples, including species that are hard to culture or of low abundance. We assessed 4,601 bacterial strains with mapped methylomes so far and showed broad applicability of mEnrich-seq. mEnrich-seq provides microbiome researchers with a versatile and cost-effective approach for selective sequencing of diverse taxa of interest.

Stuxnet fine-tunes Notch dose during development using a functional Polycomb response element.

Evolutionarily conserved Notch signaling is highly sensitive to changes in Notch receptor dose caused by intrinsic and environmental fluctuations. It is well known that epigenetic regulation responds dynamically to genetic, cellular and environmental stresses. However, it is unclear whether the Notch receptor dose is directly regulated at the epigenetic level. Here, by studying the role of the upstream epigenetic regulator Stuxnet (Stx) in Drosophila developmental signaling, we find that Stx promotes Notch receptor mRNA expression by counteracting the activity of Polycomb repressive complex 1 (PRC1). In addition, we provide evidence that Notch is a direct PRC1 target by identifying and validating in vivo the only bona fide Polycomb response element (PRE) among the seven Polycomb group (PcG)-binding sites revealed by DamID-seq and ChIP-seq analysis. Importantly, in situ deletion of this PRE results in increased Notch expression and phenotypes resembling Notch hyperactivation in cell fate specification. These results not only underscore the importance of epigenetic regulation in fine-tuning the Notch activity dose, but also the need to assess the physiological significance of omics-based PcG binding in development.

Multiomics analysis reveals the molecular mechanisms underlying virulence in Rhizoctonia and jasmonic acid-mediated resistance in Tartary buckwheat (Fagopyrum tataricum).

Rhizoctonia solani is a devastating soil-borne pathogen that seriously threatens the cultivation of economically important crops. Multiple strains with a very broad host range have been identified, but only 1 (AG1-IA, which causes rice sheath blight disease) has been examined in detail. Here, we analyzed AG4-HGI 3 originally isolated from Tartary buckwheat (Fagopyrum tataricum), but with a host range comparable to AG1-IA. Genome comparison reveals abundant pathogenicity genes in this strain. We used multiomic approaches to improve the efficiency of screening for disease resistance genes. Transcriptomes of the plant-fungi interaction identified differentially expressed genes associated with virulence in Rhizoctonia and resistance in Tartary buckwheat. Integration with jasmonate-mediated transcriptome and metabolome changes revealed a negative regulator of jasmonate signaling, cytochrome P450 (FtCYP94C1), as increasing disease resistance probably via accumulation of resistance-related flavonoids. The integration of resistance data for 320 Tartary buckwheat accessions identified a gene homolog to aspartic proteinase (FtASP), with peak expression following R. solani inoculation. FtASP exhibits no proteinase activity but functions as an antibacterial peptide that slows fungal growth. This work reveals a potential mechanism behind pathogen virulence and host resistance, which should accelerate the molecular breeding of resistant varieties in economically essential crops.

The evolutionary history of 2,658 cancers.

Cancer develops through a process of somatic evolution1,2. Sequencing data from a single biopsy represent a snapshot of this process that can reveal the timing of specific genomic aberrations and the changing influence of mutational processes3. Here, by whole-genome sequencing analysis of 2,658 cancers as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA)4, we reconstruct the life history and evolution of mutational processes and driver mutation sequences of 38 types of cancer. Early oncogenesis is characterized by mutations in a constrained set of driver genes, and specific copy number gains, such as trisomy 7 in glioblastoma and isochromosome 17q in medulloblastoma. The mutational spectrum changes significantly throughout tumour evolution in 40% of samples. A nearly fourfold diversification of driver genes and increased genomic instability are features of later stages. Copy number alterations often occur in mitotic crises, and lead to simultaneous gains of chromosomal segments. Timing analyses suggest that driver mutations often precede diagnosis by many years, if not decades. Together, these results determine the evolutionary trajectories of cancer, and highlight opportunities for early cancer detection.

Conserved intronic secondary structures with concealed branch sites regulate alternative splicing of poison exons.

Alternative splicing (AS) generates multiple RNA isoforms and increases the complexities of transcriptomes and proteomes. However, it remains unclear how RNA structures contribute to AS regulation. Here, we systematically search transcriptomes for secondary structures with concealed branch sites (BSs) in the alternatively spliced introns and predict thousands of them from six organisms, of which many are evolutionarily conserved. Intriguingly, a highly conserved stem-loop structure with concealed BSs is found in animal SF3B3 genes and colocalizes with a downstream poison exon (PE). Destabilization of this structure allows increased usage of the BSs and results in enhanced PE inclusion in human and Drosophila cells, leading to decreased expression of SF3B3. This structure is experimentally validated using an in-cell SHAPE-MaP assay. Through RNA interference screens of 28 RNA-binding proteins, we find that this stem-loop structure is sensitive to U2 factors. Furthermore, we find that SF3B3 also facilitates DNA repair and protects genome stability by enhancing interaction between ERCC6/CSB and arrested RNA polymerase II. Importantly, both Drosophila and human cells with the secondary structure mutated by genome editing exhibit altered DNA repair in vivo. This study provides a novel and common mechanism for AS regulation of PEs and reveals a physiological function of SF3B3 in DNA repair.

Deletions of singular U1 snRNA gene significantly interfere with transcription and 3'-end mRNA formation.

Small nuclear RNAs (snRNAs) are structural and functional cores of the spliceosome. In metazoan genomes, each snRNA has multiple copies/variants, up to hundreds in mammals. However, the expressions and functions of each copy/variant in one organism have not been systematically studied. Focus on U1 snRNA genes, we investigated all five copies in Drosophila melanogaster using two series of constructed strains. Analyses of transgenic flies that each have a U1 promoter-driven gfp revealed that U1:21D is the major and ubiquitously expressed copy, and the other four copies have specificities in developmental stages and tissues. Mutant strains that each have a precisely deleted copy of U1-gene exhibited various extents of defects in fly morphology or mobility, especially deletion of U1:82Eb. Interestingly, splicing was changed at limited levels in the deletion strains, while large amounts of differentially-expressed genes and alternative polyadenylation events were identified, showing preferences in the down-regulation of genes with 1-2 introns and selection of proximal sites for 3'-end polyadenylation. In vitro assays suggested that Drosophila U1 variants pulled down fewer SmD2 proteins compared to the canonical U1. This study demonstrates that all five U1-genes in Drosophila have physiological functions in development and play regulatory roles in transcription and 3'-end formation.

Single-cell bisulfite-free 5mC and 5hmC sequencing with high sensitivity and scalability.

Existing single-cell bisulfite-based DNA methylation analysis is limited by low DNA recovery, and the measurement of 5hmC at single-base resolution remains challenging. Here, we present a bisulfite-free single-cell whole-genome 5mC and 5hmC profiling technique, named Cabernet, which can characterize 5mC and 5hmC at single-base resolution with high genomic coverage. Cabernet utilizes Tn5 transposome for DNA fragmentation, which enables the discrimination between different alleles for measuring hemi-methylation status. Using Cabernet, we revealed the 5mC, hemi-5mC and 5hmC dynamics during early mouse embryo development, uncovering genomic regions exclusively governed by active or passive demethylation. We show that hemi-methylation status can be used to distinguish between pre- and post-replication cells, enabling more efficient cell grouping when integrated with 5mC profiles. The property of Tn5 naturally enables Cabernet to achieve high-throughput single-cell methylome profiling, where we probed mouse cortical neurons and embryonic day 7.5 (E7.5) embryos, and constructed the library for thousands of single cells at high efficiency, demonstrating its potential for analyzing complex tissues at substantially low cost. Together, we present a way of high-throughput methylome and hydroxymethylome detection at single-cell resolution, enabling efficient analysis of the epigenetic status of biological systems with complicated nature such as neurons and cancer cells.

Fine-Regional Role of the Claustrum in Anxiety and Higher Sensitivity to Cocaine in Adolescent Cocaine-Exposed Male Mice during Adulthood.

Adolescent cocaine exposure (ACE) induces anxiety and higher sensitivity to substances abuse during adulthood. Here, we show that the claustrum is crucial for controlling these psychiatric problems in male mice. In anxiety-like behavioral tests, the CaMKII-positive neurons in the median portion of the claustrum (MClaustrum) were triggered, and local suppression of these neurons reduced the anxiety-like behavior in ACE mice during adulthood. In contrast, the CaMKII-positive neurons in the anterior portion of the claustrum (AClaustrum) were more activated in response to subthreshold dose of cocaine induced conditioned place preference (CPP), and local suppression of these neurons blocked the acquisition of cocaine CPP in ACE mice during adulthood. Our findings for the first time identified the fine-regional role of the claustrum in regulating the anxiety and susceptibility to cocaine in ACE mice during adulthood, extending our understanding of the claustrum in substance use disorder.

Activating Lobule VI PCTH+-Med Pathway in Cerebellum Blocks the Acquisition of Methamphetamine Conditioned Place Preference in Mice.

Cerebellum has been implicated in drug addiction; however, its underlying cellular populations and neuronal circuitry remain largely unknown. In the current study, we identified a neural pathway from tyrosine hydroxylase (TH)-positive Purkinje cells (PCTH+) in cerebellar lobule VI to calcium/calmodulin-dependent protein kinase II (CaMKII)-positive glutamatergic neurons in the medial cerebellar nucleus (MedCaMKII), forming the lobule VI PCTH+-MedCaMKII pathway in male mice. In naive male mice, inhibition of PCTH+ neurons activated Med neurons. During conditioned place preference (CPP) training, exposure to methamphetamine (METH) inhibited lobule VI PCTH+ neurons while excited MedCaMKII neurons in mice. Silencing MedCaMKII using a tetanus toxin light chain (tettox) suppressed the acquisition of METH CPP in mice but resulted in motor coordination deficits in naive mice. In contrast, activating lobule VI PCTH+ terminals within Med inhibited the activity of Med neurons and subsequently blocked the acquisition of METH CPP in mice without affecting motor coordination, locomotor activity, and sucrose reinforcements in naive mice. Our findings identified a novel lobule VI PCTH+-MedCaMKII pathway within the cerebellum and explored its role in mediating the acquisition of METH-preferred behaviors.

mPFC DUSP1 mediates adolescent cocaine exposure-induced higher sensitivity to drug in adulthood.

Adolescent cocaine abuse increases the risk for developing addiction in later life, but the underlying molecular mechanism remains poorly understood. Here, we establish adolescent cocaine-exposed (ACE) male mouse models. A subthreshold dose of cocaine (sdC) treatment, insufficient to produce conditioned place preference (CPP) in adolescent mice, induces CPP in ACE mice during adulthood, along with more activated CaMKII-positive neurons, higher dual specificity protein kinase phosphatase-1 (Dusp1) mRNA, lower DUSP1 activity, and lower DUSP1 expression in CaMKII-positive neurons in the medial prefrontal cortex (mPFC). Overexpressing DUSP1 in CaMKII-positive neurons suppresses neuron activity and blocks sdC-induced CPP in ACE mice during adulthood. On the contrary, depleting DUSP1 in CaMKII-positive neurons activates more neurons and further enhances sdC-induced behavior in ACE mice during adulthood. Also, ERK1/2 might be a downstream signal of DUSP1 in the process. Our findings reveal a role of mPFC DUSP1 in ACE-induced higher sensitivity to the drug in adult mice. DUSP1 might be a potential pharmacological target to predict or treat the susceptibility to addictive drugs caused by adolescent substance use.

Loss of renal tubular G9a benefits acute kidney injury by lowering focal lipid accumulation via CES1.

Surgery-induced renal ischemia and reperfusion (I/R) injury and nephrotoxic drugs like cisplatin can cause acute kidney injury (AKI), for which there is no effective therapy. Lipid accumulation is evident following AKI in renal tubules although the mechanisms and pathological effects are unclear. Here, we report that Ehmt2-encoded histone methyltransferase G9a is upregulated in patients and mouse kidneys after AKI. Renal tubular specific knockout of G9a (Ehmt2Ksp ) or pharmacological inhibition of G9a alleviates lipid accumulation associated with AKI. Mechanistically, G9a suppresses transcription of the lipolytic enzyme Ces1; moreover, G9a and farnesoid X receptor (FXR) competitively bind to the same promoter regions of Ces1. Ces1 is consistently observed to be downregulated in the kidney of AKI patients. Pharmacological inhibition of Ces1 increases lipid accumulation, exacerbates renal I/R-injury and eliminates the beneficial effects on AKI observed in Ehmt2Ksp mice. Furthermore, lipid-lowering atorvastatin and an FXR agonist alleviate AKI by activating Ces1 and reducing renal lipid accumulation. Together, our results reveal a G9a/FXR-Ces1 axis that affects the AKI outcome via regulating renal lipid accumulation.

Sex-lethal regulates back-splicing and generation of the sex-differentially expressed circular RNAs.

Conversely to canonical splicing, back-splicing connects the upstream 3' splice site (SS) with a downstream 5'SS and generates exonic circular RNAs (circRNAs) that are widely identified and have regulatory functions in eukaryotic gene expression. However, sex-specific back-splicing in Drosophila has not been investigated and its regulation remains unclear. Here, we performed multiple RNA analyses of a variety sex-specific Drosophila samples and identified over ten thousand circular RNAs, in which hundreds are sex-differentially and -specifically back-spliced. Intriguingly, we found that expression of SXL, an RNA-binding protein encoded by Sex-lethal (Sxl), the master Drosophila sex-determination gene that is only spliced into functional proteins in females, promoted back-splicing of many female-differential circRNAs in the male S2 cells, whereas expression of a SXL mutant (SXLRRM) did not promote those events. Using a monoclonal antibody, we further obtained the transcriptome-wide RNA-binding sites of SXL through PAR-CLIP. After splicing assay of mini-genes with mutations in the SXL-binding sites, we revealed that SXL-binding on flanking exons and introns of pre-mRNAs facilitates back-splicing, whereas SXL-binding on the circRNA exons inhibits back-splicing. This study provides strong evidence that SXL has a regulatory role in back-splicing to generate sex-specific and -differential circRNAs, as well as in the initiation of sex-determination cascade through canonical forward-splicing.

Inferior Interfacial Superconductivity in 1 UC FeSe/SrVO3/SrTiO3 with Screened Interfacial Electron-Phonon Coupling.

Single-unit cell (1 UC) FeSe interfaced with TiOx or FeOx exhibits significantly enhanced superconductivity compared to that of bulk FeSe, with interfacial electron-phonon coupling (EPC) playing a crucial role. However, the reduced dimensionality in 1 UC FeSe, which may drive superconducting fluctuations, complicates our understanding of the enhancement mechanisms. We construct a new superconducting interface, 1 UC FeSe/SrVO3/SrTiO3. Here, the itinerant electrons of highly metallic SrVO3 films can screen all high-energy Fuchs-Kliewer phonons, including those of SrTiO3, making it the first FeSe/oxide system with screened interfacial EPC while maintaining the 1 UC FeSe thickness. Despite comparable doping levels, the heavily electron-doped 1 UC FeSe/SrVO3 exhibits a pairing temperature (Tg ∼ 48 K) lower than those of FeSe/SrTiO3 and FeSe/LaFeO3. Our findings disentangle the contributions of interfacial EPC from dimensionality in terms of enhancing Tg in FeSe/oxide interfaces, underscoring the critical importance of interfacial EPC. This FeSe/VOx interface also provides a platform for studying interfacial superconductivity.

Specific Inhibition of Interpeduncular Nucleus GABAergic Neurons Alleviates Anxiety-Like Behaviors in Male Mice after Prolonged Abstinence from Methamphetamine.

Anxiety is one of the most common withdrawal symptoms of methamphetamine (METH) abuse, which further drives relapse to drugs. Interpeduncular nucleus (IPN) has been implicated in anxiety-like behaviors and addiction, yet its role in METH-abstinence-induced anxiety remains unknown. Here, we found that prolonged abstinence from METH enhanced anxiety-like behaviors in male mice, accompanied by more excited IPN GABAergic neurons, as indicated by the increased c-fos expression and the enhanced neuronal excitability by electrophysiological recording in the GABAergic neurons. Using the designer receptors exclusively activated by designer drugs method, specific inhibition of IPN GABAergic neurons rescued the aberrant neuronal excitation of IPN GABAergic neurons and efficiently reduced anxiety-like behaviors, whereas it did not induce depression-like behaviors in male mice after prolonged abstinence from METH. These findings reveal that IPN GABAergic neurons should be a promising brain target to alleviate late withdrawal symptoms from METH with few side effects.SIGNIFICANCE STATEMENT Prolonged abstinence from METH triggers IPN GABAergic neurons and ultimately increases anxiety in male mice. Suppressing IPN GABAergic neurons rescues METH abstinence-induced aberrant neuronal excitation of IPN GABAergic neurons and efficiently reduces anxiety in mice.

Mechanistic Insights into Surfactant-Modulated Electrode-Electrolyte Interface for Steering H2O2 Electrosynthesis.

Electrocatalytic reactions taking place at the electrified electrode-electrolyte interface involve processes of proton-coupled electron transfer. Interfacial protons are delivered to the electrode surface via a H2O-dominated hydrogen-bond network. Less efforts are made to regulate the interfacial proton transfer from the perspective of interfacial hydrogen-bond network. Here, we present quaternary ammonium salt cationic surfactants as electrolyte additives for enhancing the H2O2 selectivity of the oxygen reduction reaction (ORR). Through in situ vibrational spectroscopy and molecular dynamics calculation, it is revealed that the surfactants are irreversibly adsorbed on the electrode surface in response to a given bias potential range, leading to the weakening of the interfacial hydrogen-bond network. This decreases interfacial proton transfer kinetics, particularly at high bias potentials, thus suppressing the 4-electron ORR pathway and achieving a highly selective 2-electron pathway toward H2O2. These results highlight the opportunity for steering H2O-involved electrochemical reactions via modulating the interfacial hydrogen-bond network.

High-Precision Measurements of 44Ca/40Ca Using Multiple-Collector Inductively Coupled Plasma Mass Spectrometry without Collision-Cell Technology.

Herein, we present a new method for determining the Ca isotopic composition of geological samples. To eliminate matrix elements from Ca, a column chromatography method was developed using a N,N,N'N' tetraoctyl-1,5-diglycolamide (TODGA) resin. The Ca isotopic compositions were measured by a multiple collector inductively coupled plasma mass spectrometry (MC-ICP-MS) without collision cell equipment, especially that direct measurement to 44Ca/40Ca can be achieved. To mitigate the interference from 40Ar during 40Ca measurement, the cold plasma technique was used to suppress the Ar+ generation, resulting in a background Ar+ intensity of <300 mV, in contrast to the conventional hot plasma conditions, which typically yield thousands of volts for Ar+ intensities. Given the potential for a concentration mismatch between the sample and bracketed standard solutions to cause an intensive shift in measured Ca isotopic compositions, a correction for the [Ca] match is needed. To avoid matrix effects arising from residue matrix elements, it is crucial to limit the concentrations below 1% of Ca for most matrix elements (including Al, Mg, K, Na, and Sr) and below 1‰ for Fe. Notably, the tolerance of residue Sr is effectively improved compared to measurements with CC-MC-ICP-MS and traditional Hot-plasma-SSB-MC-ICP-MS methods with the conventional hot plasma technique, thereby lowering the complexity of column chemistry. The measured δ44/40Ca, δ44/42Ca, and ε40Ca values for eight reference materials agree well with previously reported values within analytical uncertainties. This method demonstrates long-term precision is better than 0.10‰ (two standard deviations) for both δ values (i.e., δ44/40Ca and δ44/42Ca). We anticipate that the proposed method will benefit the growth of the Ca isotope data set and foster an increase in the application of Ca isotope in Earth science studies.

GRAS gene family in rye (Secale cereale L.): genome-wide identification, phylogeny, evolutionary expansion and expression analyses.

BACKGROUND: The GRAS transcription factor family plays a crucial role in various biological processes in different plants, such as tissue development, fruit maturation, and environmental stress. However, the GRAS family in rye has not been systematically analyzed yet. RESULTS: In this study, 67 GRAS genes in S. cereale were identified and named based on the chromosomal location. The gene structures, conserved motifs, cis-acting elements, gene replications, and expression patterns were further analyzed. These 67 ScGRAS members are divided into 13 subfamilies. All members include the LHR I, VHIID, LHR II, PFYRE, and SAW domains, and some nonpolar hydrophobic amino acid residues may undergo cross-substitution in the VHIID region. Interested, tandem duplications may have a more important contribution, which distinguishes them from other monocotyledonous plants. To further investigate the evolutionary relationship of the GRAS family, we constructed six comparative genomic maps of homologous genes between rye and different representative monocotyledonous and dicotyledonous plants. The response characteristics of 19 ScGRAS members from different subfamilies to different tissues, grains at filling stages, and different abiotic stresses of rye were systematically analyzed. Paclobutrazol, a triazole-based plant growth regulator, controls plant tissue and grain development by inhibiting gibberellic acid (GA) biosynthesis through the regulation of DELLA proteins. Exogenous spraying of paclobutrazol significantly reduced the plant height but was beneficial for increasing the weight of 1000 grains of rye. Treatment with paclobutrazol, significantly reduced gibberellin levels in grain in the filling period, caused significant alteration in the expression of the DELLA subfamily gene members. Furthermore, our findings with respect to genes, ScGRAS46 and ScGRAS60, suggest that these two family members could be further used for functional characterization studies in basic research and in breeding programmes for crop improvement. CONCLUSIONS: We identified 67 ScGRAS genes in rye and further analysed the evolution and expression patterns of the encoded proteins. This study will be helpful for further analysing the functional characteristics of ScGRAS genes.

Multi-omics identification of a key glycosyl hydrolase gene FtGH1 involved in rutin hydrolysis in Tartary buckwheat (Fagopyrum tataricum).

Rutin, a flavonoid rich in buckwheat, is important for human health and plant resistance to external stresses. The hydrolysis of rutin to quercetin underlies the bitter taste of Tartary buckwheat. In order to identify rutin hydrolysis genes, a 200 genotypes mini-core Tartary buckwheat germplasm resource was re-sequenced with 30-fold coverage depth. By combining the content of the intermediate metabolites of rutin metabolism with genome resequencing data, metabolite genome-wide association analyses (GWAS) eventually identified a glycosyl hydrolase gene FtGH1, which could hydrolyse rutin to quercetin. This function was validated both in Tartary buckwheat overexpression hairy roots and in vitro enzyme activity assays. Mutation of the two key active sites, which were determined by molecular docking and experimentally verified via overexpression in hairy roots and transient expression in tobacco leaves, exhibited abnormal subcellular localization, suggesting functional changes. Sequence analysis revealed that mutation of the FtGH1 promoter in accessions of two haplotypes might be necessary for enzymatic activity. Co-expression analysis and GWAS revealed that FtbHLH165 not only repressed FtGH1 expression, but also increased seed length. This work reveals a potential mechanism behind rutin metabolism, which should provide both theoretical support in the study of flavonoid metabolism and in the molecular breeding of Tartary buckwheat.

A model of hybrid speciation process drawn from three new poplar species originating from distant hybridization between sections.

Although numerous studies have been conducted on hybrid speciation, our understanding of this process remains limited. Through an 18-year systematic investigation of all taxa of Populus on the Qinghai-Tibet Plateau, we discovered three new taxa with clear characteristics of sect. Leucoides. Further evidence was gathered from morphology, whole-genome bioinformatics, biogeography, and breeding to demonstrate synthetically that they all originated from distant hybridization between sect. Leucoides and sect. Tacamahaca. P. gonggaensis originated from the hybridization of P. lasiocarpa with P. cathayana, P. butuoensis from the hybridization of P. wilsonii with P. szechuanica, and P. dafengensis from the hybridization of P. lasiocarpa with P. szechuanica. Due to heterosis, the three hybrid taxa possess greater ecological adaptability than their ancestral species. We propose a hybrid speciation process model that incorporates orthogonal, reverse, and backcrossing events. This model can adequately explain some crucial evolutionary concerns, such as the nuclear-cytoplasmic conflict on phylogeny and the extinction of ancestral species within the distribution range of hybrid species.