Subgenome-aware analyses reveal the genomic consequences of ancient allopolyploid hybridizations throughout the cotton family.

Malvaceae comprise some 4,225 species in 243 genera and nine subfamilies and include economically important species, such as cacao, cotton, durian, and jute, with cotton an important model system for studying the domestication of polyploids. Here, we use chromosome-level genome assemblies from representatives of five or six subfamilies (depending on the placement of Ochroma) to differentiate coexisting subgenomes and their evolution during the family's deep history. The results reveal that the allohexaploid Helicteroideae partially derive from an allotetraploid Sterculioideae and also form a component of the allodecaploid Bombacoideae and Malvoideae. The ancestral Malvaceae karyotype consists of 11 protochromosomes. Four subfamilies share a unique reciprocal chromosome translocation, and two other subfamilies share a chromosome fusion. DNA alignments of single-copy nuclear genes do not yield the same relationships as inferred from chromosome structural traits, probably because of genes originating from different ancestral subgenomes. These results illustrate how chromosome-structural data can unravel the evolutionary history of groups with ancient hybrid genomes.

Genomic structural variation contributes to evolved changes in gene expression in high-altitude Tibetan sheep.

Tibetan sheep were introduced to the Qinghai Tibet plateau roughly 3,000 B.P., making this species a good model for investigating genetic mechanisms of high-altitude adaptation over a relatively short timescale. Here, we characterize genomic structural variants (SVs) that distinguish Tibetan sheep from closely related, low-altitude Hu sheep, and we examine associated changes in tissue-specific gene expression. We document differentiation between the two sheep breeds in frequencies of SVs associated with genes involved in cardiac function and circulation. In Tibetan sheep, we identified high-frequency SVs in a total of 462 genes, including EPAS1, PAPSS2, and PTPRD. Single-cell RNA-Seq data and luciferase reporter assays revealed that the SVs had cis-acting effects on the expression levels of these three genes in specific tissues and cell types. In Tibetan sheep, we identified a high-frequency chromosomal inversion that exhibited modified chromatin architectures relative to the noninverted allele that predominates in Hu sheep. The inversion harbors several genes with altered expression patterns related to heart protection, brown adipocyte proliferation, angiogenesis, and DNA repair. These findings indicate that SVs represent an important source of genetic variation in gene expression and may have contributed to high-altitude adaptation in Tibetan sheep.

Scattered differentiation of unlinked loci across the genome underlines ecological divergence of the selfing grass Brachypodium stacei.

Ecological divergence without geographic isolation, as an early speciation process that may lead finally to reproductive isolation through natural selection, remains a captivating topic in evolutionary biology. However, the pattern of genetic divergence underlying this process across the genome may vary between species and mating systems. Here, we present evidence that Brachypodium stacei, an annual and highly selfing grass model species, has undergone sympatric ecological divergence without geographic isolation. Genomic, transcriptomic, and metabolomic analyses together with lab experiments mimicking the two opposite environmental conditions suggest that diploid B. stacei populations have diverged sympatrically in two slopes characterized by distinct biomes at Evolution Canyon I (ECI), Mount Carmel, Israel. Despite ongoing gene flow, primarily facilitated by seed dispersal, the level of gene flow has progressively decreased over time. This local adaptation involves the scattered divergence of many unlinked loci across the total genome that include both coding genes and noncoding regions. Additionally, we have identified significant differential expressions of genes related to the ABA signaling pathway and contrasting metabolome composition between the arid- vs. forest-adapted B. stacei populations in ECI. These results suggest that multiple small loci involved in environmental responses act additively to account for ecological adaptations by this selfing species in contrasting environments.

Uncovering the enigmatic evolution of bears in greater depth: The hybrid origin of the Asiatic black bear.

Bears are fascinating mammals because of their complex pattern of speciation and rapid evolution of distinct phenotypes. Interspecific hybridization has been common and has shaped the complex evolutionary history of bears. In this study, based on the largest population-level genomic dataset to date involving all Ursinae species and recently developed methods for detecting hybrid speciation, we provide explicit evidence for the hybrid origin of Asiatic black bears, which arose through historical hybridization between the ancestor of polar bear/brown bear/American black bears and the ancestor of sun bear/sloth bears. This was inferred to have occurred soon after the divergence of the two parental lineages in Eurasia due to climate-driven population expansion and dispersal. In addition, we found that the intermediate body size of this hybrid species arose from its combination of relevant genes derived from two parental lineages of contrasting sizes. This and alternate fixation of numerous other loci that had diverged between parental lineages may have initiated the reproductive isolation of the Asiatic black bear from its two parents. Our study sheds further light on the evolutionary history of bears and documents the importance of hybridization in new species formation and phenotypic evolution in mammals.

The m6A modification mediated-lncRNA POU6F2-AS1 reprograms fatty acid metabolism and facilitates the growth of colorectal cancer via upregulation of FASN.

BACKGROUND: Long noncoding RNAs (lncRNAs) have emerged as key players in tumorigenesis and tumour progression. However, the biological functions and potential mechanisms of lncRNAs in colorectal cancer (CRC) are unclear. METHODS: The novel lncRNA POU6F2-AS1 was identified through bioinformatics analysis, and its expression in CRC patients was verified via qRT-PCR and FISH. In vitro and in vivo experiments, such as BODIPY staining, Oil Red O staining, triglyceride (TAG) assays, and liquid chromatography mass spectrometry (LC-MS) were subsequently performed with CRC specimens and cells to determine the clinical significance, and functional roles of POU6F2-AS1. Biotinylated RNA pull-down, RIP, Me-RIP, ChIP, and patient-derived organoid (PDO) culture assays were performed to confirm the underlying mechanism of POU6F2-AS1. RESULTS: The lncRNA POU6F2-AS1 is markedly upregulated in CRC and associated with adverse clinicopathological features and poor overall survival in CRC patients. Functionally, POU6F2-AS1 promotes the growth and lipogenesis of CRC cells both in vitro and in vivo. Mechanistically, METTL3-induced m6A modification is involved in the upregulation of POU6F2-AS1. Furthermore, upregulated POU6F2-AS1 could tether YBX1 to the FASN promoter to induce transcriptional activation, thus facilitating the growth and lipogenesis of CRC cells. CONCLUSIONS: Our data revealed that the upregulation of POU6F2-AS1 plays a critical role in CRC fatty acid metabolism and might provide a novel promising biomarker and therapeutic target for CRC.

Subgenomic Stability of Progenitor Genomes During Repeated Allotetraploid Origins of the Same Grass Brachypodium hybridum.

Both homeologous exchanges and homeologous expression bias are generally found in most allopolyploid species. Whether homeologous exchanges and homeologous expression bias differ between repeated allopolyploid speciation events from the same progenitor species remains unknown. Here, we detected a third independent and recent allotetraploid origin for the model grass Brachypodium hybridum. Our homeologous exchange with replacement analyses indicated the absence of significant homeologous exchanges in any of the three types of wild allotetraploids, supporting the integrity of their progenitor subgenomes and the immediate creation of the amphidiploids. Further homeologous expression bias tests did not uncover significant subgenomic dominance in different tissues and conditions of the allotetraploids. This suggests a balanced expression of homeologs under similar or dissimilar ecological conditions in their natural habitats. We observed that the density of transposons around genes was not associated with the initial establishment of subgenome dominance; rather, this feature is inherited from the progenitor genome. We found that drought response genes were highly induced in the two subgenomes, likely contributing to the local adaptation of this species to arid habitats in the third allotetraploid event. These findings provide evidence for the consistency of subgenomic stability of parental genomes across multiple allopolyploidization events that led to the same species at different periods. Our study emphasizes the importance of selecting closely related progenitor species genomes to accurately assess homeologous exchange with replacement in allopolyploids, thereby avoiding the detection of false homeologous exchanges when using less related progenitor species genomes.

Ethylene-mediated stomatal responses to dehydration and rehydration in seed plants.

Ethylene, a plant hormone that significantly influences both plant growth and response to stress, plays a well-established role in stress signaling. However, its impact on stomatal opening and closure during dehydration and rehydration remains relatively unexplored and is still debated. Exogenous ethylene has been proven to induce stomatal closure through a series of signaling pathways, including the accumulation of reactive oxygen species (ROS), subsequent synthesis of nitric oxide (NO) and hydrogen sulfide (H2S), and SLOW ANION CHANNEL-ASSOCIATED 1 (SLAC1) activation. Thus, it has been suggested that ethylene might function to induce stomatal closure synergistically with abscisic acid (ABA). Furthermore, it has also been shown that increased ethylene can inhibit ABA- and jasmonic acid (JA)-induced stomatal closure, thus hindering drought-induced closure during dehydration. Simultaneously, other stresses, such as chilling, ozone pollution and K+ deficiency, inhibit drought and ABA-induced stomatal closure through an ethylene synthesis dependent way. However, ethylene has been shown to take on an opposing role during rehydration, preventing stomatal opening in the absence of ABA through its own signaling pathway. These findings offer novel insights into the function of ethylene in stomatal regulation during dehydration and rehydration, gaining a better understanding of the mechanisms underlying ethylene-induced stomatal movement in seed plants.

The ubiquitin E3 ligase SR1 modulates the submergence response by degrading phosphorylated WRKY33 in Arabidopsis.

Oxygen deprivation caused by flooding activates acclimation responses to stress and restricts plant growth. After experiencing flooding stress, plants must restore normal growth; however, which genes are dynamically and precisely controlled by flooding stress remains largely unknown. Here, we show that the Arabidopsis thaliana ubiquitin E3 ligase SUBMERGENCE RESISTANT1 (SR1) regulates the stability of the transcription factor WRKY33 to modulate the submergence response. SR1 physically interacts with WRKY33 in vivo and in vitro and controls its ubiquitination and proteasomal degradation. Both the sr1 mutant and WRKY33 overexpressors exhibited enhanced submergence tolerance and enhanced expression of hypoxia-responsive genes. Genetic experiments showed that WRKY33 functions downstream of SR1 during the submergence response. Submergence induced the phosphorylation of WRKY33, which enhanced the activation of RAP2.2, a positive regulator of hypoxia-response genes. Phosphorylated WRKY33 and RAP2.2 were degraded by SR1 and the N-degron pathway during reoxygenation, respectively. Taken together, our findings reveal that the on-and-off module SR1-WRKY33-RAP2.2 is connected to the well-known N-degron pathway to regulate acclimation to submergence in Arabidopsis. These two different but related modulation cascades precisely balance submergence acclimation with normal plant growth.

Phylogenomics reveals patterns of ancient hybridization and differential diversification that contribute to phylogenetic conflict in willows, poplars, and close relatives.

Despite the economic, ecological, and scientific importance of the genera Salix L. (willows) and Populus L. (poplars, cottonwoods, and aspens) Salicaceae, we know little about the sources of differences in species diversity between the genera and of the phylogenetic conflict that often confounds estimating phylogenetic trees. Salix subgenera and sections, in particular, have been difficult to classify, with one recent attempt termed a "spectacular failure" due to a speculated radiation of the subgenera Vetrix and Chamaetia. Here, we use targeted sequence capture to understand the evolutionary history of this portion of the Salicaceae plant family. Our phylogenetic hypothesis was based on 787 gene regions and identified extensive phylogenetic conflict among genes. Our analysis supported some previously described subgeneric relationships and confirmed the polyphyly of others. Using an fbranch analysis, we identified several cases of hybridization in deep branches of the phylogeny, which likely contributed to discordance among gene trees. In addition, we identified a rapid increase in diversification rate near the origination of the Vetrix-Chamaetia clade in Salix. This region of the tree coincided with several nodes that lacked strong statistical support, indicating a possible increase in incomplete lineage sorting due to rapid diversification. The extraordinary level of both recent and ancient hybridization in both Salix and Populus have played important roles in the diversification and diversity in these two genera.

Copper-Catalyzed Decarboxylation Cross-Coupling Cascade Reaction for Synthesis of Fused Dihydro-benzoxazinones.

A protocol for a tandem copper-catalyzed intermolecular decarboxylation cross-coupling cascade between o-bromobenzoic acids and proline or piperic acid has been disclosed. The developed protocol allows access to a variety of synthetically useful fused benzoxazinones scaffolds with high efficiency and good functional group compatibility. A mechanistically sequential approach for the decarboxylation and dehydration coupling process was presented.

Silver-Catalyzed Radical Umpolung Cross-Coupling of Silyl Enol Ethers with Activated Methylene Compounds: Access to Diverse Tricarbonyl Derivatives.

A silver-catalyzed protocol for the intermolecular radical umpolung cross-coupling protocol of silyl enol ethers with activated methylene compounds is disclosed. The protocol exhibits excellent functional group tolerance, enabling the expedient preparation of a variety of tricarbonyl compounds. Preliminary mechanistic investigations suggest that the reaction proceeds through a process involving free radicals in which silver oxide has a dual role, acting as both a catalyst and a base.

Genome evolution of the psammophyte Pugionium for desert adaptation and further speciation.

Deserts exert strong selection pressures on plants, but the underlying genomic drivers of ecological adaptation and subsequent speciation remain largely unknown. Here, we generated de novo genome assemblies and conducted population genomic analyses of the psammophytic genus Pugionium (Brassicaceae). Our results indicated that this bispecific genus had undergone an allopolyploid event, and the two parental genomes were derived from two ancestral lineages with different chromosome numbers and structures. The postpolyploid expansion of gene families related to abiotic stress responses and lignin biosynthesis facilitated environmental adaptations of the genus to desert habitats. Population genomic analyses of both species further revealed their recent divergence with continuous gene flow, and the most divergent regions were found to be centered on three highly structurally reshuffled chromosomes. Genes under selection in these regions, which were mainly located in one of the two subgenomes, contributed greatly to the interspecific divergence in microhabitat adaptation.

Haplotype-resolved genomes of two buckwheat crops provide insights into their contrasted rutin concentrations and reproductive systems.

BACKGROUND: Two widely cultivated annual buckwheat crops, Fagopyrum esculentum and F. tataricum, differ from each other in both rutin concentration and reproductive system. However, the underlying genetic mechanisms remain poorly elucidated. RESULTS: Here, we report the first haplotype-resolved chromosome-level genome assemblies of the two species. Two haplotype genomes of F. esculentum were assembled as 1.23 and 1.19 Gb with N50 = 9.8 and 12.4 Mb, respectively; the two haplotype genomes of F. tataricum were 453.7 and 446.2 Mb with N50 = 50 and 30 Mb, respectively. We further annotated protein-coding genes of each haplotype genome based on available gene sets and 48 newly sequenced transcriptomes. We found that more repetitive sequences, especially expansion of long terminal repeat retrotransposons (LTR-RTs), contributed to the large genome size of F. esculentum. Based on the well-annotated sequences, gene expressions, and luciferase experiments, we identified the sequence mutations of the promoter regions of two key genes that are likely to have greatly contributed to the high rutin concentration and selfing reproduction in F. tartaricum. CONCLUSIONS: Our results highlight the importance of high-quality genomes to identify genetic mutations underlying phenotypic differences between closely related species. F. tataricum may have been experienced stronger selection than F. esculentum through choosing these two non-coding alleles for the desired cultivation traits. These findings further suggest that genetic manipulation of the non-coding promoter regions could be widely employed for breeding buckwheat and other crops.

ACBP4-WRKY70-RAP2.12 module positively regulates submergence-induced hypoxia response in Arabidopsis thaliana.

ACYL-CoA-BINDING PROTEINs (ACBPs) play crucial regulatory roles during plant response to hypoxia, but their molecular mechanisms remain poorly understood. Our study reveals that ACBP4 serves as a positive regulator of the plant hypoxia response by interacting with WRKY70, influencing its nucleocytoplasmic shuttling in Arabidopsis thaliana. Furthermore, we demonstrate the direct binding of WRKY70 to the ACBP4 promoter, resulting in its upregulation and suggesting a positive feedback loop. Additionally, we pinpointed a phosphorylation site at Ser638 of ACBP4, which enhances submergence tolerance, potentially by facilitating WRKY70's nuclear shuttling. Surprisingly, a natural variation in this phosphorylation site of ACBP4 allowed A. thaliana to adapt to humid conditions during its historical demographic expansion. We further observed that both phosphorylated ACBP4 and oleoyl-CoA can impede the interaction between ACBP4 and WRKY70, thus promoting WRKY70's nuclear translocation. Finally, we found that the overexpression of orthologous BnaC5.ACBP4 and BnaA7.WRKY70 in Brassica napus increases submergence tolerance, indicating their functional similarity across genera. In summary, our research not only sheds light on the functional significance of the ACBP4 gene in hypoxia response, but also underscores its potential utility in breeding flooding-tolerant oilseed rape varieties.

Demographic History and Natural Selection Shape Patterns of Deleterious Mutation Load and Barriers to Introgression across Populus Genome.

Hybridization and resulting introgression are important processes shaping the tree of life and appear to be far more common than previously thought. However, how the genome evolution was shaped by various genetic and evolutionary forces after hybridization remains unresolved. Here we used whole-genome resequencing data of 227 individuals from multiple widespread Populus species to characterize their contemporary patterns of hybridization and to quantify genomic signatures of past introgression. We observe a high frequency of contemporary hybridization and confirm that multiple previously ambiguous species are in fact F1 hybrids. Seven species were identified, which experienced different demographic histories that resulted in strikingly varied efficacy of selection and burdens of deleterious mutations. Frequent past introgression has been found to be a pervasive feature throughout the speciation of these Populus species. The retained introgressed regions, more generally, tend to contain reduced genetic load and to be located in regions of high recombination. We also find that in pairs of species with substantial differences in effective population size, introgressed regions are inferred to have undergone selective sweeps at greater than expected frequencies in the species with lower effective population size, suggesting that introgression likely have higher potential to provide beneficial variation for species with small populations. Our results, therefore, illustrate that demography and recombination have interplayed with both positive and negative selection in determining the genomic evolution after hybridization.

The C2 H2 -type zinc finger transcription factor OSIC1 positively regulates stomatal closure under osmotic stress in poplar.

Salt and drought impair plant osmotic homeostasis and greatly limit plant growth and development. Plants decrease stomatal aperture to reduce water loss and maintain osmotic homeostasis, leading to improved stress tolerance. Herein, we identified the C2 H2 transcription factor gene OSMOTIC STRESS INDUCED C2 H2 1 (OSIC1) from Populus alba var. pyramidalis to be induced by salt, drought, polyethylene glycol 6000 (PEG6000) and abscisic acid (ABA). Overexpression of OSIC1 conferred transgenic poplar more tolerance to high salinity, drought and PEG6000 treatment by reducing stomatal aperture, while its mutant generated by the CRISPR/Cas9 system showed the opposite phenotype. Furthermore, OSIC1 directly up-regulates PalCuAOζ in vitro and in vivo, encoding a copper-containing polyamine oxidase, to enhance H2 O2 accumulation in guard cells and thus modulates stomatal closure when stresses occur. Additionally, ABA-, drought- and salt-induced PalMPK3 phosphorylates OSIC1 to increase its transcriptional activity to PalCuAOζ. This regulation of OSIC1 at the transcriptional and protein levels guarantees rapid stomatal closure when poplar responds to osmotic stress. Our results revealed a novel transcriptional regulatory mechanism of H2 O2 production in guard cells mediated by the OSIC1-PalCuAOζ module. These findings deepen our understanding of how perennial woody plants, like poplar, respond to osmotic stress caused by salt and drought and provide potential targets for breeding.

Multi-omics reveal differentiation and maintenance of dimorphic flowers in an alpine plant on the Qinghai-Tibet Plateau.

Dimorphic flowers growing on a single individual plant play a critical role in extreme adaption and reproductive assurance in plants and have high ecological and evolutionary significance. However, the omics bases underlying such a differentiation and maintenance remain largely unknown. We aimed to investigate this through genomic, transcriptome and metabolomic analyses of dimorphic flowers in an alpine biennial, Sinoswertia tetraptera (Gentianaceae). A high-quality chromosome-level genome sequence (903 Mb) was first assembled for S. tetraptera with 31,359 protein-coding genes annotated. Two rounds of recent independent whole-genome duplication (WGD) were revealed. Numerous genes from the recent species-specific WGD were found to be differentially expressed in the two types of flowers, and this may have helped contribute to the origin of this innovative trait. The genes with contrasting expressions between flowers were related to biosynthesis of hormones, floral pigments (carotenoids and flavonoids) and iridoid compounds, which are involved in both flower development and colour. Metabolomic analyses similarly suggested differential concentrations of these chemicals in the two types of flowers. The expression interactions between multiple genes may together lead to contrasting morphology and chemical concentration and open versus closed pollination of the dimorphic flowers in this species for reproductive assurance.

Out of the Qinghai-Tibet plateau: Genomic biogeography of the alpine monospecific genus Megadenia (Biscutelleae, Brassicaceae).

Numerous high-elevation alpine plants of the Qinghai-Tibet Plateau (QTP) also have disjunct distribution in adjacent low-altitude mountains. The out-of-QTP versus into-the-QTP hypothesis of alpine plants provide strong evidence for the highly disputed assumption of the massive ice sheet developed in the central plateau during the Last Glacial Maximum (LGM). In this study, we sequenced the genomes of most known populations of Megadenia, a monospecific alpine genus of Brassicaceae distributed primarily in the QTP, though rarely found in adjacent low-elevation mountains of north China and Russia (NC-R). All sequenced samples clustered into four geographic genetic groups: one pair was in the QTP and another was in NC-R. The latter pair is nested within the former, and these findings support the out-of-QTP hypothesis. Dating the four genetic groups and niche distribution suggested that Megadenia migrated out of the QTP to adjacent regions during the LGM. The NC-R group showed a decrease in the effective population sizes. In addition, the genes with high genetic divergences in the QTP group were mainly involved in habitat adaptations during low-altitude colonization. These findings reject the hypothesis of development massive ice sheets, and support glacial survival of alpine plants within, as well as further migration out of, the QTP.

Drivers of genomic landscapes of differentiation across a Populus divergence gradient.

Speciation, the continuous process by which new species form, is often investigated by looking at the variation of nucleotide diversity and differentiation across the genome (hereafter genomic landscapes). A key challenge lies in how to determine the main evolutionary forces at play shaping these patterns. One promising strategy, albeit little used to date, is to comparatively investigate these genomic landscapes as progression through time by using a series of species pairs along a divergence gradient. Here, we resequenced 201 whole-genomes from eight closely related Populus species, with pairs of species at different stages along the divergence gradient to learn more about speciation processes. Using population structure and ancestry analyses, we document extensive introgression between some species pairs, especially those with parapatric distributions. We further investigate genomic landscapes, focusing on within-species (i.e. nucleotide diversity and recombination rate) and among-species (i.e. relative and absolute divergence) summary statistics of diversity and divergence. We observe relatively conserved patterns of genomic divergence across species pairs. Independent of the stage across the divergence gradient, we find support for signatures of linked selection (i.e. the interaction between natural selection and genetic linkage) in shaping these genomic landscapes, along with gene flow and standing genetic variation. We highlight the importance of investigating genomic patterns on multiple species across a divergence gradient and discuss prospects to better understand the evolutionary forces shaping the genomic landscapes of diversity and differentiation.

PtoWRKY40 interacts with PtoPHR1-LIKE3 while regulating the phosphate starvation response in poplar.

Plants usually suffer from phosphorus starvation because of the low inorganic phosphate (Pi) status of most soils. To cope with this, plants have evolved an adaptive phosphate starvation response (PSR) which involves both developmental and metabolic changes regulated mainly by PHOSPHATE STARVATION RESPONSE1 (PHR1) and its homologs. Here, we elucidated how perennial woody plants, such as poplars (Populus spp.), respond to low-Pi stress. We first performed RNA-seq analysis of low-Pi-treated poplars and identified PtoWRKY40 is rapidly downregulated and protein degraded after stress. Overexpressing and knocking-down PtoWRKY40 downregulated and upregulated the expression of Pi starvation signaling genes, respectively, such as PHOSPHATE TRANSPORTER1 (PHT1)-type genes and PURPLE ACID PHOSPHATASE genes. PtoWRKY40 bound to the W box in the promoter of several PtoPHT1s and repressed their expression. Moreover, PtoWRKY40 interacted with PtoPHR1-LIKE3 (PtoPHL3), a PHR1 homolog in poplar, to inhibit the latter binding to the P1BS element and thus reduced PtoPHT1s' transcription under Pi-sufficient conditions. However, Pi deficiency decreased PtoWRKY40 abundance and therefore released its inhibition on PHT1s. In conclusion, we have uncovered a PSR mechanism mediated by PtoWRKY40 and PtoPHL3 which regulates Pi content in poplars, deepening our understanding of how poplars adapt to diverse Pi conditions and regulate appropriate responses to maintain Pi homeostasis.