The evolutionary significance of polyploidy

Patterns of genomic variation reveal a single evolutionary origin of the wild allotetraploid Mimulus sookensis

Polyploidy occurs across the tree of life and is especially common in plants. Because newly formed cytotypes are often incompatible with their progenitors, polyploidy is also said to trigger "instantaneous" speciation. If a polyploid can self-fertilize or reproduce asexually, it is even possible for one individual to produce an entirely new lineage, but how often this scenario occurs is unclear. Here, we investigate the evolutionary history of the wild allotetraploid Mimulus sookensis, which was formed through hybridization between self-compatible, diploid species in the Mimulus guttatus complex. We generate a chromosome-scale reference assembly for M. sookensis and define its distinct subgenomes. Despite previous reports suggesting multiple origins of this highly selfing polyploid, we discover patterns of population genomic variation that provide unambiguous support for a single origin. One M. sookensis subgenome is clearly derived from the selfer Mimulus nasutus, which organellar variation suggests is the maternal progenitor. The ancestor of the other subgenome is less certain, but it shares variation with both Mimulus decorus and M. guttatus, two outcrossing diploids with geographic ranges that overlap broadly with M. sookensis. This study establishes M. sookensis as an example of instantaneous speciation, likely facilitated by the polyploid's predisposition to self-fertilize.

Complex Polyploids: Origins, Genomic Composition, and Role of Introgressed Alleles

Introgression allows polyploid species to acquire new genomic content from diploid progenitors or from other unrelated diploid or polyploid lineages, contributing to genetic diversity and facilitating adaptive allele discovery. In some cases, high levels of introgression elicit the replacement of large numbers of alleles inherited from the polyploid's ancestral species, profoundly reshaping the polyploid's genomic composition. In such complex polyploids, it is often difficult to determine which taxa were the progenitor species and which taxa provided additional introgressive blocks through subsequent hybridization. Here, we use population-level genomic data to reconstruct the phylogenetic history of Betula pubescens (downy birch), a tetraploid species often assumed to be of allopolyploid origin and which is known to hybridize with at least four other birch species. This was achieved by modeling polyploidization and introgression events under the multispecies coalescent and then using an approximate Bayesian computation rejection algorithm to evaluate and compare competing polyploidization models. We provide evidence that B. pubescens is the outcome of an autoploid genome doubling event in the common ancestor of B. pendula and its extant sister species, B. platyphylla, that took place approximately 178,000-188,000 generations ago. Extensive hybridization with B. pendula, B. nana, and B. humilis followed in the aftermath of autopolyploidization, with the relative contribution of each of these species to the B. pubescens genome varying markedly across the species' range. Functional analysis of B. pubescens loci containing alleles introgressed from B. nana identified multiple genes involved in climate adaptation, while loci containing alleles derived from B. humilis revealed several genes involved in the regulation of meiotic stability and pollen viability in plant species.

Subfunctionalisation and self-repression of duplicated E1 homologues finetunes soybean flowering and adaptation

Soybean is a photoperiod-sensitive staple crop. Its photoperiodic flowering has major consequences for latitudinal adaptation and grain yield. Here, we identify and characterise a flowering locus named Time of flower 4b (Tof4b), which encodes E1-Like b (E1Lb), a homologue of the key soybean floral repressor E1. Tof4b protein physically associates with the promoters of two FLOWERING LOCUS T (FT) genes to repress their transcription and delay flowering to impart soybean adaptation to high latitudes. Three E1 homologues undergo subfunctionalisation and show differential subcellular localisation. Moreover, they all possess self-repression capability and each suppresses the two homologous counterparts. Subfunctionalisation and the transcriptional regulation of E1 genes collectively finetune flowering time and high-latitude adaptation in soybean. We propose a model for the functional fate of the three E1 genes after the soybean whole-genome duplication events, refine the molecular mechanisms underlying high-latitude adaption, and provide a potential molecular-breeding resource.

Species-tree topology impacts the inference of ancient whole-genome duplications across the angiosperm phylogeny

PREMISE

The history of angiosperms is marked by repeated rounds of ancient whole-genome duplications (WGDs). Here we used state-of-the-art methods to provide an up-to-date view of the distribution of WGDs in the history of angiosperms that considers both uncertainty introduced by different WGD inference methods and different underlying species-tree hypotheses.

METHODS

We used the distribution synonymous divergences (Ks) of paralogs and orthologs from transcriptomic and genomic data to infer and place WGDs across two hypothesized angiosperm phylogenies. We further tested these WGD hypotheses with syntenic inferences and Bayesian models of duplicate gene gain and loss.

RESULTS

The predicted number of WGDs in the history of angiosperms (~170) based on the current taxon sampling is largely similar across different inference methods, but varies in the precise placement of WGDs on the phylogeny. Ks-based methods often yield alternative hypothesized WGD placements due to variation in substitution rates among lineages. Phylogenetic models of duplicate gene gain and loss are more robust to topological variation. However, errors in species-tree inference can still produce spurious WGD hypotheses, regardless of method used.

CONCLUSIONS

Here we showed that different WGD inference methods largely agree on an average of 3.5 WGD in the history of individual angiosperm species. However, the precise placement of WGDs on the phylogeny is subject to the WGD inference method and tree topology. As researchers continue to test hypotheses regarding the impacts ancient WGDs have on angiosperm evolution, it is important to consider the uncertainty of the phylogeny as well as WGD inference methods.

Amphicarpic development in Cardamine chenopodiifolia

Amphicarpy is an unusual trait where two fruit types develop on the same plant: one above and the other belowground. This trait is not found in conventional model species. Therefore, its development and molecular genetics remain under-studied. Here, we establish the allooctoploid Cardamine chenopodiifolia as an emerging experimental system to study amphicarpy. We characterized C. chenopodiifolia development, focusing on differences in morphology and cell wall histochemistry between above- and belowground fruit. We generated a reference transcriptome with PacBio full-length transcript sequencing and analysed differential gene expression between above- and belowground fruit valves. Cardamine chenopodiifolia has two contrasting modes of seed dispersal. The main shoot fails to bolt and initiates floral primordia that grow underground where they self-pollinate and set seed. By contrast, axillary shoots bolt and develop exploding seed pods aboveground. Morphological differences between aerial explosive fruit and subterranean nonexplosive fruit were reflected in a large number of differentially regulated genes involved in photosynthesis, secondary cell wall formation and defence responses. Tools established in C. chenopodiifolia, such as a reference transcriptome, draft genome assembly and stable plant transformation, pave the way to study amphicarpy and trait evolution via allopolyploidy.

Ploidy variation induces butterfly effect on chromatin topology in wheat

Polyploidy is a prominent driver of plant diversification, accompanied with dramatic chromosomal rearrangement and epigenetic changes that affect gene expression. How chromatin interactions within and between subgenomes adapt to ploidy transition remains poorly understood. We generate open chromatin interaction maps for natural hexaploid wheat (AABBDD), extracted tetraploid wheat (AABB), diploid wheat progenitor Aegilops tauschii (DD) and resynthesized hexaploid wheat (RHW, AABBDD). Thousands of intra- and interchromosomal loops are de novo established or disappeared in AB subgenomes after separation of D subgenome, in which 37-95% of novel loops are lost again in RHW after merger of D genome. Interestingly, more than half of novel loops are formed by cascade reactions that are triggered by disruption of chromatin interaction between AB and D subgenomes. The interaction repressed genes in RHW relative to DD are expression suppressed, resulting in more balanced expression of the three homoeologs in RHW. The interaction levels of cascade anchors are decreased step-by-step. Leading single nucleotide polymorphisms of yield- and plant architecture-related quantitative trait locus are significantly enriched in cascade anchors. The expression of 116 genes interacted with these anchors are significantly correlated with the corresponding traits. Our findings reveal trans-regulation of intrachromosomal loops by interchromosomal interactions during genome merger and separation in polyploid species.

Influence of polyploidy on morphology and distribution of the Cypress Spurge (Euphorbia cyparissias, Euphorbiaceae)

Polyploidy can cause differences in phenotypic and physiological traits among different cytotypes of the same species. Polyploids may have larger organs or occupy different ecological niches than their diploid counterparts, therefore they are hypothesized to have larger distributions or prosper in stressful environments, such as higher elevations. The Cypress spurge (Euphorbia cyparissias L.; Euphorbiaceae) is a widespread European heteroploid species including di- (2x), tetra- (4x) and hexaploid (6x) cytotypes. We tested the hypotheses that polyploids are more widespread and more abundant at higher elevations and have larger organs than their diploid ancestors in the case of E. cyparissias. We also analysed whether genome downsizing had occurred after polyploidisation. We conducted a comprehensive geographic sampling of 617 populations of E. cyparissias throughout Europe. We estimated their relative genome size using flow cytometry and inferred ploidy level of each population. We scored 13 morphological traits of vegetative and seed characters and performed statistical analyses. The study indicates that polyploidisation facilitated colonisation of new areas in E. cyparissias, where the tetraploids are most widespread, whereas the diploids are limited to putative Pleistocene refugia, mostly in southern Europe. On the other hand, the three ploidies do not differ in their elevational distribution. Although some quantitative morphological traits exhibited an increasing trend with increasing ploidy, most traits did not differ significantly among the three ploidies, and there was no overall phenotypic differentiation among them. Given that individuals of different ploidies thrive in similar habitats across the same elevations, we suggest that ecological segregation following polyploidisation is a more important trigger for morphological differentiation than polyploidisation itself in autopolyploid plants. The study demonstrates that polyploidisation can be crucial for the colonisation of new areas and for range expansion, but it does not necessarily influence elevational distribution nor confer a different phenotype.

Clonal growth characteristics and diversity patterns of different Clintonia udensis (Liliaceae) diploid and tetraploid cytotypes in the Hualongshan Mountains

Polyploidization plays an important role in plant evolution and biodiversity. However, intraspecific polyploidy compared to interspecific polyploidy received less attention. Clintonia udensis (Liliaceae) possess diploid (2n = 2x = 14) and autotetraploid (2n = 4x = 28) cytotypes. In the Hualongshan Mountains, the autotetraploids grew on the northern slope, while the diploids grew on the southern slopes. The clonal growth characteristics and clonal architecture were measured and analyzed by field observations and morphological methods. The diversity level and differentiation patterns for two different cytotypes were investigated using SSR markers. The results showed that the clonal growth parameters, such as the bud numbers of each rhizome node and the ratio of rhizome branches in the autotetraploids were higher than those in the diploids. Both the diploids and autotetraploids appeared phalanx clonal architectures with short internodes between ramets. However, the ramets or genets of the diploids had a relatively scattered distribution, while those of the autotetraploids were relatively clumping. The diploids and autotetraploids all allocated more biomass to their vegetative growth. The diploids had a higher allocation to reproductive organs than that of autotetraploids, which indicated that the tetraploids invested more resources in clonal reproduction than diploids. The clone diversity and genetic diversity of the autotetraploids were higher than that of the diploids. Significant genetic differentiation between two different cytotypes was observed (P < 0.01). During establishment and evolution, C. udensis autotetraploids employed more clumping phalanx clonal architecture and exhibited more genetic variation than the diploids.

Shared single copy genes are generally reliable for inferring phylogenetic relationships among polyploid taxa

Polyploidy, or whole-genome duplication, is expected to confound the inference of species trees with phylogenetic methods for two reasons. First, the presence of retained duplicated genes requires the reconciliation of the inferred gene trees to a proposed species tree. Second, even if the analyses are restricted to shared single copy genes, the occurrence of reciprocal gene loss, where the surviving genes in different species are paralogs from the polyploidy rather than orthologs, will mean that such genes will not have evolved under the corresponding species tree and may not produce gene trees that allow inference of that species tree. Here we analyze three different ancient polyploidy events, using synteny-based inferences of orthology and paralogy to infer gene trees from nearly 17,000 sets of homologous genes. We find that the simple use of single copy genes from polyploid organisms provides reasonably robust phylogenetic signals, despite the presence of reciprocal gene losses. Such gene trees are also most often in accord with the inferred species relationships inferred from maximum likelihood models of gene loss after polyploidy: a completely distinct phylogenetic signal present in these genomes. As seen in other studies, however, we find that methods for inferring phylogenetic confidence yield high support values even in cases where the underlying data suggest meaningful conflict in the phylogenetic signals.

Ecological correlates of population genetics in Linum suffruticosum, an heterostylous polyploid and taxonomic complex endemic to the Western Mediterranean Basin

Linum suffruticosum s.l. is a taxonomic complex widespread in the Western Mediterranean basin. The complex is characterized by a high phenotypic and cytogenetic diversity, and by a unique three-dimensional heterostyly system that makes it an obligate outcrosser. We studied the patterns of genetic diversity and structure of populations throughout the entire distribution of L. suffruticosum s.l. with microsatellite markers. We analysed their relationships with various biological and ecological variables, including the morph ratio and sex organ reciprocity of populations measured with a novel multi-dimensional method. Populations consistently showed an approximate 1:1 morph ratio with high sex organ reciprocity and high genetic diversity. We found high genetic differentiation of populations, showing a pattern of isolation by distance. The Rif mountains in NW Africa were the most important genetic barrier. The taxonomic treatment within the group was not related to the genetic differentiation of populations, but to their environmental differentiation. Genetic diversity was unrelated to latitude, elevation, population size, niche suitability or breeding system. However, there was a clear influence of ploidy level on the genetic diversity of populations, and a seeming centre-periphery pattern in its distribution. Our results suggest that polyploidization events, high outcrossing rates, isolation by distance and important geographical barriers to gene flow have played major roles in the microevolutionary history of this species complex.

Haplotype-resolved genomes of octoploid species in Phyllanthaceae family reveal a critical role for polyploidization and hybridization in speciation

The Phyllanthaceae family comprises a diverse range of plants with medicinal, edible, and ornamental value, extensively cultivated worldwide. Polyploid species commonly occur in Phyllanthaceae. Due to the rather complex genomes and evolutionary histories, their speciation process has been still lacking in research. In this study, we generated chromosome-scale haplotype-resolved genomes of two octoploid species (Phyllanthus emblica and Sauropus spatulifolius) in Phyllanthaceae family. Combined with our previously reported one tetraploid (Sauropus androgynus) and one diploid species (Phyllanthus cochinchinensis) from the same family, we explored their speciation history. The three polyploid species were all identified as allopolyploids with subgenome A/B. Each of their two distinct subgenome groups from various species was uncovered to independently share a common diploid ancestor (Ancestor-AA and Ancestor-BB). Via different evolutionary routes, comprising various scenarios of bifurcating divergence, allopolyploidization (hybrid polyploidization), and autopolyploidization, they finally evolved to the current tetraploid S. androgynus, and octoploid S. spatulifolius and P. emblica, respectively. We further discuss the variations in copy number of alleles and the potential impacts within the two octoploids. In addition, we also investigated the fluctuation of metabolites with medical values and identified the key factor in its biosynthesis process in octoploids species. Our study reconstructed the evolutionary history of these Phyllanthaceae species, highlighting the critical roles of polyploidization and hybridization in their speciation processes. The high-quality genomes of the two octoploid species provide valuable genomic resources for further research of evolution and functional genomics.

Non-B-form DNA is associated with centromere stability in newly-formed polyploid wheat

Non-B-form DNA differs from the classic B-DNA double helix structure and plays a crucial regulatory role in replication and transcription. However, the role of non-B-form DNA in centromeres, especially in polyploid wheat, remains elusive. Here, we systematically analyzed seven non-B-form DNA motif profiles (A-phased DNA repeat, direct repeat, G-quadruplex, inverted repeat, mirror repeat, short tandem repeat, and Z-DNA) in hexaploid wheat. We found that three of these non-B-form DNA motifs were enriched at centromeric regions, especially at the CENH3-binding sites, suggesting that non-B-form DNA may create a favorable loading environment for the CENH3 nucleosome. To investigate the dynamics of centromeric non-B form DNA during the alloploidization process, we analyzed DNA secondary structure using CENH3 ChIP-seq data from newly formed allotetraploid wheat and its two diploid ancestors. We found that newly formed allotetraploid wheat formed more non-B-form DNA in centromeric regions compared with their parents, suggesting that non-B-form DNA is related to the localization of the centromeric regions in newly formed wheat. Furthermore, non-B-form DNA enriched in the centromeric regions was found to preferentially form on young LTR retrotransposons, explaining CENH3's tendency to bind to younger LTR. Collectively, our study describes the landscape of non-B-form DNA in the wheat genome, and sheds light on its potential role in the evolution of polyploid centromeres.

Deciphering recursive polyploidization in Lamiales and reconstructing their chromosome evolutionary trajectories

Lamiales is an order of core eudicots with abundant diversity, and many Lamiales plants have important medicinal and ornamental values. Here, we comparatively reanalyzed 11 Lamiales species with well-assembled genome sequences and found evidence that Lamiales plants, in addition to a hexaploidization or whole-genome triplication (WGT) shared by core eudicots, experienced further polyploidization events, establishing new groups in the order. Notably, we identified a whole-genome duplication (WGD) occurred just before the split of Scrophulariaceae from the other Lamiales families, such as Acanthaceae, Bignoniaceae, and Lamiaceae, suggesting its likely being the causal reason for the establishment and fast divergence of these families. We also found that a WGT occurred ∼68 to 78 million years ago (Mya), near the split of Oleaceae from the other Lamiales families, implying that it may have caused their fast divergence and the establishment of the Oleaceae family. Then, by exploring and distinguishing intra- and intergenomic chromosomal homology due to recursive polyploidization and speciation, respectively, we inferred that the Lamiales ancestral cell karyotype had 11 proto-chromosomes. We reconstructed the evolutionary trajectories from these proto-chromosomes to form the extant chromosomes in each Lamiales plant under study. We must note that most of the inferred 11 proto-chromosomes, duplicated during a WGD thereafter, have been well preserved in jacaranda (Jacaranda mimosifolia) genome, showing the credibility of the present inference implementing a telomere-centric chromosome repatterning model. These efforts are important to understand genome repatterning after recursive polyploidization, especially shedding light on the origin of new plant groups and angiosperm cell karyotype evolution.

Impact of whole-genome duplications on structural variant evolution in Cochlearia

Polyploidy, the result of whole-genome duplication (WGD), is a major driver of eukaryote evolution. Yet WGDs are hugely disruptive mutations, and we still lack a clear understanding of their fitness consequences. Here, we study whether WGDs result in greater diversity of genomic structural variants (SVs) and how they influence evolutionary dynamics in a plant genus, Cochlearia (Brassicaceae). By using long-read sequencing and a graph-based pangenome, we find both negative and positive interactions between WGDs and SVs. Masking of recessive mutations due to WGDs leads to a progressive accumulation of deleterious SVs across four ploidal levels (from diploids to octoploids), likely reducing the adaptive potential of polyploid populations. However, we also discover putative benefits arising from SV accumulation, as more ploidy-specific SVs harbor signals of local adaptation in polyploids than in diploids. Together, our results suggest that SVs play diverse and contrasting roles in the evolutionary trajectories of young polyploids.

The complex hexaploid oil-Camellia genome traces back its phylogenomic history and multi-omics analysis of Camellia oil biosynthesis

Oil-Camellia (Camellia oleifera), belonging to the Theaceae family Camellia, is an important woody edible oil tree species. The Camellia oil in its mature seed kernels, mainly consists of more than 90% unsaturated fatty acids, tea polyphenols, flavonoids, squalene and other active substances, which is one of the best quality edible vegetable oils in the world. However, genetic research and molecular breeding on oil-Camellia are challenging due to its complex genetic background. Here, we successfully report a chromosome-scale genome assembly for a hexaploid oil-Camellia cultivar Changlin40. This assembly contains 8.80 Gb genomic sequences with scaffold N50 of 180.0 Mb and 45 pseudochromosomes comprising 15 homologous groups with three members each, which contain 135 868 genes with an average length of 3936 bp. Referring to the diploid genome, intragenomic and intergenomic comparisons of synteny indicate homologous chromosomal similarity and changes. Moreover, comparative and evolutionary analyses reveal three rounds of whole-genome duplication (WGD) events, as well as the possible diversification of hexaploid Changlin40 with diploid occurred approximately 9.06 million years ago (MYA). Furthermore, through the combination of genomics, transcriptomics and metabolomics approaches, a complex regulatory network was constructed and allows to identify potential key structural genes (SAD, FAD2 and FAD3) and transcription factors (AP2 and C2H2) that regulate the metabolism of Camellia oil, especially for unsaturated fatty acids biosynthesis. Overall, the genomic resource generated from this study has great potential to accelerate the research for the molecular biology and genetic improvement of hexaploid oil-Camellia, as well as to understand polyploid genome evolution.

Autotetraploidy of rice does not potentiate the tolerance to drought stress in the seedling stage

Polyploid is considered an advantage that has evolved to be more environmentally adaptable than its diploid. To understand if doubled chromosome of diploid rice can improve drought tolerance, we evaluated the diploid (2X) and autotetraploid (4X) plants of three indica and three japonica varieties. Drought stress in the plastic bucket of four-leaf stage revealed that the drought tolerance of 4X plants was lower than that of its diploid donor plants. The assay of photosynthetic rate of all varieties showed that all 4X varieties had lower rates than their diploid donors. The capacity for reactive oxygen species production and scavenging varied among different 2X and 4X varieties. Further, transcriptomic analysis of 2X and 4X plants of four varieties under normal and drought condition showed that the wide variation of gene expression was caused by difference of varieties, not by chromosome ploidy. However, weighted gene co-expression network analysis (WGCNA) revealed that the severe interference of photosynthesis-related genes in tetraploid plants under drought stress is the primary reason for the decrease of drought tolerance in autotetraploid lines. Consistently, new transcripts analysis in autotetraploid revealed that the gene transcription related with mitochondrion and plastid of cell component was influenced most significantly. The results indicated that chromosome doubling of diploid rice weakened their drought tolerance, primarily due to disorder of photosynthesis-related genes in tetraploid plants under drought stress. Maintain tetraploid drought tolerance through chromosome doubling breeding in rice needs to start with the selection of parental varieties and more efforts.

Genome-wide identification of starch phosphorylase gene family in Rosa chinensis and expression in response to abiotic stress

Chinese rose (Rosa chinensis) is an important ornamental plant, with economic, cultural, and symbolic significance. During the application of outdoor greening, adverse environments such as high temperature and drought are often encountered, which affect its application scope and ornamental quality. The starch phosphorylase (Pho) gene family participate in the synthesis and decomposition of starch, not only related to plant energy metabolism, but also plays an important role in plant stress resistance. The role of Pho in combating salinity and high temperature stress in R. chinensis remains unknown. In this work, 4 Phos from R. chinensis were detected with Pfam number of Pho (PF00343.23) and predicted by homolog-based prediction (HBP). The Phos are characterized by sequence lengths of 821 to 997 bp, and the proteins are predicted to subcellularly located in the plastid and cytoplasm. The regulatory regions of the Phos contain abundant stress and phytohormone-responsive cis-acting elements. Based on transcriptome analysis, the Phos were found to respond to abiotic stress factors such as drought, salinity, high temperature, and plant phytohormone of jasmonic acid and salicylic acid. The response of Phos to abiotic stress factors such as salinity and high temperature was confirmed by qRT-PCR analysis. To evaluate the genetic characteristics of Phos, a total of 69 Phos from 17 species were analyzed and then classified into 3 groups in phylogenetic tree. The collinearity analysis of Phos in R. chinensis and other species was conducted for the first time. This work provides a view of evolution for the Pho gene family and indicates that Phos play an important role in abiotic stress response of R. chinensis.

Harnessing cold adaptation for postglacial colonisation: Galactinol synthase expression and raffinose accumulation in a polyploid and its progenitors

Allotetraploid white clover (Trifolium repens) formed during the last glaciation through hybridisation of two European diploid progenitors from restricted niches: one coastal, the other alpine. Here, we examine which hybridisation-derived molecular events may have underpinned white clover's postglacial niche expansion. We compared the transcriptomic frost responses of white clovers (an inbred line and an alpine-adapted ecotype), extant descendants of its progenitor species and a resynthesised white clover neopolyploid to identify genes that were exclusively frost-induced in the alpine progenitor and its derived subgenomes. From these analyses we identified galactinol synthase, the rate-limiting enzyme in biosynthesis of the cryoprotectant raffinose, and found that the extant descendants of the alpine progenitor as well as the neopolyploid white clover rapidly accumulated significantly more galactinol and raffinose than the coastal progenitor under cold stress. The frost-induced galactinol synthase expression and rapid raffinose accumulation derived from the alpine progenitor likely provided an advantage during early postglacial colonisation for white clover compared to its coastal progenitor.

Species delimitation 4.0: integrative taxonomy meets artificial intelligence

Although species are central units for biological research, recent findings in genomics are raising awareness that what we call species can be ill-founded entities due to solely morphology-based, regional species descriptions. This particularly applies to groups characterized by intricate evolutionary processes such as hybridization, polyploidy, or asexuality. Here, challenges of current integrative taxonomy (genetics/genomics + morphology + ecology, etc.) become apparent: different favored species concepts, lack of universal characters/markers, missing appropriate analytical tools for intricate evolutionary processes, and highly subjective ranking and fusion of datasets. Now, integrative taxonomy combined with artificial intelligence under a unified species concept can enable automated feature learning and data integration, and thus reduce subjectivity in species delimitation. This approach will likely accelerate revising and unraveling eukaryotic biodiversity.

Utilizing a comparative approach to assess genome evolution during diploidization in Artemisia tridentata, a keystone species of western North America

PREMISE

Polyploidization is often followed by diploidization. Diploidization is generally studied using synthetic polyploid lines and/or crop plants, but rarely using extant diploids or nonmodel plants such as Artemisia tridentata. This threatened western North American keystone species has a large genome compared to congeneric Artemisia species; dominated by diploid and tetraploid cytotypes, with multiple origins of tetraploids with genome size reduction.

METHODS

The genome of an A. tridentata sample was resequenced to study genome evolution and compared to that of A. annua, a diploid congener. Three diploid genomes of A. tridentata were compared to test for multiple diploidization events.

RESULTS

The A. tridentata genome had many chromosomal rearrangements relative to that of A. annua, while large-scale synteny of A. tridentata chromosome 3 and A. annua chromosome 4 was conserved. The three A. tridentata genomes had similar sizes (4.19-4.2 Gbp), heterozygosity (2.24-2.25%), and sequence (98.73-99.15% similarity) across scaffolds, and in k-mer analyses, similar patterns of diploid heterozygous k-mers (AB = 41%, 47%, and 47%), triploid heterozygous k-mers (AAB = 18-21%), and tetraploid k-mers (AABB = 13-17%). Biallelic SNPs were evenly distributed across scaffolds for all individuals. Comparisons of transposable element (TE) content revealed differential enrichment of TE clades.

CONCLUSIONS

Our findings suggest population-level TE differentiation after a shared polyploidization-to-diploidization event(s) and exemplify the complex processes of genome evolution. This research approached provides new resources for exploration of abiotic stress response, especially the roles of TEs in response pathways.

The evolution history of an allotetraploid mangrove tree analysed with a new tool Allo4D

Mangrove species are broadly classified as true mangroves and mangrove associates. The latter are amphibious plants that can survive in the intertidal zone and reproduce naturally in terrestrial environments. Their widespread distribution and extensive adaptability make them ideal research materials for exploring adaptive evolution. In this study, we de novo assembled two genomes of mangrove associates (the allotetraploid Barringtonia racemosa (2n = 4x = 52) and diploid Barringtonia asiatica (2n = 2x = 26)) to investigate the role of allopolyploidy in the evolutionary history of mangrove species. We developed a new allotetraploid-dividing tool Allo4D to distinguish between allotetraploid scaffold-scale subgenomes and verified its accuracy and reliability using real and simulated data. According to the two subgenomes of allotetraploid B. racemosa divided using Allo4D, the allopolyploidization event was estimated to have occurred approximately one million years ago (Mya). We found that B. racemosa, B. asiatica, and Diospyros lotus shared a whole genome duplication (WGD) event during the K-Pg (Cretaceous-Paleozoic) period. K-Pg WGD and recent allopolyploidization events contributed to the speciation of B. racemosa and its adaptation to coastal habitats. We found that genes in the glucosinolates (GSLs) pathway, an essential pathway in response to various biotic and abiotic stresses, expanded rapidly in B. racemosa during polyploidization. In summary, this study provides a typical example of the adaptation of allopolyploid plants to extreme environmental conditions. The newly developed tool, Allo4D, can effectively divide allotetraploid subgenomes and explore the evolutionary history of polyploid plants, especially for species whose ancestors are unknown or extinct.

Genomic prediction of synthetic hexaploid wheat upon tetraploid durum and diploid Aegilops parental pools

Bread wheat (Triticum aestivum L.) is a globally important food crop, which was domesticated about 8-10,000 years ago. Bread wheat is an allopolyploid, and it evolved from two hybridization events of three species. To widen the genetic base in breeding, bread wheat has been re-synthesized by crossing durum wheat (Triticum turgidum ssp. durum) and goat grass (Aegilops tauschii Coss), leading to so-called synthetic hexaploid wheat (SHW). We applied the quantitative genetics tools of "hybrid prediction"-originally developed for the prediction of wheat hybrids generated from different heterotic groups - to a situation of allopolyploidization. Our use-case predicts the phenotypes of SHW for three quantitatively inherited global wheat diseases, namely tan spot (TS), septoria nodorum blotch (SNB), and spot blotch (SB). Our results revealed prediction abilities comparable to studies in 'traditional' elite or hybrid wheat. Prediction abilities were highest using a marker model and performing random cross-validation, predicting the performance of untested SHW (0.483 for SB to 0.730 for TS). When testing parents not necessarily used in SHW, combination prediction abilities were slightly lower (0.378 for SB to 0.718 for TS), yet still promising. Despite the limited phenotypic data, our results provide a general example for predictive models targeting an allopolyploidization event and a method that can guide the use of genetic resources available in gene banks.

Surviving the Storm: The Role of Poly- and Depolyploidization in Tissues and Tumors

Polyploidization and depolyploidization are critical processes in the normal development and tissue homeostasis of diploid organisms. Recent investigations have revealed that polyaneuploid cancer cells (PACCs) exploit this ploidy variation as a survival strategy against anticancer treatment and for the repopulation of tumors. Unscheduled polyploidization and chromosomal instability in PACCs enhance malignancy and treatment resistance. However, their inability to undergo mitosis causes catastrophic cellular death in most PACCs. Adaptive ploid reversal mechanisms, such as multipolar mitosis, centrosome clustering, meiosis-like division, and amitosis, counteract this lethal outcome and drive cancer relapse. The purpose of this work is to focus on PACCs induced by cytotoxic therapy, highlighting the latest discoveries in ploidy dynamics in physiological and pathological contexts. Specifically, by emphasizing the role of "poly-depolyploidization" in tumor progression, the aim is to identify novel therapeutic targets or paradigms for combating diseases associated with aberrant ploidies.

Assess different fermentation characteristics of 54 lager yeasts based on group classification

Saccharomyces pastorianus, hybrids of Saccharomyces cerevisiae and Saccharomyces eubayanus, were generally regarded as authentic lager beer yeasts. In recent years, with more new findings of other Saccharomyces genus hybrids, yeasts used in lager beer brewing have been proved much more complicated than previous cognition. In this study, we analyzed the different fermentation characteristics of 54 yeast strains used for lager brewing in normal and very high gravity brewing based on group classification. The difference between Group Ⅰ and Group Ⅱ lager yeasts were more striking in very high gravity brewing. However, during our research progress, we realized that some yeasts used in this study were actually hybrids of S. cerevisiae and Saccharomyces kudriavzevii. Features of these hybrids could be beneficial to very high gravity brewing. We further discussed about the mechanism behind their outstanding characteristics and the reason why group classification methods of lager beer yeasts had limitations. Hybridization in yeasts is constantly getting richer. Lager yeasts could have more possibilities based on better understandings of their genetic background and roles of other Saccharomyces genus hybrids. Their heterosis shed light on innovation in brewing and other diverse fermentation industries.

Upgraded durian genome reveals the role of chromosome reshuffling during ancestral karyotype evolution, lignin biosynthesis regulation, and stress tolerance

Durian (Durio zibethinus) is a tropical fruit that has a unique flavor and aroma. It occupies a significant phylogenetic position within the Malvaceae family. Extant core-eudicot plants are reported to share seven ancestral karyotypes that have undergone reshuffling, resulting in an abundant genomic diversity. However, the ancestral karyotypes of the Malvaceae family, as well as the evolution trajectory leading to the 28 chromosomes in durian, remain poorly understood. Here, we report the high-quality assembly of the durian genome with comprehensive comparative genomic analyses. By analyzing the collinear blocks between cacao and durian, we inferred 11 Malvaceae ancestral karyotypes. These blocks were present in a single-copy form in cacao and mainly in triplicates in durian, possibly resulting from a recent whole genome triplication (WGT) event that led to hexaploidization of the durian genome around 20 (17-24) million years ago. A large proportion of the duplicated genes in durian, such as those involved in the lignin biosynthesis module for phenylpropane biosynthesis, are derived directly from whole genome duplication, which makes it an important force in reshaping its genomic architecture. Transcriptome studies have revealed that genes involved in feruloyl-CoA formations were highly preferentially expressed in fruit peels, indicating that the thorns produced on durian fruit may comprise guaiacyl and syringyl lignins. Among all the analyzed transcription factors (TFs), members of the heat shock factor family (HSF) were the most significantly upregulated under heat stress. All subfamilies of genes encoding heat shock proteins (HSPs) in the durian genome appear to have undergone expansion. The potential interactions between HSF Dzi05.397 and HSPs were examined and experimentally verified. Our study provides a high-quality durian genome and reveals the reshuffling mechanism of ancestral Malvaceae chromosomes to produce the durian genome. We also provide insights into the mechanism underlying lignin biosynthesis and heat stress tolerance.

Reproduction-associated pathways in females of gibel carp (Carassius gibelio) shed light on the molecular mechanisms of the coexistence of asexual and sexual reproduction

Gibel carp (Carassius gibelio) is a cyprinid fish that originated in eastern Eurasia and is considered as invasive in European freshwater ecosystems. The populations of gibel carp in Europe are mostly composed of asexually reproducing triploid females (i.e., reproducing by gynogenesis) and sexually reproducing diploid females and males. Although some cases of coexisting sexual and asexual reproductive forms are known in vertebrates, the molecular mechanisms maintaining such coexistence are still in question. Both reproduction modes are supposed to exhibit evolutionary and ecological advantages and disadvantages. To better understand the coexistence of these two reproduction strategies, we performed transcriptome profile analysis of gonad tissues (ovaries) and studied the differentially expressed reproduction-associated genes in sexual and asexual females. We used high-throughput RNA sequencing to generate transcriptomic profiles of gonadal tissues of triploid asexual females and males, diploid sexual males and females of gibel carp, as well as diploid individuals from two closely-related species, C. auratus and Cyprinus carpio. Using SNP clustering, we showed the close similarity of C. gibelio and C. auratus with a basal position of C. carpio to both Carassius species. Using transcriptome profile analyses, we showed that many genes and pathways are involved in both gynogenetic and sexual reproduction in C. gibelio; however, we also found that 1500 genes, including 100 genes involved in cell cycle control, meiosis, oogenesis, embryogenesis, fertilization, steroid hormone signaling, and biosynthesis were differently expressed in the ovaries of asexual and sexual females. We suggest that the overall downregulation of reproduction-associated pathways in asexual females, and their maintenance in sexual ones, allows the populations of C. gibelio to combine the evolutionary and ecological advantages of the two reproductive strategies. However, we showed that many sexual-reproduction-related genes are maintained and expressed in asexual females, suggesting that gynogenetic gibel carp retains the genetic toolkits for meiosis and sexual reproduction. These findings shed new light on the evolution of this asexual and sexual complex.

Phylogenomics and morphology of Celmisiinae (Asteraceae: Astereae): Taxonomic and evolutionary implications

The tribe Astereae (Asteraceae) includes 36 subtribes and 252 genera, and is distributed worldwide in temperate and tropical regions. One of the subtribes, Celmisiinae Saldivia, has been recently circumscribed to include six genera and ca. 160 species, and is restricted to eastern Australia, New Zealand, and New Guinea. The species show an impressive range of growth habit, from small herbs and ericoid subshrubs to medium-sized trees. They live in a wide range of habitats and are often dominant in subalpine and alpine vegetation. Despite the well-supported circumscription of Celmisiinae, uncertainties have remained about their internal relationships and classification at genus and species levels. This study exploited recent advances in high-throughput sequencing to build a robust multi-gene phylogeny for the subtribe Celmisiinae. The target enrichment Angiosperms353 bait set and the hybpiper-nf and paragone-nf pipelines were used to retrieve, infer, and assemble orthologous loci from 75 taxa representing all the main putative clades within the subtribe. Because of the diploidised ploidy level in Celmisiinae, as well as missing data in the assemblies, uncertainty remains surrounding the inference of orthology detection. However, based on a variety of gene-family sets, coalescent and concatenation-based phylogenetic reconstructions recovered similar topologies. Paralogy and missing data in the gene-families caused some problems, but the estimated phylogenies were well-supported and well-resolved. The phylogenomic evidence supported Celmisiinae and three main clades: the Pleurophyllum clade (Pleurophyllum, Macrolearia and Damnamenia), mostly in the New Zealand Subantarctic Islands, Celmisia of mainland New Zealand and Australia, and Shawia (including 'Olearia pro parte' and Pachystegia) of New Zealand, Australia and New Guinea. The results presented here add to the accumulating support for the Angiosperms353 bait set as an efficient method for documenting plant diversity.

Pan-evolutionary and regulatory genome architecture delineated by an integrated macro- and microsynteny approach

The forthcoming massive genome data generated by the Earth BioGenome Project will open up a new era of comparative genomics, for which genome synteny analysis provides an important framework. Profiling genome synteny represents an essential step in elucidating genome architecture, regulatory blocks/elements and their evolutionary history. Here we describe PanSyn, ( https://github.com/yhw320/PanSyn ), the most comprehensive and up-to-date genome synteny pipeline, providing step-by-step instructions and application examples to demonstrate its usage. PanSyn inherits both basic and advanced functions from existing popular tools, offering a user-friendly, highly customized approach for genome macrosynteny analysis and integrated pan-evolutionary and regulatory analysis of genome architecture, which are not yet available in public synteny software or tools. The advantages of PanSyn include: (i) advanced microsynteny analysis by functional profiling of microsynteny genes and associated regulatory elements; (ii) comprehensive macrosynteny analysis, including the inference of karyotype evolution from ancestors to extant species; and (iii) functional integration of microsynteny and macrosynteny for pan-evolutionary profiling of genome architecture and regulatory blocks, as well as integration with external functional genomics datasets from three- or four-dimensional genome and ENCODE projects. PanSyn requires basic knowledge of the Linux environment and Perl programming language and the ability to access a computer cluster, especially for large-scale genomic comparisons. Our protocol can be easily implemented by a competent graduate student or postdoc and takes several days to weeks to execute for dozens to hundreds of genomes. PanSyn provides yet the most comprehensive and powerful tool for integrated evolutionary and functional genomics.

Interspecific transfer of genetic information through polyploid bridges

Hybridization blurs species boundaries and leads to intertwined lineages resulting in reticulate evolution. Polyploidy, the outcome of whole genome duplication (WGD), has more recently been implicated in promoting and facilitating hybridization between polyploid species, potentially leading to adaptive introgression. However, because polyploid lineages are usually ephemeral states in the evolutionary history of life it is unclear whether WGD-potentiated hybridization has any appreciable effect on their diploid counterparts. Here, we develop a model of cytotype dynamics within mixed-ploidy populations to demonstrate that polyploidy can in fact serve as a bridge for gene flow between diploid lineages, where introgression is fully or partially hampered by the species barrier. Polyploid bridges emerge in the presence of triploid organisms, which despite critically low levels of fitness, can still allow the transfer of alleles between diploid states of independently evolving mixed-ploidy species. Notably, while marked genetic divergence prevents polyploid-mediated interspecific gene flow, we show that increased recombination rates can offset these evolutionary constraints, allowing a more efficient sorting of alleles at higher-ploidy levels before introgression into diploid gene pools. Additionally, we derive an analytical approximation for the rate of gene flow at the tetraploid level necessary to supersede introgression between diploids with nonzero introgression rates, which is especially relevant for plant species complexes, where interspecific gene flow is ubiquitous. Altogether, our results illustrate the potential impact of polyploid bridges on the (re)distribution of genetic material across ecological communities during evolution, representing a potential force behind reticulation.

A new chromosome-scale genome of wild Brassica oleracea provides insights into the domestication of Brassica crops

The cultivated diploid Brassica oleracea is an important vegetable crop, but the genetic basis of its domestication remains largely unclear in the absence of high-quality reference genomes of wild B. oleracea. Here, we report the first chromosome-level assembly of the wild Brassica oleracea L. W03 genome (total genome size, 630.7 Mb; scaffold N50, 64.6 Mb). Using the newly assembled W03 genome, we constructed a gene-based B. oleracea pangenome and identified 29 744 core genes, 23 306 dispensable genes, and 1896 private genes. We re-sequenced 53 accessions, representing six potential wild B. oleracea progenitor species. The results of the population genomic analysis showed that the wild B. oleracea populations had the highest level of diversity and represents the most closely related population to modern-day horticultural B. oleracea. In addition, the WUSCHEL gene was found to play a decisive role in domestication and to be involved in cauliflower and broccoli curd formation. We also illustrate the loss of disease-resistance genes during selection for domestication. Our results provide new insights into the domestication of B. oleracea and will facilitate the future genetic improvement of Brassica crops.

ATAD3 Proteins: Unique Mitochondrial Proteins Essential for Life in Diverse Eukaryotic Lineages

ATPase family AAA domain-containing 3 (ATAD3) proteins are unique mitochondrial proteins that arose deep in the eukaryotic lineage but that are surprisingly absent in Fungi and Amoebozoa. These ∼600-amino acid proteins are anchored in the inner mitochondrial membrane and are essential in metazoans and Arabidopsis thaliana. ATAD3s comprise a C-terminal ATPases Associated with a variety of cellular Activities (AAA+) matrix domain and an ATAD3_N domain, which is located primarily in the inner membrane space but potentially extends to the cytosol to interact with the ER. Sequence and structural alignments indicate that ATAD3 proteins are most similar to classic chaperone unfoldases in the AAA+ family, suggesting that they operate in mitochondrial protein quality control. A. thaliana has four ATAD3 genes in two distinct clades that appear first in the seed plants, and both clades are essential for viability. The four genes are generally coordinately expressed, and transcripts are highest in growing apices and imbibed seeds. Plants with disrupted ATAD3 have reduced growth, aberrant mitochondrial morphology, diffuse nucleoids and reduced oxidative phosphorylation complex I. These and other pleiotropic phenotypes are also observed in ATAD3 mutants in metazoans. Here, we discuss the distribution of ATAD3 proteins as they have evolved in the plant kingdom, their unique structure, what we know about their function in plants and the challenges in determining their essential roles in mitochondria.

The genomes of seven economic Caesalpinioideae trees provide insights into polyploidization history and secondary metabolite biosynthesis

The Caesalpinioideae subfamily contains many well-known trees that are important for economic sustainability and human health, but a lack of genomic resources has hindered their breeding and utilization. Here, we present chromosome-level reference genomes for the two food and industrial trees Gleditsia sinensis (921 Mb) and Biancaea sappan (872 Mb), the three shade and ornamental trees Albizia julibrissin (705 Mb), Delonix regia (580 Mb), and Acacia confusa (566 Mb), and the two pioneer and hedgerow trees Leucaena leucocephala (1338 Mb) and Mimosa bimucronata (641 Mb). Phylogenetic inference shows that the mimosoid clade has a much higher evolutionary rate than the other clades of Caesalpinioideae. Macrosynteny comparison suggests that the fusion and breakage of an unstable chromosome are responsible for the difference in basic chromosome number (13 or 14) for Caesalpinioideae. After an ancient whole-genome duplication (WGD) shared by all Caesalpinioideae species (CWGD, ∼72.0 million years ago [MYA]), there were two recent successive WGD events, LWGD-1 (16.2-19.5 MYA) and LWGD-2 (7.1-9.5 MYA), in L. leucocephala. Thereafter, ∼40% gene loss and genome-size contraction have occurred during the diploidization process in L. leucocephala. To investigate secondary metabolites, we identified all gene copies involved in mimosine metabolism in these species and found that the abundance of mimosine biosynthesis genes in L. leucocephala largely explains its high mimosine production. We also identified the set of all potential genes involved in triterpenoid saponin biosynthesis in G. sinensis, which is more complete than that based on previous transcriptome-derived unigenes. Our results and genomic resources will facilitate biological studies of Caesalpinioideae and promote the utilization of valuable secondary metabolites.

Whole-genome Sequencing Reveals Autooctoploidy in Chinese Sturgeon and Its Evolutionary Trajectories

The order Acipenseriformes, which includes sturgeons and paddlefishes, represents "living fossils" with complex genomes that are good models for understanding whole-genome duplication (WGD) and ploidy evolution in fishes. Here, we sequenced and assembled the first high-quality chromosome-level genome for the complex octoploid Acipenser sinensis (Chinese sturgeon), a critically endangered species that also represents a poorly understood ploidy group in Acipenseriformes. Our results show that A. sinensis is a complex autooctoploid species containing four kinds of octovalents (8n), a hexavalent (6n), two tetravalents (4n), and a divalent (2n). An analysis taking into account delayed rediploidization reveals that the octoploid genome composition of Chinese sturgeon results from two rounds of homologous WGDs, and further provides insights into the timing of its ploidy evolution. This study provides the first octoploid genome resource of Acipenseriformes for understanding ploidy compositions and evolutionary trajectories of polyploid fishes.

Evolutionary Dynamics of Chromatin Structure and Duplicate Gene Expression in Diploid and Allopolyploid Cotton

Polyploidy is a prominent mechanism of plant speciation and adaptation, yet the mechanistic understandings of duplicated gene regulation remain elusive. Chromatin structure dynamics are suggested to govern gene regulatory control. Here, we characterized genome-wide nucleosome organization and chromatin accessibility in allotetraploid cotton, Gossypium hirsutum (AADD, 2n = 4X = 52), relative to its two diploid parents (AA or DD genome) and their synthetic diploid hybrid (AD), using DNS-seq. The larger A-genome exhibited wider average nucleosome spacing in diploids, and this intergenomic difference diminished in the allopolyploid but not hybrid. Allopolyploidization also exhibited increased accessibility at promoters genome-wide and synchronized cis-regulatory motifs between subgenomes. A prominent cis-acting control was inferred for chromatin dynamics and demonstrated by transposable element removal from promoters. Linking accessibility to gene expression patterns, we found distinct regulatory effects for hybridization and later allopolyploid stages, including nuanced establishment of homoeolog expression bias and expression level dominance. Histone gene expression and nucleosome organization are coordinated through chromatin accessibility. Our study demonstrates the capability to track high-resolution chromatin structure dynamics and reveals their role in the evolution of cis-regulatory landscapes and duplicate gene expression in polyploids, illuminating regulatory ties to subgenomic asymmetry and dominance.

wgd v2: a suite of tools to uncover and date ancient polyploidy and whole-genome duplication

MOTIVATION

Major improvements in sequencing technologies and genome sequence assembly have led to a huge increase in the number of available genome sequences. In turn, these genome sequences form an invaluable source for evolutionary, ecological, and comparative studies. One kind of analysis that has become routine is the search for traces of ancient polyploidy, particularly for plant genomes, where whole-genome duplication (WGD) is rampant.

RESULTS

Here, we present a major update of a previously developed tool wgd, namely wgd v2, to look for remnants of ancient polyploidy, or WGD. We implemented novel and improved previously developed tools to (a) construct KS age distributions for the whole-paranome (collection of all duplicated genes in a genome), (b) unravel intragenomic and intergenomic collinearity resulting from WGDs, (c) fit mixture models to age distributions of gene duplicates, (d) correct substitution rate variation for phylogenetic placement of WGDs, and (e) date ancient WGDs via phylogenetic dating of WGD-retained gene duplicates. The applicability and feasibility of wgd v2 for the identification and the relative and absolute dating of ancient WGDs is demonstrated using different plant genomes.

AVAILABILITY AND IMPLEMENTATION

wgd v2 is open source and available at https://github.com/heche-psb/wgd.

Polyploidy in Cancer: Causal Mechanisms, Cancer-Specific Consequences, and Emerging Treatments

Drug resistance is the major determinant for metastatic disease and fatalities, across all cancers. Depending on the tissue of origin and the therapeutic course, a variety of biological mechanisms can support and sustain drug resistance. Although genetic mutations and gene silencing through epigenetic mechanisms are major culprits in targeted therapy, drug efflux and polyploidization are more global mechanisms that prevail in a broad range of pathologies, in response to a variety of treatments. There is an unmet need to identify patients at risk for polyploidy, understand the mechanisms underlying polyploidization, and to develop strategies to predict, limit, and reverse polyploidy thus enhancing efficacy of standard-of-care therapy that improve better outcomes. This literature review provides an overview of polyploidy in cancer and offers perspective on patient monitoring and actionable therapy.

Evolutionary divergence of subgenomes in common carp provides insights into speciation and allopolyploid success

Hybridization and polyploidization have made great contributions to speciation, heterosis, and agricultural production within plants, but there is still limited understanding and utilization in animals. Subgenome structure and expression reorganization and cooperation post hybridization and polyploidization are essential for speciation and allopolyploid success. However, the mechanisms have not yet been comprehensively assessed in animals. Here, we produced a high-fidelity reference genome sequence for common carp, a typical allotetraploid fish species cultured worldwide. This genome enabled in-depth analysis of the evolution of subgenome architecture and expression responses. Most genes were expressed with subgenome biases, with a trend of transition from the expression of subgenome A during the early stages to that of subgenome B during the late stages of embryonic development. While subgenome A evolved more rapidly, subgenome B contributed to a greater level of expression during development and under stressful conditions. Stable dominant patterns for homoeologous gene pairs both during development and under thermal stress suggest a potential fixed heterosis in the allotetraploid genome. Preferentially expressing either copy of a homoeologous gene at higher levels to confer development and response to stress indicates the dominant effect of heterosis. The plasticity of subgenomes and their shifting of dominant expression during early development, and in response to stressful conditions, provide novel insights into the molecular basis of the successful speciation, evolution, and heterosis of the allotetraploid common carp.

Convergent and/or parallel evolution of RNA-binding proteins in angiosperms after polyploidization

Increasing studies suggest that the biased retention of stress-related transcription factors (TFs) after whole-genome duplications (WGDs) could rewire gene transcriptional networks, facilitating plant adaptation to challenging environments. However, the role of posttranscriptional factors (e.g. RNA-binding proteins, RBPs) following WGDs has been largely ignored. Uncovering thousands of RBPs in 21 representative angiosperm species, we integrate genomic, transcriptomic, regulatomic, and paleotemperature datasets to unravel their evolutionary trajectories and roles in adapting to challenging environments. We reveal functional enrichments of RBP genes in stress responses and identify their convergent retention across diverse angiosperms from independent WGDs, coinciding with global cooling periods. Numerous RBP duplicates derived from WGDs are then identified as cold-induced. A significant overlap of 29 orthogroups between WGD-derived and cold-induced RBP genes across diverse angiosperms highlights a correlation between WGD and cold stress. Notably, we unveil an orthogroup (Glycine-rich RNA-binding Proteins 7/8, GRP7/8) and relevant TF duplicates (CCA1/LHY, RVE4/8, CBF2/4, etc.), co-retained in different angiosperms post-WGDs. Finally, we illustrate their roles in rewiring circadian and cold-regulatory networks at both transcriptional and posttranscriptional levels during global cooling. Altogether, we underline the adaptive evolution of RBPs in angiosperms after WGDs during global cooling, improving our understanding of plants surviving periods of environmental turmoil.

Somatic mutations that affect early genetic progression and immune microenvironment in gastric carcinoma

Gastric carcinoma (GC) is a high heterogeneity and malignant tumor with a poor prognosis. The current implementation of immunotherapy in GC is limited due to the insufficient exploration of immune-related mutations and speculated early mutation events. Therefore, we performed whole-exome sequencing on 40 patients with GC to explore their genetic characteristics, shedding light on the order of genetic events, somatic mutations impacting the immune microenvironment, and potential biomarkers for immunotherapy. Regarding genetic events, TP53 disruptions were identified as frequent and early events in GC progression, often occurring alongside other gene mutations. The mutations occurring in GANS, SMAD4, and POLE were early independent events. Patients harboring CSMD3, FAT4, FLG, KMT2C, LRP1B, MUC5B, MUC16, PLEC, RNF43, SYNE1, TP53, TTN, XIRP2, and ZFHX4 mutations tended to have decreased B cells, T cells, macrophage, neutrophil, and dendritic cells infiltration, except for the ARID1A gene mutations. We also found patients with microsatellite instability-high tumors had higher homologous recombination deficiency (HRD) scores. HRD showed a positive correlation with tumor mutational burden, which might serve as indirect evidence supporting the potential of HRD as a biomarker for GC. These findings highlighted GC's high heterogeneity and complexity and provided valuable insights into the somatic mutations that affect early genetic progression and immune microenvironment.

De novo genome assembly and population genomics of a shrub tree Barthea barthei (Hance) krass provide insights into the adaptive color variations

Flower color is a classic example of an ecologically important trait under selection in plants. Understanding the genetic mechanisms underlying shifts in flower color can provide key insights into ecological speciation. In this study, we investigated the genetic basis of flower color divergence in Barthea barthei, a shrub tree species exhibiting natural variation in flower color. We assembled a high-quality genome assembly for B. barthei with a contig N50 of 2.39 Mb and a scaffold N50 of 16.21 Mb. The assembly was annotated with 46,430 protein-coding genes and 1,560 non-coding RNAs. Genome synteny analysis revealed two recent tetraploidization events in B. barthei, estimated to have occurred at approximately 17 and 63 million years ago. These tetraploidization events resulted in massive duplicated gene content, with over 70% of genes retained in collinear blocks. Gene family members of the core regulators of the MBW complex were significantly expanded in B. barthei compared to Arabidopsis, suggesting that these duplications may have provided raw genetic material for the evolution of novel regulatory interactions and the diversification of anthocyanin pigmentation. Transcriptome profiling of B. barthei flowers revealed differential expression of 9 transcription factors related to anthocyanin biosynthesis between the two ecotypes. Six of these differentially expressed transcription factors were identified as high-confidence candidates for adaptive evolution based on positive selection signals. This study provides insights into the genetic basis of flower color divergence and the evolutionary mechanisms underlying ecological adaptation in plants.

Cytonuclear interplay in auto- and allopolyploids: a multifaceted perspective from the Festuca-Lolium complex

Restoring cytonuclear stoichiometry is necessary after whole-genome duplication (WGD) and interspecific/intergeneric hybridization in plants. We investigated this phenomenon in auto- and allopolyploids of the Festuca-Lolium complex providing insights into the mechanisms governing cytonuclear interactions in early polyploid and hybrid generations. Our study examined the main processes potentially involved in restoring the cytonuclear balance after WGD comparing diploids and new and well-established autopolyploids. We uncovered that both the number of chloroplasts and the number of chloroplast genome copies were significantly higher in the newly established autopolyploids and grew further in more established autopolyploids. The increase in the copy number of the chloroplast genome exceeded the rise in the number of chloroplasts and fully compensated for the doubling of the nuclear genome. In addition, changes in nuclear and organelle gene expression were insignificant. Allopolyploid Festuca × Lolium hybrids displayed potential structural conflicts in parental protein variants within the cytonuclear complexes. While biased maternal allele expression has been observed in numerous hybrids, our results suggest that its role in cytonuclear stabilization in the Festuca × Lolium hybrids is limited. This study provides insights into the restoration of the cytonuclear stoichiometry, yet it emphasizes the need for future research to explore post-transcriptional regulation and its impact on cytonuclear gene expression stoichiometry. Our findings may enhance the understanding of polyploid plant evolution, with broader implications for the study of cytonuclear interactions in diverse biological contexts.

Genome-wide characterization and expression profiling of the HD-ZIP gene family in Acoraceae under salinity and cold stress

The Homeodomain-Leucine Zipper (HD-ZIP) transcription factors play a pivotal role in governing various aspects of plant growth, development, and responses to abiotic stress. Despite the well-established importance of HD-ZIPs in many plants, their functions in Acoraceae, the basal lineage of monocots, remain largely unexplored. Using recently published whole-genome data, we identified 137 putative HD-ZIPs in two Acoraceae species, Acorus gramineus and Acorus calamus. These HD-ZIP genes were further classified into four subfamilies (I, II, III, IV) based on phylogenetic and conserved motif analyses, showcasing notable variations in exon-intron patterns among different subfamilies. Two microRNAs, miR165/166, were found to specifically target HD-ZIP III genes with highly conserved binding sites. Most cis-acting elements identified in the promoter regions of Acoraceae HD-ZIPs are involved in modulating light and phytohormone responsiveness. Furthermore, our study revealed an independent duplication event in Ac. calamus and a one-to-multiple correspondence between HD-ZIP genes of Ac. calamus and Ac. gramineus. Expression profiles obtained from qRT-PCR demonstrated that HD-ZIP I genes are strongly induced by salinity stress, while HD-ZIP II members have contrasting stress responses in two species. HD-ZIP III and IV genes show greater sensitivity in stress-bearing roots. Taken together, these findings contribute valuable insights into the roles of HD-ZIP genes in stress adaptation and plant resilience in basal monocots, illuminating their multifaceted roles in plant growth, development, and response to abiotic stress.

Evolutionary analysis of gene ages across TADs associates chromatin topology with whole-genome duplications

Topologically associated domains (TADs) are interaction subnetworks of chromosomal regions in 3D genomes. TAD boundaries frequently coincide with genome breaks while boundary deletion is under negative selection, suggesting that TADs may facilitate genome rearrangements and evolution. We show that genes co-localize by evolutionary age in humans and mice, resulting in TADs having different proportions of younger and older genes. We observe a major transition in the age co-localization patterns between the genes born during vertebrate whole-genome duplications (WGDs) or before and those born afterward. We also find that genes recently duplicated in primates and rodents are more frequently essential when they are located in old-enriched TADs and interact with genes that last duplicated during the WGD. Therefore, the evolutionary relevance of recent genes may increase when located in TADs with established regulatory networks. Our data suggest that TADs could play a role in organizing ancestral functions and evolutionary novelty.

Deciphering the differential expression patterns of yield-related negative regulators in hexaploid wheat cultivars and hybrids at different growth stages

BACKGROUND

Hexaploid bread wheat underwent a series of polyploidization events through interspecific hybridizations that conferred adaptive plasticity and resulted in duplication and neofunctionalization of major agronomic genes. The genetic architecture of polyploid wheat not only confers adaptive plasticity but also offers huge genetic diversity. However, the contribution of different gene copies (homeologs) encoded from different subgenomes (A, B, D) at different growth stages remained unexplored.

METHODS

In this study, hybrid of elite cultivars of wheat were developed via reciprocal crosses (cytoplasm swapping) and phenotypically evaluated. We assessed differential expression profiles of yield-related negative regulators in these cultivars and their F1 hybrids and identified various cis-regulatory signatures by employing bioinformatics tools. Furthermore, the preferential expression patterns of the syntenic triads encoded from A, B, and D subgenomes were assessed to decipher their functional redundancy at six different growth stages.

RESULTS

Hybrid progenies showed better heterosis such as up to 17% increase in the average number of grains and up to 50% increase in average thousand grains weight as compared to mid-parents. Based on the expression profiling, our results indicated significant dynamic transcriptional expression patterns, portraying the different homeolog-dominance at the same stage in the different cultivars and their hybrids. Albeit belonging to same syntenic triads, a dynamic trend was observed in the regulatory signatures of these genes that might be influencing their expression profiles.

CONCLUSION

These findings can substantially contribute and provide insights for the selective introduction of better cultivars into traditional and hybrid breeding programs which can be harnessed for the improvement of future wheat.

Whole-genome duplication in the Multicellularity Long Term Evolution Experiment

Whole-genome duplication (WGD) is widespread across eukaryotes and can promote adaptive evolution1-4. However, given the instability of newly-formed polyploid genomes5-7, understanding how WGDs arise in a population, persist, and underpin adaptations remains a challenge. Using our ongoing Multicellularity Long Term Evolution Experiment (MuLTEE)8, we show that diploid snowflake yeast (Saccharomyces cerevisiae) under selection for larger multicellular size rapidly undergo spontaneous WGD. From its origin within the first 50 days of the experiment, tetraploids persist for the next 950 days (nearly 5,000 generations, the current leading edge of our experiment) in ten replicate populations, despite being genomically unstable. Using synthetic reconstruction, biophysical modeling, and counter-selection experiments, we found that tetraploidy evolved because it confers immediate fitness benefits in this environment, by producing larger, longer cells that yield larger clusters. The same selective benefit also maintained tetraploidy over long evolutionary timescales, inhibiting the reversion to diploidy that is typically seen in laboratory evolution experiments. Once established, tetraploidy facilitated novel genetic routes for adaptation, playing a key role in the evolution of macroscopic multicellular size via the origin of evolutionarily conserved aneuploidy. These results provide unique empirical insights into the evolutionary dynamics and impacts of WGD, showing how it can initially arise due to its immediate adaptive benefits, be maintained by selection, and fuel long-term innovations by creating additional dimensions of heritable genetic variation.

The genetic consequences of range expansion and its influence on diploidization in polyploids

Despite newly formed polyploids being subjected to myriad fitness consequences, the relative prevalence of polyploidy both contemporarily and in ancestral branches of the tree of life suggests alternative advantages that outweigh these consequences. One proposed advantage is that polyploids may more easily colonize novel habitats such as deglaciated areas. However, previous research conducted in diploids suggests that range expansion comes with a fitness cost as deleterious mutations may fix rapidly on the expansion front. Here, we interrogate the potential consequences of expansion in polyploids by conducting spatially explicit forward-in-time simulations to investigate how ploidy and inheritance patterns impact the relative ability of polyploids to expand their range. We show that under realistic dominance models, autopolyploids suffer greater fitness reductions than diploids as a result of range expansion due to the fixation of increased mutational load that is masked in the range core. Alternatively, the disomic inheritance of allopolyploids provides a shield to this fixation resulting in minimal fitness consequences. In light of this advantage provided by disomy, we investigate how range expansion may influence cytogenetic diploidization through the reversion to disomy in autotetraploids. We show that under a wide range of parameters investigated for two models of diploidization, disomy frequently evolves more rapidly on the expansion front than in the range core, and that this dynamic inheritance model has additional effects on fitness. Together our results point to a complex interaction between dominance, ploidy, inheritance, and recombination on fitness as a population spreads across a geographic range.

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

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

Recurrent allopolyploidy and its implications for conservation in vascular plants: a commentary on 'Population genomics of the Isoetes appalachiana (Isoetaceae) complex supports a "diploids-first" approach to conservation'

This article comments on:

David Wickell, Jacob Landis, Elizabeth Zimmer and Fay-Wei Li, Population genomics of the Isoetes appalachiana (Isoetaceae) complex supports a ‘diploids-first’ approach to conservation, Annals of Botany, Volume 133, Issue 2, 01 February 2024, Pages 261–272, https://doi.org/10.1093/aob/mcad180