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.

Integrated genomic, transcriptomic, and metabolomic analyses of Ilex hylonoma provide insights into the triterpenoid saponin biosynthesis.

Ilex is known for its rich content of secondary metabolites, particularly triterpenoid saponins. These compounds hold significant value in natural remedies and herbal medicine. However, the molecular mechanisms responsible for triterpenoid biosynthesis in plants of this genus remain largely unexplored. In this study, we successfully generated the first chromosome-scale genome of Ilex hylonoma. The assembly, comprising 20 anchored chromosomes, has an N50 contig size of 2.13 Mb and a scaffold size of 33.68 Mb. Comparative genome analyses with two other congeners with available chromosome-level genomes suggested that an end-to-end chromosome fusion event likely contributed to the reduction in chromosome number from n = 20 to n = 19 within this genus. By integrating transcriptomic and metabolomic data, we identified the gene expression patterns and metabolite profiles of I. hylonoma across three commonly utilized medicinal tissues. We subsequently pinpointed candidate genes involved in the regulation of triterpenoid saponin biosynthesis, including CYP450 genes, UGT genes, and associated transcription factors. Furthermore, yeast heterologous expression analysis revealed that ihyl08363 catalyzed the conversion of ß-amyrin into oleanolic acid, while ihyl04303 catalyzed the C-2α hydroxylation of oleanolic acid to produce maslinic acid. This integrated analysis provides valuable insights into the biosynthesis of important triterpenoid saponins in medicinal Ilex plants.

Geography, altitude, agriculture and hypoxia.

Reduced oxygen availability (hypoxia) represents a key plant abiotic stress in natural and agricultural systems, but conversely it is also an important component of normal growth and development. We review recent advances that demonstrate how genetic adaptations associated with hypoxia impact the known plant oxygen sensing mechanism through the PLANT CYSTEINE OXIDASE (PCO) N-degron pathway. Only three protein substrates of this pathway have been identified, and all adaptations identified to date are associated with the most important of these, the group VII ETHYLENE RESPONSE FACTOR transcription factors (ERFVIIs). We discuss how geography, altitude, and agriculture have all shaped molecular responses to hypoxia, and how these responses have emerged at different taxonomic levels through the evolution of land plants. Understanding how ecological and agricultural genetic variation acts positively to enhance hypoxia tolerance will provide novel tools and concepts to improve the performance of crops in the face of increasing extreme flooding events.

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.

Divergent Flowering Time Responses to Increasing Temperatures Are Associated With Transcriptome Plasticity and Epigenetic Modification Differences at FLC Promoter Region of Arabidopsis thaliana.

Understanding the genetic, and transcriptomic changes that drive the phenotypic plasticity of fitness traits is a central question in evolutionary biology. In this study, we utilised 152 natural Swedish Arabidopsis thaliana accessions with re-sequenced genomes, transcriptomes and methylomes and measured flowering times (FTs) under two temperature conditions (10°C and 16°C) to address this question. We revealed that the northern accessions exhibited advanced flowering in response to decreased temperature, whereas the southern accessions delayed their flowering, indicating a divergent flowering response. This contrast in flowering responses was associated with the isothermality of their native ranges, which potentially enables the northern accessions to complete their life cycle more rapidly in years with shorter growth seasons. At the transcriptome level, we observed extensive rewiring of gene co-expression networks, with the expression of 25 core genes being associated with the mean FT and its plastic variation. Notably, variations in FLC expression sensitivity between northern and southern accessions were found to be associated with the divergence FT response. Further analysis suggests that FLC expression sensitivity is associated with differences in CG, CHG and CHH methylation at the promoter region. Overall, our study revealed the association between transcriptome plasticity and flowering time plasticity among different accessions, providing evidence for its relevance in ecological adaptation. These findings offer deeper insights into the genetics of rapid responses to environmental changes and ecological adaptation.

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.