Functional Divergence between Subgenomes and Gene Pairs after Whole Genome Duplications

Genome-Wide Identification of the Sulfate Transporter Gene Family Reveals That BolSULTR2;1 Regulates Plant Resistance to Alternaria brassicicola Through the Modulation of Glutathione Biosynthesis in Broccoli

Sulfate transporters (SULTRs) are key players that regulate sulfur acquisition and distribution within plants, thereby influencing cellular redox hemostasis under pathogen attacks, such as Alternaria brassicicola (Ab). In this study, a total of 23 BolSULTR (Brassica oleracea SULTR) genes were identified from the Brassica genome. These BolSULTRs are distributed across nine chromosomes, with all collinear BolSULTR gene pairs undergoing purifying selections. Phylogenetic analysis reveals that the SULTR family is evolutionarily conserved among plant kingdoms. qRT-PCR analysis demonstrated that the expression of BolSULTRs varies across different plant organs and is modulated by hormonal signals. Furthermore, transcriptome analysis identified several BolSULTRs whose expression levels were depressed in Ab-challenged leaves in broccoli. Among them, the BolSULTR2;1 gene emerged as a key player in the plant's response to Ab. Virus-induced gene silencing (VIGS) of BolSULTR2;1s resulted in elevated glutathione (GSH) levels and enhanced tolerance to Ab. Taken together, these findings underscore the role of BolSULTR2;1 in maintaining redox homeostasis and enhancing plant disease resistance, suggesting its potential as a target for genome editing to develop broccoli varieties with improved pathogen tolerance.

A stepwise guide for pangenome development in crop plants: an alfalfa (Medicago sativa) case study

BACKGROUND

The concept of pangenomics and the importance of structural variants is gaining recognition within the plant genomics community. Due to advancements in sequencing and computational technology, it has become feasible to sequence the entire genome of numerous individuals of a single species at a reasonable cost. Pangenomes have been constructed for many major diploid crops, including rice, maize, soybean, sorghum, pearl millet, peas, sunflower, grapes, and mustards. However, pangenomes for polyploid species are relatively scarce and are available in only few crops including wheat, cotton, rapeseed, and potatoes.

MAIN BODY

In this review, we explore the various methods used in crop pangenome development, discussing the challenges and implications of these techniques based on insights from published pangenome studies. We offer a systematic guide and discuss the tools available for constructing a pangenome and conducting downstream analyses. Alfalfa, a highly heterozygous, cross pollinated and autotetraploid forage crop species, is used as an example to discuss the concerns and challenges offered by polyploid crop species. We conducted a comparative analysis using linear and graph-based methods by constructing an alfalfa graph pangenome using three publicly available genome assemblies. To illustrate the intricacies captured by pangenome graphs for a complex crop genome, we used five different gene sequences and aligned them against the three graph-based pangenomes. The comparison of the three graph pangenome methods reveals notable variations in the genomic variation captured by each pipeline.

CONCLUSION

Pangenome resources are proving invaluable by offering insights into core and dispensable genes, novel gene discovery, and genome-wide patterns of variation. Developing user-friendly online portals for linear pangenome visualization has made these resources accessible to the broader scientific and breeding community. However, challenges remain with graph-based pangenomes including compatibility with other tools, extraction of sequence for regions of interest, and visualization of genetic variation captured in pangenome graphs. These issues necessitate further refinement of tools and pipelines to effectively address the complexities of polyploid, highly heterozygous, and cross-pollinated species.

Dosage sensitivity shapes balanced expression and gene longevity of homoeologs after whole-genome duplications in angiosperms

Following whole-genome duplication (WGD), duplicate gene pairs (homoeologs) can evolve varying degrees of expression divergence. However, the determinants influencing these relative expression level differences (RFPKM) between homoeologs remain elusive. In this study, we analyzed the RFPKM between homoeologs in 3 angiosperms, Nymphaea colorata, Nelumbo nucifera, and Acorus tatarinowii, all having undergone a single WGD since the origin of angiosperms. Our results show significant positive correlations in RFPKM of homoeologs among tissues within the same species, and among orthologs across these 3 species, indicating convergent expression balance/bias between homoeologous gene copies following independent WGDs. We linked RFPKM between homoeologs to gene attributes associated with dosage-balance constraints, such as protein-protein interactions, lethal-phenotype scores in Arabidopsis (Arabidopsis thaliana) orthologs, domain numbers, and expression breadth. Notably, homoeologs with lower RFPKM often had more interactions and higher lethal-phenotype scores, indicating selective pressures favoring balanced expression. Also, homoeologs with lower RFPKM were more likely to be retained after WGDs in angiosperms. Within Nelumbo, greater RFPKM between homoeologs correlated with increased cis- and trans-regulatory differentiation between species, highlighting the ongoing escalation of gene expression divergence. We further found that expression degeneration in 1 copy of homoeologs is inclined toward nonfunctionalization. Our research highlights the importance of balanced expression, shaped by dosage-balance constraints, in the evolutionary retention of homoeologs in plants.

Genome evolution of the ancient hexaploid Platanus × acerifolia (London planetree)

Whole-genome duplication (WGD; i.e., polyploidy) and chromosomal rearrangement (i.e., genome shuffling) significantly influence genome structure and organization. Many polyploids show extensive genome shuffling relative to their pre-WGD ancestors. No reference genome is currently available for Platanaceae (Proteales), one of the sister groups to the core eudicots. Moreover, Platanus × acerifolia (London planetree; Platanaceae) is a widely used street tree. Given the pivotal phylogenetic position of Platanus and its 2-y flowering transition, understanding its flowering-time regulatory mechanism has significant evolutionary implications; however, the impact of Platanus genome evolution on flowering-time genes remains unknown. Here, we assembled a high-quality, chromosome-level reference genome for P. × acerifolia using a phylogeny-based subgenome phasing method. Comparative genomic analyses revealed that P. × acerifolia (2n = 42) is an ancient hexaploid with three subgenomes resulting from two sequential WGD events; Platanus does not seem to share any WGD with other Proteales or with core eudicots. Each P. × acerifolia subgenome is highly similar in structure and content to the reconstructed pre-WGD ancestral eudicot genome without chromosomal rearrangements. The P. × acerifolia genome exhibits karyotypic stasis and gene sub-/neo-functionalization and lacks subgenome dominance. The copy number of flowering-time genes in P. × acerifolia has undergone an expansion compared to other noncore eudicots, mainly via the WGD events. Sub-/neo-functionalization of duplicated genes provided the genetic basis underlying the unique flowering-time regulation in P. × acerifolia. The P. × acerifolia reference genome will greatly expand understanding of the evolution of genome organization, genetic diversity, and flowering-time regulation in angiosperms.

Functional Differentiation of the Duplicated Gene BrrCIPK9 in Turnip (Brassica rapa var. rapa)

Gene duplication is a key biological process in the evolutionary history of plants and an important driving force for the diversification of genomic and genetic systems. Interactions between the calcium sensor calcineurin B-like protein (CBL) and its target, CBL-interacting protein kinase (CIPK), play important roles in the plant's response to various environmental stresses. As a food crop with important economic and research value, turnip (Brassica rapa var. rapa) has been well adapted to the environment of the Tibetan Plateau and become a traditional crop in the region. The BrrCIPK9 gene in turnip has not been characterized. In this study, two duplicated genes, BrrCIPK9.1 and BrrCIPK9.2, were screened from the turnip genome. Based on the phylogenetic analysis, BrrCIPK9.1 and BrrCIPK9.2 were found located in different sub-branches on the phylogenetic tree. Real-time fluorescence quantitative PCR analyses revealed their differential expression levels between the leaves and roots and in response to various stress treatments. The differences in their interactions with BrrCBLs were also revealed by yeast two-hybrid analyses. The results indicate that BrrCIPK9.1 and BrrCIPK9.2 have undergone Asparagine-alanine-phenylalanine (NAF) site divergence during turnip evolution, which has resulted in functional differences between them. Furthermore, BrrCIPK9.1 responded to high-pH (pH 8.5) stress, while BrrCIPK9.2 retained its ancestral function (low K+), thus providing further evidence of their functional divergence. These functional divergence genes facilitate turnip's good adaptation to the extreme environment of the Tibetan Plateau. In summary, the results of this study reveal the characteristics of the duplicated BrrCIPK9 genes and provide a basis for further functional studies of BrrCBLs-BrrCIPKs in turnip.

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.

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.

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.

Origin, evolution, and diversification of the wall-associated kinase gene family in plants

KEY MESSAGE

The study of the origin, evolution, and diversification of the wall-associated kinase gene family in plants facilitates their functional investigations in the future. Wall-associated kinases (WAKs) make up one subfamily of receptor-like kinases (RLKs), and function directly in plant cell elongation and responses to biotic and abiotic stresses. The biological functions of WAKs have been extensively characterized in angiosperms; however, the origin and evolutionary history of the WAK family in green plants remain unclear. Here, we performed a comprehensive analysis of the WAK family to reveal its origin, evolution, and diversification in green plants. In total, 1061 WAK genes were identified in 37 species from unicellular algae to multicellular plants, and the results showed that WAK genes probably originated before bryophyte differentiation and were widely distributed in land plants, especially angiosperms. The phylogeny indicated that the land plant WAKs gave rise to five clades and underwent lineage-specific expansion after species differentiation. Cis-acting elements and expression patterns analyses of WAK genes in Arabidopsis and rice demonstrated the functional diversity of WAK genes in these two species. Many gene gains and losses have occurred in angiosperms, leading to an increase in the number of gene copies. The evolutionary trajectory of the WAK family during polyploidization was uncovered using Gossypium species. Our results provide insights into the evolution of WAK genes in green plants, facilitating their functional investigations in the future.

Telomere-to-telomere pear (Pyrus pyrifolia) reference genome reveals segmental and whole genome duplication driving genome evolution

Previously released pear genomes contain a plethora of gaps and unanchored genetic regions. Here, we report a telomere-to-telomere (T2T) gap-free genome for the red-skinned pear, 'Yunhong No. 1' (YH1; Pyrus pyrifolia), which is mainly cultivated in Yunnan Province (southwest China), the pear's primary region of origin. The YH1 genome is 501.20 Mb long with a contig N50 length of 29.26 Mb. All 17 chromosomes were assembled to the T2T level with 34 characterized telomeres. The 17 centromeres were predicted and mainly consist of centromeric-specific monomers (CEN198) and long terminal repeat (LTR) Gypsy elements (≥74.73%). By filling all unclosed gaps, the integrity of YH1 is markedly improved over previous P. pyrifolia genomes ('Cuiguan' and 'Nijisseiki'). A total of 1531 segmental duplication (SD) driven duplicated genes were identified and enriched in stress response pathways. Intrachromosomal SDs drove the expansion of disease resistance genes, suggesting the potential of SDs in adaptive pear evolution. A large proportion of duplicated gene pairs exhibit dosage effects or sub-/neo-functionalization, which may affect agronomic traits like stone cell content, sugar content, and fruit skin russet. Furthermore, as core regulators of anthocyanin biosynthesis, we found that MYB10 and MYB114 underwent various gene duplication events. Multiple copies of MYB10 and MYB114 displayed obvious dosage effects, indicating role differentiation in the formation of red-skinned pear fruit. In summary, the T2T gap-free pear genome provides invaluable resources for genome evolution and functional genomics.

A near-complete genome assembly of the allotetrapolyploid Cenchrus fungigraminus (JUJUNCAO) provides insights into its evolution and C4 photosynthesis

JUJUNCAO (Cenchrus fungigraminus; 2n = 4x = 28) is a Cenchrus grass with the highest biomass production among cultivated plants, and it can be used for mushroom cultivation, animal feed, and biofuel production. Here, we report a nearly complete genome assembly of JUJUNCAO and reveal that JUJUNCAO is an allopolyploid that originated ∼2.7 million years ago (mya). Its genome consists of two subgenomes, and subgenome A shares high collinear synteny with pearl millet. We also investigated the genome evolution of JUJUNCAO and suggest that the ancestral karyotype of Cenchrus split into the A and B ancestral karyotypes of JUJUNCAO. Comparative transcriptome and DNA methylome analyses revealed functional divergence of homeologous gene pairs between the two subgenomes, which was a further indication of asymmetric DNA methylation. The three types of centromeric repeat in the JUJUNCAO genome (CEN137, CEN148, and CEN156) may have evolved independently within each subgenome, with some introgressions of CEN156 from the B to the A subgenome. We investigated the photosynthetic characteristics of JUJUNCAO, revealing its typical C4 Kranz anatomy and high photosynthetic efficiency. NADP-ME and PEPCK appear to cooperate in the major C4 decarboxylation reaction of JUJUNCAO, which is different from other C4 photosynthetic subtypes and may contribute to its high photosynthetic efficiency and biomass yield. Taken together, our results provide insights into the highly efficient photosynthetic mechanism of JUJUNCAO and provide a valuable reference genome for future genetic and evolutionary studies, as well as genetic improvement of Cenchrus grasses.

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.

Genome-wide analysis of NBS-LRR genes revealed contribution of disease resistance from Saccharum spontaneum to modern sugarcane cultivar

INTRODUCTION

During plant evolution, nucleotide-binding sites (NBS) and leucine-rich repeat (LRR) genes have made significant contributions to plant disease resistance. With many high-quality plant genomes sequenced, identification and comprehensive analyses of NBS-LRR genes at whole genome level are of great importance to understand and utilize them.

METHODS

In this study, we identified the NBS-LRR genes of 23 representative species at whole genome level, and researches on NBS-LRR genes of four monocotyledonous grass species, Saccharum spontaneum, Saccharum officinarum, Sorghum bicolor and Miscanthus sinensis, were focused.

RESULTS AND DISCUSSION

We found that whole genome duplication, gene expansion, and allele loss could be factors affecting the number of NBS-LRR genes in the species, and whole genome duplication is likely to be the main cause of the number of NBS-LRR genes in sugarcane. Meanwhile, we also found a progressive trend of positive selection on NBS-LRR genes. These studies further elucidated the evolutionary pattern of NBS-LRR genes in plants. Transcriptome data from multiple sugarcane diseases revealed that more differentially expressed NBS-LRR genes were derived from S. spontaneum than from S. officinarum in modern sugarcane cultivars, and the proportion was significantly higher than the expected. This finding reveals that S. spontaneum has a greater contribution to disease resistance for modern sugarcane cultivars. In addition, we observed allelespecific expression of seven NBS-LRR genes under leaf scald, and 125 NBS-LRR genes responding to multiple diseases were identified. Finally, we built a plant NBS-LRR gene database to facilitate subsequent analysis and use of NBSLRR genes obtained here. In conclusion, this study complemented and completed the research of plant NBS-LRR genes, and discussed how NBS-LRR genes responding to sugarcane diseases, which provided a guide and genetic resources for further research and utilization of NBS-LRR genes.

Comparative gene retention analysis in barley, wild emmer, and bread wheat pangenome lines reveals factors affecting gene retention following gene duplication

BACKGROUND

Gene duplication is a prevalent phenomenon and a major driving force underlying genome evolution. The process leading to the fixation of gene duplicates following duplication is critical to understand how genome evolves but remains fragmentally understood. Most previous studies on gene retention are based on gene duplicate analyses in single reference genome. No population-based comparative gene retention analysis has been performed to date.

RESULTS

Taking advantage of recently published genomic data in Triticeae, we dissected a divergent homogentisate phytyltransferase (HPT2) lineage caught in the middle stage of gene fixation following duplication. The presence/absence of HPT2 in barley (diploid), wild emmer (tetraploid), and bread wheat (hexaploid) pangenome lines appears to be associated with gene dosage constraint and environmental adaption. Based on these observations, we adopted a phylogeny-based orthology inference approach and performed comparative gene retention analyses across barley, wild emmer, and bread wheat. This led to the identification of 326 HPT2-pattern-like genes at whole genome scale, representing a pool of gene duplicates in the middle stage of gene fixation. Majority of these HPT2-pattern-like genes were identified as small-scale duplicates, such as dispersed, tandem, and proximal duplications. Natural selection analyses showed that HPT2-pattern-like genes have experienced relaxed selection pressure, which is generally accompanied with partial positive selection and transcriptional divergence. Functional enrichment analyses showed that HPT2-pattern-like genes are over-represented with molecular-binding and defense response functions, supporting the potential role of environmental adaption during gene retention. We also observed that gene duplicates from larger gene family are more likely to be lost, implying a gene dosage constraint effect. Further comparative gene retention analysis in barley and bread wheat pangenome lines revealed combined effects of species-specific selection and gene dosage constraint.

CONCLUSIONS

Comparative gene retention analyses at the population level support gene dosage constraint, environmental adaption, and species-specific selection as three factors that may affect gene retention following gene duplication. Our findings shed light on the evolutionary process leading to the retention of newly formed gene duplicates and will greatly improve our understanding on genome evolution via duplication.

A MYB Transcription Factor Atlas Provides Insights into the Evolution of Environmental Adaptations in Plants

The MYB transcription factor superfamily includes key regulators of plant development and responses to environmental changes. The diversity of lifestyles and morphological characteristics exhibited by plants are potentially associated with the genomic dynamics of the MYB superfamily. With the release of the plant genomes, a comprehensive phylogenomic analysis of the MYB superfamily across Viridiplantae is allowed. The present study performed phylogenetic, phylogenomic, syntenic, horizontal gene transfer, and neo/sub-functionalization analysis of the MYB superfamily to explore the evolutionary contributions of MYB members to species diversification, trait formation, and environmental adaptation in 437 different plant species. We identified major changes in copy number variation and genomic context within subclades across lineages. Multiple MYB subclades showed highly conserved copy number patterns and synteny across flowering plants, whereas others were more dynamic and showed lineage-specific patterns. As examples of lineage-specific morphological divergence, we hypothesize that the gain of a MYB orthogroup associated with flower development and environmental responses and an orthogroup associated with auxin and wax biosynthesis in angiosperms were correlated with the emergence of flowering plants, unbiased neo-/sub-functionalization of gene duplicates contributed to environmental adaptation, and species-specific neo-/sub-functionalization contributed to phenotype divergence between species. Transposable element insertion in promoter regions may have facilitated the sub-/neo-functionalization of MYB genes and likely played a tissue-specific role contributing to sub-/neo-functionalization in plant root tissues. This study provides new insights into the evolutionary divergence of the MYB superfamily across major flowering and non-flowering lineages and emphasizes the need for lineage-/tissue-specific characterization to further understand trait variability and environmental adaptation.

Flavan-3-ol Galloylation-Related Functional Gene Cluster and the Functional Diversification of SCPL Paralogs in Camellia sp

The high accumulation of galloylated flavan-3-ols in Camellia sp. is a noteworthy phenomenon. We identified a flavan-3-ol galloylation-related functional gene cluster in tannin-rich plant Camellia sp., which included UGT84A22 and SCPL-AT gene clusters. We investigated the possible correlation between the accumulation of metabolites and the expression of SCPL-ATs and UGT84A22. The results revealed that C. sinensis, C. ptilophylla, and C. oleifera accumulated galloylated cis-flavan-3-ols (EGCG), galloylated trans-flavan-3-ols (GCG), and hydrolyzed tannins, respectively; however, C. nitidissima did not accumulate any galloylated compounds. C. nitidissima exhibited no expression of SCPL-AT or UGT84A22, whereas the other three species of Camellia exhibited various expression patterns. This indicated that the functions of the paralogs of SCPL-AT vary. Enzymatic analysis revealed that SCPL5 was neofunctionalized as a noncatalytic chaperone paralog, a type of chaerone-like protein, associating with flavan-3-ol galloylation; moreover, CsSCPL4 was subfunctionalized in association with the galloylation of cis- and trans-flavan-3-ols. In C. nitidissima, an SCPL4 homolog was noted with mutations in two cysteine residues forming a disulfide bond, which suggested that this homolog was defunctionalized. The findings of this study improve our understanding of the functional diversification of SCPL paralogs in Camellia sp.

Plant pan-genomics and its applications

Plant genomes are so highly diverse that a substantial proportion of genomic sequences are not shared among individuals. The variable DNA sequences, along with the conserved core sequences, compose the more sophisticated pan-genome that represents the collection of all non-redundant DNA in a species. With rapid progress in genome sequencing technologies, pan-genome research in plants is now accelerating. Here we review recent advances in plant pan-genomics, including major driving forces of structural variations that constitute the variable sequences, methodological innovations for representing the pan-genome, and major successes in constructing plant pan-genomes. We also summarize recent efforts toward decoding the remaining dark matter in telomere-to-telomere or gapless plant genomes. These new genome resources, which have remarkable advantages over numerous previously assembled less-than-perfect genomes, are expected to become new references for genetic studies and plant breeding.

Potential allopolyploid origin of Ericales revealed with gene-tree reconciliation

Few incidents of ancient allopolyploidization (polyploidization by hybridization or merging diverged genomes) were previously revealed, although there is significant evidence for the accumulation of whole genome duplications (WGD) in plants. Here, we focused on Ericales, one of the largest and most diverse angiosperm orders with significant ornamental and economic value. Through integrating 24 high-quality whole genome data selected from ~ 200 Superasterids genomes/species and an algorithm of topology-based gene-tree reconciliation, we explored the evolutionary history of in Ericales with ancient complex. We unraveled the allopolyploid origin of Ericales and detected extensive lineage-specific gene loss following the polyploidization. Our study provided a new hypothesis regarding the origin of Ericales and revealed an instructive perspective of gene loss as a pervasive source of genetic variation and adaptive phenotypic diversity in Ericales.

Transcription factor retention through multiple polyploidization steps in wheat

Whole-genome duplication is widespread in plant evolutionary history and is followed by nonrandom gene loss to return to a diploid state. Across multiple angiosperm species, the retained genes tend to be dosage-sensitive regulatory genes such as transcription factors, yet data for younger polyploid species is sparse. Here, we analyzed the retention, expression, and genetic variation in transcription factors in the recent allohexaploid bread wheat (Triticum aestivum L.). By comparing diploid, tetraploid, and hexaploid wheat, we found that, following each of two hybridization and whole-genome duplication events, the proportion of transcription factors in the genome increased. Transcription factors were preferentially retained over other genes as homoeologous groups in tetraploid and hexaploid wheat. Across cultivars, transcription factor homoeologs contained fewer deleterious missense mutations than nontranscription factors, suggesting that transcription factors are maintained as three functional homoeologs in hexaploid wheat populations. Transcription factor homoeologs were more strongly coexpressed than nontranscription factors, indicating conservation of function between homoeologs. We found that the B3, MADS-M-type, and NAC transcription factor families were less likely to have three homoeologs present than other families, which was associated with low expression levels and high levels of tandem duplication. Together, our results show that transcription factors are preferentially retained in polyploid wheat genomes although there is variation between families. Knocking out one transcription factor homoeolog to alter gene dosage, using TILLING or CRISPR, could generate new phenotypes for wheat breeding.

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.

A cluster of transcripts identifies a transition stage initiating leafy head growth in heading morphotypes of Brassica

Leaf heading is an important and economically valuable horticultural trait in many vegetables. The formation of a leafy head is a specialized leaf morphogenesis characterized by the emergence of the enlarged incurving leaves. However, the transcriptional regulation mechanisms underlying the transition to leaf heading remain unclear. We carried out large-scale time-series transcriptome assays covering the major vegetative growth phases of two headingBrassica crops, Chinese cabbage and cabbage, with the non-heading morphotype Taicai as the control. A regulatory transition stage that initiated the heading process is identified, accompanied by a developmental switch from rosette leaf to heading leaf in Chinese cabbages. This transition did not exist in the non-heading control. Moreover, we reveal that the heading transition stage is also conserved in the cabbage clade. Chinese cabbage acquired through domestication a leafy head independently from the origins of heading in other cabbages; phylogenetics supports that the ancestor of all cabbages is non-heading. The launch of the transition stage is closely associated with the ambient temperature. In addition, examination of the biological activities in the transition stage identified the ethylene pathway as particularly active, and we hypothesize that this pathway was targeted for selection for domestication to form the heading trait specifically in Chinese cabbage. In conclusion, our findings on the transcriptome transition that initiated the leaf heading in Chinese cabbage and cabbage provide a new perspective for future studies of leafy head crops.

A High-Quality Haplotype-Resolved Genome of Common Bermudagrass (Cynodon dactylon L.) Provides Insights Into Polyploid Genome Stability and Prostrate Growth

Common bermudagrass (Cynodon dactylon L.) is an important perennial warm-season turfgrass species with great economic value. However, the reference genome is still deficient in C. dactylon, which severely impedes basic studies and breeding studies. In this study, a high-quality haplotype-resolved genome of C. dactylon cultivar Yangjiang was successfully assembled using a combination of multiple sequencing strategies. The assembled genome is approximately 1.01 Gb in size and is comprised of 36 pseudo chromosomes belonging to four haplotypes. In total, 76,879 protein-coding genes and 529,092 repeat sequences were annotated in the assembled genome. Evolution analysis indicated that C. dactylon underwent two rounds of whole-genome duplication events, whereas syntenic and transcriptome analysis revealed that global subgenome dominance was absent among the four haplotypes. Genome-wide gene family analyses further indicated that homologous recombination-regulating genes and tiller-angle-regulating genes all showed an adaptive evolution in C. dactylon, providing insights into genome-scale regulation of polyploid genome stability and prostrate growth. These results not only facilitate a better understanding of the complex genome composition and unique plant architectural characteristics of common bermudagrass, but also offer a valuable resource for comparative genome analyses of turfgrasses and other plant species.

A Dual-Successive-Screen Model at Pollen/Stigma and Pollen Tube/Ovary Explaining Paradoxical Self-Incompatibility Diagnosis in the Olive Tree—An Interpretative Update of the Literature

The ‘pollen test’ and ‘fruit set test’ following controlled crossing combinations of parents are the most commonly used methods for pollination incompatibility studies in Olea europaea L. Self-incompatibility (SI), with diagnoses based on the pollen test and pollen germination, indicating self-compatibility, is not always followed by fruit set in this species. To solve this dispute, we have reconciled all observations into a new model. Mismatches between field and laboratory data and between methods are resolved by the dual-successive-screen model (DSSM) supposing two different loci for the expression of the two SI mechanisms. Pollen/stigma is controlled by diallelic SI, or DSI, inferring two G1 and G2 compatibility/incompatibility (C/I) groups for varieties, whereas pollen tubes in ovaries are controlled by poly-allelic SI or PASI with twenty C/I groups. To explain the selfing of varieties, we have suggested that some determinants in the pollen tube and stigma are unstable and degrade (DS-D for degradation of S-determinant) after three to five days, enabling some pollen tubes to avoid being rejected, hence reaching ovules. DSI and PASI in the DSSM and DS-D mechanisms, plus the andromonoecy of the olive tree, complexify SI studies. Inferences from DSSM and DS-D mechanisms in olive orchard practice are detailed.

Characterization and expression analysis of wall-associated kinase (WAK) and WAK-like family in cotton

The wall-associated kinases (WAKs) and WAK-like kinases (WAKLs) form a group of receptor-like kinases (RLKs) with extracellular domains tightly linked to the cell wall. The WAKs/WAKLs have been known to be involved in plant growth, development, and stress responses. However, the functions of WAKs/WAKLs are less well known in cotton. In this study, 58, 66, and 99 WAK/WAKL genes were identified in Gossypium arboreum, G. raimondii, and G. hirsutum, respectively. Phylogenetic analysis showed they were classified into five groups, with two groups specific to cotton. Collinearity analysis revealed that segmental and tandem duplications resulted in expansion of the WAK/WAKL gene family in cotton. Moreover, the Ka/Ks ratios indicated this family was exposed to purifying selection pressure during evolution. The structures of the GhWAK/WAKL genes and encoded proteins suggested the functions of WAKs/WAKLs in cotton were conserved. Transient expression of four WAK/WAKL-GFP fusion constructs in Arabidopsis protoplasts indicated that they were localized on the plasma membrane. The cis-elements in the GhWAK/WAKL promoters were responsive to multiple phytohormones and abiotic stresses. Expression profiling showed that GhWAK/WAKL genes were induced by various abiotic stresses. This study provides insights into the evolution of WAK/WAKL genes and presents fundamental information for further analysis in cotton.

Genome-Wide Identification of Soybean ABC Transporters Relate to Aluminum Toxicity

ATP-binding cassette (ABC) transporter proteins are a gene super-family in plants and play vital roles in growth, development, and response to abiotic and biotic stresses. The ABC transporters have been identified in crop plants such as rice and buckwheat, but little is known about them in soybean. Soybean is an important oil crop and is one of the five major crops in the world. In this study, 255 ABC genes that putatively encode ABC transporters were identified from soybean through bioinformatics and then categorized into eight subfamilies, including 7 ABCAs, 52 ABCBs, 48 ABCCs, 5 ABCDs, 1 ABCEs, 10 ABCFs, 111 ABCGs, and 21 ABCIs. Their phylogenetic relationships, gene structure, and gene expression profiles were characterized. Segmental duplication was the main reason for the expansion of the GmABC genes. Ka/Ks analysis suggested that intense purifying selection was accompanied by the evolution of GmABC genes. The genome-wide collinearity of soybean with other species showed that GmABCs were relatively conserved and that collinear ABCs between species may have originated from the same ancestor. Gene expression analysis of GmABCs revealed the distinct expression pattern in different tissues and diverse developmental stages. The candidate genes GmABCB23, GmABCB25, GmABCB48, GmABCB52, GmABCI1, GmABCI5, and GmABCI13 were responsive to Al toxicity. This work on the GmABC gene family provides useful information for future studies on ABC transporters in soybean and potential targets for the cultivation of new germplasm resources of aluminum-tolerant soybean.

Genome-Wide Analysis of Coding and Non-coding RNA Reveals a Conserved miR164-NAC-mRNA Regulatory Pathway for Disease Defense in Populus

MicroRNAs (miRNAs) contribute to plant defense responses by increasing the overall genetic diversity; however, their origins and functional importance in plant defense remain unclear. Here, we employed Illumina sequencing technology to assess how miRNA and messenger RNA (mRNA) populations vary in the Chinese white poplar (Populus tomentosa) during a leaf black spot fungus (Marssonina brunnea) infection. We sampled RNAs from infective leaves at conidia germinated stage [12 h post-inoculation (hpi)], infective vesicles stage (24 hpi), and intercellular infective hyphae stage (48 hpi), three essential stages associated with plant colonization and biotrophic growth in M. brunnea fungi. In total, 8,938 conserved miRNA-target gene pairs and 3,901 Populus-specific miRNA-target gene pairs were detected. The result showed that Populus-specific miRNAs (66%) were more involved in the regulation of the disease resistance genes. By contrast, conserved miRNAs (>80%) target more whole-genome duplication (WGD)-derived transcription factors (TFs). Among the 1,023 WGD-derived TF pairs, 44.9% TF pairs had only one paralog being targeted by a miRNA that could be due to either gain or loss of a miRNA binding site after the WGD. A conserved hierarchical regulatory network combining promoter analyses and hierarchical clustering approach uncovered a miR164–NAM, ATAF, and CUC (NAC) transcription factor–mRNA regulatory module that has potential in Marssonina defense responses. Furthermore, analyses of the locations of miRNA precursor sequences reveal that pseudogenes and transposon contributed a certain proportion (∼30%) of the miRNA origin. Together, these observations provide evolutionary insights into the origin and potential roles of miRNAs in plant defense and functional innovation.

The contributions from the progenitor genomes of the mesopolyploid Brassiceae are evolutionarily distinct but functionally compatible

The members of the tribe Brassiceae share a whole-genome triplication (WGT), and one proposed model for its formation is a two-step pair of hybridizations producing hexaploid descendants. However, evidence for this model is incomplete, and the evolutionary and functional constraints that drove evolution after the hexaploidy are even less understood. Here, we report a new genome sequence of Crambe hispanica, a species sister to most sequenced Brassiceae. Using this new genome and three others that share the hexaploidy, we traced the history of gene loss after the WGT using the Polyploidy Orthology Inference Tool (POInT). We confirm the two-step formation model and infer that there was a significant temporal gap between those two allopolyploidizations, with about a third of the gene losses from the first two subgenomes occurring before the arrival of the third. We also, for the 90,000 individual genes in our study, make parental subgenome assignments, inferring, with measured uncertainty, from which of the progenitor genomes of the allohexaploidy each gene derives. We further show that each subgenome has a statistically distinguishable rate of homoeolog losses. There is little indication of functional distinction between the three subgenomes: the individual subgenomes show no patterns of functional enrichment, no excess of shared protein-protein or metabolic interactions between their members, and no biases in their likelihood of having experienced a recent selective sweep. We propose a "mix and match" model of allopolyploidy, in which subgenome origin drives homoeolog loss propensities but where genes from different subgenomes function together without difficulty.

Evolution of AITR family genes in cotton and their functions in abiotic stress tolerance

Abiotic stresses are major environmental factors inhibiting plant growth and development. AITRs (ABA-induced transcription repressors) are a novel family of transcription factors regulating ABA (abscisic acid) signalling and plant responses to abiotic stresses in Arabidopsis. However, the composition and evolution history of AITRs and their roles in the cotton genus are largely unknown. A total of 12 putative AITRs genes were identified in cultivated tetraploid cotton, Gossypium hirsutum. Phylogenetic analysis of GhAITRs in these tetraploid cottons and their closely related species implicate ancient genome-wide duplication occurring after speciation of Gossypium, and Theobroma could generate duplicates of GhAITRs. Duplicated GhAITRs were stably inherited following diploid speciation and further allotetraploidy in Gossypium. Homologous GhAITRs shared common expression patterns in response to ABA, drought and salinity treatments, and drought tolerance induced in transgenic Arabidopsis plants expressing GhAITR-A1. Together, our findings reveal that duplicates in the GhAITRs gene family were achieved by whole genome duplication rather than three individual duplication events, and that GhAITRs function as transcription repressors and are involved in the regulation of plant responses to ABA and drought stress. These results provide insights towards the improvement of abiotic stress tolerance in cotton using GhAITRs.

Whole-Genome Sequence of Synthesized Allopolyploids in Cucumis Reveals Insights into the Genome Evolution of Allopolyploidization

The importance of allopolyploidy in plant evolution has been widely recognized. The genetic changes triggered by allopolyploidy, however, are not yet fully understood due to inconsistent phenomena reported across diverse species. The construction of synthetic polyploids offers a controlled approach to systematically reveal genomic changes that occur during the process of polyploidy. This study reports the first fully sequenced synthetic allopolyploid constructed from a cross between Cucumis sativus and C. hystrix, with high-quality assembly. The two subgenomes are confidently partitioned and the C. sativus-originated subgenome predominates over the C. hystrix-originated subgenome, retaining more sequences and showing higher homeologous gene expression. Most of the genomic changes emerge immediately after interspecific hybridization. Analysis of a series of genome sequences from several generations (S0, S4-S13) of C. ×hytivus confirms that genomic changes occurred in the very first generations, subsequently slowing down as the process of diploidization is initiated. The duplicated genome of the allopolyploid with double genes from both parents broadens the genetic base of C. ×hytivus, resulting in enhanced phenotypic plasticity. This study provides novel insights into plant polyploid genome evolution and demonstrates a promising strategy for the development of a wide array of novel plant species and varieties through artificial polyploidization.

Gene Duplication Accelerates the Pace of Protein Gain and Loss from Plant Organelles

Organelle biogenesis and function is dependent on the concerted action of both organellar-encoded (if present) and nuclear-encoded proteins. Differences between homologous organelles across the Plant Kingdom arise, in part, as a result of differences in the cohort of nuclear-encoded proteins that are targeted to them. However, neither the rate at which differences in protein targeting accumulate nor the evolutionary consequences of these changes are known. Using phylogenomic approaches coupled to ancestral state estimation, we show that the plant organellar proteome has diversified in proportion with molecular sequence evolution such that the proteomes of plant chloroplasts and mitochondria lose or gain on average 3.6 proteins per million years. We further demonstrate that changes in organellar protein targeting are associated with an increase in the rate of molecular sequence evolution and that such changes predominantly occur in genes with regulatory rather than metabolic functions. Finally, we show that gain and loss of protein target signals occurs at a higher rate following gene duplication, revealing that gene and genome duplication are a key facilitator of plant organelle evolution.

Conservation and trans-regulation of histone modification in the A and B subgenomes of polyploid wheat during domestication and ploidy transition

BACKGROUND

Polyploidy has played a prominent role in the evolution of plants and many other eukaryotic lineages. However, how polyploid genomes adapt to the abrupt presence of two or more sets of chromosomes via genome regulation remains poorly understood. Here, we analyzed genome-wide histone modification and gene expression profiles in relation to domestication and ploidy transition in the A and B subgenomes of polyploid wheat.

RESULTS

We found that epigenetic modification patterns by two typical euchromatin histone markers, H3K4me3 and H3K27me3, for the great majority of homoeologous triad genes in A and B subgenomes were highly conserved between wild and domesticated tetraploid wheats and remained stable in the process of ploidy transitions from hexaploid to extracted tetraploid and then back to resynthesized hexaploid. However, a subset of genes was differentially modified during tetraploid and hexaploid wheat domestication and in response to ploidy transitions, and these genes were enriched for particular gene ontology (GO) terms. The extracted tetraploid wheat manifested higher overall histone modification levels than its hexaploid donor, and which were reversible and restored to normal levels in the resynthesized hexaploid. Further, while H3K4me3 marks were distally distributed along each chromosome and significantly correlated with subgenome expression as expected, H3K27me3 marks showed only a weak distal bias and did not show a significant correlation with gene expression.

CONCLUSIONS

Our results reveal overall high stability of histone modification patterns in the A and B subgenomes of polyploid wheat during domestication and in the process of ploidy transitions. However, modification levels of a subset of functionally relevant genes in the A and B genomes were trans-regulated by the D genome in hexaploid wheat.

Genome-wide identification and characterization of olfactory receptor genes in common carp (Cyprinus carpio)

The environment contains a large extent of chemical information, which could be detected as olfactory sense. Olfactory in vertebrates plays important roles on many aspects during life time, including localizing prey or food, avoiding predators, mating behavior and social communication. Considering the essential role of olfactory receptors in the specific recognition of diverse stimuli, understanding the evolutionary dynamics of olfactory receptors in teleost means a lot, especially in the allotetraploid common carp, who has undergone the fourth whole-genome duplication event. Here, we identified the whole set of olfactory receptor genes in representative teleosts and found a significant contraction in common carp when compared with other teleosts. Odorant receptor genes (OR) occupy the most among four groups of olfactory receptors, including 33 functional genes and 16 pseudogenes. Furthermore, 6 trace amine-associated receptor (TAAR) genes (including 1 pseudogene), 7 odorant-related-A receptor genes, and 10 olfactory C family receptor genes (including 3 pseudogenes) were identified in common carp. Phylogenetic and motif analysis were performed to illustrate the phylogenetic relationship and structural conservation of teleost olfactory receptors. Selection pressure analysis suggested that olfactory receptor groups in common carp were all under relaxed purifying-selection. Additionally, gene expression divergences for olfactory receptor genes were investigated during embryonic development stages of common carp. We aim to determine the abundance of common carp olfactory receptor genes, explore the evolutionary fate and expression dynamics, and provide some genomic clues for the evolution of polyploid olfactory after whole-genome duplication and for future studies of teleost olfactory.

Genome evolution in Oryza allopolyploids of various ages: Insights into the process of diploidization

The prevalence and recurrence of whole-genome duplication in plants and its major role in evolution have been well recognized. Despite great efforts, many aspects of genome evolution, particularly the temporal progression of genomic responses to allopolyploidy and the underlying mechanisms, remain poorly understood. The rice genus Oryza consists of both recently formed and older allopolyploid species, representing an attractive system for studying the genome evolution after allopolyploidy. In this study, through screening BAC libraries and sequencing and annotating the targeted BAC clones, we generated orthologous genomic sequences surrounding the DEP1 locus, a major grain yield QTL in cultivated rice, from four Oryza polyploids of various ages and their likely diploid genome donors or close relatives. Based on sequenced DEP1 region and published data from three other genomic regions, we investigated the temporal evolutionary dynamics of four polyploid genomes at both genetic and expression levels. In the recently formed BBCC polyploid, Oryza minuta, genome dominance was not observed and its short-term responses to allopolyploidy are mainly manifested as a high proportion of homoeologous gene pairs showing unequal expression. This could partly be explained by parental legacy, rewiring of divergent regulatory networks and epigenetic modulation. Moreover, we detected an ongoing diploidization process in this genus, and suggest that the expression divergence driven by changes of selective constraint probably plays a big role in the long-term diploidization. These findings add novel insights into our understanding of genome evolution after allopolyploidy, and could facilitate crop improvements through hybridization and polyploidization.

Unbiased Subgenome Evolution in Allotetraploid Species of Ephedra and Its Implications for the Evolution of Large Genomes in Gymnosperms

The evolutionary dynamics of polyploid genomes and consequences of polyploidy have been studied extensively in angiosperms but very rarely in gymnosperms. The gymnospermous genus Ephedra is characterized by a high frequency of polyploidy, and thus provides an ideal system to investigate the evolutionary mode of allopolyploid genomes and test whether subgenome dominance has occurred in gymnosperms. Here, we sequenced transcriptomes of two allotetraploid species of Ephedra and their putative diploid progenitors, identified expressed homeologs, and analyzed alternative splicing and homeolog expression based on PacBio Iso-Seq and Illumina RNA-seq data. We found that the two subgenomes of the allotetraploids had similar numbers of expressed homeologs, similar percentages of homeologs with dominant expression, and approximately equal numbers of isoforms with alternative splicing, showing an unbiased subgenome evolution as in a few polyploid angiosperms, with a divergence of the two subgenomes at ∼8 Ma. In addition, the nuclear DNA content of the allotetraploid species is almost equal to the sum of two putative progenitors, suggesting limited genome restructuring after allotetraploid speciation. The allopolyploid species of Ephedra might have undergone slow diploidization, and the unbiased subgenome evolution implies that the formation of large genomes in gymnosperms could be attributed to even and slow fractionation following polyploidization.

Expression Partitioning of Duplicate Genes at Single Cell Resolution in Arabidopsis Roots

Gene duplication is a key evolutionary phenomenon, prevalent in all organisms but particularly so in plants, where whole genome duplication (WGD; polyploidy) is a major force in genome evolution. Much effort has been expended in attempting to understand the evolution of duplicate genes, addressing such questions as why some paralog pairs rapidly return to single copy status whereas, in other pairs, both paralogs are retained and may diverge in expression pattern or function. The effect of a gene - its site of expression and thus the initial locus of its function - occurs at the level of a cell comprising a single cell type at a given state of the cell's development. Using Arabidopsis thaliana single cell transcriptomic data we categorized patterns of expression for 11,470 duplicate gene pairs across 36 cell clusters comprising nine cell types and their developmental states. Among these 11,470 pairs, 10,187 (88.8%) had at least one copy expressed in at least one of the 36 cell clusters. Pairs produced by WGD more often had both paralogs expressed in root cells than did pairs produced by small scale duplications. Three quarters of gene pairs expressed in the 36 cell clusters (7,608/10,187) showed extreme expression bias in at least one cluster, including 352 cases of reciprocal bias, a pattern consistent with expression subfunctionalization. More than twice as many pairs showed reciprocal expression bias between cell states than between cell types or between roots and leaves. A group of 33 gene pairs with reciprocal expression bias showed evidence of concerted divergence of gene networks in stele vs. epidermis. Pairs with both paralogs expressed without bias were less likely to have paralogs with divergent mutant phenotypes; such bias-free pairs showed evidence of preservation by maintenance of dosage balance. Overall, we found considerable evidence of shifts in gene expression following duplication, including in >80% of pairs encoding 7,653 genes expressed ubiquitously in all root cell types and states for which we inferred the polarity of change.

Maize BIG GRAIN1 homolog overexpression increases maize grain yield

The Zea Mays BIG GRAIN 1 HOMOLOG 1 (ZM-BG1H1) was ectopically expressed in maize. Elite commercial hybrid germplasm was yield tested in diverse field environment locations representing commercial models. Yield was measured in 101 tests across all 4 events, 26 locations over 2 years, for an average yield gain of 355 kg/ha (5.65 bu/ac) above control, with 83% tests broadly showing yield gains (range +2272 kg/ha to -1240 kg/ha), with seven tests gaining more than one metric ton per hectare. Plant and ear height were slightly elevated, and ear and tassel flowering time were delayed one day, but ASI was unchanged, and these traits did not correlate to yield gain. ZM-BG1H1 overexpression is associated with increased ear kernel row number and total ear kernel number and mass, but individual kernels trended slightly smaller and less dense. The ZM-BG1H1 protein is detected in the plasma membrane like rice OS-BG1. Five predominant native ZM-BG1H1 alleles exhibit little structural and expression variation compared to the large increased expression conferred by these ectopic alleles.

Expansion of MIR482/2118 by a class-II transposable element in cotton

Some plant microRNA (miRNA) families contain multiple members generating identical or highly similar mature miRNA variants. Mechanisms underlying the expansion of miRNA families remain elusive, although tandem and/or segmental duplications have been proposed. In this study of two tetraploid cottons, Gossypium hirsutum and Gossypium barbadense, and their extant diploid progenitors, Gossypium arboreum and Gossypium raimondii, we investigated the gain and loss of members of the miR482/2118 superfamily, which modulates the expression of nucleotide-binding site leucine-rich repeat (NBS-LRR) disease resistance genes. We found significant expansion of MIR482/2118d in G. barbadense, G. hirsutum and G. raimondii, but not in G. arboreum. Several newly expanded MIR482/2118d loci have mutated to produce different miR482/2118 variants with altered target-gene specificity. Based on detailed analysis of sequences flanking these MIR482/2118 loci, we found that this expansion of MIR482/2118d and its derivatives resulted from an initial capture of an MIR482/2118d by a class-II DNA transposable element (TE) in G. raimondii prior to the tetraploidization event, followed by transposition to new genomic locations in G. barbadense, G. hirsutum and G. raimondii. The 'GosTE' involved in the capture and proliferation of MIR482/2118d and its derivatives belongs to the PIF/Harbinger superfamily, generating a 3-bp target site duplication upon insertion at new locations. All orthologous MIR482/2118 loci in the two diploids were retained in the two tetraploids, but mutation(s) in miR482/2118 were observed across all four species as well as in different cultivars of both G. barbadense and G. hirsutum, suggesting a dynamic co-evolution of miR482/2118 and its NBS-LRR targets. Our results provide fresh insights into the mechanisms contributing to MIRNA proliferation and enrich our knowledge on TEs.

Distinct Expression and Methylation Patterns for Genes with Different Fates following a Single Whole-Genome Duplication in Flowering Plants

For most sequenced flowering plants, multiple whole-genome duplications (WGDs) are found. Duplicated genes following WGD often have different fates that can quickly disappear again, be retained for long(er) periods, or subsequently undergo small-scale duplications. However, how different expression, epigenetic regulation, and functional constraints are associated with these different gene fates following a WGD still requires further investigation due to successive WGDs in angiosperms complicating the gene trajectories. In this study, we investigate lotus (Nelumbo nucifera), an angiosperm with a single WGD during the K-pg boundary. Based on improved intraspecific-synteny identification by a chromosome-level assembly, transcriptome, and bisulfite sequencing, we explore not only the fundamental distinctions in genomic features, expression, and methylation patterns of genes with different fates after a WGD but also the factors that shape post-WGD expression divergence and expression bias between duplicates. We found that after a WGD genes that returned to single copies show the highest levels and breadth of expression, gene body methylation, and intron numbers, whereas the long-retained duplicates exhibit the highest degrees of protein-protein interactions and protein lengths and the lowest methylation in gene flanking regions. For those long-retained duplicate pairs, the degree of expression divergence correlates with their sequence divergence, degree in protein-protein interactions, and expression level, whereas their biases in expression level reflecting subgenome dominance are associated with the bias of subgenome fractionation. Overall, our study on the paleopolyploid nature of lotus highlights the impact of different functional constraints on gene fate and duplicate divergence following a single WGD in plant.

In Support of Winge's Theory of "Hybridization Followed by Chromosome Doubling"

Polyploidy-or chromosome doubling-plays a significant role in plant speciation and evolution. Much of the existing evidence indicates that fusion of unreduced (or 2n) gametes is the major pathway responsible for polyploid formation. In the early 1900s, a theory was put forward that the mechanism of "hybridization followed by chromosome doubling" would enable the survival and development of the hybrid zygote by providing each chromosome with a homolog with which to pair. However, to date there is only scant empirical evidence supporting this theory. In our previous study, interspecific-interploid crosses between the tetraploid Hylocereus megalanthus, as the female parent, and the diploid H. undatus, as the male parent, yielded only allopentaploids, allohexaploids, and 5x-and 6x-aneuploids instead of the expected allotriploids. No viable hybrids were obtained from the reciprocal cross. Since H. undatus underwent normal meiosis with regular pairing in the pollen mother cells and only reduced pollen grains were observed, the allohexaploids obtained supported the concept of "chromosome doubling." In this work, we report ploidy level, fruit morphology, and pollen viability and diameter in a group of putative hybrids obtained from an embryo rescue procedure following controlled H. megalanthus × H. undatus crosses, with the aim to elucidate, for the first time, the timing and developmental stage of the chromosome doubling. As in our previous report, no triploids were obtained, but tetraploids, pentaploids, hexaploids, and 5x- and 6x-aneuploids were found in the regenerated plants. The tetraploids exhibited the morphological features of the maternal parent and could not be considered true hybrids. Based on our previous studies, we can assume that the pentaploids were a result of a fertilization event between one unreduced (2n) female gamete from the tetraploid H. megalanthus and a normal (n) haploid male gamete from H. undatus. All the allohexaploids obtained from the embryo rescue technique where those that regenerated from fertilized ovules 10 days after pollination (at the pro-embryo stage), showing that the chromosome doubling event occurred at a very early development stage, i.e., at the zygote stage or shortly after zygote formation. These allohexaploids thus constitute empirical evidence of "hybridization followed by chromosome doubling."

The differential expression patterns of paralogs in response to stresses indicate expression and sequence divergences

BACKGROUND

Theoretically, paralogous genes generated through whole genome duplications should share identical expression levels due to their identical sequences and chromatin environments. However, functional divergences and expression differences have arisen due to selective pressures throughout evolution. A comprehensive investigation of the expression patterns of paralogous gene pairs in response to various stresses and a study of correlations between the expression levels and sequence divergences of the paralogs are needed.

RESULTS

In this study, we analyzed the expression patterns of paralogous genes under different types of stress and investigated the correlations between the expression levels and sequence divergences of the paralogs. We analyzed the differential expression patterns of the paralogs under four different types of stress (drought, cold, infection, and herbivory) and classified them into three main types according to their expression patterns. We then further analyzed the differential expression patterns under various degrees of stress and constructed corresponding co-expression networks of differentially expressed paralogs and transcription factors. Finally, we investigated the correlations between the expression levels and sequence divergences of the paralogs and identified positive correlations between expression level and sequence divergence. With regard to sequence divergence, we identified correlations between selective pressures and phylogenetic relationships.

CONCLUSIONS

These results shed light on differential expression patterns of paralogs in response to environmental stresses and are helpful for understanding the relationships between expression levels and sequences divergences.

Homoeologous gene expression and co-expression network analyses and evolutionary inference in allopolyploids

Polyploidy is a widespread phenomenon throughout eukaryotes. Due to the coexistence of duplicated genomes, polyploids offer unique challenges for estimating gene expression levels, which is essential for understanding the massive and various forms of transcriptomic responses accompanying polyploidy. Although previous studies have explored the bioinformatics of polyploid transcriptomic profiling, the causes and consequences of inaccurate quantification of transcripts from duplicated gene copies have not been addressed. Using transcriptomic data from the cotton genus (Gossypium) as an example, we present an analytical workflow to evaluate a variety of bioinformatic method choices at different stages of RNA-seq analysis, from homoeolog expression quantification to downstream analysis used to infer key phenomena of polyploid expression evolution. In general, EAGLE-RC and GSNAP-PolyCat outperform other quantification pipelines tested, and their derived expression dataset best represents the expected homoeolog expression and co-expression divergence. The performance of co-expression network analysis was less affected by homoeolog quantification than by network construction methods, where weighted networks outperformed binary networks. By examining the extent and consequences of homoeolog read ambiguity, we illuminate the potential artifacts that may affect our understanding of duplicate gene expression, including an overestimation of homoeolog co-regulation and the incorrect inference of subgenome asymmetry in network topology. Taken together, our work points to a set of reasonable practices that we hope are broadly applicable to the evolutionary exploration of polyploids.

Homeologs of Brassica SOC1, a central regulator of flowering time, are differentially regulated due to partitioning of evolutionarily conserved transcription factor binding sites in promoters

Evolution of Brassica genome post-polyploidization reveals asymmetrical genome fractionation and copy number variation. Herein, we describe the impact of promoter divergence among SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1 (SOC1) homeologs on expression and function in Brassica spp. SOC1, a regulated floral pathway integrator, is conserved as 3 redundant homeologs in diploid Brassicas. Even with high sequence identity within coding regions (92.8-100%), the spatio-temporal expression patterns of 9 SOC1 homologs in B. juncea and B. nigra indicates regulatory divergence. While LF and MF2 SOC1 homeologs are upregulated during floral transition, MF1 is barely expressed. Also, MF2 homeolog levels do not decline post-flowering, unlike LF. To investigate the underlying source of divergence, we analyzed the sequence and phylogeny of all reported (22) and isolated (21) upstream regions of Brassica SOC1. Full length upstream regions (4712-19189 bp) reveal 5 ubiquitously conserved ancestral Blocks, harboring binding sites of 18 TFs (TFBSs) characterized in Arabidopsis thaliana. The orthologs of these TFBSs are differentially conserved among Brassica SOC1 homeologs, imparting expression divergence. No crucial TFBSs are exclusively lost from LF_SOC1 promoter, while MF1_SOC1 has lost NF-Y binding site crucial for SOC1 activation by CONSTANS. MF2_SOC1 homeologs have lost important TFBSs (SEP3, AP1 and SMZ), responsible for SOC1 repression post-flowering. BjuAALF_SOC1 promoter (proximal 2 kb) shows ubiquitous reporter expression in B. juncea cv. Varuna transgenics, while BjuAAMF1_SOC1 promoter shows absence of reporter expression, validating the impact of TFBS divergence. Conservation of the original primary protein sequence is discovered in B. rapa homeologs (46) of 18 TFs. Co-regulation pattern of these TFs appeared similar for B. rapa LF and MF2 SOC1 homeologs; MF1 shows significant variation. Strong regulatory association is recorded for AP1, AP2, SEP3, FLC and CONSTANS/NF-Y, highlighting their importance in homeolog-specific SOC1 regulation. Correlation of B. juncea AP1, AP2 and FLC expression with SOC1 homeologs also complies with the TFBS differences. We thus conclude that redundant SOC1 loci contribute differentially to cumulative expression of SOC1 due to divergent selection of ancestral TFBSs.

Tissue-specific gene expression and protein abundance patterns are associated with fractionation bias in maize

BACKGROUND

Maize experienced a whole-genome duplication event approximately 5 to 12 million years ago. Because this event occurred after speciation from sorghum, the pre-duplication subgenomes can be partially reconstructed by mapping syntenic regions to the sorghum chromosomes. During evolution, maize has had uneven gene loss between each ancient subgenome. Fractionation and divergence between these genomes continue today, constantly changing genetic make-up and phenotypes and influencing agronomic traits.

RESULTS

Here we regenerate the subgenome reconstructions for the most recent maize reference genome assembly. Based on both expression and abundance data for homeologous gene pairs across multiple tissues, we observed functional divergence of genes across subgenomes. Although the genes in the larger maize subgenome are often expressing more highly than their homeologs in the smaller subgenome, we observed cases where homeolog expression dominance switches in different tissues. We demonstrate for the first time that protein abundances are higher in the larger subgenome, but they also show tissue-specific dominance, a pattern similar to RNA expression dominance. We also find that pollen expression is uniquely decoupled from protein abundance.

CONCLUSION

Our study shows that the larger subgenome has a greater range of functional assignments and that there is a relative lack of overlap between the subgenomes in terms of gene functions than would be suggested by similar patterns of gene expression and protein abundance. Our study also revealed that some reactions are catalyzed uniquely by the larger and smaller subgenomes. The tissue-specific, nonequivalent expression-level dominance pattern observed here implies a change in regulatory control which favors differentiated selective pressure on the retained duplicates leading to eventual change in gene functions.

Generation of Transgenic Self-Incompatible Arabidopsis thaliana Shows a Genus-Specific Preference for Self-Incompatibility Genes

Brassicaceae species employ both self-compatibility and self-incompatibility systems to regulate post-pollination events. Arabidopsis halleri is strictly self-incompatible, while the closely related Arabidopsis thaliana has transitioned to self-compatibility with the loss of functional S-locus genes during evolution. The downstream signaling protein, ARC1, is also required for the self-incompatibility response in some Arabidopsis and Brassica species, and its gene is deleted in the A. thaliana genome. In this study, we attempted to reconstitute the SCR-SRK-ARC1 signaling pathway to restore self-incompatibility in A. thaliana using genes from A. halleri and B. napus, respectively. Several of the transgenic A. thaliana lines expressing the A. halleriSCR₁₃-SRK₁₃-ARC1 transgenes displayed self-incompatibility, while all the transgenic A. thaliana lines expressing the B. napusSCR₁-SRK₁-ARC1 transgenes failed to show any self-pollen rejection. Furthermore, our results showed that the intensity of the self-incompatibility response in transgenic A. thaliana plants was not associated with the expression levels of the transgenes. Thus, this suggests that there are differences between the Arabidopsis and Brassica self-incompatibility signaling pathways, which perhaps points to the existence of other factors downstream of B. napusSRK that are absent in Arabidopsis species.

Comparative genome/transcriptome analysis probes Boraginales' phylogenetic position, WGDs in Boraginales, and key enzyme genes in the alkannin/shikonin core pathway

Boraginales (the forget-me-not order) is a core group within the lamiids clade. However, until now, no genome from Boraginales has been reported, and published transcriptomes are also rare. Here, we report the first Boraginales species de novo genome (i.e. Echium plantagineum genome) and seven other Boraginales species transcriptomes to probe three issues: (i) Boraginales' phylogenetic position within the lamiids clade; (ii) potential whole genome duplications (WGDs) in Boraginales; and (iii) candidate key enzyme genes in the alkannin/shikonin core pathway. The results showed that: (i) Boraginales was most probably closer to the Solanales/Gentianales clade than the Lamiales clade, at least based on the single-copy orthologous genes from genome/transcriptome data; (ii) after the gamma (γ) event, Boraginaceae (classified into the Boraginales I clade) probably underwent at least two rounds of WGD, whereas Heliotropiaceae and Ehretiaceae (classified into the Boraginales II clade) probably underwent only one round of WGD; and (iii) several candidate key enzyme genes in the alkannin/shikonin core pathway were inferred, e.g. genes corresponding to geranyl cyclase, naphthol hydroxylase and O-acyl transferase.

Evolutionary history and structure of nuclear matrix constituent proteins, the plant analogues of lamins

Nuclear matrix constituent proteins (NMCPs), the structural components of the plant lamina, are considered to be the analogues of lamins in plants based on numerous structural and functional similarities. Current phylogenetic knowledge suggests that, in contrast to lamins, which are widely distributed in eukaryotes, NMCPs are taxonomically restricted to Streptophyta. At present, most information about NMCPs comes from angiosperms, and virtually no data are available from more ancestral groups. In angiosperms, the NMCP family comprises two phylogenetic groups, NMCP1 and NMCP2, which evolved from the NMCP1 and NMCP2 progenitor genes. Based on sequence conservation and the presence of NMCP-specific domains, we determined the structure and number of NMCP genes present in different Streptophyta clades. We analysed 91 species of embryophytes and report additional NMCP sequences from mosses, liverworts, clubmosses, horsetail, ferns, gymnosperms, and Charophyta algae. Our results confirm an origin of NMCPs in Charophyta (the earliest diverging group of Streptophyta), resolve the number and structure of NMCPs in the different clades, and propose the emergence of additional NMCP homologues by whole-genome duplication events. Immunofluorescence microscopy demonstrated localization of a basal NMCP from the moss Physcomitrella patens at the nuclear envelope, suggesting a functional conservation for basal and more evolved NMCPs.

Paleogenomics: reconstruction of plant evolutionary trajectories from modern and ancient DNA

How contemporary plant genomes originated and evolved is a fascinating question. One approach uses reference genomes from extant species to reconstruct the sequence and structure of their common ancestors over deep timescales. A second approach focuses on the direct identification of genomic changes at a shorter timescale by sequencing ancient DNA preserved in subfossil remains. Merged within the nascent field of paleogenomics, these complementary approaches provide insights into the evolutionary forces that shaped the organization and regulation of modern genomes and open novel perspectives in fostering genetic gain in breeding programs and establishing tools to predict future population changes in response to anthropogenic pressure and global warming.

Identification and Characterization of the EXO70 Gene Family in Polyploid Wheat and Related Species

The EXO70 gene family is involved in different biological processes in plants, ranging from plant polar growth to plant immunity. To date, analysis of the EXO70 gene family has been limited in Triticeae species, e.g., hexaploidy Triticum aestivum and its ancestral/related species. By in silico analysis of multiple Triticeae sequence databases, a total of 200 EXO70 members were identified. By homologue cloning approaches, 15 full-length cDNA of EXO70s were cloned from diploid Haynaldia villosa. Phylogenetic relationship analysis of 215 EXO70 members classified them into three groups (EXO70.1, EXO70.2, and EXO70.3) and nine subgroups (EXO70A to EXO70I). The distribution of most EXO70 genes among different species/sub-genomes were collinear, implying their orthologous relationship. The EXO70A subgroup has the most introns (at least five introns), while the remaining seven subgroups have only one intron on average. The expression profiling of EXO70 genes from wheat revealed that 40 wheat EXO70 genes were expressed in at least one tissue (leaf, stem, or root), of which 25 wheat EXO70 genes were in response to at least one biotic stress (stripe rust or powdery mildew) or abiotic stress (drought or heat). Subcellular localization analysis showed that ten EXO70-V proteins had distinct plasma membrane localization, EXO70I1-V showed a distinctive spotted pattern on the membrane. The 15 EXO70-V genes were differentially expressed in three tissue. Apart from EXO70D2-V, the remaining EXO70-V genes were in response to at least one stress (flg22, chitin, powdery mildew, drought, NaCl, heat, or cold) or phytohormones (salicylic acid, methyl jasmonate, ethephon, or abscisic acid) and hydrogen peroxide treatments. This research provides a genome-wide glimpse of the Triticeae EXO70 gene family and those up- or downregulated genes require further validation of their biological roles in response to biotic/abiotic stresses.

Natural variation in Brassica FT homeologs influences multiple agronomic traits including flowering time, silique shape, oil profile, stomatal morphology and plant height in B. juncea

Natural structural variants of regulatory proteins causing quantitative phenotypic consequences have not been reported in plants. Herein, we show that 28 natural structural variants of FT homeologs, isolated from 6 species of Brassica, differ with respect to amino-acid substitutions in regions critical for interactions with FD and represent two evolutionarily distinct categories. Analysis of structural models of selected candidates from Brassica juncea (BjuFT_AAMF1) and Brassica napus (BnaFT_CCLF) predicted stronger binding between BjuFT and Arabidopsis thaliana FD. Over-expression of BjuFT and BnaFT in wild type and ft-10 mutant backgrounds of Arabidopsis validated higher potency of BjuFT in triggering floral transition. Analysis of gain-of-function and artificial miRNA mediated silenced lines of B. juncea implicated Brassica FT in multiple agronomic traits beyond flowering, consistent with a pleiotropic effect. Several dependent and independent traits such as lateral branching, silique shape, seed size, oil-profile, stomatal morphology and plant height were found altered in mutant lines. Enhanced FT levels caused early flowering, which in turn was positively correlated to a higher proportion of desirable fatty acids (PUFA). However, higher FT levels also resulted in altered silique shape and reduced seed size, suggesting trait trade-offs. Modulation of FT levels for achieving optimal balance of trait values and parsing pair-wise interactions among a reportoire of regulatory protein homeologs in polyploid genomes are indeed future areas of crop research.