The evolutionary significance of polyploidy

Genetic Redundancy in Rye Shows in a Variety of Ways

Fifty years ago Susumu Ohno formulated the famous C-value paradox, which states that there is no correlation between the physical sizes of the genome, i.e., the amount of DNA, and the complexity of the organism, and highlighted the problem of genome redundancy. DNA that does not have a positive effect on the fitness of organisms has been characterized as "junk or selfish DNA". The controversial concept of junk DNA remains viable. Rye is a convenient subject for yet another test of the correctness and scientific significance of this concept. The genome of cultivated rye, Secale cereale L., is considered one of the largest among species of the tribe Triticeae and thus it tops the average angiosperm genome and the genomes of its closest evolutionary neighbors, such as species of barley, Hordeum (by approximately 30-35%), and diploid wheat species, Triticum (approximately 25%). The review provides an analysis of the structural organization of various regions of rye chromosomes with a description of the molecular mechanisms contributing to their size increase during evolution and the classes of DNA sequences involved in these processes. The history of the development of the concept of eukaryotic genome redundancy is traced and the current state of this problem is discussed.

Phenotypic variation of a new synthetic allotetraploid Arabidopsis kamchatica enhanced in natural environment

The phenotypic variation of vegetative organs and reproductive organs of newly synthesized and natural Arabidopsis kamchatica genotypes was investigated in both a controlled environment and a natural environment in an experimental garden. When we compared the variation of their leaf shape as a vegetative organ, the synthetic A. kamchatica individuals grown in the garden showed larger variation compared with the individuals incubated in a growth chamber, suggesting enhanced phenotypic variation in a natural fluctuating environment. In contrast, the natural A. kamchatica genotypes did not show significant change in variation by growth condition. The phenotypic variation of floral organs by growth condition was much smaller in both synthetic and natural A. kamchatica genotypes, and the difference in variation width between the growth chamber and the garden was not significant in each genotype as well as among genotypes. The higher phenotypic variation in synthetic leaf may imply flexible transcriptomic regulation of a newly synthesized polyploid compared with a natural polyploid.

Developing a CRISPR System in Nongenetic Model Polyploids

The genetic consequences following polyploidy (i.e., whole-genome duplication; WGD) vary greatly across organisms and through time since polyploidization. At the gene level in allopolyploids, changes include loss/retention of both parental gene copies, function/expression divergence between the two parental copies, and silencing of one parental copy. Functional studies of genes with different retention patterns contribute to a better understanding of the genetic factors underlying the success of polyploids. Most research on gene functions to date focuses on a few well-established genetic models or crops. However, many species that best exemplify the polyploidy process are nongenetic models; the lack of an efficient genome editing system hinders functional studies in these systems. In this chapter, we discuss the considerations of developing CRISPR, a robust and efficient genome editing system, in polyploid plants that are not genetic models. We use diploid and polyploid Tragopogon (Asteraceae) as examples of a well-studied evolutionary model system for which abundant genetic and genomic resources are lacking. Using this system, we provide our protocols for sgRNA design, plasmid construction, a useful protoplast transient assay, and a plant transformation method we developed for this system. We also provide suggestions for possible modifications to these protocols to help promote successful application to other non-models. With the rapid applications of CRISPR in plant sciences, the broad adaptation of CRISPR in studies of the evolutionary significance of WGD holds enormous potential. We hope our studies and methods developed for polyploid Tragopogon will provide a guideline for establishing a CRISPR system in other nongenetic model polyploids of evolutionary or other interest.

Sequencing and Assembly of Polyploid Genomes

Polyploidy has been observed throughout major eukaryotic clades and has played a vital role in the evolution of angiosperms. Recent polyploidizations often result in highly complex genome structures, posing challenges to genome assembly and phasing. Recent advances in sequencing technologies and genome assembly algorithms have enabled high-quality, near-complete chromosome-level assemblies of polyploid genomes. Advances in novel sequencing technologies include highly accurate single-molecule sequencing with HiFi reads, chromosome conformation capture with Hi-C technique, and linked reads sequencing. Additionally, new computational approaches have also significantly improved the precision and reliability of polyploid genome assembly and phasing, such as HiCanu, hifiasm, ALLHiC, and PolyGembler. Herein, we review recently published polyploid genomes and compare the various sequencing, assembly, and phasing approaches that are utilized in these genome studies. Finally, we anticipate that accurate and telomere-to-telomere chromosome-level assembly of polyploid genomes could ultimately become a routine procedure in the near future.

New insights into the genome and transmission of the microsporidian pathogen Nosema muscidifuracis

Introduction

Nosema is a diverse genus of unicellular microsporidian parasites of insects and other arthropods. Nosema muscidifuracis infects parasitoid wasp species of Muscidifurax zaraptor and M. raptor (Hymenoptera: Pteromalidae), causing ~50% reduction in longevity and ~90% reduction in fecundity.

Methods and Results

Here, we report the first assembly of the N. muscidifuracis genome (14,397,169 bp in 28 contigs) of high continuity (contig N50 544.3 Kb) and completeness (BUSCO score 97.0%). A total of 2,782 protein-coding genes were annotated, with 66.2% of the genes having two copies and 24.0% of genes having three copies. These duplicated genes are highly similar, with a sequence identity of 99.3%. The complex pattern suggests extensive gene duplications and rearrangements across the genome. We annotated 57 rDNA loci, which are highly GC-rich (37%) in a GC-poor genome (25% genome average). Nosema-specific qPCR primer sets were designed based on 18S rDNA annotation as a diagnostic tool to determine its titer in host samples. We discovered high Nosema titers in Nosema-cured M. raptor and M. zaraptor using heat treatment in 2017 and 2019, suggesting that the remedy did not completely eliminate the Nosema infection. Cytogenetic analyses revealed heavy infections of N. muscidifuracis within the ovaries of M. raptor and M. zaraptor, consistent with the titer determined by qPCR and suggesting a heritable component of infection and per ovum vertical transmission.

Discussion

The parasitoids-Nosema system is laboratory tractable and, therefore, can serve as a model to inform future genome manipulations of Nosema-host system for investigations of Nosemosis.

Where whole-genome duplication is most beneficial: Adaptation of mangroves to a wide salinity range between land and sea

Whole-genome duplication (WGD) is believed to increase the chance of adaptation to a new environment. This conjecture may apply particularly well to new environments that are not only different but also more variable than ancestral habitats. One such prominent environment is the interface between land and sea, which has been invaded by woody plants, collectively referred as mangroves, multiple times. Here, we use two distantly related mangrove species (Avicennia marina and Rhizophora apiculata) to explore the effects of WGD on the adaptive process. We found that a high proportion of duplicated genes retained after WGD have acquired derived differential expression in response to salt gradient treatment. The WGD duplicates differentially expressed in at least one copy usually (>90%) diverge from their paralogues' expression profiles. Furthermore, both species evolved in parallel to have one paralogue expressed at a high level in both fresh water and hypersaline conditions but at a lower level at medium salinity. The pattern contrasts with the conventional view of monotone increase/decrease as salinity increases. Differentially expressed copies have thus probably acquired a new role in salinity tolerance. Our results indicate that the WGD duplicates may have evolved to function collaboratively in coping with different salinity levels, rather than specializing in the intermediate salinity optimal for mangrove plants. In conclusion, WGD and the retained duplicates appear to be an effective solution for adaptation to new and unstable environments.

Analysis of 5S rDNA Genomic Organization Through the RepeatExplorer2 Pipeline: A Simplified Protocol

The ribosomal RNA genes (rDNA) are universal genome components with a housekeeping function, given the crucial role of ribosomal RNA in the synthesis of ribosomes and thus for life-on-Earth. Therefore, their genomic organization is of considerable interest for biologists, in general. Ribosomal RNA genes have also been largely used to establish phylogenetic relationships, and to identify allopolyploid or homoploid hybridization.Here, we demonstrate how high-throughput sequencing data, through graph clustering implemented in RepeatExplorer2 pipeline ( https://repeatexplorer-elixir.cerit-sc.cz/galaxy/ ), can be helpful to decipher the genomic organization of 5S rRNA genes. We show that the linear shapes of cluster graphs are reminiscent to the linked organization of 5S and 35S rDNA (L-type arrangement) while the circular graphs correspond to their separate arrangement (S-type). We further present a simplified protocol based on the paper by (Garcia et al., Front Plant Sci 11:41, 2020) about the use of graph clustering of 5S rDNA homoeologs (S-type) to identify hybridization events in the species history. We found that the graph complexity (i.e., graph circularity in this case) is related to ploidy and genome complexity, with diploids typically showing circular-shaped graphs while allopolyploids and other interspecific hybrids display more complex graphs, with usually two or more interconnected loops representing intergenic spacers. When a three-genomic comparative clustering analysis from a given hybrid (homoploid/allopolyploid) and its putative progenitor species (diploids) is performed, it is possible to identify the corresponding homoeologous 5S rRNA gene families, and to elucidate the contribution of each putative parental genome to the 5S rDNA pool of the hybrid. Thus, the analysis of 5S rDNA cluster graphs by RepeatExplorer, together with information coming from other sources (e.g., morphology, cytogenetics) is a complementary approach for the determination of allopolyploid or homoploid hybridization and even ancient introgression events.

Molecular evidence for adaptive evolution of drought tolerance in wild cereals

The considerable drought tolerance of wild cereal crop progenitors has diminished during domestication in the pursuit of higher productivity. Regaining this trait in cereal crops is essential for global food security but requires novel genetic insight. Here, we assessed the molecular evidence for natural variation of drought tolerance in wild barley (Hordeum spontaneum), wild emmer wheat (Triticum dicoccoides), and Brachypodium species collected from dry and moist habitats at Evolution Canyon, Israel (ECI). We report that prevailing moist vs dry conditions have differentially shaped the stomatal and photosynthetic traits of these wild cereals in their respective habitats. We present the genomic and transcriptomic evidence accounting for differences, including co-expression gene modules, correlated with physiological traits, and selective sweeps, driven by the xeric site conditions on the African Slope (AS) at ECI. Co-expression gene module 'circadian rhythm' was linked to significant drought-induced delay in flowering time in Brachypodium stacei genotypes. African Slope-specific differentially expressed genes are important in barley drought tolerance, verified by silencing Disease-Related Nonspecific Lipid Transfer 1 (DRN1), Nonphotochemical Quenching 4 (NPQ4), and Brassinosteroid-Responsive Ring-H1 (BRH1). Our results provide new genetic information for the breeding of resilient wheat and barley in a changing global climate with increasingly frequent drought events.

How plants conquered land: evolution of terrestrial adaptation

The transition of plants from water to land is considered one of the most significant events in the evolution of life on Earth. The colonization of land by plants, accompanied by their morphological, physiological and developmental changes, resulted in plant biodiversity. Besides significantly influencing oxygen levels in the air and on land, plants manufacture organic matter from CO₂ and water with the help of sunlight, paving the way for the diversification of nonplant lineages ranging from microscopic organisms to animals. Land plants regulate the climate by adjusting total biomass and energy flow. At the genetic level, these innovations are achieved through the rearrangement of pre-existing genetic information. Advances in genome sequencing technology are revamping our understanding of plant evolution. This study highlights the morphological and genomic innovations that allow plants to integrate life on Earth.

The genome of the king protea, Protea cynaroides

The king protea (Protea cynaroides), an early-diverging eudicot, is the most iconic species from the Megadiverse Cape Floristic Region, and the national flower of South Africa. Perhaps best known for its iconic flower head, Protea is a key genus for the South African horticulture industry and cut-flower market. Ecologically, the genus and the family Proteaceae are important models for radiation and adaptation, particularly to soils with limited phosphorus bio-availability. Here, we present a high-quality chromosome-scale assembly of the P. cynaroides genome as the first representative of the fynbos biome. We reveal an ancestral whole-genome duplication event that occurred in the Proteaceae around the late Cretaceous that preceded the divergence of all crown groups within the family and its extant diversity in all Southern continents. The relatively stable genome structure of P. cynaroides is invaluable for comparative studies and for unveiling paleopolyploidy in other groups, such as the distantly related sister group Ranunculales. Comparative genomics in sequenced genomes of the Proteales shows loss of key arbuscular mycorrhizal symbiosis genes likely ancestral to the family, and possibly the order. The P. cynaroides genome empowers new research in plant diversification, horticulture and adaptation, particularly to nutrient-poor soils.

POInT: Modeling Polyploidy in the Era of Ubiquitous Genomics

Thirteen years ago, we described an evolutionary modeling tool that could resolve the orthology relationships among the homologous genomic regions created by a whole-genome duplication. This tool, which we subsequently named POInT (the Polyploid Orthology Inference Tool), was originally only useful for studying a genome duplication known from bakers' yeast and its relatives. Now, with hundreds of genome sequences that contain the relicts of ancient polyploidy available, POInT can be used to study dozens of different polyploidies, asking both questions about the history of individual events and about the commonalities and differences seen between those events. In this chapter, I give a brief history of the development of POInT as an illustration of the interconnected nature of computational biology research. I then further describe how POInT operates and some of the strengths and drawbacks of its structure. I close with a few examples of discoveries we have made using it.

Ecological niches in the polyploid complex Linum suffruticosum s.l

Introduction

The high frequency of polyploidy in the evolutionary history of many plant groups occurring in the Mediterranean region is likely a consequence of its dynamic paleogeographic and climatic history. Polyploids frequently have distinct characteristics that allow them to overcome the minority cytotype exclusion. Such traits may enable polyploid individuals to grow in habitats different from their parentals and/or expand to new areas, leading to spatial segregation. Therefore, the successful establishment of polyploid lineages has long been associated with niche divergence or niche partitioning and the ability of polyploids to cope with different, often more stressful, conditions. In this study, we aimed to explore the role of environmental variables associated with the current distribution patterns of cytotypes within the polyploid complex Linum suffruticosum s.l..

Methods

The distribution and environmental niches of the five main cytotypes of Linum suffruticosum s.l. (diploids, tetraploids, hexaploids, octoploids and decaploids) were studied across its distribution range. Realized environmental niche of each cytotype was determined using niche modelling tools, such as maximum entropy modelling and niche equivalency and similarity tests.

Results

Differences in the environmental conditions of L. suffruticosum s.l. cytotypes were observed, with polyploids being associated with habitats of increased drought and soil pH, narrower temperature ranges and decreased soil water and cation exchange capacities. Diploids present the widest environmental niche, and polyploids occupy part of the diploid niche. Although some polyploids have equivalent potential ecological niches, cytotypes do not co-occur in nature. Additionally, the ecological niche of this polyploid complex is different between continents, with North African habitats being characterised by differences in soil texture, higher pH, and low cation exchange capacity, precipitation and soil water capacity and higher temperatures than habitats in southwest Europe.

Discussion

The different ecological conditions played a role in the distribution of cytotypes, but the mosaic distribution could not be entirely explained by the environmental variables included in this study. Other factors, such as reproductive isolation and competitive interactions among cytotypes, could further explain the current diversity and distribution patterns in white flax. This study provides relevant data on the niche requirements of each cytotype for further competition and reciprocal transplant experiments. further competition and reciprocal transplant experiments.

Phylotranscriptomics of liverworts: revisiting the backbone phylogeny and ancestral gene duplications

BACKGROUND AND AIMS

With some 7300 extant species, liverworts (Marchantiophyta) represent one of the major land plant lineages. The backbone relationships, such as the phylogenetic position of Ptilidiales, and the occurrence and timing of whole-genome duplications, are still contentious.

METHODS

Based on analyses of the newly generated transcriptome data for 38 liverworts and complemented with those publicly available, we reconstructed the evolutionary history of liverworts and inferred gene duplication events along the 55 taxon liverwort species tree.

KEY RESULTS

Our phylogenomic study provided an ordinal-level liverwort nuclear phylogeny and identified extensive gene tree conflicts and cyto-nuclear incongruences. Gene duplication analyses based on integrated phylogenomics and Ks distributions indicated no evidence of whole-genome duplication events along the backbone phylogeny of liverworts.

CONCLUSIONS

With a broadened sampling of liverwort transcriptomes, we re-evaluated the backbone phylogeny of liverworts, and provided evidence for ancient hybridizations followed by incomplete lineage sorting that shaped the deep evolutionary history of liverworts. The lack of whole-genome duplication during the deep evolution of liverworts indicates that liverworts might represent one of the few major embryophyte lineages whose evolution was not driven by whole-genome duplications.

syntenet: an R/Bioconductor package for the inference and analysis of synteny networks

SUMMARY

Interpreting and visualizing synteny relationships across several genomes is a challenging task. We previously proposed a network-based approach for better visualization and interpretation of large-scale microsynteny analyses. Here, we present syntenet, an R package to infer and analyze synteny networks from whole-genome protein sequence data. The package offers a simple and complete framework, including data preprocessing, synteny detection and network inference, network clustering and phylogenomic profiling, and microsynteny-based phylogeny inference. Graphical functions are also available to create publication-ready plots. Synteny networks inferred with syntenet can highlight taxon-specific gene clusters that likely contributed to the evolution of important traits, and microsynteny-based phylogenies can help resolve phylogenetic relationships under debate.

AVAILABILITY AND IMPLEMENTATION

syntenet is available on Bioconductor (https://bioconductor.org/packages/syntenet), and the source code is available on a GitHub repository (https://github.com/almeidasilvaf/syntenet).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Comparative Analysis of Physiological, Hormonal and Transcriptomic Responses Reveal Mechanisms of Saline-Alkali Tolerance in Autotetraploid Rice (Oryza sativa L.)

Saline-alkali soil has posed challenges to the growth of agricultural crops, while polyploidy often show greater adaptability in diverse and extreme environments including saline-alkali stress, but its defense mechanisms in rice remain elusive. Herein, we explored the mechanisms of enhanced saline-alkali tolerance of autotetraploid rice 93-11T relative to diploid rice 93-11D, based on physiological, hormonal and transcriptomic profilings. Physiologically, the enhanced saline-alkali tolerance in 93-11T was manifested in higher soluble sugar accumulation and stronger superoxide dismutase (SOD) and peroxidase (POD) activities in leaves during 24 h after saline-alkali shock. Furthermore, various hormone levels in leaves of 93-11T altered greatly, such as the negative correlation between salicylic acid (SA) and the other four hormones changed to positive correlation due to polyploidy. Global transcriptome profiling revealed that the upregulated differentially expressed genes (DEGs) in leaves and roots of 93-11T were more abundant than that in 93-11D, and there were more DEGs in roots than in leaves under saline-alkali stress. Genes related to phytohormone signal transduction of auxin (AUX) and SA in roots, lignin biosynthesis in leaves or roots, and wax biosynthesis in leaves were obviously upregulated in 93-11T compared with 93-11D under saline-alkali condition. Collectively, 93-11T subjected to saline-alkali stress possibly possesses higher osmotic regulation ability due to cuticular wax synthesis, stronger negative regulation of reactive oxygen species (ROS) production by increasing the SA levels and maintaining relative lower levels of IAA, and higher antioxidant capacity by increasing activities of SOD and POD, as well as lignin biosynthesis. Our research provides new insights for exploring the mechanisms of saline-alkali tolerance in polyploid rice and discovering new gene targets for rice genetic improvement.

A theory of rapid evolutionary change explaining the de novo appearance of megakaryocytes and platelets in mammals

Platelets are found only in mammals. Uniquely, they have a log Gaussian volume distribution and are produced from megakaryocytes, large cells that have polyploid nuclei. In this Hypothesis, we propose that a possible explanation for the origin of megakaryocytes and platelets is that, ∼220 million years ago, an inheritable change occurred in a mammalian ancestor that caused the haemostatic cell line of the animal to become polyploid. This inheritable change occurred specifically in the genetic programme of the cell lineage from which the haemostatic cell originated and led, because of increase in cell size, to its fragmentation into cytoplasmic particles (platelets) in the pulmonary circulatory system, as found in modern mammals. We hypothesize that these fragments originating from the new large haemostatic polyploid cells proved to be more efficient at stopping bleeding, and, therefore, the progeny of this ancestor prospered through natural selection. We also propose experimental strategies that could provide evidence to support this hypothesis.

The Current Developments in Medicinal Plant Genomics Enabled the Diversification of Secondary Metabolites' Biosynthesis

Medicinal plants produce important substrates for their adaptation and defenses against environmental factors and, at the same time, are used for traditional medicine and industrial additives. Plants have relatively little in the way of secondary metabolites via biosynthesis. Recently, the whole-genome sequencing of medicinal plants and the identification of secondary metabolite production were revolutionized by the rapid development and cheap cost of sequencing technology. Advances in functional genomics, such as transcriptomics, proteomics, and metabolomics, pave the way for discoveries in secondary metabolites and related key genes. The multi-omics approaches can offer tremendous insight into the variety, distribution, and development of biosynthetic gene clusters (BGCs). Although many reviews have reported on the plant and medicinal plant genome, chemistry, and pharmacology, there is no review giving a comprehensive report about the medicinal plant genome and multi-omics approaches to study the biosynthesis pathway of secondary metabolites. Here, we introduce the medicinal plant genome and the application of multi-omics tools for identifying genes related to the biosynthesis pathway of secondary metabolites. Moreover, we explore comparative genomics and polyploidy for gene family analysis in medicinal plants. This study promotes medicinal plant genomics, which contributes to the biosynthesis and screening of plant substrates and plant-based drugs and prompts the research efficiency of traditional medicine.

Homoeologous exchange enables rapid evolution of tolerance to salinity and hyper-osmotic stresses in a synthetic allotetraploid wheat

The link between polyploidy and enhanced adaptation to environmental stresses could be a result of polyploidy itself harbouring higher tolerance to adverse conditions, or polyploidy possessing higher evolvability than diploids under stress conditions. Natural polyploids are inherently unsuitable to disentangle these two possibilities. Using selfed progenies of a synthetic allotetraploid wheat AT3 (AADD) along with its diploid parents, Triticum urartu TMU38 (AA) and Aegilops tauschii TQ27 (DD), we addressed the foregoing issue under abiotic salinity and hyper-osmotic (drought-like) stress. Under short duration of both stresses, euploid plants of AT3 showed intermediate tolerance of diploid parents; under life-long duration of both stresses, tolerant individuals to either stress emerged from selfed progenies of AT3, but not from comparable-sized diploid parent populations. Tolerance to both stresses were conditioned by the same two homoeologous exchanges (HEs; 2DS/2AS and 3DL/3AL), and at least one HE needed to be at the homozygous state. Transcriptomic analyses revealed that hyper-up-regulation of within-HE stress responsive genes of the A sub-genome origin is likely responsible for the dual-stress tolerant phenotypes. Our results suggest that HE-mediated inter-sub-genome rearrangements can be an important mechanism leading to adaptive evolution in allopolyploids as well as a promising target for genetic manipulation in crop improvement.

WGDI: A user-friendly toolkit for evolutionary analyses of whole-genome duplications and ancestral karyotypes

Evidence of whole-genome duplications (WGDs) and subsequent karyotype changes has been detected in most major lineages of living organisms on Earth. To clarify the complex resulting multi-layered patterns of gene collinearity in genome analyses, there is a need for convenient and accurate toolkits. To meet this need, we developed WGDI (Whole-Genome Duplication Integrated analysis), a Python-based command-line tool that facilitates comprehensive analysis of recursive polyploidization events and cross-species genome alignments. WGDI supports three main workflows (polyploid inference, hierarchical inference of genomic homology, and ancestral chromosome karyotyping) that can improve the detection of WGD and characterization of WGD-related events based on high-quality chromosome-level genomes. Significantly, it can extract complete synteny blocks and facilitate reconstruction of detailed karyotype evolution. This toolkit is freely available at GitHub (https://github.com/SunPengChuan/wgdi). As an example of its application, WGDI convincingly clarified karyotype evolution in Aquilegia coerulea and Vitis vinifera following WGDs and rejected the hypothesis that Aquilegia contributed as a parental lineage to the allopolyploid origin of core dicots.

Genome-Wide Characterization and Expression Profiling of HD-Zip Genes in ABA-Mediated Processes in Fragaria vesca

Members of homeodomain-leucine zipper (HD-Zip) transcription factors can play their roles by modulating abscisic acid (ABA) signaling in Arabidopsis. So far, our knowledge of the functions of HD-Zips in woodland strawberries (Fragaria vesca), a model plant for studying ABA-mediated fruit ripening, is limited. Here, we identified a total of 31 HD-Zip genes (FveHDZ1-31) in F. vesca, and classified them into four subfamilies (I to IV). Promoter analyses show that the ABA-responsive element, ABRE, is prevalent in the promoters of subfamily I and II FveHDZs. RT-qPCR results demonstrate that 10 of the 14 investigated FveHDZs were consistently >1.5-fold up-regulated or down-regulated in expression in response to exogenous ABA, dehydration, and ABA-induced senescence in leaves. Five of the six consistently up-regulated genes are from subfamily I and II. Thereinto, FveHDZ4, and 20 also exhibited significantly enhanced expression along with increased ABA content during fruit ripening. In yeast one-hybrid assays, FveHDZ4 proteins could bind the promoter of an ABA signaling gene FvePP2C6. Collectively, our results strongly support that the FveHDZs, particularly those from subfamilies I and II, are involved in the ABA-mediated processes in F. vesca, providing a basis for further functional characterization of the HD-Zips in strawberry and other plants.

Inbreeding depression in polyploid species: a meta-analysis

Whole-genome duplication is a common mutation in eukaryotes with far-reaching phenotypic effects, the resulting morphological and fitness consequences and how they affect the survival of polyploid lineages are intensively studied. Another important factor may also determine the probability of establishment and success of polyploid lineages: inbreeding depression. Inbreeding depression is expected to play an important role in the establishment of neopolyploid lineages, their capacity to colonize new environments, and in the simultaneous evolution of ploidy and other life-history traits such as self-fertilization. Both theoretically and empirically, there is no consensus on the consequences of polyploidy on inbreeding depression. In this meta-analysis, we investigated the effect of polyploidy on the evolution of inbreeding depression, by performing a meta-analysis within angiosperm species. The main results of our study are that the consequences of polyploidy on inbreeding depression are complex and depend on the time since polyploidization. We found that young polyploid lineages have a much lower amount of inbreeding depression than their diploid relatives and their established counterparts. Natural polyploid lineages are intermediate and have a higher amount of inbreeding depression than synthetic neopolyploids, and a smaller amount than diploids, suggesting that the negative effect of polyploidy on inbreeding depression decreases with time since polyploidization.

Tree2GD: a phylogenomic method to detect large-scale gene duplication events

MOTIVATION

Whole-genome duplication events have long been discovered throughout the evolution of eukaryotes, contributing to genome complexity and biodiversity and leaving traces in the descending organisms. Therefore, an accurate and rapid phylogenomic method is needed to identify the retained duplicated genes on various lineages across the target taxonomy.

RESULTS

Here, we present Tree2GD, an integrated method to identify large-scale gene duplication events by automatically perform multiple procedures, including sequence alignment, recognition of homolog, gene tree/species tree reconciliation, Ks distribution of gene duplicates and synteny analyses. Application of Tree2GD on 2 datasets, 12 metazoan genomes and 68 angiosperms, successfully identifies all reported whole-genome duplication events exhibited by these species, showing effectiveness and efficiency of Tree2GD on phylogenomic analyses of large-scale gene duplications.

AVAILABILITY AND IMPLEMENTATION

Tree2GD is written in Python and C++ and is available at https://github.com/Dee-chen/Tree2gd.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A common whole-genome paleotetraploidization in Cucurbitales

Cucurbitales are an important order of flowering plants known for encompassing edible plants of economic and medicinal value and numerous ornamental plants of horticultural value. By reanalyzing the genomes of two representative families (Cucurbitaceae and Begoniaceae) in Cucurbitales, we found that the previously identified Cucurbitaceae common paleotetraploidization that occurred shortly after the core-eudicot-common hexaploidization event is shared by Cucurbitales, including Begoniaceae. We built a multigenome alignment framework for Cucurbitales by identifying orthologs and paralogs and systematically redating key evolutionary events in Cucurbitales. Notably, characterizing the gene retention levels and genomic fractionation patterns between subgenomes generated from different polyploidizations in Cucurbitales suggested the autopolyploid nature of the Begoniaceae common tetraploidization and the allopolyploid nature of the Cucurbitales common tetraploidization and the Cucurbita-specific tetraploidization. Moreover, we constructed the ancestral Cucurbitales karyotype comprising 17 proto-chromosomes, confirming that the most recent common ancestor of Cucurbitaceae contained 15 proto-chromosomes and rejecting the previous hypothesis for an ancestral Cucurbitaceae karyotype with 12 proto-chromosomes. In addition, we found that the polyploidization and tandem duplication events promoted the expansion of gene families involved in the cucurbitacin biosynthesis pathway; however, gene loss and chromosomal rearrangements likely limited the expansion of these gene families.

Genome-wide identification and analysis of LEA_2 gene family in alfalfa (Medicago sativa L.) under aluminum stress

Late embryonic development abundant proteins (LEAs) are a large family of proteins commonly existing in plants. LEA_2 is the largest subfamily in the LEA, it plays an important role in plant resistance to abiotic stress. In order to explore the characteristics of LEA_2 gene family members in alfalfa (Medicago sativa L.), 155 members of LEA_2 (MsLEA_2) family were identified from alfalfa genome. Bioinformatics analysis was conducted from the aspects of phylogenetic relationship, chromosome distribution, chromosome colinearity, physical and chemical properties, motif composition, exon-intron structure, cis-element and so on. Expression profiles of MsLEA_2 gene were obtained based on Real-time fluorescent quantitative PCR (qRT-PCR) analysis and previous RNA-seq data under aluminum (Al) stress. Bioinformatics results were shown that the MsLEA_2 genes are distributed on all 32 chromosomes. Among them, 85 genes were present in the gene clusters, accounting for 54.83%, and chromosome Chr7.3 carries the largest number of MsLEA_2 (19 LEA_2 genes on Chr7.3). Chr7.3 has a unique structure of MsLEA_2 distribution, which reveals a possible special role of Chr7.3 in ensuring the function of MsLEA_2. Transcriptional structure analysis revealed that the number of exons in each gene varies from 1 to 3, and introns varies from 0 to 2. Cis-element analysis identified that the promoter region of MsLEA_2 is rich in ABRE, MBS, LTR, and MeJARE, indicating MsLEA_2 has stress resistance potential under abiotic stress. RNA-seq data and qRT-PCR analyses showed that most of the MsLEA_2 members were up-regulated when alfalfa exposed to Al stress. This study revealed that phylogenetic relationship and possible function of LEA_ 2 gene in alfalfa, which were helpful for the functional analysis of LEA_ 2 proteins in the future and provided a new theoretical basis for improving Al tolerance of alfalfa.

Salmonidae Genome: Features, Evolutionary and Phylogenetic Characteristics

The salmon family is one of the most iconic and economically important fish families, primarily possessing meat of excellent taste as well as irreplaceable nutritional and biological value. One of the most common and, therefore, highly significant members of this family, the Atlantic salmon (Salmo salar L.), was not without reason one of the first fish species for which a high-quality reference genome assembly was produced and published. Genomic advancements are becoming increasingly essential in both the genetic enhancement of farmed salmon and the conservation of wild salmon stocks. The salmon genome has also played a significant role in influencing our comprehension of the evolutionary and functional ramifications of the ancestral whole-genome duplication event shared by all Salmonidae species. Here we provide an overview of the current state of research on the genomics and phylogeny of the various most studied subfamilies, genera, and individual salmonid species, focusing on those studies that aim to advance our understanding of salmonid ecology, physiology, and evolution, particularly for the purpose of improving aquaculture production. This review should make potential researchers pay attention to the current state of research on the salmonid genome, which should potentially attract interest in this important problem, and hence the application of new technologies (such as genome editing) in uncovering the genetic and evolutionary features of salmoniforms that underlie functional variation in traits of commercial and scientific importance.

A draft genome of the medicinal plant Cremastra appendiculata (D. Don) provides insights into the colchicine biosynthetic pathway

Cremastra appendiculata (D. Don) Makino is a rare terrestrial orchid with a high market value as an ornamental and Chinese traditional medicinal herb with a wide range of pharmacological properties. The pseudobulbs of C. appendiculata are one of the primary sources of the famous traditional Chinese medicine "Shancigu", which has been clinically used for treating many diseases, especially, as the main component to treat gout. The lack of genetic research and genome data restricts the modern development and clinical use of C. appendiculata. Here, we report a 2.3 Gb chromosome-level genome of C. appendiculata. We identify a series of candidates of 35 candidate genes responsible for colchicine biosynthesis, among which O-methyltransferase (OMT) gene exhibits an important role in colchicine biosynthesis. Co-expression analysis reveal purple and green-yellow module have close relationships with pseudobulb parts and comprise most of the colchicine pathway genes. Overall, our genome data and the candidate genes reported here set the foundation to decipher the colchicine biosynthesis pathways in medicinal plants.

Low-affinity SPL binding sites contribute to subgenome expression divergence in allohexaploid wheat

Expression divergence caused by genetic variation and crosstalks among subgenomes of the allohexaploid bread wheat (Triticum aestivum. L., BBAADD) is hypothesized to increase its adaptability and/or plasticity. However, the molecular basis of expression divergence remains unclear. Squamosa promoter-binding protein-like (SPL) transcription factors are critical for a wide array of biological processes. In this study, we constructed expression regulatory networks by combining DAP-seq for 40 SPLs, ATAC-seq, and RNA-seq. Our findings indicate that a group of low-affinity SPL binding regions (SBRs) were targeted by diverse SPLs and caused different sequence preferences around the core GTAC motif. The SBRs including the low-affinity ones are evolutionarily conserved, enriched GWAS signals related to important agricultural traits. However, those SBRs are highly diversified among the cis-regulatory regions (CREs) of syntenic genes, with less than 8% SBRs coexisting in triad genes, suggesting that CRE variations are critical for subgenome differentiations. Knocking out of TaSPL7A/B/D and TaSPL15A/B/D subfamily further proved that both high- and low-affinity SBRs played critical roles in the differential expression of genes regulating tiller number and spike sizes. Our results have provided baseline data for downstream networks of SPLs and wheat improvements and revealed that CRE variations are critical sources for subgenome divergence in the allohexaploid wheat.

Frequent spontaneous structural rearrangements promote rapid genome diversification in a Brassica napus F1 generation

In a cross between two homozygous Brassica napus plants of synthetic and natural origin, we demonstrate that novel structural genome variants from the synthetic parent cause immediate genome diversification among F1 offspring. Long read sequencing in twelve F1 sister plants revealed five large-scale structural rearrangements where both parents carried different homozygous alleles but the heterozygous F1 genomes were not identical heterozygotes as expected. Such spontaneous rearrangements were part of homoeologous exchanges or segmental deletions and were identified in different, individual F1 plants. The variants caused deletions, gene copy-number variations, diverging methylation patterns and other structural changes in large numbers of genes and may have been causal for unexpected phenotypic variation between individual F1 sister plants, for example strong divergence of plant height and leaf area. This example supports the hypothesis that spontaneous de novo structural rearrangements after de novo polyploidization can rapidly overcome intense allopolyploidization bottlenecks to re-expand crops genetic diversity for ecogeographical expansion and human selection. The findings imply that natural genome restructuring in allopolyploid plants from interspecific hybridization, a common approach in plant breeding, can have a considerably more drastic impact on genetic diversity in agricultural ecosystems than extremely precise, biotechnological genome modifications.

Tetraploidy accelerates adaptation under drug selection in a fungal pathogen

Baseline ploidy significantly impacts evolutionary trajectories and, specifically, tetraploidy is associated with higher rates of adaptation relative to haploidy and diploidy. While the majority of experimental evolution studies investigating ploidy use the budding yeast Saccharomyces cerivisiae, the fungal pathogen Candida albicans is a powerful system to investigate ploidy dynamics, particularly in the context of acquiring antifungal drug resistance. C. albicans laboratory and clinical strains are predominantly diploid, but have been isolated as haploid and polyploid. Here, we evolved diploid and tetraploid C. albicans for ~60 days in the antifungal drug caspofungin. Tetraploid-evolved lines adapted faster than diploid-evolved lines and reached higher levels of caspofungin resistance. While diploid-evolved lines generally maintained their initial genome size, tetraploid-evolved lines rapidly underwent genome-size reductions and did so prior to caspofungin adaptation. While clinical resistance was largely due to mutations in FKS1, these mutations were caused by substitutions in diploid, and indels in tetraploid isolates. Furthermore, fitness costs in the absence of drug selection were significantly less in tetraploid-evolved lines compared to the diploid-evolved lines. Taken together, this work supports a model of adaptation in which the tetraploid state is transient but its ability to rapidly transition ploidy states improves adaptive outcomes and may drive drug resistance in fungal pathogens.

Convergent evolution of AP2/ERF III and IX subfamilies through recurrent polyploidization and tandem duplication during eudicot adaptation to paleoenvironmental changes

Whole-genome duplication (WGD or polyploidization) has been suggested as a genetic contributor to angiosperm adaptation to environmental changes. However, many eudicot lineages did not undergo recent WGD (R-WGD) around and/or after the Cretaceous-Paleogene (K-Pg) boundary, times of severe environmental changes; how those plants survived has been largely ignored. Here, we collected 22 plants from major branches of the eudicot phylogeny and classified them into two groups according to the occurrence or absence of R-WGD: 12 R-WGD-containing plants (R-WGD-Y) and 10 R-WGD-lacking plants (R-WGD-N). Subsequently, we identified 496 gene-rich families in R-WGD-Y and revealed that members of the AP2/ERF transcription factor family were convergently over-retained after R-WGDs and showed exceptional cold stimulation. The evolutionary trajectories of the AP2/ERF family were then compared between R-WGD-Y and R-WGD-N to reveal convergent expansions of the AP2/ERF III and IX subfamilies through recurrent independent WGDs and tandem duplications (TDs) after the radiation of the plants. The expansions showed coincident enrichments in- times around and/or after the K-Pg boundary, when global cooling was a major environmental stressor. Consequently, convergent expansions and co-retentions of AP2/ERF III C-repeat binding factor (CBF) duplicates and their regulons in different eudicot lineages contributed to the rewiring of cold-specific regulatory networks. Moreover, promoter analysis of cold-responsive AP2/ERF genes revealed an underlying cis-regulatory code (G-box: CACGTG). We propose a seesaw model of WGDs and TDs in the convergent expansion of AP2/ERF III and IX genes that has contributed to eudicot adaptation during paleoenvironmental changes, and we suggest that TD may be a reciprocal/alternative mechanism for genetic innovation in plants that lack WGD.

Transposable elements orchestrate subgenome-convergent and -divergent transcription in common wheat

The success of common wheat as a global staple crop was largely attributed to its genomic diversity and redundancy due to the merge of different genomes, giving rise to the major question how subgenome-divergent and -convergent transcription is mediated and harmonized in a single cell. Here, we create a catalog of genome-wide transcription factor-binding sites (TFBSs) to assemble a common wheat regulatory network on an unprecedented scale. A significant proportion of subgenome-divergent TFBSs are derived from differential expansions of particular transposable elements (TEs) in diploid progenitors, which contribute to subgenome-divergent transcription. Whereas subgenome-convergent transcription is associated with balanced TF binding at loci derived from TE expansions before diploid divergence. These TFBSs have retained in parallel during evolution of each diploid, despite extensive unbalanced turnover of the flanking TEs. Thus, the differential evolutionary selection of paleo- and neo-TEs contribute to subgenome-convergent and -divergent regulation in common wheat, highlighting the influence of TE repertory plasticity on transcriptional plasticity in polyploid.

Evolution of isoform-level gene expression patterns across tissues during lotus species divergence

Both gene duplication and alternative splicing (AS) drive the functional diversity of gene products in plants, yet the relative contributions of the two key mechanisms to the evolution of gene function are largely unclear. Here, we studied AS in two closely related lotus plants, Nelumbo lutea and Nelumbo nucifera, and the outgroup Arabidopsis thaliana, for both single-copy and duplicated genes. We show that most splicing events evolved rapidly between orthologs and that the origin of lineage-specific splice variants or isoforms contributed to gene functional changes during species divergence within Nelumbo. Single-copy genes contain more isoforms, have more AS events conserved across species, and show more complex tissue-dependent expression patterns than their duplicated counterparts. This suggests that expression divergence through isoforms is a mechanism to extend the expression breadth of genes with low copy numbers. As compared to isoforms of local, small-scale duplicates, isoforms of whole-genome duplicates are less conserved and display a less conserved tissue bias, pointing towards their contribution to subfunctionalization. Through comparative analysis of isoform expression networks, we identified orthologous genes of which the expression of at least some of their isoforms displays a conserved tissue bias across species, indicating a strong selection pressure for maintaining a stable expression pattern of these isoforms. Overall, our study shows that both AS and gene duplication contributed to the diversity of gene function during the evolution of lotus.

Genome-wide chromatin accessibility analysis unveils open chromatin convergent evolution during polyploidization in cotton

Allopolyploidization, resulting in divergent genomes in the same cell, is believed to trigger a “genome shock”, leading to broad genetic and epigenetic changes. However, little is understood about chromatin and gene-expression dynamics as underlying driving forces during allopolyploidization. Here, we examined the genome-wide DNase I-hypersensitive site (DHS) and its variations in domesticated allotetraploid cotton ( Gossypium hirsutum and Gossypium barbadense , AADD) and its extant AA ( Gossypium arboreum ) and DD ( Gossypium raimondii ) progenitors. We observed distinct DHS distributions between G. arboreum and G. raimondii . In contrast, the DHSs of the two subgenomes of G. hirsutum and G. barbadense showed a convergent distribution. This convergent distribution of DHS was also present in the wild allotetraploids Gossypium darwinii and G. hirsutum var. yucatanense , but absent from a resynthesized hybrid of G. arboreum and G. raimondii , suggesting that it may be a common feature in polyploids, and not a consequence of domestication after polyploidization. We revealed that putative cis -regulatory elements (CREs) derived from polyploidization-related DHSs were dominated by several families, including Dof, ERF48, and BPC1. Strikingly, 56.6% of polyploidization-related DHSs were derived from transposable elements (TEs). Moreover, we observed positive correlations between DHS accessibility and the histone marks H3K4me3, H3K27me3, H3K36me3, H3K27ac, and H3K9ac, indicating that coordinated interplay among histone modifications, TEs, and CREs drives the DHS landscape dynamics under polyploidization. Collectively, these findings advance our understanding of the regulatory architecture in plants and underscore the complexity of regulome evolution during polyploidization.

Temporal stability of spatial cytotype structure in mixed-ploidy populations of Centaurea stoebe

Spatial segregation of cytotypes reduces the negative effect of frequency-dependent mating on the fitness of minority cytotype(s) and thus allows its establishment and coexistence with the majority cytotype in mixed-ploidy populations. Despite its evolutionary importance, the stability of spatial segregation is largely unknown. Furthermore, closely related sympatric cytotypes that differ in their life histories might exhibit contrasting spatial dynamics over time. We studied the temporal stability of spatial structure at a secondary contact zone of co-occurring monocarpic diploids and polycarpic tetraploids of Centaurea stoebe, whose tetraploid cytotype has undergone a rapid range expansion in Europe and became invasive in North America. Eleven years after the initial screening, we re-assessed the microspatial distribution of diploids and tetraploids and their affinities to varying vegetation-cover density in three mixed-ploidy populations in Central Europe. We found that overall, spatial patterns and frequencies of both cytotypes in all sites were very similar over time, with one exception. At one site, in one previously purely 2x patch, diploids completely disappeared due to intensive succession by shrubby vegetation. The remaining spatial patterns, however, showed the same cytotype clumping and higher frequency of 2x despite subtle changes in vegetation-cover densities. In contrast to the expected expansion of polycarpic tetraploids having higher colonization ability when compared to diploids, the tetraploids remained confined to their former microsites and showed no spatial expansion. Spatial patterns of coexisting diploids and tetraploids, which exhibit contrasting life histories, did not change over more than a decade. Such temporal stability is likely caused by relatively stable habitat conditions and very limited seed dispersal. Our results thus imply that in the absence of a disturbance regime connected with frequent human- or animal-mediated seed dispersal, spatial patterns may be very stable over time, thus contributing to the long-term coexistence of cytotypes.

Genome-wide identification, evolutionary analysis, and antimicrobial activity prediction of CC chemokines in allotetraploid common carp, Cyprinus carpio

Chemokines are a group of secreted small molecules which are essential for cell migration in physiological and pathological conditions by binding to specific chemokine receptors. They are structurally classified into five groups, namely CXC, CC, CX3C, XC and CX. CC chemokine group is the largest one among them. In this study, we identified and characterized 61 CC chemokines from allotetraploid common carp (Cyprinus carpio). The sequence analyses showed that the majority of CC chemokines had an N-terminal signal peptide, and an SCY domain, and all CC chemokines were located in the extracellular region. Phylogenetic, evolutionary and syntenic analyses confirmed that CC chemokines were annotated as 11 different types (CCL19, CCL20, CCL25, CCL27, CCL32, CCL33, CCL34, CCL35, CCL36, CCL39, and CCL44), which exhibited unique gene arrangement pattern and chromosomal location respectively. Furthermore, genome synteny analyses between common carp and four representative teleost species indicated expansion of common carp CC chemokines resulted from the whole genome duplication (WGD) event. Additionally, the continuous evolution of gene CCL25s in teleost afforded a novel viewpoint to explain the WGD event in teleost. Then, we predicted the three-dimensional structures and probable function regions of common carp CC chemokines. All the CC chemokines core structures were constituted of an N-loop, a three-stranded β-sheet, and a C-terminal helix. Finally, 43 CC chemokines were predicted to have probable general antimicrobial activity. Their tertiary structures, cationic and amphiphilic physicochemical property supported the viewpoint. To verify the prediction, six recombinant CCL19s proteins were prepared and the antibacterial activity against Escherichia coli and Aeromonas hydrophila were verified. The results supported our prediction that rCCL19a.1s (rCCL19a.1_a, rCCL19a.1_b) and rCCL19bs (rCCL19b_a, rCCL19b_b), especially rCCL19bs, exhibited extremely significant inhibition to the growth of both E. coli and A. hydrophila. On the contrary, two rCCL19a.2s had no significant inhibitory effect. These studies suggested that CC chemokines were essential in immune system evolution and not monofunctional during pathogen infection.

Innovations and stepwise evolution of CBFs/DREB1s and their regulatory networks in angiosperms

The C-repeat binding factors/dehydration-responsive element binding protein 1s (CBFs/DREB1s) have been identified as major regulators of cold acclimation in many angiosperm plants. However, their origin and evolutionary process associated to cold responsiveness are still lacking. By integrating multi-omics data of genomes, transcriptomes, and CBFs/DREB1s genome-wide binding profiles, we unveil the origin and evolution of CBFs/DREB1s and their regulatory network. Gene collinearity and phylogeny analyses show that CBF/DREB1 is an innovation evolved from tandem duplication-derived DREB III gene. A subsequent event of ε-whole genome duplication led to two CBF/DREB1 archetypes (Clades I and II) in ancient angiosperms. In contrast to cold-insensitivity of Clade I and their parent DREB III genes, Clade II evolved a further innovation in cold-sensitive response and was stepwise expanded in eudicots and monocots by independent duplications. In geological time, the duplication events were mainly enriched around the Cretaceous-Paleogene (K-Pg) boundary and/or in the Late Cenozoic Ice Age, when the global average temperature significantly decreased. Consequently, the duplicated CBF/DREB1 genes contributed to the rewiring of CBFs/DREB1s-regulatory network for cold tolerance. Altogether, our results highlight an origin and convergent evolution of CBFs/DREB1s and their regulatory network probably for angiosperms adaptation to global cooling.

The genomes of chicory, endive, great burdock and yacon provide insights into Asteraceae palaeo-polyploidization history and plant inulin production

Inulin is an important reserve polysaccharide in Asteraceae plants, and is also widely used as a sweetener, a source of dietary fibre and prebiotic. Nevertheless, a lack of genomic resources for inulin-producing plants has hindered extensive studies on inulin metabolism and regulation. Here, we present chromosome-level reference genomes for four inulin-producing plants: chicory (Cichorium intybus), endive (Cichorium endivia), great burdock (Arctium lappa) and yacon (Smallanthus sonchifolius), with assembled genome sizes of 1.28, 0.89, 1.73 and 2.72 Gb, respectively. We found that the chicory, endive and great burdock genomes were shaped by whole genome triplication (WGT-1), and the yacon genome was shaped by WGT-1 and two subsequent whole genome duplications (WGD-2 and WGD-3). A yacon unique whole genome duplication (WGD-3) occurred 5.6-5.8 million years ago. Our results also showed the genome size difference between chicory and endive is largely due to LTR retrotransposons, and rejected a previous hypothesis that chicory is an ancestor of endive. Furthermore, we identified fructan-active-enzyme and transcription-factor genes, and found there is one copy in chicory, endive and great burdock but two copies in yacon for most of these genes, except for the 1-FEH II gene which is significantly expanded in chicory. Interestingly, inulin synthesis genes 1-SST and 1-FFT are located close to each other, as are the degradation genes 1-FEH I and 1-FEH II. Finally, we predicted protein structures for 1-FFT genes to explore the mechanism determining inulin chain length.

Accurate Chromosome Identification in the Prunus Subgenus Cerasus (Prunus pseudocerasus) and its Relatives by Oligo-FISH

A precise, rapid and straightforward approach to chromosome identification is fundamental for cytogenetics studies. However, the identification of individual chromosomes was not previously possible for Chinese cherry or other Prunus species due to the small size and similar morphology of their chromosomes. To address this issue, we designed a pool of oligonucleotides distributed across specific pseudochromosome regions of Chinese cherry. This oligonucleotide pool was amplified through multiplex PCR with specific internal primers to produce probes that could recognize specific chromosomes. External primers modified with red and green fluorescence tags could produce unique signal barcoding patterns to identify each chromosome concomitantly. The same oligonucleotide pool could also discriminate all chromosomes in other Prunus species. Additionally, the 5S/45S rDNA probes and the oligo pool were applied in two sequential rounds of fluorescence in situ hybridization (FISH) localized to chromosomes and showed different distribution patterns among Prunus species. At the same time, comparative karyotype analysis revealed high conservation among P. pseudocerasus, P. avium, and P. persica. Together, these findings establish this oligonucleotide pool as the most effective tool for chromosome identification and the analysis of genome organization and evolution in the genus Prunus.

A gain-of-function allele of a DREB transcription factor gene ameliorates drought tolerance in wheat

Drought is a major environmental factor limiting wheat production worldwide. However, the genetic components underlying wheat drought tolerance are largely unknown. Here, we identify a DREB transcription factor gene (TaDTG6-B) by genome-wide association study that is tightly associated with drought tolerance in wheat. Candidate gene association analysis revealed that a 26-bp deletion in the TaDTG6-B coding region induces a gain-of-function for TaDTG6-BDel574, which exhibits stronger transcriptional activation, protein interactions, and binding activity to dehydration-responsive elements (DRE)/CRT cis-elements than the TaDTG6-BIn574 encoded by the allele lacking the deletion, thus conferring greater drought tolerance in wheat seedlings harboring this variant. Knockdown of TaDTG6-BDel574 transcripts attenuated drought tolerance in transgenic wheat, whereas its overexpression resulted in enhanced drought tolerance without accompanying phenotypic abnormalities. Furthermore, the introgression of the TaDTG6-BDel574 elite allele into drought-sensitive cultivars improved their drought tolerance, thus providing a valuable genetic resource for wheat breeding. We also identified 268 putative target genes that are directly bound and transcriptionally regulated by TaDTG6-BDel574. Further analysis showed that TaDTG6-BDel574 positively regulates TaPIF1 transcription to enhance wheat drought tolerance. These results describe the genetic basis and accompanying mechanism driving phenotypic variation in wheat drought tolerance, and provide a novel genetic resource for crop breeding programs.

Glutamate 1-semialdehyde aminotransferase is connected to GluTR by GluTR-binding protein and contributes to the rate-limiting step of 5-aminolevulinic acid synthesis

Tetrapyrroles play fundamental roles in crucial processes including photosynthesis, respiration, and catalysis. In plants, 5-aminolevulinic acid (ALA) is the common precursor of tetrapyrroles. ALA is synthesized from activated glutamate by the enzymes glutamyl-tRNA reductase (GluTR) and glutamate-1-semialdehyde aminotransferase (GSAAT). ALA synthesis is recognized as the rate-limiting step in this pathway. We aimed to explore the contribution of GSAAT to the control of ALA synthesis and the formation of a protein complex with GluTR. In Arabidopsis thaliana, two genes encode GSAAT isoforms: GSA1 and GSA2. A comparison of two GSA knockout mutants with the wild-type revealed the correlation of reduced GSAAT activity and ALA-synthesizing capacity in leaves with lower chlorophyll content. Growth and green pigmentation were more severely impaired in gsa2 than in gsa1, indicating the predominant role of GSAAT2 in ALA synthesis. Interestingly, GluTR accumulated to higher levels in gsa2 than in the wild-type and was mainly associated with the plastid membrane. We propose that the GSAAT content modulates the amount of soluble GluTR available for ALA synthesis. Several different biochemical approaches revealed the GSAAT-GluTR interaction through the assistance of GluTR-binding protein (GBP). A modeled structure of the tripartite protein complex indicated that GBP mediates the stable association of GluTR and GSAAT for adequate ALA synthesis.

Transcriptional network underpinning ploidy-related elevated leaf potassium in neo-tetraploids

Whole-genome duplication generates a tetraploid from a diploid. Newly created tetraploids (neo-tetraploids) of Arabidopsis (Arabidopsis thaliana) have elevated leaf potassium (K), compared to their diploid progenitor. Micro-grafting has previously established that this elevated leaf K is driven by processes within the root. Here, mutational analysis revealed that the K+-uptake transporters K+ TRANSPORTER 1 (AKT1) and HIGH AFFINITY K+ TRANSPORTER 5 (HAK5) are not necessary for the difference in leaf K caused by whole-genome duplication. However, the endodermis and salt overly sensitive and abscisic acid-related signaling were necessary for the elevated leaf K in neo-tetraploids. Contrasting the root transcriptomes of neo-tetraploid and diploid wild-type and mutants that suppress the neo-tetraploid elevated leaf K phenotype allowed us to identify a core set of 92 differentially expressed genes associated with the difference in leaf K between neo-tetraploids and their diploid progenitor. This core set of genes connected whole-genome duplication with the difference in leaf K between neo-tetraploids and their diploid progenitors. The set of genes is enriched in functions such as cell wall and Casparian strip development and ion transport in the endodermis, root hairs, and procambium. This gene set provides tools to test the intriguing idea of recreating the physiological effects of whole-genome duplication within a diploid genome.

Hybridization order is not the driving factor behind biases in duplicate gene losses among the hexaploid Solanaceae

We model the post-hexaploidy evolution of four genomes from the Solanaceae, a group of flowering plants comprising tomatoes, potatoes and their relatives. The hexaploidy that these genomes descend from occurred through two sequential allopolyploidy events and was marked by the unequal losses of duplicated genes from the different progenitor subgenomes. In contrast with the hexaploid Brassiceae (broccoli and its relatives), where the subgenome with the most surviving genes arrived last in the hexaploidy, among the Solanaceae the most preserved subgenome descends from one of the original two tetraploid progenitors. In fact, the last-arriving subgenome in these plants actually has the fewest surviving genes in the modern genomes. We explore whether the distribution of repetitive elements (REs) in these genomes can explain the biases in gene losses, but while the signals we find are broadly consistent with a role for high RE density in driving gene losses, the REs turn over so quickly that little signal of the RE condition at the time of paleopolyploidy is extant in the modern genomes.

Sublinear scaling of the cellular proteome with ploidy

Ploidy changes are frequent in nature and contribute to evolution, functional specialization and tumorigenesis. Analysis of model organisms of different ploidies revealed that increased ploidy leads to an increase in cell and nuclear volume, reduced proliferation, metabolic changes, lower fitness, and increased genomic instability, but the underlying mechanisms remain poorly understood. To investigate how gene expression changes with cellular ploidy, we analyzed isogenic series of budding yeasts from 1N to 4N. We show that mRNA and protein abundance scales allometrically with ploidy, with tetraploid cells showing only threefold increase in protein abundance compared to haploids. This ploidy-dependent sublinear scaling occurs via decreased rRNA and ribosomal protein abundance and reduced translation. We demonstrate that the activity of Tor1 is reduced with increasing ploidy, which leads to diminished rRNA gene repression via a Tor1-Sch9-Tup1 signaling pathway. mTORC1 and S6K activity are also reduced in human tetraploid cells and the concomitant increase of the Tup1 homolog Tle1 downregulates the rDNA transcription. Our results suggest that the mTORC1-Sch9/S6K-Tup1/TLE1 pathway ensures proteome remodeling in response to increased ploidy.

Alteration of synonymous codon usage bias accompanies polyploidization in wheat

The diploidization of polyploid genomes is accompanied by genomic variation, including synonymous nucleotide substitutions that may lead to synonymous codon usage bias (SCUB). SCUB can mirror the evolutionary specialization of plants, but its effect on the formation of polyploidies is not well documented. We explored this issue here with hexaploid wheat and its progenitors. Synonymous codons (SCs) ending in either cytosine (NNC) or guanidine (NNG) were more frequent than those ending in either adenosine (NNA) or thymine (NNT), and the preference for NNC/G codons followed the increase in genome ploidy. The ratios between NNC/G and NNA/T codons gradually decreased in genes with more introns, and the difference in these ratios between wheat and its progenitors diminished with increasing ploidy. SCUB frequencies were heterogeneous among exons, and the bias preferred to NNA/T in more internal exons, especially for genes with more exons; while the preference did not appear to associate with ploidy. The SCUB alteration of the progenitors was different during the formation of hexaploid wheat, so that SCUB was the homogeneous among A, B and D subgenomes. DNA methylation-mediated conversion from cytosine to thymine weakened following the increase of genome ploidy, coinciding with the stronger bias for NNC/G SCs in the genome as a function of ploidy, suggesting that SCUB contribute to the epigenetic variation in hexaploid wheat. The patterns in SCUB mirrored the formation of hexaploid wheat, which provides new insight into genome shock-induced genetic variation during polyploidization. SCs representing non-neutral synonymous mutations can be used for genetic dissection and improvement of agricultural traits of wheat and other polyploidies.

Multiple origins, one evolutionary trajectory: gradual evolution characterizes distinct lineages of allotetraploid Brachypodium

The "genomic shock" hypothesis posits that unusual challenges to genome integrity such as whole genome duplication may induce chaotic genome restructuring. Decades of research on polyploid genomes have revealed that this is often, but not always the case. While some polyploids show major chromosomal rearrangements and derepression of transposable elements in the immediate aftermath of whole genome duplication, others do not. Nonetheless, all polyploids show gradual diploidization over evolutionary time. To evaluate these hypotheses, we produced a chromosome-scale reference genome for the natural allotetraploid grass Brachypodium hybridum, accession "Bhyb26." We compared 2 independently derived accessions of B. hybridum and their deeply diverged diploid progenitor species Brachypodium stacei and Brachypodium distachyon. The 2 B. hybridum lineages provide a natural timecourse in genome evolution because one formed 1.4 million years ago, and the other formed 140 thousand years ago. The genome of the older lineage reveals signs of gradual post-whole genome duplication genome evolution including minor gene loss and genome rearrangement that are missing from the younger lineage. In neither B. hybridum lineage do we find signs of homeologous recombination or pronounced transposable element activation, though we find evidence supporting steady post-whole genome duplication transposable element activity in the older lineage. Gene loss in the older lineage was slightly biased toward 1 subgenome, but genome dominance was not observed at the transcriptomic level. We propose that relaxed selection, rather than an abrupt genomic shock, drives evolutionary novelty in B. hybridum, and that the progenitor species' similarity in transposable element load may account for the subtlety of the observed genome dominance.

Systematic and functional analysis of non-specific lipid transfer protein family genes in sugarcane under Xanthomonas albilineans infection and salicylic acid treatment

Plant non-specific lipid transfer proteins (nsLTPs) are small basic proteins that play a significant regulatory role in a wide range of physiological processes. To date, no genome-wide survey and expression analysis of this gene family in sugarcane has been performed. In this study we identified the nsLTP gene family in Saccharum spontaneum and carried out expression profiling of nsLTPs in two sugarcane cultivars (Saccharum spp.) that have different resistance to leaf scald caused by Xanthomonas albilineans (Xa) infection. The effect of stress related to exogenous salicylic acid (SA) treatment was also examined. At a genome-wide level, S. spontaneum AP85-441 had 71 SsnsLTP genes including 66 alleles. Tandem (9 gene pairs) and segmental (36 gene pairs) duplication events contributed to SsnsLTP gene family expansion. Five SsnsLTP proteins were predicted to interact with five other proteins. Expression of ShnsLTPI.8/10/Gb.1 genes was significantly upregulated in LCP85-384 (resistant cultivar), but downregulated in ROC20 (susceptible cultivar), suggesting that these genes play a positive regulatory role in response of sugarcane to Xa infection. Conversely, ShnsLTPGa.4/Ge.3 appears to act as a negative regulator in response Xa infection. The majority (16/17) of tested genes were positively induced in LCP85-384 72 h after SA treatment. In both cultivars, but particularly in LCP85-384, ShnsLTPIV.3/VIII.1 genes were upregulated at all time-points, suggesting that the two genes might act as positive regulators under SA stress. Meanwhile, both cultivars showed downregulated ShnsLTPGb.1 gene expression, indicating its potential negative role in SA treatment responses. Notably, the ShnsLTPGb.1 gene had contrasting effects, with positive regulation of gene expression in response to Xa infection and negative regulation induced by SA stress. Together, our results provide valuable information for elucidating the function of ShnsLTP family members under two stressors and identified novel gene sources for development of sugarcane that are tolerant of environmental stimuli.

A high-quality Buxus austro-yunnanensis (Buxales) genome provides new insights into karyotype evolution in early eudicots

BACKGROUND

Eudicots are the most diverse group of flowering plants that compromise five well-defined lineages: core eudicots, Ranunculales, Proteales, Trochodendrales, and Buxales. However, the phylogenetic relationships between these five lineages and their chromosomal evolutions remain unclear, and a lack of high-quality genome analyses for Buxales has hindered many efforts to address this knowledge gap.

RESULTS

Here, we present a high-quality chromosome-level genome of Buxus austro-yunnanensis (Buxales). Our phylogenomic analyses revealed that Buxales and Trochodendrales are genetically similar and classified as sisters. Additionally, both are sisters to the core eudicots, while Ranunculales was found to be the first lineage to diverge from these groups. Incomplete lineage sorting and hybridization were identified as the main contributors to phylogenetic discordance (34.33%) between the lineages. In fact, B. austro-yunnanensis underwent only one whole-genome duplication event, and collinear gene phylogeny analyses suggested that separate independent polyploidizations occurred in the five eudicot lineages. Using representative genomes from these five lineages, we reconstructed the ancestral eudicot karyotype (AEK) and generated a nearly gapless karyotype projection for each eudicot species. Within core eudicots, we recovered one common chromosome fusion event in asterids and malvids, respectively. Further, we also found that the previously reported fused AEKs in Aquilegia (Ranunculales) and Vitis (core eudicots) have different fusion positions, which indicates that these two species have different karyotype evolution histories.

CONCLUSIONS

Based on our phylogenomic and karyotype evolution analyses, we revealed the likely relationships and evolutionary histories of early eudicots. Ultimately, our study expands genomic resources for early-diverging eudicots.