Differentiation of the maize subgenomes by genome dominance and both ancient and ongoing gene loss

Systematic approach of polyploidy as an evolutionary genetic and genomic phenomenon in horticultural crops

Polyploidy is thought to be an evolutionary and systematic mechanism for gene flow and phenotypic advancement in flowering plants. It is a natural phenomenon that promotes diversity by creating new permutations enhancing the prime potentials as compared to progenitors. Two different pathways have been recognized in studying polyploidy in nature; mitotic or somatic chromosome doubling and cytogenetics variation. Secondly, the vital influence of being polyploid is its heritable property (unreduced reproductive cells) formed during first and second-division restitution (FDR & SDR). Different approaches either chemical (Colchicine, Oryzalin, Caffeine, Trifuralin, or phosphoric amides) or gaseous i.e. Nitrous oxide have been deliberated as strong polyploidy causing agents. A wide range of cytogenetic practices like chromosomes study, ploidy, genome analysis, and plant morphology and anatomy have been studied in different plant species. Flow cytometry for ploidy and chromosome analysis through fluorescence and genomic in situ hybridization (FISH & GISH) are the basic methods to evaluate heredity substances sampled from leaves and roots. Many horticultural crops have been developed successfully and released commercially for consumption. Moreover, some deep detailed studies are needed to check the strong relationship between unique morphological features and genetic makeup concerning genes and hormonal expression in a strong approach.

Genome of Halimeda opuntia reveals differentiation of subgenomes and molecular bases of multinucleation and calcification in algae

Algae mostly occur either as unicellular (microalgae) or multicellular (macroalgae) species, both being uninucleate. There are important exceptions, however, as some unicellular algae are multinucleate and macroscopic, some of which inhabit tropical seas and contribute to biocalcification and coral reef robustness. The evolutionary mechanisms and ecological significance of multinucleation and associated traits (e.g., rapid wound healing) are poorly understood. Here, we report the genome of Halimeda opuntia, a giant multinucleate unicellular chlorophyte characterized by interutricular calcification. We achieve a high-quality genome assembly that shows segregation into four subgenomes, with evidence for polyploidization concomitant with historical sea level and climate changes. We further find myosin VIII missing in H. opuntia and three other unicellular multinucleate chlorophytes, suggesting a potential mechanism that may underpin multinucleation. Genome analysis provides clues about how the unicellular alga could survive fragmentation and regenerate, as well as potential signatures for extracellular calcification and the coupling of calcification with photosynthesis. In addition, proteomic alkalinity shifts were found to potentially confer plasticity of H. opuntia to ocean acidification (OA). Our study provides crucial genetic information necessary for understanding multinucleation, cell regeneration, plasticity to OA, and different modes of calcification in algae and other organisms, which has important implications in reef conservation and bioengineering.

The Zea mays PeptideAtlas: A New Maize Community Resource

This study presents the Maize PeptideAtlas resource (www.peptideatlas.org/builds/maize) to help solve questions about the maize proteome. Publicly available raw tandem mass spectrometry (MS/MS) data for maize collected from ProteomeXchange were reanalyzed through a uniform processing and metadata annotation pipeline. These data are from a wide range of genetic backgrounds and many sample types and experimental conditions. The protein search space included different maize genome annotations for the B73 inbred line from MaizeGDB, UniProtKB, NCBI RefSeq, and for the W22 inbred line. 445 million MS/MS spectra were searched, of which 120 million were matched to 0.37 million distinct peptides. Peptides were matched to 66.2% of proteins in the most recent B73 nuclear genome annotation. Furthermore, most conserved plastid- and mitochondrial-encoded proteins (NCBI RefSeq annotations) were identified. Peptides and proteins identified in the other B73 genome annotations will improve maize genome annotation. We also illustrate the high-confidence detection of unique W22 proteins. N-terminal acetylation, phosphorylation, ubiquitination, and three lysine acylations (K-acetyl, K-malonyl, and K-hydroxyisobutyryl) were identified and can be inspected through a PTM viewer in PeptideAtlas. All matched MS/MS-derived peptide data are linked to spectral, technical, and biological metadata. This new PeptideAtlas is integrated in MaizeGDB with a peptide track in JBrowse.

Divergent Roles of the Auxin Response Factors in Lemongrass (Cymbopogon flexuosus (Nees ex Steud.) W. Watson) during Plant Growth

Auxin Response Factors (ARFs) make up a plant-specific transcription factor family that mainly couples perception of the phytohormone, auxin, and gene expression programs and plays an important and multi-faceted role during plant growth and development. Lemongrass (Cymbopogon flexuosus) is a representative Cymbopogon species widely used in gardening, beverages, fragrances, traditional medicine, and heavy metal phytoremediation. Biomass yield is an important trait for several agro-economic purposes of lemongrass, such as landscaping, essential oil production, and phytoremediation. Therefore, we performed gene mining of CfARFs and identified 26 and 27 CfARF-encoding genes in each of the haplotype genomes of lemongrass, respectively. Phylogenetic and domain architecture analyses showed that CfARFs can be divided into four groups, among which groups 1, 2, and 3 correspond to activator, repressor, and ETTN-like ARFs, respectively. To identify the CfARFs that may play major roles during the growth of lemongrass plants, RNA-seq was performed on three tissues (leaf, stem, and root) and four developmental stages (3-leaf, 4-leaf, 5-leaf. and mature stages). The expression profiling of CfARFs identified several highly expressed activator and repressor CfARFs and three CfARFs (CfARF3, 18, and 35) with gradually increased levels during leaf growth. Haplotype-resolved transcriptome analysis revealed that biallelic expression dominance is frequent among CfARFs and contributes to their gene expression patterns. In addition, co-expression network analysis identified the modules enriched with CfARFs. By establishing orthologous relationships among CfARFs, sorghum ARFs, and maize ARFs, we showed that CfARFs were mainly expanded by whole-genome duplications, and that the duplicated CfARFs might have been divergent due to differential expression and variations in domains and motifs. Our work provides a detailed catalog of CfARFs in lemongrass, representing a first step toward characterizing CfARF functions, and may be useful in molecular breeding to enhance lemongrass plant growth.

Prioritizing Maize Metabolic Gene Regulators through Multi-Omic Network Integration

Elucidating gene regulatory networks is a major area of study within plant systems biology. Phenotypic traits are intricately linked to specific gene expression profiles. These expression patterns arise primarily from regulatory connections between sets of transcription factors (TFs) and their target genes. Here, we integrated 46 co-expression networks, 283 protein-DNA interaction (PDI) assays, and 16 million SNPs used to identify expression quantitative trait loci (eQTL) to construct TF-target networks. In total, we analyzed ∼4.6M interactions to generate four distinct types of TF-target networks: co-expression, PDI, trans -eQTL, and cis -eQTL combined with PDIs. To functionally annotate TFs based on their target genes, we implemented three different network integration strategies. We evaluated the effectiveness of each strategy through TF loss-of function mutant inspection and random network analyses. The multi-network integration allowed us to identify transcriptional regulators of several biological processes. Using the topological properties of the fully integrated network, we identified potential functionally redundant TF paralogs. Our findings retrieved functions previously documented for numerous TFs and revealed novel functions that are crucial for informing the design of future experiments. The approach here-described lays the foundation for the integration of multi-omic datasets in maize and other plant systems.

GRAPHICAL ABSTRACT

Homoeologs in Allopolyploids: Navigating Redundancy as Both an Evolutionary Opportunity and a Technical Challenge-A Transcriptomics Perspective

Allopolyploidy in plants involves the merging of two or more distinct parental genomes into a single nucleus, a significant evolutionary process in the plant kingdom. Transcriptomic analysis provides invaluable insights into allopolyploid plants by elucidating the fate of duplicated genes, revealing evolutionary novelties and uncovering their environmental adaptations. By examining gene expression profiles, scientists can discern how duplicated genes have evolved to acquire new functions or regulatory roles. This process often leads to the development of novel traits and adaptive strategies that allopolyploid plants leverage to thrive in diverse ecological niches. Understanding these molecular mechanisms not only enhances our appreciation of the genetic complexity underlying allopolyploidy but also underscores their importance in agriculture and ecosystem resilience. However, transcriptome profiling is challenging due to genomic redundancy, which is further complicated by the presence of multiple chromosomes sets and the variations among homoeologs and allelic genes. Prior to transcriptome analysis, sub-genome phasing and homoeology inference are essential for obtaining a comprehensive view of gene expression. This review aims to clarify the terminology in this field, identify the most challenging aspects of transcriptome analysis, explain their inherent difficulties, and suggest reliable analytic strategies. Furthermore, bulk RNA-seq is highlighted as a primary method for studying allopolyploid gene expression, focusing on critical steps like read mapping and normalization in differential gene expression analysis. This approach effectively captures gene expression from both parental genomes, facilitating a comprehensive analysis of their combined profiles. Its sensitivity in detecting low-abundance transcripts allows for subtle differences between parental genomes to be identified, crucial for understanding regulatory dynamics and gene expression balance in allopolyploids.

Evolution and subfunctionalization of CIPK6 homologous genes in regulating cotton drought resistance

The occurrence of whole-genome duplication or polyploidy may promote plant adaptability to harsh environments. Here, we clarify the evolutionary relationship of eight GhCIPK6 homologous genes in upland cotton (Gossypium hirsutum). Gene expression and interaction analyses indicate that GhCIPK6 homologous genes show significant functional changes after polyploidy. Among these, GhCIPK6D1 and GhCIPK6D3 are significantly up-regulated by drought stress. Functional studies reveal that high GhCIPK6D1 expression promotes cotton drought sensitivity, while GhCIPK6D3 expression promotes drought tolerance, indicating clear functional differentiation. Genetic and biochemical analyses confirm the synergistic negative and positive regulation of cotton drought resistance through GhCBL1A1-GhCIPK6D1 and GhCBL2A1-GhCIPK6D3, respectively, to regulate stomatal movement by controlling the directional flow of K+ in guard cells. These results reveal differentiated roles of GhCIPK6 homologous genes in response to drought stress in upland cotton following polyploidy. The work provides a different perspective for exploring the functionalization and subfunctionalization of duplicated genes in response to polyploidization.

Doubled haploid technology and synthetic apomixis: Recent advances and applications in future crop breeding

Doubled haploid (DH) technology and synthetic apomixis approaches can considerably shorten breeding cycles and enhance breeding efficiency. Compared with traditional breeding methods, DH technology offers the advantage of rapidly generating inbred lines, while synthetic apomixis can effectively fix hybrid vigor. In this review, we focus on (i) recent advances in identifying and characterizing genes responsible for haploid induction (HI), (ii) the molecular mechanisms of HI, (iii) spontaneous haploid genome doubling, and (iv) crop synthetic apomixis. We also discuss the challenges and potential solutions for future crop breeding programs utilizing DH technology and synthetic apomixis. Finally, we provide our perspectives about how to integrate DH and synthetic apomixis for precision breeding and de novo domestication.

Allopolyploid subgenome identification and implications for evolutionary analysis

Whole-genome duplications (WGDs) are widespread genomic events in eukaryotes that are hypothesized to contribute to the evolutionary success of many lineages, including flowering plants, Saccharomyces yeast, and vertebrates. WGDs generally can be classified into autopolyploids (ploidy increase descended from one species) or allopolyploids (ploidy increase descended from multiple species). Assignment of allopolyploid progenitor species (called subgenomes in the polyploid) is important to understanding the biology and evolution of polyploids, including the asymmetric subgenome evolution following hybridization (biased fractionation). Here, I review the different methodologies used to identify the ancestors of allopolyploid subgenomes, discuss the advantages and disadvantages of these methods, and outline the implications of how these methods affect the subsequent evolutionary analysis of these genomes.

Genetic Diversity and Population Structure of a Large USDA Sesame Collection

Sesame, Sesamum indicum L., is one of the oldest domesticated crops used for its oil and protein in many parts of the world. To build genomic resources for sesame that could be used to improve sesame productivity and responses to stresses, a USDA sesame germplasm collection of 501 accessions originating from 36 countries was used in this study. The panel was genotyped using genotyping-by-sequencing (GBS) technology to explore its genetic diversity and population structure and the relatedness among its accessions. A total of 24,735 high-quality single-nucleotide polymorphism (SNP) markers were identified over the 13 chromosomes. The marker density was 1900 SNP per chromosome, with an average polymorphism information content (PIC) value of 0.267. The marker polymorphisms and heterozygosity estimators indicated the usefulness of the identified SNPs to be used in future genetic studies and breeding activities. The population structure, principal components analysis (PCA), and unrooted neighbor-joining phylogenetic tree analyses classified two distinct subpopulations, indicating a wide genetic diversity within the USDA sesame collection. Analysis of molecular variance (AMOVA) revealed that 29.5% of the variation in this population was due to subpopulations, while 57.5% of the variation was due to variation among the accessions within the subpopulations. These results showed the degree of differentiation between the two subpopulations as well as within each subpopulation. The high fixation index (FST) between the distinguished subpopulations indicates a wide genetic diversity and high genetic differentiation among and within the identified subpopulations. The linkage disequilibrium (LD) pattern averaged 161 Kbp for the whole sesame genome, while the LD decay ranged from 168 Kbp at chromosome LG09 to 123 Kbp in chromosome LG05. These findings could explain the complications of linkage drag among the traits during selections. The selected accessions and genotyped SNPs provide tools to enhance genetic gain in sesame breeding programs through molecular approaches.

The ancestral karyotype of the Heliantheae Alliance, herbicide resistance, and human allergens: Insights from the genomes of common and giant ragweed

Ambrosia artemisiifolia and Ambrosia trifida (Asteraceae) are important pest species and the two greatest sources of aeroallergens globally. Here, we took advantage of a hybrid to simplify genome assembly and present chromosome-level assemblies for both species. These assemblies show high levels of completeness with Benchmarking Universal Single-Copy Ortholog (BUSCO) scores of 94.5% for A. artemisiifolia and 96.1% for A. trifida and long terminal repeat (LTR) Assembly Index values of 26.6 and 23.6, respectively. The genomes were annotated using RNA data identifying 41,642 genes in A. artemisiifolia and 50,203 in A. trifida. More than half of the genome is composed of repetitive elements, with 62% in A. artemisiifolia and 69% in A. trifida. Single copies of herbicide resistance-associated genes PPX2L, HPPD, and ALS were found, while two copies of the EPSPS gene were identified; this latter observation may reveal a possible mechanism of resistance to the herbicide glyphosate. Ten of the 12 main allergenicity genes were also localized, some forming clusters with several copies, especially in A. artemisiifolia. The evolution of genome structure has differed among these two species. The genome of A. trifida has undergone greater rearrangement, possibly the result of chromoplexy. In contrast, the genome of A. artemisiifolia retains a structure that makes the allotetraploidization of the most recent common ancestor of the Heliantheae Alliance the clearest feature of its genome. When compared to other Heliantheae Alliance species, this allowed us to reconstruct the common ancestor's karyotype-a key step for furthering of our understanding of the evolution and diversification of this economically and allergenically important group.

The additive function of YIGE2 and YIGE1 in regulating maize ear length

Ear length (EL) is a key trait that greatly contributes to yield in maize. Although dozens of EL quantitative trait loci have been mapped, very few causal genes have been cloned, and the molecular mechanisms remain largely unknown. Our previous study showed that YIGE1 is involved in sugar and auxin pathways to regulate ear inflorescence meristem (IM) development and thus affects EL in maize. Here, we reveal that YIGE2, the paralog of YIGE1, regulates maize ear development and EL through auxin pathway. Knockout of YIGE2 causes a significant decrease of auxin level, IM length, floret number, EL, and grain yield. yige1 yige2 double mutants had even shorter IM and ears implying that these two genes redundantly regulate IM development and EL. The genes controlling auxin levels are differential expressed in yige1 yige2 double mutants, leading to lower auxin level. These results elucidated the critical role of YIGE2 and the redundancy between YIGE2 and YIGE1 in maize ear development, providing a new genetic resource for maize yield improvement.

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

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

Genome organization and botanical diversity

The rich diversity of angiosperms, both the planet's dominant flora and the cornerstone of agriculture, is integrally intertwined with a distinctive evolutionary history. Here, we explore the interplay between angiosperm genome organization and botanical diversity, empowered by genomic approaches ranging from genetic linkage mapping to analysis of gene regulation. Commonality in the genetic hardware of plants has enabled robust comparative genomics that has provided a broad picture of angiosperm evolution and implicated both general processes and specific elements in contributing to botanical diversity. We argue that the hardware of plant genomes-both in content and in dynamics-has been shaped by selection for rather substantial differences in gene regulation between plants and animals such as maize and human, organisms of comparable genome size and gene number. Their distinctive genome content and dynamics may reflect in part the indeterminate development of plants that puts strikingly different demands on gene regulation than in animals. Repeated polyploidization of plant genomes and multiplication of individual genes together with extensive rearrangement and differential retention provide rich raw material for selection of morphological and/or physiological variations conferring fitness in specific niches, whether natural or artificial. These findings exemplify the burgeoning information available to employ in increasing knowledge of plant biology and in modifying selected plants to better meet human needs.

The genome and population genomics of allopolyploid Coffea arabica reveal the diversification history of modern coffee cultivars

Coffea arabica, an allotetraploid hybrid of Coffea eugenioides and Coffea canephora, is the source of approximately 60% of coffee products worldwide, and its cultivated accessions have undergone several population bottlenecks. We present chromosome-level assemblies of a di-haploid C. arabica accession and modern representatives of its diploid progenitors, C. eugenioides and C. canephora. The three species exhibit largely conserved genome structures between diploid parents and descendant subgenomes, with no obvious global subgenome dominance. We find evidence for a founding polyploidy event 350,000-610,000 years ago, followed by several pre-domestication bottlenecks, resulting in narrow genetic variation. A split between wild accessions and cultivar progenitors occurred ~30.5 thousand years ago, followed by a period of migration between the two populations. Analysis of modern varieties, including lines historically introgressed with C. canephora, highlights their breeding histories and loci that may contribute to pathogen resistance, laying the groundwork for future genomics-based breeding of C. arabica.

Dynamics of accessible chromatin regions and subgenome dominance in octoploid strawberry

Subgenome dominance has been reported in diverse allopolyploid species, where genes from one subgenome are preferentially retained and are more highly expressed than those from other subgenome(s). However, the molecular mechanisms responsible for subgenome dominance remain poorly understood. Here, we develop genome-wide map of accessible chromatin regions (ACRs) in cultivated strawberry (2n = 8x = 56, with A, B, C, D subgenomes). Each ACR is identified as an MNase hypersensitive site (MHS). We discover that the dominant subgenome A contains a greater number of total MHSs and MHS per gene than the submissive B/C/D subgenomes. Subgenome A suffers fewer losses of MHS-related DNA sequences and fewer MHS fragmentations caused by insertions of transposable elements. We also discover that genes and MHSs related to stress response have been preferentially retained in subgenome A. We conclude that preservation of genes and their cognate ACRs, especially those related to stress responses, play a major role in the establishment of subgenome dominance in octoploid strawberry.

Cytonuclear Interactions and Subgenome Dominance Shape the Evolution of Organelle-Targeted Genes in the Brassica Triangle of U

The interaction and coevolution between nuclear and cytoplasmic genomes are one of the fundamental hallmarks of eukaryotic genome evolution and, 2 billion yr later, are still major contributors to the formation of new species. Although many studies have investigated the role of cytonuclear interactions following allopolyploidization, the relative magnitude of the effect of subgenome dominance versus cytonuclear interaction on genome evolution remains unclear. The Brassica triangle of U features 3 diploid species that together have formed 3 separate allotetraploid species on similar evolutionary timescales, providing an ideal system for understanding the contribution of the cytoplasmic donor to hybrid polyploid. Here, we investigated the evolutionary pattern of organelle-targeted genes in Brassica carinata (BBCC) and 2 varieties of Brassica juncea (AABB) at the whole-genome level, with particular focus on cytonuclear enzyme complexes. We found partial evidence that plastid-targeted genes experience selection to match plastid genomes, but no obvious corresponding signal in mitochondria-targeted genes from these 2 separately formed allopolyploids. Interestingly, selection acting on plastid genomes always reduced the retention rate of plastid-targeted genes encoded by the B subgenome, regardless of whether the Brassica nigra (BB) subgenome was contributed by the paternal or maternal progenitor. More broadly, this study illustrates the distinct selective pressures experienced by plastid- and mitochondria-targeted genes, despite a shared pattern of inheritance and natural history. Our study also highlights an important role for subgenome dominance in allopolyploid genome evolution, even in genes whose function depends on separately inherited molecules.

New whole-genome alignment tools are needed for tapping into plant diversity

Genome alignment is one of the most foundational methods for genome sequence studies. With rapid advances in sequencing and assembly technologies, these newly assembled genomes present challenges for alignment tools to meet the increased complexity and scale. Plant genome alignment is technologically challenging because of frequent whole-genome duplications (WGDs) as well as chromosome rearrangements and fractionation, high nucleotide diversity, widespread structural variation, and high transposable element (TE) activity causing large proportions of repeat elements. We summarize classical pairwise and multiple genome alignment (MGA) methods, and highlight techniques that are widely used or are being developed by the plant research community. We also outline the remaining challenges for precise genome alignment and the interpretation of alignment results in plants.

Identification and characterization of the Remorin gene family in Saccharum and the involvement of ScREM1.5e-1/-2 in SCMV infection on sugarcane

Introduction

Remorins (REMs) are plant-specific membrane-associated proteins that play important roles in plant-pathogen interactions and environmental adaptations. Group I REMs are extensively involved in virus infection. However, little is known about the REM gene family in sugarcane (Saccharum spp. hyrid), the most important sugar and energy crop around world.

Methods

Comparative genomics were employed to analyze the REM gene family in Saccharum spontaneum. Transcriptomics or RT-qPCR were used to analyze their expression files in different development stages or tissues under different treatments. Yeast two hybrid, bimolecular fluorescence complementation and co-immunoprecipitation assays were applied to investigate the protein interaction.

Results

In this study, 65 REMs were identified from Saccharum spontaneum genome and classified into six groups based on phylogenetic tree analysis. These REMs contain multiple cis-elements associated with growth, development, hormone and stress response. Expression profiling revealed that among different SsREMs with variable expression levels in different developmental stages or different tissues. A pair of alleles, ScREM1.5e-1/-2, were isolated from the sugarcane cultivar ROC22. ScREM1.5e-1/-2 were highly expressed in leaves, with the former expressed at significantly higher levels than the latter. Their expression was induced by treatment with H2O2, ABA, ethylene, brassinosteroid, SA or MeJA, and varied upon Sugarcane mosaic virus (SCMV) infection. ScREM1.5e-1 was localized to the plasma membrane (PM), while ScREM1.5e-2 was localized to the cytoplasm or nucleus. ScREM1.5e-1/-2 can self-interact and interact with each other, and interact with VPgs from SCMV, Sorghum mosaic virus, or Sugarcane streak mosaic virus. The interactions with VPgs relocated ScREM1.5e-1 from the PM to the cytoplasm.

Discussion

These results reveal the origin, distribution and evolution of the REM gene family in sugarcane and may shed light on engineering sugarcane resistance against sugarcane mosaic pathogens.

Polyploid speciation in Zea (Poaceae): cytogenetic insights

MAIN CONCLUSION

The analysis of meiotic pairing affinities and genomic formulae in species and hybrids of Zea allowed us to speculate an evolutionary model to recreate the ancient polyploidization of maize and allied species. The meiotic pairing affinities and the genomic formulae analysis in Zea species and hybrids obtained in new and previous crosses, together with the molecular data known in the genus, allowed us to speculate an evolutionary model to attempt to recreate the ancient polyploidization process of Zea species. We propose that x = 5 semispecies are the ancestors of all modern species of the genus. The complex evolutionary process that originated the different taxa could be included hybridization between sympatric diploid ancestral semispecies (2n = 10) and recurrent duplication of the hybrid chromosome number, resulting in distinct auto- and allopolyploids. After the merger and doubling of independent genomes would have undergone cytological and genetical diploidization, implying revolutionary changes in genome organization and genic balance processes. Based on the meiotic behaviour of the 2n = 30 hybrids, that showed homoeology between the A subgenomes of all parental species, we propose that this subgenome A would be pivotal in all the species and would have conserved the rDNA sequences and the pairing regulator locus (PrZ). In the hypothetical model postulated here, the ancestral semispecies with the pivotal subgenome A would have had a wide geographic distribution, co-occurring and hybridizing with the semispecies harbouring B subgenomes, thus enabling sympatric speciation.

Chromosomal-scale genomes of two Rosa species provide insights into genome evolution and ascorbate accumulation

Rosa roxburghii and Rosa sterilis, two species belonging to the Rosaceae family, are widespread in the southwest of China. These species have gained recognition for their remarkable abundance of ascorbate in their fresh fruits, making them an ideal vitamin C resource. In this study, we generated two high-quality chromosome-scale genome assemblies for R. roxburghii and R. sterilis, with genome sizes of 504 and 981.2 Mb, respectively. Notably, we present a haplotype-resolved, chromosome-scale assembly for diploid R. sterilis. Our results indicated that R. sterilis originated from the hybridization of R. roxburghii and R. longicuspis. Genome analysis revealed the absence of recent whole-genome duplications in both species and identified a series of duplicated genes that possibly contributing to the accumulation of flavonoids. We identified two genes in the ascorbate synthesis pathway, GGP and GalLDH, that show signs of positive selection, along with high expression levels of GDP-d-mannose 3', 5'-epimerase (GME) and GDP-l-galactose phosphorylase (GGP) during fruit development. Furthermore, through co-expression network analysis, we identified key hub genes (MYB5 and bZIP) that likely regulate genes in the ascorbate synthesis pathway, promoting ascorbate biosynthesis. Additionally, we observed the expansion of terpene synthase genes in these two species and tissue expression patterns, suggesting their involvement in terpenoid biosynthesis. Our research provides valuable insights into genome evolution and the molecular basis of the high concentration of ascorbate in these two Rosa species.

Genome dosage alteration caused by chromosome pyramiding and shuffling effects on karyotypic heterogeneity, reproductive diversity, and phenotypic variation in Zea-Tripsacum allopolyploids

KEY MESSAGE

We developed an array of Zea-Tripsacum tri-hybrid allopolyploids with multiple ploidies. We unveiled that changes in genome dosage due to the chromosomes pyramiding and shuffling of three species effects karyotypic heterogeneity, reproductive diversity, and phenotypic variation in Zea-Tripsacum allopolyploids. Polyploidy, or whole genome duplication, has played a major role in evolution and speciation. The genomic consequences of polyploidy have been extensively studied in many plants; however, the extent of chromosomal variation, genome dosage, phenotypic diversity, and heterosis in allopolyploids derived from multiple species remains largely unknown. To address this question, we synthesized an allohexaploid involving Zea mays, Tripsacum dactyloides, and Z. perennis by chromosomal pyramiding. Subsequently, an allooctoploid and an allopentaploid were obtained by hybridization of the allohexaploid with Z. perennis. Moreover, we constructed three populations with different ploidy by chromosomal shuffling (allopentaploid × Z. perennis, allohexaploid × Z. perennis, and allooctoploid × Z. perennis). We have observed 3 types of sexual reproductive modes and 2 types of asexual reproduction modes in the tri-species hybrids, including 2n gamete fusion (2n + n), haploid gamete fusion (n + n), polyspermy fertilization (n + n + n) or 2n gamete fusion (n + 2n), haploid gametophyte apomixis, and asexual reproduction. The tri-hybrids library presents extremely rich karyotype heterogeneity. Chromosomal compensation appears to exist between maize and Z. perennis. A rise in the ploidy of the trihybrids was linked to a higher frequency of chromosomal translocation. Variation in the degree of phenotypic diversity observed in different segregating populations suggested that genome dosage effects phenotypic manifestation. These findings not only broaden our understanding of the mechanisms of polyploid formation and reproductive diversity but also provide a novel insight into genome pyramiding and shuffling driven genome dosage effects and phenotypic diversity.

Separating phases of allopolyploid evolution with resynthesized and natural Capsella bursa-pastoris

Allopolyploidization is a frequent evolutionary transition in plants that combines whole-genome duplication (WGD) and interspecific hybridization. The genome of an allopolyploid species results from initial interactions between parental genomes and long-term evolution. Distinguishing the contributions of these two phases is essential to understanding the evolutionary trajectory of allopolyploid species. Here, we compared phenotypic and transcriptomic changes in natural and resynthesized Capsella allotetraploids with their diploid parental species. We focused on phenotypic traits associated with the selfing syndrome and on transcription-level phenomena such as expression-level dominance (ELD), transgressive expression (TRE), and homoeolog expression bias (HEB). We found that selfing syndrome, high pollen, and seed quality in natural allotetraploids likely resulted from long-term evolution. Similarly, TRE and most down-regulated ELD were only found in natural allopolyploids. Natural allotetraploids also had more ELD toward the self-fertilizing parental species than resynthesized allotetraploids, mirroring the establishment of the selfing syndrome. However, short-term changes mattered, and 40% of the cases of ELD in natural allotetraploids were already observed in resynthesized allotetraploids. Resynthesized allotetraploids showed striking variation of HEB among chromosomes and individuals. Homoeologous synapsis was its primary source and may still be a source of genetic variation in natural allotetraploids. In conclusion, both short- and long-term mechanisms contributed to transcriptomic and phenotypic changes in natural allotetraploids. However, the initial gene expression changes were largely reshaped during long-term evolution leading to further morphological changes.

High quality genomes produced from single MinION flow cells clarify polyploid and demographic histories of critically endangered Fraxinus (ash) species

With populations of threatened and endangered species declining worldwide, efforts are being made to generate high quality genomic records of these species before they are lost forever. Here, we demonstrate that data from single Oxford Nanopore Technologies (ONT) MinION flow cells can, even in the absence of highly accurate short DNA-read polishing, produce high quality de novo plant genome assemblies adequate for downstream analyses, such as synteny and ploidy evaluations, paleodemographic analyses, and phylogenomics. This study focuses on three North American ash tree species in the genus Fraxinus (Oleaceae) that were recently added to the International Union for Conservation of Nature (IUCN) Red List as critically endangered. Our results support a hexaploidy event at the base of the Oleaceae as well as a subsequent whole genome duplication shared by Syringa, Osmanthus, Olea, and Fraxinus. Finally, we demonstrate the use of ONT long-read sequencing data to reveal patterns in demographic history.

The genetic and molecular basis of haploinsufficiency in flowering plants

In diploid organisms, haploinsufficiency can be defined as the requirement for more than one fully functional copy of a gene. In contrast to most genes, whose loss-of-function alleles are recessive, loss-of-function alleles of haploinsufficient genes are dominant. However, forward and reverse genetic screens are biased toward obtaining recessive, loss-of-function mutations, and therefore, dominant mutations of all types are underrepresented in mutant collections. Despite this underrepresentation, haploinsufficient loci have intriguing implications for studies of genome evolution, gene dosage, stability of protein complexes, genetic redundancy, and gene expression. Here we review examples of haploinsufficiency in flowering plants and describe the underlying molecular mechanisms and evolutionary forces driving haploinsufficiency. Finally, we discuss the masking of haploinsufficiency by genetic redundancy, a widespread phenomenon among angiosperms.

Phased gap-free genome assembly of octoploid cultivated strawberry illustrates the genetic and epigenetic divergence among subgenomes

The genetic and epigenetic mechanisms underlying the coexistence and coordination of the four diverged subgenomes (ABCD) in octoploid strawberries (Fragaria × ananassa) remains poorly understood. In this study, we have assembled a haplotype-phased gap-free octoploid genome for the strawberry, which allowed us to uncover the sequence, structure, and epigenetic divergences among the subgenomes. The diploid progenitors of the octoploid strawberry, apart from subgenome A (Fragaria vesca), have been a subject of public controversy. Phylogenomic analyses revealed a close relationship between diploid species Fragaria iinumae and subgenomes B, C, and D. Subgenome A, closely related to F. vesca, retains the highest number of genes, exhibits the lowest content of transposable elements (TEs), experiences the strongest purifying selection, shows the lowest DNA methylation levels, and displays the highest expression level compared to the other three subgenomes. Transcriptome and DNA methylome analyses revealed that subgenome A-biased genes were enriched in fruit development biological processes. In contrast, although subgenomes B, C, and D contain equivalent amounts of repetitive sequences, they exhibit diverged methylation levels, particularly for TEs located near genes. Taken together, our findings provide valuable insights into the evolutionary patterns of subgenome structure, divergence and epigenetic dynamics in octoploid strawberries, which could be utilized in strawberry genetics and breeding research.

Gene expression bias between the subgenomes of allopolyploid hybrids is an emergent property of the kinetics of expression

Hybridization coupled to polyploidy, or allopolyploidy, has dramatically shaped the evolution of flowering plants, teleost fishes, and other lineages. Studies of recently formed allopolyploid plants have shown that the two subgenomes that merged to form that new allopolyploid do not generally express their genes equally. Instead, one of the two subgenomes expresses its paralogs more highly on average. Meanwhile, older allopolyploidy events tend to show biases in duplicate losses, with one of the two subgenomes retaining more genes than the other. Since reduced expression is a pathway to duplicate loss, understanding the origins of expression biases may help explain the origins of biased losses. Because we expect gene expression levels to experience stabilizing selection, our conceptual frameworks for how allopolyploid organisms form tend to assume that the new allopolyploid will show balanced expression between its subgenomes. It is then necessary to invoke phenomena such as differences in the suppression of repetitive elements to explain the observed expression imbalances. Here we show that, even for phenotypically identical diploid progenitors, the inherent kinetics of gene expression give rise to biases between the expression levels of the progenitor genes in the hybrid. Some of these biases are expected to be gene-specific and not give rise to global differences in progenitor gene expression. However, particularly in the case of allopolyploids formed from progenitors with different genome sizes, global expression biases favoring one subgenome are expected immediately on formation. Hence, expression biases are arguably the expectation upon allopolyploid formation rather than a phenomenon needing explanation. In the future, a deeper understanding of the kinetics of allopolyploidy may allow us to better understand both biases in duplicate losses and hybrid vigor.

GRASSY TILLERS1 (GT1) and SIX-ROWED SPIKE1 (VRS1) homologs share conserved roles in growth repression

Crop engineering and de novo domestication using gene editing are new frontiers in agriculture. However, outside of well-studied crops and model systems, prioritizing engineering targets remains challenging. Evolution can guide us, revealing genes with deeply conserved roles that have repeatedly been selected in the evolution of plant form. Homologs of the transcription factor genes GRASSY TILLERS1 (GT1) and SIX-ROWED SPIKE1 (VRS1) have repeatedly been targets of selection in domestication and evolution, where they repress growth in many developmental contexts. This suggests a conserved role for these genes in regulating growth repression. To test this, we determined the roles of GT1 and VRS1 homologs in maize (Zea mays) and the distantly related grass brachypodium (Brachypodium distachyon) using gene editing and mutant analysis. In maize, gt1; vrs1-like1 (vrl1) mutants have derepressed growth of floral organs. In addition, gt1; vrl1 mutants bore more ears and more branches, indicating broad roles in growth repression. In brachypodium, Bdgt1; Bdvrl1 mutants have more branches, spikelets, and flowers than wild-type plants, indicating conserved roles for GT1 and VRS1 homologs in growth suppression over ca. 59 My of grass evolution. Importantly, many of these traits influence crop productivity. Notably, maize GT1 can suppress growth in arabidopsis (Arabidopsis thaliana) floral organs, despite ca. 160 My of evolution separating the grasses and arabidopsis. Thus, GT1 and VRS1 maintain their potency as growth regulators across vast timescales and in distinct developmental contexts. This work highlights the power of evolution to inform gene editing in crop improvement.

Maternal dominance contributes to subgenome differentiation in allopolyploid fishes

Teleost fishes, which are the largest and most diverse group of living vertebrates, have a rich history of ancient and recent polyploidy. Previous studies of allotetraploid common carp and goldfish (cyprinids) reported a dominant subgenome, which is more expressed and exhibits biased gene retention. However, the underlying mechanisms contributing to observed 'subgenome dominance' remains poorly understood. Here we report high-quality genomes of twenty-one cyprinids to investigate the origin and subsequent subgenome evolution patterns following three independent allopolyploidy events. We identify the closest extant relatives of the diploid progenitor species, investigate genetic and epigenetic differences among subgenomes, and conclude that observed subgenome dominance patterns are likely due to a combination of maternal dominance and transposable element densities in each polyploid. These findings provide an important foundation to understanding subgenome dominance patterns observed in teleost fishes, and ultimately the role of polyploidy in contributing to evolutionary innovations.

A genome assembly for Orinus kokonorica provides insights into the origin, adaptive evolution and further diversification of two closely related grass genera

Ancient whole-genome duplication (WGD) or polyploidization is prevalent in plants and has played a crucial role in plant adaptation. However, the underlying genomic basis of ecological adaptation and subsequent diversification after WGD are still poorly understood in most plants. Here, we report a chromosome-scale genome assembly for the genus Orinus (Orinus kokonorica as representative) and preform comparative genomics with its closely related genus Cleistogenes (Cleistogenes songorica as representative), both belonging to a newly named subtribe Orininae of the grass subfamily Chloridoideae. The two genera may share one paleo-allotetraploidy event before 10 million years ago, and the two subgenomes of O. kokonorica display neither fractionation bias nor global homoeolog expression dominance. We find substantial genome rearrangements and extensive structural variations (SVs) between the two species. With comparative transcriptomics, we demonstrate that functional innovations of orthologous genes may have played an important role in promoting adaptive evolution and diversification of the two genera after polyploidization. In addition, copy number variations and extensive SVs between orthologs of flower and rhizome related genes may contribute to the morphological differences between the two genera. Our results provide new insights into the adaptive evolution and subsequent diversification of the two genera after polyploidization.

Evolutionary analysis of conserved non-coding elements subsequent to whole-genome duplication in opium poppy

Whole-genome duplication (WGD) leads to the duplication of both coding and non-coding sequences within an organism's genome, providing an abundant supply of genetic material that can drive evolution, ultimately contributing to plant adaptation and speciation. Although non-coding sequences contain numerous regulatory elements, they have been understudied compared to coding sequences. In order to address this gap, we explored the evolutionary patterns of regulatory sequences, coding sequences and transcriptomes using conserved non-coding elements (CNEs) as regulatory element proxies following the recent WGD event in opium poppy (Papaver somniferum). Our results showed similar evolutionary patterns in subgenomes of regulatory and coding sequences. Specifically, the biased or unbiased retention of coding sequences reflected the same pattern as retention levels in regulatory sequences. Further, the divergence of gene expression patterns mediated by regulatory element variations occurred at a more rapid pace than that of gene coding sequences. However, gene losses were purportedly dependent on relaxed selection pressure in coding sequences. Specifically, the rapid evolution of tissue-specific benzylisoquinoline alkaloid production in P. somniferum was associated with regulatory element changes. The origin of a novel stem-specific ACR, which utilized ancestral cis-elements as templates, is likely to be linked to the evolutionary trajectory behind the transition of the PSMT1-CYP719A21 cluster from high levels of expression solely in P. rhoeas root tissue to its elevated expression in P. somniferum stem tissue. Our findings demonstrate that rapid regulatory element evolution can contribute to the emergence of new phenotypes and provide valuable insights into the high evolvability of regulatory elements.

Subgenome phasing for complex allopolyploidy: case-based benchmarking and recommendations

Accurate subgenome phasing is crucial for understanding the origin, evolution and adaptive potential of polyploid genomes. SubPhaser and WGDI software are two common methodologies for subgenome phasing in allopolyploids, particularly in scenarios lacking known diploid progenitors. Triggered by a recent debate over the subgenomic origins of the cultivated octoploid strawberry, we examined four well-documented complex allopolyploidy cases as benchmarks, to evaluate and compare the accuracy of the two software. Our analysis demonstrates that the subgenomic structure phased by both software is in line with prior research, effectively tracing complex allopolyploid evolutionary trajectories despite the limitations of each software. Furthermore, using these validated methodologies, we revisited the controversial issue regarding the progenitors of the octoploid strawberry. The results of both methodologies reaffirm Fragaria vesca and Fragaria iinumae as progenitors of the octoploid strawberry. Finally, we propose recommendations for enhancing the accuracy of subgenome phasing in future studies, recognizing the potential of integrated tools for advanced complex allopolyploidy research and offering a new roadmap for robust subgenome-based phylogenetic analysis.

Cold stress induces differential gene expression of retained homeologs in Camelina sativa cv Suneson

Camelina sativa (L.) Crantz, a member of the Brassicaceae, has potential as a biofuel feedstock which is attributable to the production of fatty acids in its seeds, its fast growth cycle, and low input requirements. While a genome assembly is available for camelina, it was generated from short sequence reads and is thus highly fragmented in nature. Using long read sequences, we generated a chromosome-scale, highly contiguous genome assembly (644,491,969 bp) for the spring biotype cultivar 'Suneson' with an N50 contig length of 12,031,512 bp and a scaffold N50 length of 32,184,682 bp. Annotation of protein-coding genes revealed 91,877 genes that encode 133,355 gene models. We identified a total of 4,467 genes that were significantly up-regulated under cold stress which were enriched in gene ontology terms associated with "response to cold" and "response to abiotic stress". Coexpression analyses revealed multiple coexpression modules that were enriched in genes differentially expressed following cold stress that had putative functions involved in stress adaptation, specifically within the plastid. With access to a highly contiguous genome assembly, comparative analyses with Arabidopsis thaliana revealed 23,625 A. thaliana genes syntenic with 45,453 Suneson genes. Of these, 24,960 Suneson genes were syntenic to 8,320 A. thaliana genes reflecting a 3 camelina homeolog to 1 Arabidopsis gene relationship and retention of all three homeologs. Some of the retained triplicated homeologs showed conserved gene expression patterns under control and cold-stressed conditions whereas other triplicated homeologs displayed diverged expression patterns revealing sub- and neo-functionalization of the homeologs at the transcription level. Access to the chromosome-scale assembly of Suneson will enable both basic and applied research efforts in the improvement of camelina as a sustainable biofuel feedstock.

Dynamics of cis-regulatory sequences and transcriptional divergence of duplicated genes in soybean

Transcriptional divergence of duplicated genes after whole genome duplication (WGD) has been described in many plant lineages and is often associated with subgenome dominance, a genome-wide mechanism. However, it is unknown what underlies the transcriptional divergence of duplicated genes in polyploid species that lack subgenome dominance. Soybean is a paleotetraploid with a WGD that occurred 5 to 13 Mya. Approximately 50% of the duplicated genes retained from this WGD exhibit transcriptional divergence. We developed accessible chromatin region (ACR) datasets from leaf, flower, and seed tissues using MNase-hypersensitivity sequencing. We validated enhancer function of several ACRs associated with known genes using CRISPR/Cas9-mediated genome editing. The ACR datasets were used to examine and correlate the transcriptional patterns of 17,111 pairs of duplicated genes in different tissues. We demonstrate that ACR dynamics are correlated with divergence of both expression level and tissue specificity of individual gene pairs. Gain or loss of flanking ACRs and mutation of cis-regulatory elements (CREs) within the ACRs can change the balance of the expression level and/or tissue specificity of the duplicated genes. Analysis of DNA sequences associated with ACRs revealed that the extensive sequence rearrangement after the WGD reshaped the CRE landscape, which appears to play a key role in the transcriptional divergence of duplicated genes in soybean. This may represent a general mechanism for transcriptional divergence of duplicated genes in polyploids that lack subgenome dominance.

The lack of negative association between TE load and subgenome dominance in synthesized Brassica allotetraploids

Polyploidization is important to the evolution of plants. Subgenome dominance is a distinct phenomenon associated with most allopolyploids. A gene on the dominant subgenome tends to express to higher RNA levels in all organs as compared to the expression of its syntenic paralogue (homoeolog). The mechanism that underlies the formation of subgenome dominance remains unknown, but there is evidence for the involvement of transposon/DNA methylation density differences nearby the genes of parents as being causal. The subgenome with lower density of transposon and methylation near genes is positively associated with subgenome dominance. Here, we generated eight generations of allotetraploid progenies from the merging of parental genomes Brassica rapa and Brassica oleracea. We found that transposon/methylation density differ near genes between the parental (rapa:oleracea) existed in the wide hybrid, persisted in the neotetraploids (the synthetic Brassica napus), but these neotetraploids expressed no expected subgenome dominance. This absence of B. rapa vs. B. oleracea subgenome dominance is particularly significant because, while there is no negative relationship between transposon/methylation level and subgenome dominance in the neotetraploids, the more ancient parental subgenomes for all Brassica did show differences in transposon/methylation densities near genes and did express, in the same samples of cells, biased gene expression diagnostic of subgenome dominance. We conclude that subgenome differences in methylated transposon near genes are not sufficient to initiate the biased gene expressions defining subgenome dominance. Our result was unexpected, and we suggest a "nuclear chimera" model to explain our data.

Insights into the Evolution of Ohnologous Sequences and Their Epigenetic Marks Post-WGD in Malus Domestica

A Whole Genome Duplication (WGD) event occurred several Ma in a Rosaceae ancestor, giving rise to the Maloideae subfamily which includes today many pome fruits such as pear (Pyrus communis) and apple (Malus domestica). This complete and well-conserved genome duplication makes the apple an organism of choice to study the early evolutionary events occurring to ohnologous chromosome fragments. In this study, we investigated gene sequence evolution and expression, transposable elements (TE) density, and DNA methylation level. Overall, we identified 16,779 ohnologous gene pairs in the apple genome, confirming the relatively recent WGD. We identified several imbalances in QTL localization among duplicated chromosomal fragments and characterized various biases in genome fractionation, gene transcription, TE densities, and DNA methylation. Our results suggest a particular chromosome dominance in this autopolyploid species, a phenomenon that displays similarities with subgenome dominance that has only been described so far in allopolyploids.

Control of cellularization, nuclear localization, and antipodal cell cluster development in maize embryo sacs

The maize female gametophyte contains four cell types: two synergids, an egg cell, a central cell, and a variable number of antipodal cells. In maize, these cells are produced after three rounds of free-nuclear divisions followed by cellularization, differentiation, and proliferation of the antipodal cells. Cellularization of the eight-nucleate syncytium produces seven cells with two polar nuclei in the central cell. Nuclear localization is tightly controlled in the embryo sac. This leads to precise allocation of the nuclei into the cells upon cellularization. Nuclear positioning within the syncytium is highly correlated with their identity after cellularization. Two mutants are described with extra polar nuclei, abnormal antipodal cell morphology, and reduced antipodal cell number, as well as frequent loss of antipodal cell marker expression. Mutations in one of these genes, indeterminate gametophyte2 encoding a MICROTUBULE ASSOCIATED PROTEIN65-3 homolog, shows a requirement for MAP65-3 in cellularization of the syncytial embryo sac as well as for normal seed development. The timing of the effects of ig2 suggests that the identity of the nuclei in the syncytial female gametophyte can be changed very late before cellularization.

Improved pearl millet genomes representing the global heterotic pool offer a framework for molecular breeding applications

High-quality reference genome assemblies, representative of global heterotic patterns, offer an ideal platform to accurately characterize and utilize genetic variation in the primary gene pool of hybrid crops. Here we report three platinum grade de-novo, near gap-free, chromosome-level reference genome assemblies from the active breeding germplasm in pearl millet with a high degree of contiguity, completeness, and accuracy. An improved Tift genome (Tift23D2B1-P1-P5) assembly has a contig N50 ~ 7,000-fold (126 Mb) compared to the previous version and better alignment in centromeric regions. Comparative genome analyses of these three lines clearly demonstrate a high level of collinearity and multiple structural variations, including inversions greater than 1 Mb. Differential genes in improved Tift genome are enriched for serine O-acetyltransferase and glycerol-3-phosphate metabolic process which play an important role in improving the nutritional quality of seed protein and disease resistance in plants, respectively. Multiple marker-trait associations are identified for a range of agronomic traits, including grain yield through genome-wide association study. Improved genome assemblies and marker resources developed in this study provide a comprehensive framework/platform for future applications such as marker-assisted selection of mono/oligogenic traits as well as whole-genome prediction and haplotype-based breeding of complex traits.

Dosage-sensitivity shapes how genes transcriptionally respond to allopolyploidy and homoeologous exchange in resynthesized Brassica napus

The gene balance hypothesis proposes that selection acts on the dosage (i.e. copy number) of genes within dosage-sensitive portions of networks, pathways, and protein complexes to maintain balanced stoichiometry of interacting proteins, because perturbations to stoichiometric balance can result in reduced fitness. This selection has been called dosage balance selection. Dosage balance selection is also hypothesized to constrain expression responses to dosage changes, making dosage-sensitive genes (those encoding members of interacting proteins) experience more similar expression changes. In allopolyploids, where whole-genome duplication involves hybridization of diverged lineages, organisms often experience homoeologous exchanges that recombine, duplicate, and delete homoeologous regions of the genome and alter the expression of homoeologous gene pairs. Although the gene balance hypothesis makes predictions about the expression response to homoeologous exchanges, they have not been empirically tested. We used genomic and transcriptomic data from 6 resynthesized, isogenic Brassica napus lines over 10 generations to identify homoeologous exchanges, analyzed expression responses, and tested for patterns of genomic imbalance. Groups of dosage-sensitive genes had less variable expression responses to homoeologous exchanges than dosage-insensitive genes, a sign that their relative dosage is constrained. This difference was absent for homoeologous pairs whose expression was biased toward the B. napus A subgenome. Finally, the expression response to homoeologous exchanges was more variable than the response to whole-genome duplication, suggesting homoeologous exchanges create genomic imbalance. These findings expand our knowledge of the impact of dosage balance selection on genome evolution and potentially connect patterns in polyploid genomes over time, from homoeolog expression bias to duplicate gene retention.

ZmGI2 regulates flowering time through multiple flower development pathways in maize

GIGANTEA (GI) encodes a component of the circadian clock core oscillator and has been identified as a regulatory pathway of the circadian rhythm and photoperiodic flowering in model plants. However, the regulatory pathway of GI affecting flowering time is unknown in maize. Here, we identified that the zmgi2 mutant flowered earlier than the wild type under long day (LD) conditions, whereas the difference in flowering time was not apparent under short day (SD) conditions. The 24 h optimal expression of the gene in the stem apex meristems (SAM) appeared at 9 h after dawn under LD conditions and at 11 h after dawn under SD conditions. DAP-Seq and RNA-Seq further revealed that ZmGI2 delays flowering by directly binding to the upstream regions of ZmVOZs, ZmZCN8 and ZmFPF1 to repress the expression of these genes and by directly binding to the upstream regions of ZmARR11, ZmDOF and ZmUBC11 to promote the expression of these genes. The genetic and biochemical evidence suggests a model for the potential role of ZmGI2 in regulating the flowering time-dependent photoperiodic pathway. This study provides novel insights into the function of ZmGIs in maize and further demonstrates their potential importance for floral transition. These results contribute to a comprehensive understanding of the molecular mechanisms and regulatory networks of GI transcription factors in regulating flowering time in maize.

Design, execution, and interpretation of plant RNA-seq analyses

Genomics has transformed our understanding of the genetic architecture of traits and the genetic variation present in plants. Here, we present a review of how RNA-seq can be performed to tackle research challenges addressed by plant sciences. We discuss the importance of experimental design in RNA-seq, including considerations for sampling and replication, to avoid pitfalls and wasted resources. Approaches for processing RNA-seq data include quality control and counting features, and we describe common approaches and variations. Though differential gene expression analysis is the most common analysis of RNA-seq data, we review multiple methods for assessing gene expression, including detecting allele-specific gene expression and building co-expression networks. With the production of more RNA-seq data, strategies for integrating these data into genetic mapping pipelines is of increased interest. Finally, special considerations for RNA-seq analysis and interpretation in plants are needed, due to the high genome complexity common across plants. By incorporating informed decisions throughout an RNA-seq experiment, we can increase the knowledge gained.

Homoeologous evolution of the allotetraploid genome of Poa annua L

BACKGROUND

Poa annua (annual bluegrass) is an allotetraploid turfgrass, an agronomically significant weed, and one of the most widely dispersed plant species on earth. Here, we report the chromosome-scale genome assemblies of P. annua's diploid progenitors, P. infirma and P. supina, and use multi-omic analyses spanning all three species to better understand P. annua's evolutionary novelty.

RESULTS

We find that the diploids diverged from their common ancestor 5.5 - 6.3 million years ago and hybridized to form P. annua ≤ 50,000 years ago. The diploid genomes are similar in chromosome structure and most notably distinguished by the divergent evolutionary histories of their transposable elements, leading to a 1.7 × difference in genome size. In allotetraploid P. annua, we find biased movement of retrotransposons from the larger (A) subgenome to the smaller (B) subgenome. We show that P. annua's B subgenome is preferentially accumulating genes and that its genes are more highly expressed. Whole-genome resequencing of several additional P. annua accessions revealed large-scale chromosomal rearrangements characterized by extensive TE-downsizing and evidence to support the Genome Balance Hypothesis.

CONCLUSIONS

The divergent evolutions of the diploid progenitors played a central role in conferring onto P. annua its remarkable phenotypic plasticity. We find that plant genes (guided by selection and drift) and transposable elements (mostly guided by host immunity) each respond to polyploidy in unique ways and that P. annua uses whole-genome duplication to purge highly parasitized heterochromatic sequences. The findings and genomic resources presented here will enable the development of homoeolog-specific markers for accelerated weed science and turfgrass breeding.

Patterns, mechanisms, and consequences of homoeologous exchange in allopolyploid angiosperms: a genomic and epigenomic perspective

Allopolyploids result from hybridization between different evolutionary lineages coupled with genome doubling. Homoeologous chromosomes (chromosomes with common shared ancestry) may undergo recombination immediately after allopolyploid formation and continue over successive generations. The outcome of this meiotic pairing behavior is dynamic and complex. Homoeologous exchanges (HEs) may lead to the formation of unbalanced gametes, reduced fertility, and selective disadvantage. By contrast, HEs could act as sources of novel evolutionary substrates, shifting the relative dosage of parental gene copies, generating novel phenotypic diversity, and helping the establishment of neo-allopolyploids. However, HE patterns vary among lineages, across generations, and even within individual genomes and chromosomes. The causes and consequences of this variation are not fully understood, though interest in this evolutionary phenomenon has increased in the last decade. Recent technological advances show promise in uncovering the mechanistic basis of HEs. Here, we describe recent observations of the common patterns among allopolyploid angiosperm lineages, underlying genomic and epigenomic features, and consequences of HEs. We identify critical research gaps and discuss future directions with far-reaching implications in understanding allopolyploid evolution and applying them to the development of important phenotypic traits of polyploid crops.

Two-step model of paleohexaploidy, ancestral genome reshuffling and plasticity of heat shock response in Asteraceae

An ancient hexaploidization event in the most but not all Asteraceae plants, may have been responsible for shaping the genomes of many horticultural, ornamental, and medicinal plants that promoting the prosperity of the largest angiosperm family on the earth. However, the duplication process of this hexaploidy, as well as the genomic and phenotypic diversity of extant Asteraceae plants caused by paleogenome reorganization, are still poorly understood. We analyzed 11 genomes from 10 genera in Asteraceae, and redated the Asteraceae common hexaploidization (ACH) event ~70.7-78.6 million years ago (Mya) and the Asteroideae specific tetraploidization (AST) event ~41.6-46.2 Mya. Moreover, we identified the genomic homologies generated from the ACH, AST and speciation events, and constructed a multiple genome alignment framework for Asteraceae. Subsequently, we revealed biased fractionations between the paleopolyploidization produced subgenomes, suggesting the ACH and AST both are allopolyplodization events. Interestingly, the paleochromosome reshuffling traces provided clear evidence for the two-step duplications of ACH event in Asteraceae. Furthermore, we reconstructed ancestral Asteraceae karyotype (AAK) that has 9 paleochromosomes, and revealed a highly flexible reshuffling of Asteraceae paleogenome. Of specific significance, we explored the genetic diversity of Heat Shock Transcription Factors (Hsfs) associated with recursive whole-genome polyploidizations, gene duplications, and paleogenome reshuffling, and revealed that the expansion of Hsfs gene families enable heat shock plasticity during the genome evolution of Asteraceae. Our study provides insights on polyploidy and paleogenome remodeling for the successful establishment of Asteraceae, and is helpful for further communication and exploration of the diversification of plant families and phenotypes.

Transposon signatures of allopolyploid genome evolution

Hybridization brings together chromosome sets from two or more distinct progenitor species. Genome duplication associated with hybridization, or allopolyploidy, allows these chromosome sets to persist as distinct subgenomes during subsequent meioses. Here, we present a general method for identifying the subgenomes of a polyploid based on shared ancestry as revealed by the genomic distribution of repetitive elements that were active in the progenitors. This subgenome-enriched transposable element signal is intrinsic to the polyploid, allowing broader applicability than other approaches that depend on the availability of sequenced diploid relatives. We develop the statistical basis of the method, demonstrate its applicability in the well-studied cases of tobacco, cotton, and Brassica napus, and apply it to several cases: allotetraploid cyprinids, allohexaploid false flax, and allooctoploid strawberry. These analyses provide insight into the origins of these polyploids, revise the subgenome identities of strawberry, and provide perspective on subgenome dominance in higher polyploids.

Polyploidy and diploidization in soybean

Polyploidy is widespread and particularly common in angiosperms. The prevalence of polyploidy in the plant suggests it as a crucial driver of diversification and speciation. The paleopolyploid soybean (Glycine max) is one of the most important crops of plant protein and oil for humans and livestock. Soybean experienced two rounds of whole genome duplication around 13 and 59 million years ago. Due to the relatively slow process of post-polyploid diploidization, most genes are present in multiple copies across the soybean genome. Growing evidence suggests that polyploidization and diploidization could cause rapid and dramatic changes in genomic structure and epigenetic modifications, including gene loss, transposon amplification, and reorganization of chromatin architecture. This review is focused on recent progresses about genetic and epigenetic changes during polyploidization and diploidization of soybean and represents the challenges and potentials for application of polyploidy in soybean breeding.

A Study of GUS Expression in Arabidopsis as a Tool for the Evaluation of Gene Evolution, Function and the Role of Expression Derived from Gene Duplication

Gene duplication played a fundamental role in eukaryote evolution and different copies of a given gene can be present in extant species, often with expressions and functions differentiated during evolution. We assume that, when such differentiation occurs in a gene copy, this may be indicated by its maintenance in all the derived species. To verify this hypothesis, we compared the histological expression domains of the three β-glucuronidase genes (AtGUS) present in Arabidopsis thaliana with the GUS evolutionary tree in angiosperms. We found that AtGUS gene expression overlaps in the shoot apex, the floral bud and the root hairs. In the root apex, AtGUS3 expression differs completely from AtGUS1 and AtGUS2, whose transcripts are present in the root cap meristem and columella, in the staminal cell niche, in the epidermis and in the proximal cortex. Conversely, AtGUS3 transcripts are limited to the old border-like cells of calyptra and those found along the protodermal cell line. The GUS evolutionary tree reveals that the two main clusters (named GUS1 and GUS3) originate from a duplication event predating angiosperm radiation. AtGUS3 belongs to the GUS3 cluster, while AtGUS1 and AtGUS2, which originate from a duplication event that occurred in an ancestor of the Brassicaceae family, are found together in the GUS1 cluster. There is another, previously undescribed cluster, called GUS4, originating from a very ancient duplication event. While the copy of GUS4 has been lost in many species, copies of GUS3 and GUS1 have been conserved in all species examined.

Different evolutionary patterns of TIR1/AFBs and AUX/IAAs and their implications for the morphogenesis of land plants

BACKGROUND

The plant hormone auxin is widely involved in plant growth, development, and morphogenesis, and the TIR1/AFB and AUX/IAA proteins are closely linked to rapid auxin response and signal transmission. However, their evolutionary history, historical patterns of expansion and contraction, and changes in interaction relationships are still unknown.

RESULTS

Here, we analyzed the gene duplications, interactions, and expression patterns of TIR1/AFBs and AUX/IAAs to understand their underlying mechanisms of evolution. The ratios of TIR1/AFBs to AUX/IAAs range from 4:2 in Physcomitrium patens to 6:29 in Arabidopsis thaliana and 3:16 in Fragaria vesca. Whole-genome duplication (WGD) and tandem duplication have contributed to the expansion of the AUX/IAA gene family, but numerous TIR1/AFB gene duplicates were lost after WGD. We further analyzed the expression profiles of TIR1/AFBs and AUX/IAAs in different tissue parts of Physcomitrium patens, Selaginella moellendorffii, Arabidopsis thaliana and Fragaria vesca, and found that TIR1/AFBs and AUX/IAAs were highly expressed in all tissues in P. patens, S. moellendorffii. In A. thaliana and F. vesca, TIR1/AFBs maintained the same expression pattern as the ancient plants with high expression in all tissue parts, while AUX/IAAs appeared tissue-specific expression. In F. vesca, 11 AUX/IAAs interacted with TIR1/AFBs with different interaction strengths, and the functional specificity of AUX/IAAs was related to their ability to bind TIR1/AFBs, thus promoting the development of specific higher plant organs. Verification of the interactions among TIR1/AFBs and AUX/IAAs in Marchantia polymorpha and F. vesca also showed that the regulation of AUX/IAA members by TIR1/AFBs became more refined over the course of plant evolution.

CONCLUSIONS

Our results indicate that specific interactions and specific gene expression patterns both contributed to the functional diversification of TIR1/AFBs and AUX/IAAs.

A chromosome-scale genome assembly of Artemisia argyi reveals unbiased subgenome evolution and key contributions of gene duplication to volatile terpenoid diversity

Artemisia argyi Lévl. et Vant., a perennial Artemisia herb with an intense fragrance, is widely used in traditional medicine in China and many other Asian countries. Here, we present a chromosome-scale genome assembly of A. argyi comprising 3.89 Gb assembled into 17 pseudochromosomes. Phylogenetic and comparative genomic analyses revealed that A. argyi underwent a recent lineage-specific whole-genome duplication (WGD) event after divergence from Artemisia annua, resulting in two subgenomes. We deciphered the diploid ancestral genome of A. argyi, and unbiased subgenome evolution was observed. The recent WGD led to a large number of duplicated genes in the A. argyi genome. Expansion of the terpene synthase (TPS) gene family through various types of gene duplication may have greatly contributed to the diversity of volatile terpenoids in A. argyi. In particular, we identified a typical germacrene D synthase gene cluster within the expanded TPS gene family. The entire biosynthetic pathways of germacrenes, (+)-borneol, and (+)-camphor were elucidated in A. argyi. In addition, partial deletion of the amorpha-4,11-diene synthase (ADS) gene and loss of function of ADS homologs may have resulted in the lack of artemisinin production in A. argyi. Our study provides new insights into the genome evolution of Artemisia and lays a foundation for further improvement of the quality of this important medicinal plant.