Genomic and expression plasticity of polyploidy

Polyploidy-mediated variations in glutamate receptor proteins linked to Fusarium wilt resistance in upland cotton

Cotton production in the US faces a serious threat from Fusarium oxysporum f. sp. vasinfectum race 4 (FOV4), a soil-borne fungus causing Fusarium wilt by infecting the roots and vascular system of susceptible cotton, leading to rapid wilting and death. Here, we investigate genetic mechanisms of resistance to FOV4 in the highly resistant upland cotton genotype "U1" using an early-generation segregating biparental population ("U1" × "CSX8308") with comprehensive genomic resources. Reference-grade genomic assemblies of the parents revealed minor structural variations between "U1" haplotypes, a high degree of collinearity at chromosome synteny and micro-synteny levels, and significant divergence from "CSX8308" with 8.9 million SNPs. QTL analysis identified significant markers on chromosomes D03 and A02 linked to reduced Fusarium wilt severity. Within these regions, two glutamate-receptor-like (GLR) genes showed structural variation and overlapped between translocated segments on A02 and D03, suggesting a rare but important reinforcing effect of parallel evolution between susceptible and resistant genotypes. Transcriptome profiles of "U1" under FOV4 infection reveal activation of calcium-binding proteins and transcription factors regulating plant hormones (ethylene, abscisic acid, jasmonic acid, and salicylic acid), along with enzymes involved in cell wall remodeling and phytoalexin production. Advancing cotton improvement depends on incorporating durable genetic disease resistance into high-yielding, high-quality cultivars.

Ecological divergence of sibling allopolyploid marsh orchids is associated with species specific plasticity and distinct fungal communities

Phenotypic plasticity, the dynamic adjustment of traits to environmental variations, is crucial for enabling species to exploit broader niches and withstand suboptimal conditions. This adaptability is particularly relevant for newly formed allopolyploids, which possess redundant gene copies and must become established in diverse environments distinct from their parents and other relatives. By evaluating gene expression and root mycobiome among two ecologically divergent sibling allopolyploid marsh orchids (Dactylorhiza majalis and D. traunsteineri) in reciprocal transplants at localities where both species are native, we aimed to understand the drivers of species persistence in the face of interspecific gene flow. Despite consistent abiotic differences characterising the alternative environments at each locality, the majority of gene expression differences between the allopolyploids appears to be plastic. Ecologically relevant processes, such as photosynthesis and transmembrane transport, include some genes that are differentially expressed between the two orchids regardless of the environment, while others change their activity plastically in one species or the other. This suggests that although plasticity helps define the specific ecological range of each sibling allopolyploid, it also mediates gene flow between them, thereby preventing differentiation. Extending our investigations to the root mycobiome, we uncover more diverse fungal communities for either species when grown in the environment with nutrient-poor soils, indicating that both abiotic and biotic factors drive the distribution of sibling marsh orchids. Altogether, our results indicate that plasticity can simultaneously promote diversification and homogenisation of lineages, influencing the establishment and persistence of recurrently formed allopolyploid species.

Refining polyploid breeding in sweet potato through allele dosage enhancement

Allele dosage plays a key role in the phenotypic variation of polyploids. Here we present a genome-wide variation map of hexaploid sweet potato that captures allele dosage information, constructed from deep sequencing of 294 hexaploid accessions. Genome-wide association studies identified quantitative trait loci with dosage effects on 23 agronomic traits. Our analyses reveal that sweet potato breeding has progressively increased the dosage of favourable alleles to enhance trait performance. Notably, the Mesoamerican gene pool has evolved towards higher dosages of favourable alleles at multiple loci, which have been increasingly introgressed into modern Chinese cultivars. We substantiated the breeding-driven dosage accumulation through transgenic validation of IbEXPA4, an expansin gene influencing tuberous root weight. In addition, we explored causative sequence variations that alter the expression of the Orange gene, which regulates flesh colour. Our findings illuminate the breeding history of sweet potato and establish a foundation for leveraging allele dosages in polyploid breeding practices.

Dosage effect genes modulate grain development in synthesized Triticum durum-Haynaldia villosa allohexaploid

Polyploidization in plants often leads to increased cell size and grain size, which may be affected by the increased genome dosage and transcription abundance. The synthesized Triticum durum (AABB)-Haynaldia villosa (VV) amphiploid (AABBVV) has significantly increased grain size, especially grain length, than the tetraploid and diploid parents. To investigate how polyploidization affects grain development at the transcriptional level, we perform transcriptome analysis using the immature seeds of T. durum, H. villosa, and the amphiploid. The dosage effect genes are contributed more by differentially expressed genes from genome V of H. villosa. The dosage effect genes overrepresent grain development-related genes. Interestingly, the vernalization gene TaVRN1 is among the positive dosage effect genes in the T. durum‒H. villosa and T. turgidum‒Ae. tauschii amphiploids. The expression levels of TaVRN1 homologs are positively correlated with the grain size and weight. The TaVRN1-B1 or TaVRN1-D1 mutation shows delayed florescence, decreased cell size, grain size, and grain yield. These data indicate that dosage effect genes could be one of the important explanations for increased grain size by regulating grain development. The identification and functional validation of dosage effect genes may facilitate the finding of valuable genes for improving wheat yield.

Strawberry phenotypic plasticity in flowering time is driven by the interaction between genetic loci and temperature

Flowering time (FT), which determines when fruits or seeds can be harvested, is subject to phenotypic plasticity, that is, the ability of a genotype to display different phenotypes in response to environmental variation. Here, we investigated how the environment affects the genetic architecture of FT in cultivated strawberry (Fragaria × ananassa) and modifies its quantitative trait locus (QTL) effects. To this end, we used a bi-parental segregating population grown for 2 years at widely divergent latitudes (five European countries) and combined climatic variables with genomic data (Affymetrix SNP array). Examination, using different phenological models, of the response of FT to photoperiod, temperature, and global radiation indicated that temperature is the main driver of FT in strawberry. We next characterized in the segregating population the phenotypic plasticity of FT by using three statistical approaches that generated plasticity parameters including reaction norm parameters. We detected 25 FT QTLs summarized as 10 unique QTLs. Mean values and plasticity parameter QTLs were co-localized in three of them, including the major 6D_M QTL whose effect is strongly modulated by temperature. The design and validation of a genetic marker for the 6D_M QTL offers great potential for breeding programs, for example selecting early-flowering strawberry varieties well adapted to different environmental conditions.

Integrative analyses of long and short-read RNA sequencing reveal the spliced isoform regulatory network of seedling growth dynamics in upland cotton

The polyploid genome of cotton has significantly increased the transcript complexity. Recent advances in full-length transcript sequencing are now widely used to characterize the complete landscape of transcriptional events. Such studies in cotton can help us to explore the genetic mechanisms of the cotton seedling growth. Through long-read single-molecule RNA sequencing, this study compared the transcriptomes of three yield contrasting genotypes of upland cotton. Our analysis identified different numbers of spliced isoforms from 31,166, 28,716, and 28,713 genes in SJ48, Z98, and DT8 cotton genotypes, respectively, most of which were novel compared to previous cotton reference transcriptomes, and showed significant differences in the number of exon structures and coding sequence length due to intron retention. Quantification of isoform expression revealed significant differences in expression in the root and leaf of each genotype. An array of key isoform target genes showed protein kinase or phosphorylation functions, and their protein interaction network contained most of the circadian oscillator proteins. Spliced isoforms from the GIGANTEA (GI) protien were differentially regulated in each genotype and might be expected to regulate translational activities, including the sequence and function of target proteins. In addition, these spliced isoforms generate diurnal expression profiles in cotton leaves, which may alter the transcriptional regulatory network of seedling growth. Silencing of the novel spliced GI isoform Gh_A02G0645_N17 significantly affected biomass traits, contributed to variable growth, and increased transcription of the early flowering pathway gene ELF in cotton. Our high-throughput hybrid sequencing results will be useful to dissect functional differences among spliced isoforms in the polyploid cotton genome.

Polyploidy in maize: from evolution to breeding

Polyploidy played an important role in the evolution of the three most important crops: wheat, maize and rice, each of them providing a unique model for studying allopolyploidy, segmental alloploidy or paleopolyploidy. However, its genetic and evolutionary role is still vague. The undelying mechanisms and consequences of polyploidy remain fundamental objectives in the study of eukaryotes. Maize is one of the underutilized crops at the polyploid level. This species has no stable natural polyploids, the existing ones being artificially obtained. From the experimental polyploid series of maize, only the tetraploid forms (4n = 40) are of interest. They are characterized by some valuable morphological, physiological and biochemical features, superior to the diploid forms from which they originated, but also by some drawbacks such as: reduced fertility, slower development, longer vegetation period, low productivity and adaptedness. Due to these barriers to using tetraploids in field production, maize tetraploids primarily found utility in scientific studies regarding genetic variability, inbreeding, heterosis and gene dosage effect. Since the first mention of a triploid maize plant to present, many scientists and schools, devoted their efforts to capitalize on the use of polyploidy in maize. Despite its common disadvantages as a crop, significant progress in developing tetraploid maize with good agronomic performance was achieved leading to registered tetraploid maize varieties. In this review we summarize and discuss the different aspects of polyploidy in maize, such as evolutionary context, methods of induction, morphology, fertility issue, inheritance patterns, gene expression and potential use.

Characteristics of duplicated gene expression and DNA methylation regulation in different tissues of allopolyploid Brassica napus

Plant polyploidization increases the complexity of epigenomes and transcriptional regulation, resulting in genome evolution and enhanced adaptability. However, few studies have been conducted on the relationship between gene expression and epigenetic modification in different plant tissues after allopolyploidization. In this study, we studied gene expression and DNA methylation modification patterns in four tissues (stems, leaves, flowers and siliques) of Brassica napusand its diploid progenitors. On this basis, the alternative splicing patterns and cis-trans regulation patterns of four tissues in B. napus and its diploid progenitors were also analyzed. It can be seen that the number of alternative splicing occurs in the B. napus is higher than that in the diploid progenitors, and the IR type increases the most during allopolyploidy. In addition, we studied the fate changes of duplicated genes after allopolyploidization in B. napus. We found that the fate of most duplicated genes is conserved, but the number of neofunctionalization and specialization is also large. The genetic fate of B. napus was classified according to five replication types (WGD, PD, DSD, TD, TRD). This study also analyzed generational transmission analysis of expression and DNA methylation patterns. Our study provides a reference for the fate differentiation of duplicated genes during allopolyploidization.

Water stress modulates terpene biosynthesis and morphophysiology at different ploidal levels in Lippia alba (Mill.) N. E. Brown (Verbenaceae)

Monoterpenes are the main component in essential oils of Lippia alba. In this species, the chemical composition of essential oils varies with genome size: citral (geraniol and neral) is dominant in diploids and tetraploids, and linalool in triploids. Because environmental stress impacts various metabolic pathways, we hypothesized that stress responses in L. alba could alter the relationship between genome size and essential oil composition. Water stress affects the flowering, production, and reproduction of plants. Here, we evaluated the effect of water stress on morphophysiology, essential oil production, and the expression of genes related to monoterpene synthesis in diploid, triploid, and tetraploid accessions of L. alba cultivated in vitro for 40 days. First, using transcriptome data, we performed de novo gene assembly and identified orthologous genes using phylogenetic and clustering-based approaches. The expression of candidate genes related to terpene biosynthesis was estimated by real-time quantitative PCR. Next, we assessed the expression of these genes under water stress conditions, whereby 1% PEG-4000 was added to MS medium. Water stress modulated L. alba morphophysiology at all ploidal levels. Gene expression and essential oil production were affected in triploid accessions. Polyploid accessions showed greater growth and metabolic tolerance under stress compared to diploids. These results confirm the complex regulation of metabolic pathways such as the production of essential oils in polyploid genomes. In addition, they highlight aspects of genotype and environment interactions, which may be important for the conservation of tropical biodiversity.

Growth, ultrastructural and physiological characteristics of Abelmoschus cytotypes under elevated ozone stress: a study on ploidy-specific responses

Tropospheric ozone (O3 ) is a significant abiotic stressor whose rising concentration negatively influences plant growth. Studies related to the differential response of Abelmoschus cytotypes to elevated O3 treatment are scarce and need further exploration to recognise the role of polyploidisation in stress tolerance. In this study, we analysed the changes in growth pattern, ultrastructure, physiology and foliar protein profile occurring under O3 stress in Abelmoschus moschatus (monoploid), Abelmoschus esculentus (diploid) and Abelmoschus caillei (triploid). Our findings showed that higher stomatal conductance in A. moschatus triggered higher O3 intake, causing damage to stomatal cells and photosynthetic pigments. Additionally, it caused a reduction in photosynthetic rates, leading to reduced plant growth, total biomass and economic yield. This O3 -induced toxicity was less in diploid and triploid cytotypes of Abelmoschus . Protein profiling by sodium dodecyl sulpate-polyacrylamide gel electrophoresis showed a significant decrease in the commonly found RuBisCO larger and smaller subunits. The decrease was more prominent in monoploid compared to diploid and triploid. This study provides crucial data for research that aim to enhance plant ability to withstand O3 induced oxidative stress. Our findings may help in developing a tolerant variety through plant breeding techniques, which will be economically more advantageous in reaching the objective of sustainable production at the high O3 levels projected under a climate change scenario.

Investigation of regulatory divergence between homoeologs in the recently formed allopolyploids, Tragopogon mirus and T. miscellus (Asteraceae)

Polyploidy is an important evolutionary process throughout eukaryotes, particularly in flowering plants. Duplicated gene pairs (homoeologs) in allopolyploids provide additional genetic resources for changes in molecular, biochemical, and physiological mechanisms that result in evolutionary novelty. Therefore, understanding how divergent genomes and their regulatory networks reconcile is vital for unraveling the role of polyploidy in plant evolution. Here, we compared the leaf transcriptomes of recently formed natural allotetraploids (Tragopogon mirus and T. miscellus) and their diploid parents (T. porrifolius X T. dubius and T. pratensis X T. dubius, respectively). Analysis of 35 400 expressed loci showed a significantly higher level of transcriptomic additivity compared to old polyploids; only 22% were non-additively expressed in the polyploids, with 5.9% exhibiting transgressive expression (lower or higher expression in the polyploids than in the diploid parents). Among approximately 7400 common orthologous regions (COREs), most loci in both allopolyploids exhibited expression patterns that were vertically inherited from their diploid parents. However, 18% and 20.3% of the loci showed novel expression bias patterns in T. mirus and T. miscellus, respectively. The expression changes of 1500 COREs were explained by cis-regulatory divergence (the condition in which the two parental subgenomes do not interact) between the diploid parents, whereas only about 423 and 461 of the gene expression changes represent trans-effects (the two parental subgenomes interact) in T. mirus and T. miscellus, respectively. The low degree of both non-additivity and trans-effects on gene expression may present the ongoing evolutionary processes of the newly formed Tragopogon polyploids (~80-90 years).

Genome resources for the elite bread wheat cultivar Aikang 58 and mining of elite homeologous haplotypes for accelerating wheat improvement

Despite recent progress in crop genomics studies, the genomic changes brought about by modern breeding selection are still poorly understood, thus hampering genomics-assisted breeding, especially in polyploid crops with compound genomes such as common wheat (Triticum aestivum). In this work, we constructed genome resources for the modern elite common wheat variety Aikang 58 (AK58). Comparative genomics between AK58 and the landrace cultivar Chinese Spring (CS) shed light on genomic changes that occurred through recent varietal improvement. We also explored subgenome diploidization and divergence in common wheat and developed a homoeologous locus-based genome-wide association study (HGWAS) approach, which was more effective than single homoeolog-based GWAS in unraveling agronomic trait-associated loci. A total of 123 major HGWAS loci were detected using a genetic population derived from AK58 and CS. Elite homoeologous haplotypes (HHs), formed by combinations of subgenomic homoeologs of the associated loci, were found in both parents and progeny, and many could substantially improve wheat yield and related traits. We built a website where users can download genome assembly sequence and annotation data for AK58, perform blast analysis, and run JBrowse. Our work enriches genome resources for wheat, provides new insights into genomic changes during modern wheat improvement, and suggests that efficient mining of elite HHs can make a substantial contribution to genomics-assisted breeding in common wheat and other polyploid crops.

Horizontal acquisition of a DNA ligase improves DNA damage tolerance in eukaryotes

Bdelloid rotifers are part of the restricted circle of multicellular animals that can withstand a wide range of genotoxic stresses at any stage of their life cycle. In this study, bdelloid rotifer Adineta vaga is used as a model to decipher the molecular basis of their extreme tolerance. Proteomic analysis shows that a specific DNA ligase, different from those usually involved in DNA repair in eukaryotes, is strongly over-represented upon ionizing radiation. A phylogenetic analysis reveals its orthology to prokaryotic DNA ligase E, and its horizontal acquisition by bdelloid rotifers and plausibly other eukaryotes. The fungus Mortierella verticillata, having a single copy of this DNA Ligase E homolog, also exhibits an increased radiation tolerance with an over-expression of this DNA ligase E following X-ray exposure. We also provide evidence that A. vaga ligase E is a major contributor of DNA breaks ligation activity, which is a common step of all important DNA repair pathways. Consistently, its heterologous expression in human cell lines significantly improves their radio-tolerance. Overall, this study highlights the potential of horizontal gene transfers in eukaryotes, and their contribution to the adaptation to extreme conditions.

Exploring the Role of TaPLC1-2B in Heat Tolerance at Seedling and Adult Stages of Wheat through Transcriptome Analysis

Heat stress is a major abiotic stress that can cause serious losses of a crop. Our previous work identified a gene involved in heat stress tolerance in wheat, TaPLC1-2B. To further investigate its mechanisms, in the present study, TaPLC1-2B RNAi-silenced transgenic wheat and the wild type were comparatively analyzed at both the seedling and adult stages, with or without heat stress, using transcriptome sequencing. A total of 15,549 differentially expressed genes (DEGs) were identified at the adult stage and 20,535 DEGs were detected at the seedling stage. After heat stress, an enrichment of pathways such as phytohormones and mitogen-activated protein kinase signaling was mainly found in the seedling stage, and pathways related to metabolism, glycerophospholipid metabolism, circadian rhythms, and ABC transporter were enriched in the adult stage. Auxin and abscisic acid were downregulated in the seedling stage and vice versa in the adult stage; and the MYB, WRKY, and no apical meristem gene families were downregulated in the seedling stage in response to heat stress and upregulated in the adult stage in response to heat stress. This study deepens our understanding of the mechanisms of TaPLC1-2B in regard to heat stress in wheat at the seedling and adult stages.

Polyploidization and genomic selection integration for grapevine breeding: a perspective

Grapevines are economically important woody perennial crops widely cultivated for their fruits that are used for making wine, grape juice, raisins, and table grapes. However, grapevine production is constantly facing challenges due to climate change and the prevalence of pests and diseases, causing yield reduction, lower fruit quality, and financial losses. To ease the burden, continuous crop improvement to develop superior grape genotypes with desirable traits is imperative. Polyploidization has emerged as a promising tool to generate genotypes with novel genetic combinations that can confer desirable traits such as enhanced organ size, improved fruit quality, and increased resistance to both biotic and abiotic stresses. While previous studies have shown high polyploid induction rates in Vitis spp., rigorous screening of genotypes among the produced polyploids to identify those exhibiting desired traits remains a major bottleneck. In this perspective, we propose the integration of the genomic selection approach with omics data to predict genotypes with desirable traits among the vast unique individuals generated through polyploidization. This integrated approach can be a powerful tool for accelerating the breeding of grapevines to develop novel and improved grapevine varieties.

Exploring the genetic landscape of nitrogen uptake in durum wheat: genome-wide characterization and expression profiling of NPF and NRT2 gene families

Nitrate uptake by plants primarily relies on two gene families: Nitrate transporter 1/peptide transporter (NPF) and Nitrate transporter 2 (NRT2). Here, we extensively characterized the NPF and NRT2 families in the durum wheat genome, revealing 211 NPF and 20 NRT2 genes. The two families share many Cis Regulatory Elements (CREs) and Transcription Factor binding sites, highlighting a partially overlapping regulatory system and suggesting a coordinated response for nitrate transport and utilization. Analyzing RNA-seq data from 9 tissues and 20 cultivars, we explored expression profiles and co-expression relationships of both gene families. We observed a strong correlation between nucleotide variation and gene expression within the NRT2 gene family, implicating a shared selection mechanism operating on both coding and regulatory regions. Furthermore, NPF genes showed highly tissue-specific expression profiles, while NRT2s were mainly divided in two co-expression modules, one expressed in roots (NAR2/NRT3 dependent) and the other induced in anthers and/ovaries during maturation. Our evidences confirmed that the majority of these genes were retained after small-scale duplication events, suggesting a neo- or sub-functionalization of many NPFs and NRT2s. Altogether, these findings indicate that the expansion of these gene families in durum wheat could provide valuable genetic variability useful to identify NUE-related and candidate genes for future breeding programs in the context of low-impact and sustainable agriculture.

Genome shock in a synthetic allotetraploid wheat invokes subgenome-partitioned gene regulation, meiotic instability, and karyotype variation

It is becoming increasingly evident that interspecific hybridization at the homoploid level or coupled with whole-genome duplication (i.e. allopolyploidization) has played a major role in biological evolution. However, the direct impacts of hybridization and allopolyploidization on genome structure and function, phenotype, and fitness remains to be fully understood. Synthetic hybrids and allopolyploids are trackable experimental systems that can be used to address this issue. In this study, we resynthesized a pair of reciprocal F1 hybrids and corresponding reciprocal allotetraploids using the two diploid progenitor species of bread wheat (Triticum aestivum, BBAADD), namely T. urartu (AA) and Aegilops tauschii (DD). By comparing phenotypes related to growth, development, and fitness, and by analysing genome expression in both hybrids and allotetraploids in relation to the parents, we found that the types and trends of karyotype variation in the immediately formed allotetraploids were correlated with both instability of meiosis and chromosome- and subgenome-biased expression. We determined clear advantages of allotetraploids over diploid F1 hybrids in several morphological traits including fitness that mirrored the tissue- and developmental stage-dependent subgenome-partitioning of the allotetraploids. The allotetraploids were meiotically unstable primarily due to homoeologous pairing that varied dramatically among the chromosomes. Nonetheless, the manifestation of organismal karyotype variation and the occurrence of meiotic irregularity were not concordant, suggesting a role of functional constraints probably imposed by subgenome- and chromosome-biased gene expression. Our results provide new insights into the direct impacts and consequences of hybridization and allopolyploidization that are relevant to evolution and likely to be informative for future crop improvement approaches using synthetic polyploids.

Genome-Wide Identification of TCP Gene Family in Dendrobium and Their Expression Patterns in Dendrobium chrysotoxum

The TCP gene family are plant-specific transcription factors that play important roles in plant growth and development. Dendrobium chrysotoxum, D. nobile, and D. huoshanense are orchids with a high ornamental value, but few studies have investigated the specific functions of TCPs in Dendrobium flower development. In this study, we used these three Dendrobium species to analyze TCPs, examining their physicochemical properties, phylogenetic relationships, gene structures, and expression profiles. A total of 50 TCPs were identified across three Dendrobium species; they were divided into two clades-Class-I (PCF subfamily) and Class-II (CIN and CYC/TB1 subfamilies)-based on their phylogenetic relationships. Our sequence logo analysis showed that almost all Dendrobium TCPs contain a conserved TCP domain, as well as the existence of fewer exons, and the cis-regulatory elements of the TCPs were mostly related to light response. In addition, our transcriptomic data and qRT-PCR results showed that DchTCP2 and DchTCP13 had a significant impact on lateral organs. Moreover, changes in the expression level of DchTCP4 suggested its important role in the phenotypic variation of floral organs. Therefore, this study provides a significant reference for the further exploration of TCP gene functions in the regulation of different floral organs in Dendrobium orchids.

Developmental Instability and Gene Dysregulation in an Extracted Tetraploid from Hexaploid Wheat

The BBAA subgenomes of hexaploid common wheat can be 'extracted' to constitute a viable and self-reproducing novel tetraploid wheat, termed extracted tetraploid wheat (ETW). Prior studies have shown ETW manifesting phenotypic abnormalities and alteration in gene expression and epigenetic modifications. No population level investigation has been conducted, leaving the issue unclear regarding whether developmental stability, an essential property evolved in all natural organisms, might have been undermined in ETW. Here, we measured variations in five morphological traits and somatic chromosomal stability in populations of ETW and of its hexaploid donor, a resynthesized hexaploid and a natural tetraploid wheat. We observed phenotypic defects in ETW. Meanwhile, we documented much greater within-population variations in ETW than in the other wheat genotypes, most probably due to disrupted developmental stability in ETW. Also, somatic structural chromosome variations were detected only in ETW. Comparative transcriptome analyses indicated that the disrupted developmental stability of ETW is likely linked to massive dysregulation of genome-wide gene expression rather than to genetic mutations. Population network analysis of gene expression implicated intrinsic connectivity among the variable traits, while gene set enrichment analysis provided possible links between dysregulated gene expression and interlaced trait variation.

Effects of Allopolyploidization and Homoeologous Chromosomal Segment Exchange on Homoeolog Expression in a Synthetic Allotetraploid Wheat under Variable Environmental Conditions

Allopolyploidy through the combination of divergent genomes into a common nucleus at doubled dosage is known as a potent genetic and evolutionary force. As a macromutation, a striking feature of allopolyploidy in comparison with other mutational processes is that 'genome shock' can be evoked, thereby generating rapid and saltational biological consequences. A major manifestation of genome shock is genome-wide gene expression rewiring, which previously remained to be fully elucidated. Here, using a large set of RNAseq-based transcriptomic data of a synthetic allotetraploid wheat (genome AADD) and its parental species, we performed in-depth analyses of changes in the genome-wide gene expression under diverse environmental conditions at the subgenome (homoeolog) level and investigated the additional effects of homoeologous chromosomal segment exchanges (abbreviated HEs). We show that allopolyploidy caused large-scale changes in gene expression that were variable across the conditions and exacerbated by both stresses and HEs. Moreover, although both subgenomes (A and D) showed clear commonality in the changes, they responded differentially under variable conditions. The subgenome- and condition-dependent differentially expressed genes were enriched for different gene ontology terms implicating different biological functions. Our results provide new insights into the direct impacts of allopolyploidy on condition-dependent changes in subgenome expression and the additional effects of HEs in nascent allopolyploidy.

Genetic factors inherited from both diploid parents interact to affect genome stability and fertility in resynthesized allotetraploid Brassica napus

Established allopolyploids are known to be genomically stable and fertile. However, in contrast, most newly resynthesized allopolyploids are infertile and meiotically unstable. Identifying the genetic factors responsible for genome stability in newly formed allopolyploid is key to understanding how 2 genomes come together to form a species. One hypothesis is that established allopolyploids may have inherited specific alleles from their diploid progenitors which conferred meiotic stability. Resynthesized Brassica napus lines are often unstable and infertile, unlike B. napus cultivars. We tested this hypothesis by characterizing 41 resynthesized B. napus lines produced by crosses between 8 Brassica rapa and 8 Brassica oleracea lines for copy number variation resulting from nonhomologous recombination events and fertility. We resequenced 8 B. rapa and 5 B. oleracea parent accessions and analyzed 19 resynthesized lines for allelic variation in a list of meiosis gene homologs. SNP genotyping was performed using the Illumina Infinium Brassica 60K array for 3 individuals per line. Self-pollinated seed set and genome stability (number of copy number variants) were significantly affected by the interaction between both B. rapa and B. oleracea parental genotypes. We identified 13 putative meiosis gene candidates which were significantly associated with frequency of copy number variants and which contained putatively harmful mutations in meiosis gene haplotypes for further investigation. Our results support the hypothesis that allelic variants inherited from parental genotypes affect genome stability and fertility in resynthesized rapeseed.

Chromosome-specific painting in Thinopyrum species using bulked oligonucleotides

Chromosome-specific painting probes were developed to identify the individual chromosomes from 1 to 7E in Thinopyrum species and detect alien genetic material of the E genome in a wheat background. The E genome of Thinopyrum is closely related to the ABD genome of wheat (Triticum aestivum L.) and harbors genes conferring beneficial traits to wheat, including high yield, disease resistance, and unique end-use quality. Species of Thinopyrum vary from diploid (2n = 2x = 14) to decaploid (2n = 10x = 70), and chromosome structural variation and differentiation have arisen during polyploidization. To investigate the variation and evolution of the E genome, we developed a complete set of E genome-specific painting probes for identification of the individual chromosomes 1E to 7E based on the genome sequences of Th. elongatum (Host) D. R. Dewey and wheat. By using these new probes in oligonucleotide-based chromosome painting, we showed that Th. bessarabicum (PI 531711, EbEb) has a close genetic relationship with diploid Th. elongatum (EeEe), with five chromosomes (1E, 2E, 3E, 6E, and 7E) maintaining complete synteny in the two species except for a reciprocal translocation between 4 and 5Eb. All 14 pairs of chromosomes of tetraploid Th. elongatum have maintained complete synteny with those of diploid Th. elongatum (Thy14), but the two sets of E genomes have diverged. This study also demonstrated that the E genome-specific painting probes are useful for rapid and effective detection of the alien genetic material of E genome in wheat-Thinopyrum derived lines.

Genome-Wide Identification and Analysis of the Plant Cysteine Oxidase (PCO) Gene Family in Brassica napus and Its Role in Abiotic Stress Response

Plant Cysteine Oxidase (PCO) is a plant O₂-sensing enzyme catalyzing the oxidation of cysteine to Cys-sulfinic acid at the N-termini of target proteins. To better understand the Brassica napus PCO gene family, PCO genes in B. napus and related species were analyzed. In this study, 20, 7 and 8 PCO genes were identified in Brassica napus, Brassica rapa and Brassica oleracea, respectively. According to phylogenetic analysis, the PCOs were divided into five groups: PCO1, PCO2, PCO3, PCO4 and PCO5. Gene organization and motif distribution analysis suggested that the PCO gene family was relatively conserved during evolution. According to the public expression data, PCO genes were expressed in different tissues at different developmental stages. Moreover, qRT-PCR data showed that most of the Bna/Bra/BoPCO5 members were expressed in leaves, roots, flowers and siliques, suggesting an important role in both vegetative and reproductive development. Expression of BnaPCO was induced by various abiotic stress, especially waterlogging stress, which was consistent with the result of cis-element analysis. In this study, the PCO gene family of Brassicaceae was analyzed for the first time, which contributes to a comprehensive understanding of the origin and evolution of PCO genes in Brassicaceae and the function of BnaPCO in abiotic stress responses.

Autopolyploids of Arabidopsis thaliana are more phenotypically plastic than their diploid progenitors

BACKGROUND AND AIMS

Polyploids are often hypothesized to have increased phenotypic plasticity compared with their diploid progenitors, but recent work suggests that the relationship between whole-genome duplication (WGD) and plasticity is not so straightforward. Impacts of WGD on plasticity are moderated by other evolutionary processes in nature, which has impeded generalizations regarding the effects of WGD alone. We assessed shifts in phenotypic plasticity and mean trait values accompanying WGD, as well as the adaptive consequences of these shifts.

METHODS

To isolate WGD effects, we compared two diploid lineages of Arabidopsis thaliana wiht corresponding autotetraploids grown across different salt and nutrient conditions in a growth chamber.

KEY RESULTS

For the few cases in which diploids and polyploids differed in plasticity, polyploids were more plastic, consistent with hypotheses that WGD increases plasticity. Under stress, increased plasticity was often adaptive (associated with higher total seed mass), but in other cases plasticity was unrelated to fitness. Mean trait values and plasticity were equally likely to be affected by WGD, but the adaptive consequences of these shifts were often context dependent or lineage specific. For example, polyploids had extended life spans, a shift that was adaptive in one polyploid lineage under amenable conditions but was maladaptive in the other lineage under stress.

CONCLUSIONS

Our work shows that increased phenotypic plasticity can result from WGD alone, independent of other evolutionary processes. We find that the effects of WGD can differ depending on the genotype of the progenitor and the environmental context. Though our experiment was limited to two genotypes of a single species, these findings support the idea that WGD can indeed increase plasticity.

Recent allopolyploidy alters Spartina microRNA expression in response to xenobiotic-induced stress

Environmental contamination by xenobiotics represents a major threat for natural ecosystems and public health. In response, xenobiotic detoxification is a fundamental trait of organisms for developmental plasticity and stress tolerance, but the underlying molecular mechanisms remain poorly understood in plants. To decipher this process, we explored the consequences of allopolyploidy on xenobiotic tolerance in the genus Spartina Schreb. Specifically, we focused on microRNAs (miRNAs) owing to their central function in the regulation of gene expression patterns, including responses to stress. Small RNA-Seq was conducted on the parents S. alterniflora and S. maritima, their F1 hybrid S. x townsendii and the allopolyploid S. anglica under phenanthrene-induced stress (phe), a model Polycyclic Aromatic Hydrocarbon (PAH) compound. Differentially expressed miRNAs in response to phe were specifically identified within species. In complement, the respective impacts of hybridization and genome doubling were detected, through changes in miRNA expression patterns between S. x townsendii, S. anglica and the parents. The results support the impact of allopolyploidy in miRNA-guided regulation of plant response to phe. In total, we identified 17 phe-responsive miRNAs in Spartina among up-regulated MIR156 and down-regulated MIR159. We also describe novel phe-responsive miRNAs as putative Spartina-specific gene expression regulators in response to stress. Functional validation using Arabidopsis (L.) Heynh. T-DNA lines inserted in homologous MIR genes was performed, and the divergence of phe-responsive miRNA regulatory networks between Arabidopsis and Spartina was discussed.

Identification of CDPKs involved in TaNOX7 mediated ROS production in wheat

As the critical sensors and decoders of calcium signal, calcium-dependent protein kinase (CDPK) has become the focus of current research, especially in plants. However, few resources are available on the properties and functions of CDPK gene family in Triticum aestivum (TaCDPK). Here, a total of 79 CDPK genes were identified in the wheat genome. These TaCDPKs could be classified into four subgroups on phylogenesis, while they may be classified into two subgroups based on their tissue and organ-spatiotemporal expression profiles or three subgroups according to their induced expression patterns. The analysis on the signal network relationships and interactions of TaCDPKs and NADPH (reduced nicotinamide adenine dinucleotide phosphate oxidases, NOXs), the key producers for reactive oxygen species (ROS), showed that there are complicated cross-talks between these two family proteins. Further experiments demonstrate that, two members of TaCDPKs, TaCDPK2/4, can interact with TaNOX7, an important member of wheat NOXs, and enhanced the TaNOX7-mediated ROS production. All the results suggest that TaCDPKs are highly expressed in wheat with distinct tissue or organ-specificity and stress-inducible diversity, and play vital roles in plant development and response to biotic and abiotic stresses by directly interacting with TaNOXs for ROS production.

Genome assemblies of Vigna reflexo-pilosa (créole bean) and its progenitors, Vigna hirtella and Vigna trinervia, revealed homoeolog expression bias and expression-level dominance in the allotetraploid

Vigna reflexo-pilosa (créole bean) is a wild legume belonging to the subgenus Ceratoropis and is widely distributed in Asia. Créole bean is the only tetraploid species in the genus Vigna, and it has been shown to derive from the hybridization of Vigna hirtella and Vigna trinervia. In this study, we combined the long-read PacBio technology with the chromatin contact mapping (Hi-C) technique to obtain a chromosome-level assembly of V. reflexo-pilosa. The final assembly contained 998,724,903 bases with an N50 length of 42,545,650 bases. Our gene prediction recovered 99.4% of the highly conserved orthologs based on the BUSCO analysis. To investigate homoeolog expression bias and expression level dominance in the tetraploid, we also sequenced and assembled the genomes of its progenitors. Overall, the majority of the homoeolog pairs (72.9%) displayed no expression bias, and among those that exhibited biased expression, 16.3% showed unbalanced homoeolog expression bias toward the V. trinervia subgenome. Moreover, 41.2% and 36.2% of the expressed gene pairs exhibited transgressive expression and expression level dominance, respectively. Interestingly, the genome-wide expression level dominance in the tetraploid was biased toward the V. trinervia subgenome. The analysis of methylation patterns also revealed that the average methylation levels in coding regions were higher in the V. hirtella subgenome than those in the V. trinervia subgenome. The genomic/transcriptomic resources for these three species are useful not only for the development of elite cultivars in Vigna breeding programs but also to researchers studying comparative genomics and investigating genomic/epigenomic changes following polyploid events.

Validation of Reference Genes in a Population of Blueberry (Vaccinium corymbosum) Plants Regenerated in Colchicine

For the first time we report the validation of reference genes in plants from a population of blueberry (Vaccinium corymbosum) clones cultured in vitro on a colchicine-supplemented medium. Nodal segment explants of the cultivar Duke were regenerated by organogenesis under different periods of colchicine 1 mg/L exposure (1, 2, 3, 5, 30 days). The clones selected for the study showed variability for phenotypic traits after 2 years of adaptation to field conditions, compared to plants of the donor genotype that were regenerated on a medium without colchicine. Vaccinium myrtillus (GAPDH) and Vaccinium macrocarpon (ATP1, NADH, RPOB and COX2) were used as reference genomes for primer design. The results show that colchicine treatments can cause genomic changes in blueberry plants. At the molecular level, exposure of plants to colchicine in early periods could promote an increase in gene expression of specific genes such as ATP1, COX2, GAPDH, MATK, NADH and RPOB. However, prolonged exposure (30 days) could decrease gene expression of the genes studied. For qPCR assays, the primers designed for ATP1, COX2, GAPDH and MATK genes showed high efficiency. In addition, the GAPDH, ATP1, NADH and COX2 genes showed high stability and could be recommended as potential reference genes for gene expression assays in Vaccinium.

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.

Variation in cytonuclear expression accommodation among allopolyploid plants

Cytonuclear coevolution is a common feature among plants, which coordinates gene expression and protein products between the nucleus and organelles. Consequently, lineage-specific differences may result in incompatibilities between the nucleus and cytoplasm in hybrid taxa. Allopolyploidy is also a common phenomenon in plant evolution. The hybrid nature of allopolyploids may result in cytonuclear incompatibilities, but the massive nuclear redundancy created during polyploidy affords additional avenues for resolving cytonuclear conflict (i.e. cytonuclear accommodation). Here we evaluate expression changes in organelle-targeted nuclear genes for 6 allopolyploid lineages that represent 4 genera (i.e. Arabidopsis, Arachis, Chenopodium, and Gossypium) and encompass a range in polyploid ages. Because incompatibilities between the nucleus and cytoplasm could potentially result in biases toward the maternal homoeolog and/or maternal expression level, we evaluate patterns of homoeolog usage, expression bias, and expression-level dominance in cytonuclear genes relative to the background of noncytonuclear expression changes and to the diploid parents. Although we find subsets of cytonuclear genes in most lineages that match our expectations of maternal preference, these observations are not consistent among either allopolyploids or categories of organelle-targeted genes. Our results indicate that cytonuclear expression evolution may be subtle and variable among genera and genes, likely reflecting a diversity of mechanisms to resolve nuclear-cytoplasmic incompatibilities in allopolyploid species.

Morphological, Transcriptome, and Hormone Analysis of Dwarfism in Tetraploids of Populus alba × P. glandulosa

Breeding for dwarfism is an important approach to improve lodging resistance. Here, we performed comparative analysis of the phenotype, transcriptome, and hormone contents between diploids and tetraploids of poplar 84K (Populus alba × P. glandulosa). Compared with diploids, the indole-3-acetic acid (IAA) and gibberellin (GA₃) contents were increased, whereas the jasmonic acid (JA) and abscisic acid (ABA) contents were decreased in tetraploids. RNA-sequencing revealed that differentially expressed genes (DEGs) in leaves of tetraploids were mainly involved in plant hormone pathways. Most DEGs associated with IAA and GA promotion of plant growth and development were downregulated, whereas most DEGs associated with ABA and JA promotion of plant senescence were upregulated. Weighted gene co-expression network analysis indicated that certain transcription factors may be involved in the regulation of genes involved in plant hormone pathways. Thus, the altered expression of some genes in the plant hormone pathways may lead to a reduction in IAA and GA contents, as well as an elevation in ABA and JA contents, resulting in the dwarfing of tetraploids. The results show that polyploidization is a complex biological process affected by multiple plant hormone signals, and it provides a foundation for further exploration of the mechanism of tetraploids dwarfing in forest trees.

Induced polyploidy deeply influences reproductive life cycles, related phytochemical features, and phytohormonal activities in blackberry species

In some cases, polyploidy is an important phenomenon in the evolution of fruit crops. Polyploidy can be used in fruit breeding programs to develop varieties with higher yields and better fruit quality, as well as better adaptation to adverse environmental conditions. In this study, three wild species of blackberry were subjected to different degrees of induced polyploidy, and the effects of which were evaluated on morphological, physiological, and phytohormonal traits. With the aim of gaining a deep insight into the generative phase of plant growth and development, different levels of induced polyploidy were evaluated on the three blackberry species, i.e., Rubus persicus Bioss. (2x, 4x, and 8x), R. caesius L. (2x and 4x), and R. hirtus Schreb. (2x and 4x). The results showed that the polyploid plants performed significantly better than their diploid counterparts in terms of morphological traits such as flower count per spike and berry weight, as well as biochemical traits such as total soluble solids in the leaves. Induced polyploidy increased berry weight and drupe count per fruit. Microscopic examinations revealed a smaller number of viable pollen in the polyploids, compared to the diploids. Electron microscopy showed that the octaploid R. persicus had larger conical cells on the flower surface, compared to the diploid R. persicus. Correlation analysis showed that the ratio of indoleacetic acid to jasmonic acid changed synergistically with the total soluble solids in the leaves during the fruit set. The ploidy level correlated significantly with the number of pistils, leaf green index, total soluble solids in the leaves, and glucose content in floral nectar. Overall, induced polyploidy allowed Rubus to develop advantageous traits that can benefit future breeding programs and expand reproductive research in blackberries.

Repeat Dynamics across Timescales: A Perspective from Sibling Allotetraploid Marsh Orchids (Dactylorhiza majalis s.l.)

To provide insights into the fate of transposable elements (TEs) across timescales in a post-polyploidization context, we comparatively investigate five sibling Dactylorhiza allotetraploids (Orchidaceae) formed independently and sequentially between 500 and 100K generations ago by unidirectional hybridization between diploids D. fuchsii and D. incarnata. Our results first reveal that the paternal D. incarnata genome shows a marked increased content of LTR retrotransposons compared to the maternal species, reflected in its larger genome size and consistent with a previously hypothesized bottleneck. With regard to the allopolyploids, in the youngest D. purpurella both genome size and TE composition appear to be largely additive with respect to parents, whereas for polyploids of intermediate ages we uncover rampant genome expansion on a magnitude of multiple entire genomes of some plants such as Arabidopsis. The oldest allopolyploids in the series are not larger than the intermediate ones. A putative tandem repeat, potentially derived from a non-autonomous miniature inverted-repeat TE (MITE) drives much of the genome dynamics in the allopolyploids. The highly dynamic MITE-like element is found in higher proportions in the maternal diploid, D. fuchsii, but is observed to increase in copy number in both subgenomes of the allopolyploids. Altogether, the fate of repeats appears strongly regulated and therefore predictable across multiple independent allopolyploidization events in this system. Apart from the MITE-like element, we consistently document a mild genomic shock following the allopolyploidizations investigated here, which may be linked to their relatively large genome sizes, possibly associated with strong selection against further genome expansions.

Parental legacy versus regulatory innovation in salt stress responsiveness of allopolyploid cotton (Gossypium) species

Polyploidy provides an opportunity for evolutionary innovation and species diversification, especially under stressful conditions. In allopolyploids, the conditional dynamics of homoeologous gene expression can be either inherited from ancestral states pre-existing in the parental diploids or novel upon polyploidization, the latter potentially permitting a wider range of phenotypic responses to stresses. To gain insight into regulatory mechanisms underlying the diversity of salt resistance in Gossypium species, we compared global transcriptomic responses to modest salinity stress in two allotetraploid (AD-genome) cotton species, Gossypium hirsutum and G. mustelinum, relative to their model diploid progenitors (A-genome and D-genome). Multivariate and pairwise analyses of salt-responsive changes revealed a profound alteration of gene expression for about one third of the transcriptome. Transcriptional responses and associated functional implications of salt acclimation varied across species, as did species-specific coexpression modules among species and ploidy levels. Salt responsiveness in both allopolyploids was strongly biased toward the D-genome progenitor. A much lower level of transgressive downregulation was observed in the more salt-tolerant G. mustelinum than in the less tolerant G. hirsutum. By disentangling inherited effects from evolved responses, we show that expression biases that are not conditional upon salt stress approximately equally reflect parental legacy and regulatory novelty upon allopolyploidization, whereas stress-responsive biases are predominantly novel, or evolved, in allopolyploids. Overall, our work suggests that allopolyploid cottons acquired a wide range of stress response flexibility relative to their diploid ancestors, most likely mediated by complex suites of duplicated genes and regulatory factors.

Chromosome-level genome assembly and characterization of Sophora Japonica

Sophora japonica is a medium-size deciduous tree belonging to Leguminosae family and famous for its high ecological, economic and medicinal value. Here, we reveal a draft genome of S. japonica, which was ∼511.49 Mb long (contig N50 size of 17.34 Mb) based on Illumina, Nanopore and Hi-C data. We reliably assembled 110 contigs into 14 chromosomes, representing 91.62% of the total genome, with an improved N50 size of 31.32 Mb based on Hi-C data. Further investigation identified 271.76 Mb (53.13%) of repetitive sequences and 31,000 protein-coding genes, of which 30,721 (99.1%) were functionally annotated. Phylogenetic analysis indicates that S. japonica separated from Arabidopsis thaliana and Glycine max ∼107.53 and 61.24 million years ago, respectively. We detected evidence of species-specific and common-legume whole-genome duplication events in S. japonica. We further found that multiple TF families (e.g. BBX and PAL) have expanded in S. japonica, which might have led to its enhanced tolerance to abiotic stress. In addition, S. japonica harbours more genes involved in the lignin and cellulose biosynthesis pathways than the other two species. Finally, population genomic analyses revealed no obvious differentiation among geographical groups and the effective population size continuously declined since 2 Ma. Our genomic data provide a powerful comparative framework to study the adaptation, evolution and active ingredients biosynthesis in S. japonica. More importantly, our high-quality S. japonica genome is important for elucidating the biosynthesis of its main bioactive components, and improving its production and/or processing.

Expression Level Dominance and Homeolog Expression Bias Upon Cold Stress in the F1 Hybrid Between the Invasive Sphagneticola trilobata and the Native S. calendulacea in South China, and Implications for Its Invasiveness

The role of hybridization is significant in biological invasion, and thermotolerance is a trait critical to range expansions. The South American Sphagneticola trilobata is now widespread in South China, threatening the native S. calendulacea by competition and hybridization. Furthermore, upon formation, their F1 hybrid can quickly replace both parents. In this study, the three taxa were used as a model to investigate the consequences of hybridization on cold tolerance, particularly the effect of subgenome dominance in the hybrid. Upon chilling treatments, physiological responses and transcriptome profiles were compared across different temperature points to understand their differential responses to cold. While both parents showed divergent responses, the hybrid’s responses showed an overall resemblance to S. calendulacea, but the contribution of homeolog expression bias to cold stress was not readily evident in the F1 hybrid possibly due to inherent bias that comes with the sampling location. Our findings provided insights into the role of gene expression in differential cold tolerance, and further contribute to predicting the invasive potential of other hybrids between S. trilobata and its congeners around the world.

Comprehensive Evolutionary Analysis of CPP Genes in Brassica napus L. and Its Two Diploid Progenitors Revealing the Potential Molecular Basis of Allopolyploid Adaptive Advantage Under Salt Stress

Allopolyploids exist widely in nature and have strong environmental adaptability. The typical allopolyploid Brassica napus L. is a widely cultivated crop, but whether it is superior to its diploid progenitors in abiotic stress resistance and the key genes that may be involved are not fully understood. Cystein-rich polycomb-like protein (CPP) genes encode critical transcription factors involved in the response of abiotic stress, including salt stress. To explore the potential molecular basis of allopolyploid adaptation to salt stress, we comprehensively analyzed the characteristics and salt stress response of the CPP genes in B. napus and its two diploid progenitors in this study. We found some molecular basis that might be associated with the adaptability of B. napus, including the expansion of the CPP gene family, the acquisition of introns by some BnCPPs, and abundant cis-acting elements upstream of BnCPPs. We found two duplication modes (whole genome duplication and transposed duplication) might be the main reasons for the expansion of CPP gene family in B. napus during allopolyploidization. CPP gene expression levels and several physiological indexes were changed in B. napus and its diploid progenitors after salt stress, suggesting that CPP genes might play important roles in the response of salt stress. We found that some BnCPPs might undergo new functionalization or subfunctionalization, and some BnCPPs also show biased expression, which might contribute to the adaptation of B. napus under saline environment. Compared with diploid progenitors, B. napus showed stronger physiological responses, and BnCPP gene expression also showed higher changes after salt stress, indicating that the allopolyploid B. napus had an adaptive advantage under salt stress. This study could provide evidence for the adaptability of polyploid and provide important clues for the study of the molecular mechanism of salt stress resistance in B. napus.

Subgenome dominance and its evolutionary implications in crop domestication and breeding

Polyploidization or whole-genome duplication (WGD) is a well-known speciation and adaptation mechanism in angiosperms, while subgenome dominance is a crucial phenomenon in allopolyploids, established following polyploidization. The dominant subgenomes contribute more to genome evolution and homoeolog expression bias, both of which confer advantages for short-term phenotypic adaptation and long-term domestication. In this review, we firstly summarize the probable mechanistic basis for subgenome dominance, including the effects of genetic [transposon, genetic incompatibility, and homoeologous exchange (HE)], epigenetic (DNA methylation and histone modification), and developmental and environmental factors on this evolutionary process. We then move to Brassica rapa, a typical allopolyploid with subgenome dominance. Polyploidization provides the B. rapa genome not only with the genomic plasticity for adapting to changeable environments, but also an abundant genetic basis for morphological variation, making it a representative species for subgenome dominance studies. According to the 'two-step theory', B. rapa experienced genome fractionation twice during WGD, in which most of the genes responding to the environmental cues and phytohormones were over-retained, enhancing subgenome dominance and consequent adaption. More than this, the pangenome of 18 B. rapa accessions with different morphotypes recently constructed provides further evidence to reveal the impacts of polyploidization and subgenome dominance on intraspecific diversification in B. rapa. Above and beyond the fundamental understanding of WGD and subgenome dominance in B. rapa and other plants, however, it remains elusive why subgenome dominance has tissue- and spatiotemporal-specific features and could shuffle between homoeologous regions of different subgenomes by environments in allopolyploids. We lastly propose acceleration of the combined application of resynthesized allopolyploids, omics technology, and genome editing tools to deepen mechanistic investigations of subgenome dominance, both genetic and epigenetic, in a variety of species and environments. We believe that the implications of genomic and genetic basis of a variety of ecologically, evolutionarily, and agriculturally interesting traits coupled with subgenome dominance will be uncovered and aid in making new discoveries and crop breeding.

Transcriptome shock in interspecific F1 allotriploid hybrids between Brassica species

Interspecific hybridization drives the evolution of angiosperms and can be used to introduce novel alleles for important traits or to activate heterosis in crop breeding. Hybridization brings together gene expression networks from two different species, potentially causing global alterations of gene expression in the F1 plants which is called 'transcriptome shock'. Here, we explored such a transcriptome shock in allotriploid Brassica hybrids. We generated interspecific F1 allotriploid hybrids between the allotetraploid species Brassica napus and three accessions of the diploid species Brassica rapa. RNA-seq of the F1 hybrids and the parental plants revealed that 26.34-30.89% of genes were differentially expressed between the parents. We also analyzed expression level dominance and homoeolog expression bias between the parents and the F1 hybrids. The expression-level dominance biases of the Ar, An, and Cn subgenomes was genotype and stage dependent, whereas significant homoeolog expression bias was observed among three subgenomes from different parents. Furthermore, more genes were involved in trans regulation than in cis regulation in allotriploid F1 hybrids. Our findings provide new insights into the transcriptomic responses of cross-species hybrids and hybrids showing heterosis, as well as a new method for promoting the breeding of desirable traits in polyploid Brassica species.

Genome-Wide Identification of Direct Targets of ZjVND7 Reveals the Putative Roles of Whole-Genome Duplication in Sour Jujube in Regulating Xylem Vessel Differentiation and Drought Tolerance

The effects of whole-genome duplication span multiple levels. Previous study reported that the autotetraploid sour jujube exhibited superior drought tolerance than diploid. However, the difference in water transport system between diploids and autotetraploids and its mechanism remain unclear. Here, we found the number of xylem vessels and parenchyma cells in autotetraploid sour jujube increased to nearly twice that of diploid sour jujube, which may be closely related to the differences in xylem vessel differentiation-related ZjVND7 targets between the two ploidy types. Although the five enriched binding motifs are different, the most reliable motif in both diploid and autotetraploid sour jujube was CTTNAAG. Additionally, ZjVND7 targeted 236 and 321 genes in diploids and autotetraploids, respectively. More identified targeted genes of ZjVND7 were annotated to xylem development, secondary wall synthesis, cell death, cell division, and DNA endoreplication in autotetraploids than in diploids. SMR1 plays distinct roles in both proliferating and differentiated cells. Under drought stress, the binding signal of ZjVND7 to ZjSMR1 was stronger in autotetraploids than in diploids, and the fold-changes in the expression of ZjVND7 and ZjSMR1 were larger in the autotetraploids than in the diploids. These results suggested that the targeted regulation of ZjVND7 on ZjSMR1 may play valuable roles in autotetraploids in the response to drought stress. We hypothesized that the binding of ZjVND7 to ZjSMR1 might play a role in cell division and transdifferentiation from parenchyma cells to vessels in the xylem. This regulation could prolong the cell cycle and regulate endoreplication in response to drought stress and abscisic acid, which may be stronger in polyploids.

Interfered chromosome pairing at high temperature promotes meiotic instability in autotetraploid Arabidopsis

Changes in environmental temperature affect multiple meiotic processes in flowering plants. Polyploid plants derived from whole-genome duplication (WGD) have enhanced genetic plasticity and tolerance to environmental stress but face challenges in organizing and segregating doubled chromosome sets. In this study, we investigated the impact of increased environmental temperature on male meiosis in autotetraploid Arabidopsis (Arabidopsis thaliana). Under low to mildly increased temperatures (5°C-28°C), irregular chromosome segregation universally occurred in synthetic autotetraploid Columbia-0 (Col-0). Similar meiotic lesions occurred in autotetraploid rice (Oryza sativa L.) and allotetraploid canola (Brassica napus cv Westar), but not in evolutionarily derived hexaploid wheat (Triticum aestivum). At extremely high temperatures, chromosome separation and tetrad formation became severely disordered due to univalent formation caused by the suppression of crossing-over. We found a strong correlation between tetravalent formation and successful chromosome pairing, both of which were negatively correlated with temperature elevation, suggesting that increased temperature interferes with crossing-over predominantly by impacting homolog pairing. We also showed that loading irregularities of axis proteins ASY1 and ASY4 co-localize on the chromosomes of the syn1 mutant and the heat-stressed diploid and autotetraploid Col-0, revealing that heat stress affects the lateral region of synaptonemal complex (SC) by impacting the stability of the chromosome axis. Moreover, we showed that chromosome axis and SC in autotetraploid Col-0 are more sensitive to increased temperature than those in diploid Arabidopsis. Taken together, our data provide evidence suggesting that WGD negatively affects the stability and thermal tolerance of meiotic recombination in newly synthetic autotetraploid Arabidopsis.

Open chromatin interaction maps reveal functional regulatory elements and chromatin architecture variations during wheat evolution

BACKGROUND

Bread wheat (Triticum aestivum) is an allohexaploid that is generated by two subsequent allopolyploidization events. The large genome size (16 Gb) and polyploid complexity impede our understanding of how regulatory elements and their interactions shape chromatin structure and gene expression in wheat. The open chromatin enrichment and network Hi-C (OCEAN-C) is a powerful antibody-independent method to detect chromatin interactions between open chromatin regions throughout the genome.

RESULTS

Here we generate open chromatin interaction maps for hexaploid wheat and its tetraploid and diploid relatives using OCEAN-C. The anchors of chromatin loops show high chromatin accessibility and are concomitant with several active histone modifications, with 67% of them interacting with multiple loci. Binding motifs of various transcription factors are significantly enriched in the hubs of open chromatin interactions (HOCIs). The genes linked by HOCIs represent higher expression level and lower coefficient expression variance than the genes linked by other loops, which suggests HOCIs may coordinate co-expression of linked genes. Thousands of interchromosomal loops are identified, while limited interchromosomal loops (0.4%) are identified between homoeologous genes in hexaploid wheat. Moreover, we find structure variations contribute to chromatin interaction divergence of homoeologs and chromatin topology changes between different wheat species. The genes with discrepant chromatin interactions show expression alteration in hexaploid wheat compared with its tetraploid and diploid relatives.

CONCLUSIONS

Our results reveal open chromatin interactions in different wheat species, which provide new insights into the role of open chromatin interactions in gene expression during the evolution of polyploid wheat.

The activation of gene expression and alternative splicing in the formation and evolution of allopolyploid Brassica napus

Allopolyploids contain two or more sets of subgenomes. To establish a compatible relationship between subgenomes, a series of gene expression changes occurred in allopolyploids. What evolutionary changes of transcripts have taken place in Brassica napus during the early establishment and subsequent evolution was a fascinating scientific question. Here, we study this issue using a set of materials (natural, resynthesized B. napus and their progenitors/parents) by long-read RNA sequencing technology. The results showed that more genes were up-regulated in resynthesized B. napus compared with its two parents, and more up-regulated expressed genes were observed in natural B. napus compared with resynthesized B. napus. The presence of up-regulation genes in organism may help it adapt to the influence of "genomic shock" and cope with natural environment. Isoforms are produced from precursor mRNAs by alternative splicing (AS) events, and more than 60% of novel isoforms were identified in all materials, which could improve the reference genome information of B. napus. We found that the isoform numbers, the number of genes potentially involved in AS and alternative polyadenylation increased in B. napus after evolution, which may involve in the adaptation of plants to natural environment. In addition, all identified isoforms were functional annotated by searching 7 databases. In general, this study could improve our overall understanding of the full-length transcriptome of B. napus, and help us recognize the significant gene expression changes and isoform abundance changes occurred in allopolyploid B. napus during evolution.

Genome-Wide Identification and Characterization of the RCI2 Gene Family in Allotetraploid Brassica napus Compared with Its Diploid Progenitors

Brassica napus and its diploid progenitors (B. rapa and B. oleracea) are suitable for studying the problems associated with polyploidization. As an important anti-stress protein, RCI2 proteins widely exist in various tissues of plants, and are crucial to plant growth, development, and stress response. In this study, the RCI2 gene family was comprehensively identified and analyzed, and 9, 9, and 24 RCI2 genes were identified in B. rapa, B. oleracea, and B. napus, respectively. Phylogenetic analysis showed that all of the identified RCI2 genes were divided into two groups, and further divided into three subgroups. Ka/Ks analysis showed that most of the identified RCI2 genes underwent a purifying selection after the duplication events. Moreover, gene structure analysis showed that the structure of RCI2 genes is largely conserved during polyploidization. The promoters of the RCI2 genes in B. napus contained more cis-acting elements, which were mainly involved in plant development and growth, plant hormone response, and stress responses. Thus, B. napus might have potential advantages in some biological aspects. In addition, the changes of RCI2 genes during polyploidization were also discussed from the aspects of gene number, gene structure, gene relative location, and gene expression, which can provide reference for future polyploidization analysis.

Genome-wide identification and characterization of the CLASP_N gene family in upland cotton (Gossypium hirsutum L.)

BACKGROUND

Cytoplasmic linker-associated proteins (CLASPs) are tubule proteins that can bind to microtubules and participate in regulating the structure and function of microtubules, which significantly affects the development and growth of plants. These proteins have been identified in Arabidopsis; however, little research has been performed in upland cotton.

METHODS

In this study, the whole genome of the CLASP_N family was analyzed to provide theoretical support for the function of this gene family in the development of upland cotton fiber. Bioinformatics was used to analyze the family characteristics of CLASP_N in upland cotton, such as member identification, sequence characteristics, conserved domain structure and coevolutionary relationships. Real-time fluorescent quantitative PCR (qRT-PCR) was used to clarify the expression pattern of the upland cotton CLASP_N gene family in cotton fiber.

RESULTS

At the genome-wide level, we identified 16 upland cotton CLASP_N genes. A chromosomal localization analysis revealed that these 16 genes were located on 13 chromosomes. The motif results showed that all CLASP_N proteins have the CLASP_N domain. Gene structure analysis showed that the structure and length of exons and introns were consistent in the subgroups. In the evolutionary analysis with other species, the gene family clearly diverged from the other species in the evolutionary process. A promoter sequence analysis showed that this gene family contains a large number of cis-acting elements related to a variety of plant hormones. qRT-PCR was used to clarify the expression pattern of the upland cotton CLASP_N gene family in cotton fiber and leaves, and Gh210800 was found to be highly expressed in the later stages of fiber development. The results of this study provide a foundation for further research on the molecular role of the CLASP_N genes in cotton fiber development.

Spatial and Temporal Transcriptomic Heredity and Asymmetry in an Artificially Constructed Allotetraploid Wheat (AADD)

Polyploidy, or whole-genome duplication (WGD), often induces dramatic changes in gene expression due to "transcriptome shock. " However, questions remain about how allopolyploidy (the merging of multiple nuclear genomes in the same nucleus) affects gene expression within and across multiple tissues and developmental stages during the initial foundation of allopolyploid plants. Here, we systematically investigated the immediate effect of allopolyploidy on gene expression variation in an artificial allopolyploidy system consisting of a constructed allotetraploid wheat (AADD genome, accession AT2) and its diploid progenitors Triticum urartu and Aegilops tauschii. We performed comprehensive RNA sequencing of 81 samples from different genotypes, tissues, and developmental stages. First, we found that intrinsic interspecific differences between the diploid parents played a major role in establishing the expression architecture of the allopolyploid. Nonetheless, allopolyploidy per se also induced dramatic and asymmetric patterns of differential gene expression between the subgenomes, and genes from the D subgenome exhibited a more drastic response. Second, analysis of homoeolog expression bias (HEB) revealed that the D subgenome exhibited significant expression bias and that de novo-generated HEB was attributed mainly to asymmetrical differential gene expression. Homoeolog-specific expression (HSE) analyses showed that the cis-only regulatory pattern was predominant in AT2, reflecting significant divergence between the parents. Co-expression network analysis revealed that homoeolog expression connectivity (HEC) was significantly correlated with sequence divergence in cis elements between subgenomes. Interestingly, allopolyploidy-induced reconstruction of network modules was also associated with different HSE patterns. Finally, a transcriptome atlas of spike development demonstrated that the phenotypic similarity of AT2 to T. urartu may be attributed to the combination of relatively stable expression of A-subgenome genes and drastic downregulation of their D-subgenome homoeologs. These findings provide a broad, multidimensional characterization of allopolyploidy-induced transcriptomic responses and suggest that allopolyploidy can have immediate and complex regulatory effects on the expression of nuclear genes.

Impact of polyploidy on plant tolerance to abiotic and biotic stresses

Polyploidy, defined as the coexistence of three or more complete sets of chromosomes in an organism's cells, is considered as a pivotal moving force in the evolutionary history of vascular plants and has played a major role in the domestication of several crops. In the last decades, improved cultivars of economically important species have been developed artificially by inducing autopolyploidy with chemical agents. Studies on diverse species have shown that the anatomical and physiological changes generated by either natural or artificial polyploidization can increase tolerance to abiotic and biotic stresses as well as disease resistance, which may positively impact on plant growth and net production. The aim of this work is to review the current literature regarding the link between plant ploidy level and tolerance to abiotic and biotic stressors, with an emphasis on the physiological and molecular mechanisms responsible for these effects, as well as their impact on the growth and development of both natural and artificially generated polyploids, during exposure to adverse environmental conditions. We focused on the analysis of those types of stressors in which more progress has been made in the knowledge of the putative morpho-physiological and/or molecular mechanisms involved, revealing both the factors in common, as well as those that need to be addressed in future research.

Genome-wide characterization and analysis of Golden 2-Like transcription factors related to leaf chlorophyll synthesis in diploid and triploid Eucalyptus urophylla

Golden 2-Like (GLK) transcription factors play a crucial role in chloroplast development and chlorophyll synthesis in many plant taxa. To date, no systematic analysis of GLK transcription factors in tree species has been conducted. In this study, 40 EgrGLK genes in the Eucalyptus grandis genome were identified and divided into seven groups based on the gene structure and motif composition. The EgrGLK genes were mapped to 11 chromosomes and the distribution of genes on chromosome was uneven. Phylogenetic analysis of GLK proteins between E. grandis and other species provided information for the high evolutionary conservation of GLK genes among different species. Prediction of cis-regulatory elements indicated that the EgrGLK genes were involved in development, light response, and hormone response. Based on the finding that the content of chlorophyll in mature leaves was the highest, and leaf chlorophyll content of triploid Eucalyptus urophylla was higher than that of the diploid control, EgrGLK expression pattern in leaves of triploid and diploid E. urophylla was examined by means of transcriptome analysis. Differential expression of EgrGLK genes in leaves of E. urophylla of different ploidies was consistent with the trend in chlorophyll content. To further explore the relationship between EgrGLK expression and chlorophyll synthesis, co-expression networks were generated, which indicated that EgrGLK genes may have a positive regulatory relationship with chlorophyll synthesis. In addition, three EgrGLK genes that may play an important role in chlorophyll synthesis were identified in the co-expression networks. And the prediction of miRNAs targeting EgrGLK genes showed that miRNAs might play an important role in the regulation of EgrGLK gene expression. This research provides valuable information for further functional characterization of GLK genes in Eucalyptus.

Genome-wide unbalanced expression bias and expression level dominance toward Brassica oleracea in artificially synthesized intergeneric hybrids of Raphanobrassica

Raphanobrassica (RrRrCrCr, 2n = 4x = 36), which is generated by distant hybridization between the maternal parent Raphanus sativus (RsRs, 2n = 2x = 18) and the paternal parent Brassica oleracea (C°C°, 2n = 2x = 18), displays intermediate silique phenotypes compared to diploid progenitors. However, the hybrid shares much more similarities in silique phenotypes with those of B. oleracea than those of R. sativus. Strikingly, the silique of Raphanobrassica is obviously split into two parts. To investigate the gene expression patterns behind these phenomena, transcriptome analysis was performed on the upper, middle, and lower sections of pods (RCsiu, RCsim, and RCsil), seeds in the upper and lower sections of siliques (RCseu and RCsel) from Raphanobrassica, whole pods (Rsi and Csi) and all seeds in the siliques (Rse and Cse) from R. sativus and B. oleracea. Transcriptome shock was observed in all five aforementioned tissues of Raphanobrassica. Genome-wide unbalanced biased expression and expression level dominance were also discovered, and both of them were toward B. oleracea in Raphanobrassica, which is consistent with the observed phenotypes. The present results reveal the global gene expression patterns of different sections of siliques of Raphanobrassica, pods, and seeds of B. oleracea and R. sativus, unraveling the tight correlation between global gene expression patterns and phenotypes of the hybrid and its parents.

Prospects and challenges of epigenomics in crop improvement

BACKGROUND

The advent of high-throughput epigenome mapping techniques has ushered in a new era of multiomics with powerful tools now available to map and record genomic output at different levels. Integrating the different components of the epigenome from these multiomics measures allows investigations of cis-regulatory elements on a genome-scale. Mapping of chromatin state, chromatin accessibility dynamics, and higher-order chromatin structure enables a new level of understanding of cell fate determination, identity and function in normal growth and development, disease resistance, and yield.

OBJECTIVE

In this paper, the recent advances in epigenomics research of rice, maize, and wheat are reviewed, and the development trends of epigenomics of major crops in the coming years are projected.

METHODS

We highlight the role of epigenomics in regulating growth and development and identifying potential distal cis-regulatory elements in three major crops, and discuss the prospects and challenges for new epigenetics-mediated breeding technologies in crop improvement.

CONCLUSION

In this review, we summarize and analyze recent epigenomic advances in three major crops epigenomics and discuss possibilities and challenges for future research in the field.