Coordinated Functional Divergence of Genes after Genome Duplication in Arabidopsis thaliana

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

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

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.

Genetic Channelization Mechanism of Four Chalcone Isomerase Homologous Genes for Synergistic Resistance to Fusarium wilt in Gossypium barbadense L

Duplication events occur very frequently during plant evolution. The genes in the duplicated pathway or network can evolve new functions through neofunctionalization and subfunctionalization. Flavonoids are secondary metabolites involved in plant development and defense. Our previous transcriptomic analysis of F6 recombinant inbred lines (RILs) and the parent lines after Fusarium oxysporum f. sp. vasinfectum (Fov) infection showed that CHI genes have important functions in cotton. However, there are few reports on the possible neofunctionalization differences of CHI family paralogous genes involved in Fusarium wilt resistance in cotton. In this study, the resistance to Fusarium wilt, expression of metabolic pathway-related genes, metabolite content, endogenous hormone content, reactive oxygen species (ROS) content and subcellular localization of four paralogous CHI family genes in cotton were investigated. The results show that the four paralogous CHI family genes may play a synergistic role in Fusarium wilt resistance. These results revealed a genetic channelization mechanism that can regulate the metabolic flux homeostasis of flavonoids under the mediation of endogenous salicylic acid (SA) and methyl jasmonate (MeJA) via the four paralogous CHI genes, thereby achieving disease resistance. Our study provides a theoretical basis for studying the evolutionary patterns of homologous plant genes and using homologous genes for molecular breeding.

Predicting gene expression divergence between single-copy orthologs in two species

Predicting gene expression divergence is integral to understanding the emergence of new biological functions and associated traits. Whereas several sophisticated methods have been developed for this task, their applications are either limited to duplicate genes or require expression data from more than two species. Thus, here we present PiXi, the first machine learning framework for predicting gene expression divergence between single-copy orthologs in two species. PiXi models gene expression evolution as an Ornstein-Uhlenbeck process, and overlays this model with multi-layer neural network, random forest, and support vector machine architectures for making predictions. It outputs the predicted class "conserved" or "diverged" for each pair of orthologs, as well as their predicted expression optima in the two species. We show that PiXi has high power and accuracy in predicting gene expression divergence between single-copy orthologs, as well as high accuracy and precision in estimating their expression optima in the two species, across a wide range of evolutionary scenarios, with the globally best performance achieved by a multi-layer neural network. Moreover, application of our best performing PiXi predictor to empirical gene expression data from single-copy orthologs residing at different loci in two species of Drosophila reveals that approximately 23% underwent expression divergence after positional relocation. Further analysis shows that several of these "diverged" genes are involved in the electron transport chain of the mitochondrial membrane, suggesting that new chromatin environments may impact energy production in Drosophila. Thus, by providing a toolkit for predicting gene expression divergence between single-copy orthologs in two species, PiXi can shed light on the origins of novel phenotypes across diverse biological processes and study systems.

Enrichment of intrinsically disordered residues in ohnologs facilitates abiotic stress resilience in Brassica rapa

Arabidopsis thaliana and Brassica rapa are in the same evolutionary lineage, although the latter experienced an additional whole genome triplication event. Therefore, it would be intriguing to investigate the traits that gene duplication imposes to mediate plant stress tolerance. Here, we noticed that B. rapa abiotic stress resistance (ASR) genes which code at least one stress responsive domain have a significantly higher number of paralogs than A. thaliana. Analysing the disordered content of the ASR genes in both species, we found that intrinsically disordered residues (IDR) are specifically enriched in whole genome duplication (WGD) derived paralogs. Subsequently, domain similarity analysis between WGD pairs of both species has revealed that majority of WGD pairs in B. rapa did not share domains with each other. Furthermore, domain enrichment analysis has shown that B. rapa paralogs contain 36 distinct stress responsive enriched domains, significantly higher than A. thaliana paralogs. Next, we performed MSA to investigate the domain conservation between orthologs and ohnologs pairs, we found that 80.13% of B. rapa ohnologs acquire new domains, depicting the fact that ohnologs play a significant role in stress-related behaviours. The average IDR content of the ohnologs enriching new domains after gene duplication in B. rapa (0.19), is also significantly higher than A. thaliana (0.04). Interestingly, we also found that all of these attributes i.e., exhibiting higher number of WGD paralogs and enhancement of IDR in ASR genes of B. rapa compared to A. thaliana is exclusive for ASR genes only. No such significant differences were observed in randomly selected non-ASR genes between the two species. Together these results provide strong support for the hypothesis that augmentation of IDR content followed by a whole genome duplication event imposes the stress resistance potentiality in B. rapa. This research will shed light on the mechanism of how B. rapa is able to successfully adapt to stress over the evolutionary timescale.

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

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

Evolutionary analysis of the LORELEI gene family in plants reveals regulatory subfunctionalization

A signaling complex comprising members of the LORELEI (LRE)-LIKE GPI-anchored protein (LLG) and Catharanthus roseus RECEPTOR-LIKE KINASE 1-LIKE (CrRLK1L) families perceive RAPID ALKALINIZATION FACTOR (RALF) peptides and regulate growth, reproduction, immunity, and stress responses in Arabidopsis (Arabidopsis thaliana). Genes encoding these proteins are members of multigene families in most angiosperms and could generate thousands of signaling complex variants. However, the links between expansion of these gene families and the functional diversification of this critical signaling complex as well as the evolutionary factors underlying the maintenance of gene duplicates remain unknown. Here, we investigated LLG gene family evolution by sampling land plant genomes and explored the function and expression of angiosperm LLGs. We found that LLG diversity within major land plant lineages is primarily due to lineage-specific duplication events, and that these duplications occurred both early in the history of these lineages and more recently. Our complementation and expression analyses showed that expression divergence (i.e. regulatory subfunctionalization), rather than functional divergence, explains the retention of LLG paralogs. Interestingly, all but one monocot and all eudicot species examined had an LLG copy with preferential expression in male reproductive tissues, while the other duplicate copies showed highest levels of expression in female or vegetative tissues. The single LLG copy in Amborella trichopoda is expressed vastly higher in male compared to in female reproductive or vegetative tissues. We propose that expression divergence plays an important role in retention of LLG duplicates in angiosperms.

De novo genes in Arachis hypogaea cv. Tifrunner: systematic identification, molecular evolution, and potential contributions to cultivated peanut

De novo genes are derived from non-coding sequences, and they can play essential roles in organisms. Cultivated peanut (Arachis hypogaea) is a major oil and protein crop derived from a cross between Arachis duranensis and Arachis ipaensis. However, few de novo genes have been documented in Arachis. Here, we identified 381 de novo genes in A. hypogaea cv. Tifrunner based on comparison with five closely related Arachis species. There are distinct differences in gene expression patterns and gene structures between conserved and de novo genes. The identified de novo genes originated from ancestral sequence regions associated with metabolic and biosynthetic processes, and they were subsequently integrated into existing regulatory networks. De novo paralogs and homoeologs were identified in A. hypogaea cv. Tifrunner. De novo paralogs and homoeologs with conserved expression have mismatching cis-acting elements under normal growth conditions. De novo genes potentially have pluripotent functions in responses to biotic stresses as well as in growth and development based on quantitative trait locus data. This work provides a foundation for future research examining gene birth processes and gene function in Arachis and related taxa.

Diversification of SEC15a and SEC15b isoforms of an exocyst subunit in seed plants is manifested in their specific roles in Arabidopsis sporophyte and male gametophyte

The exocyst complex is an octameric evolutionarily conserved tethering complex engaged in the regulation of polarized secretion in eukaryotic cells. Here, we focus on the systematic comparison of two isoforms of the SEC15 exocyst subunit, SEC15a and SEC15b. We infer that SEC15 gene duplication and diversification occurred in the common ancestor of seed plants (Spermatophytes). In Arabidopsis, SEC15a represents the main SEC15 isoform in the male gametophyte, and localizes to the pollen tube tip at the plasma membrane. Although pollen tubes of sec15a mutants are impaired, sporophytes show no phenotypic deviations. Conversely, SEC15b is the dominant isoform in the sporophyte and localizes to the plasma membrane in root and leaf cells. Loss-of-function sec15b mutants exhibit retarded elongation of hypocotyls and root hairs, a loss of apical dominance, dwarfed plant stature and reduced seed coat mucilage formation. Surprisingly, the sec15b mutants also exhibit compromised pollen tube elongation in vitro, despite its very low expression in pollen, suggesting a non-redundant role for the SEC15b isoform there. In pollen tubes, SEC15b localizes to distinct cytoplasmic structures. Reciprocally to this, SEC15a also functions in the sporophyte, where it accumulates at plasmodesmata. Importantly, although overexpressed SEC15a could fully complement the sec15b phenotypic deviations in the sporophyte, the pollen-specific overexpression of SEC15b was unable to fully compensate for the loss of SEC15a function in pollen. We conclude that the SEC15a and SEC15b isoforms evolved in seed plants, with SEC15a functioning mostly in pollen and SEC15b functioning mostly in the sporophyte.

Adaptive Variation in Homolog Number Within Transcript Families Promotes Expression Divergence in Reef-Building Coral

Gene expression, especially in multi-species experiments, is used to gain insight into the genetic basis of how organisms adapt and respond to changing environments. However, evolutionary processes which can influence gene expression patterns between species such as the presence of paralogs which arise from gene duplication events are rarely accounted for. Paralogous transcripts can alter the transcriptional output of a gene and thus exclusion of these transcripts can obscure important biological differences between species. To address this issue, we investigated how differences in transcript family size is associated with divergent gene expression patterns in five species of Caribbean reef-building corals. We demonstrate that transcript families that are rapidly evolving in terms of size have increased levels of expression divergence. Additionally, these rapidly evolving transcript families are enriched for multiple biological processes, with genes involved in the coral innate immune system demonstrating pronounced variation in homolog number between species. Overall, this investigation demonstrates the importance of incorporating paralogous transcripts when comparing gene expression across species by influencing both transcriptional output and the number of transcripts within biological processes. As this investigation was based on transcriptome assemblies, additional insights into the relationship between gene duplications and expression patterns will likely emergence once more genome assemblies are available for study.

Chromosome-scale assembly and evolution of the tetraploid Salvia splendens (Lamiaceae) genome

Polyploidization plays a key role in plant evolution, but the forces driving the fate of homoeologs in polyploid genomes, i.e., paralogs resulting from a whole-genome duplication (WGD) event, remain to be elucidated. Here, we present a chromosome-scale genome assembly of tetraploid scarlet sage (Salvia splendens), one of the most diverse ornamental plants. We found evidence for three WGD events following an older WGD event shared by most eudicots (the γ event). A comprehensive, spatiotemporal, genome-wide analysis of homoeologs from the most recent WGD unveiled expression asymmetries, which could be associated with genomic rearrangements, transposable element proximity discrepancies, coding sequence variation, selection pressure, and transcription factor binding site differences. The observed differences between homoeologs may reflect the first step toward sub- and/or neofunctionalization. This assembly provides a powerful tool for understanding WGD and gene and genome evolution and is useful in developing functional genomics and genetic engineering strategies for scarlet sage and other Lamiaceae species.

A small peptide inhibits siRNA amplification in plants by mediating autophagic degradation of SGS3/RDR6 bodies

Selective autophagy mediates specific degradation of unwanted cytoplasmic components to maintain cellular homeostasis. The suppressor of gene silencing 3 (SGS3) and RNA-dependent RNA polymerase 6 (RDR6)-formed bodies (SGS3/RDR6 bodies) are essential for siRNA amplification in planta. However, whether autophagy receptors regulate selective turnover of SGS3/RDR6 bodies is unknown. By analyzing the transcriptomic response to virus infection in Arabidopsis, we identified a virus-induced small peptide 1 (VISP1) composed of 71 amino acids, which harbor a ubiquitin-interacting motif that mediates interaction with autophagy-related protein 8. Overexpression of VISP1 induced selective autophagy and compromised antiviral immunity by inhibiting SGS3/RDR6-dependent viral siRNA amplification, whereas visp1 mutants exhibited opposite effects. Biochemistry assays demonstrate that VISP1 interacted with SGS3 and mediated autophagic degradation of SGS3/RDR6 bodies. Further analyses revealed that overexpression of VISP1, mimicking the sgs3 mutant, impaired biogenesis of endogenous trans-acting siRNAs and up-regulated their targets. Collectively, we propose that VISP1 is a small peptide receptor functioning in the crosstalk between selective autophagy and RNA silencing.

Molecular mechanisms involved in functional macroevolution of plant transcription factors

Transcription factors (TFs) are key components of the transcriptional regulation machinery. In plants, they accompanied the evolution from unicellular aquatic algae to complex flowering plants that dominate the land environment. The adaptations of the body plan and physiological responses required changes in the biological functions of TFs. Some ancestral gene regulatory networks are highly conserved, while others evolved more recently and only exist in particular lineages. The recent emergence of novel model organisms provided the opportunity for comparative studies, producing new insights to infer these evolutionary trajectories. In this review, we comprehensively revisit the recent literature on TFs of nonseed plants and algae, focusing on the molecular mechanisms driving their functional evolution. We discuss the particular contribution of changes in DNA-binding specificity, protein-protein interactions and cis-regulatory elements to gene regulatory networks. Current advances have shown that these evolutionary processes were shaped by changes in TF expression pattern, not through great innovation in TF protein sequences. We propose that the role of TFs associated with environmental and developmental regulation was unevenly conserved during land plant evolution.

Comparative regulomics supports pervasive selection on gene dosage following whole genome duplication

BACKGROUND

Whole genome duplication (WGD) events have played a major role in eukaryotic genome evolution, but the consequence of these extreme events in adaptive genome evolution is still not well understood. To address this knowledge gap, we used a comparative phylogenetic model and transcriptomic data from seven species to infer selection on gene expression in duplicated genes (ohnologs) following the salmonid WGD 80-100 million years ago.

RESULTS

We find rare cases of tissue-specific expression evolution but pervasive expression evolution affecting many tissues, reflecting strong selection on maintenance of genome stability following genome doubling. Ohnolog expression levels have evolved mostly asymmetrically, by diverting one ohnolog copy down a path towards lower expression and possible pseudogenization. Loss of expression in one ohnolog is significantly associated with transposable element insertions in promoters and likely driven by selection on gene dosage including selection on stoichiometric balance. We also find symmetric expression shifts, and these are associated with genes under strong evolutionary constraints such as ribosome subunit genes. This possibly reflects selection operating to achieve a gene dose reduction while avoiding accumulation of "toxic mutations". Mechanistically, ohnolog regulatory divergence is dictated by the number of bound transcription factors in promoters, with transposable elements being one likely source of novel binding sites driving tissue-specific gains in expression.

CONCLUSIONS

Our results imply pervasive adaptive expression evolution following WGD to overcome the immediate challenges posed by genome doubling and to exploit the long-term genetic opportunities for novel phenotype evolution.

Polyploidy: an evolutionary and ecological force in stressful times

Polyploidy has been hypothesized to be both an evolutionary dead-end and a source for evolutionary innovation and species diversification. Although polyploid organisms, especially plants, abound, the apparent nonrandom long-term establishment of genome duplications suggests a link with environmental conditions. Whole-genome duplications seem to correlate with periods of extinction or global change, while polyploids often thrive in harsh or disturbed environments. Evidence is also accumulating that biotic interactions, for instance, with pathogens or mutualists, affect polyploids differently than nonpolyploids. Here, we review recent findings and insights on the effect of both abiotic and biotic stress on polyploids versus nonpolyploids and propose that stress response in general is an important and even determining factor in the establishment and success of polyploidy.

LAZY Gene Family in Plant Gravitropism

Adapting to the omnipresent gravitational field was a fundamental basis driving the flourishing of terrestrial plants on the Earth. Plants have evolved a remarkable capability that not only allows them to live and develop within the Earth's gravity field, but it also enables them to use the gravity vector to guide the growth of roots and shoots, in a process known as gravitropism. Triggered by gravistimulation, plant gravitropism is a highly complex, multistep process that requires many organelles and players to function in an intricate coordinated way. Although this process has been studied for several 100 years, much remains unclear, particularly the early events that trigger the relocation of the auxin efflux carrier PIN-FORMED (PIN) proteins, which presumably leads to the asymmetrical redistribution of auxin. In the past decade, the LAZY gene family has been identified as a crucial player that ensures the proper redistribution of auxin and a normal tropic response for both roots and shoots upon gravistimulation. LAZY proteins appear to be participating in the early steps of gravity signaling, as the mutation of LAZY genes consistently leads to altered auxin redistribution in multiple plant species. The identification and characterization of the LAZY gene family have significantly advanced our understanding of plant gravitropism, and opened new frontiers of investigation into the novel molecular details of the early events of gravitropism. Here we review current knowledge of the LAZY gene family and the mechanism modulated by LAZY proteins for controlling both roots and shoots gravitropism. We also discuss the evolutionary significance and conservation of the LAZY gene family in plants.

EXO70A2 Is Critical for Exocyst Complex Function in Pollen Development

Pollen development, pollen grain germination, and pollen tube elongation are crucial biological processes in angiosperm plants that need precise regulation to deliver sperm cells to ovules for fertilization. Highly polarized secretion at a growing pollen tube tip requires the exocyst tethering complex responsible for specific targeting of secretory vesicles to the plasma membrane. Here, we demonstrate that Arabidopsis (Arabidopsis thaliana) EXO70A2 (At5g52340) is the main exocyst EXO70 isoform in the male gametophyte, governing the conventional secretory function of the exocyst, analogous to EXO70A1 (At5g03540) in the sporophyte. Our analysis of a CRISPR-generated exo70a2 mutant revealed that EXO70A2 is essential for efficient pollen maturation, pollen grain germination, and pollen tube growth. GFP:EXO70A2 was localized to the nucleus and cytoplasm in developing pollen grains and later to the apical domain in growing pollen tube tips characterized by intensive exocytosis. Moreover, EXO70A2 could substitute for EXO70A1 function in the sporophyte, but not vice versa, indicating partial functional redundancy of these two closely related isoforms and higher specificity of EXO70A2 for pollen development-related processes. Phylogenetic analysis revealed that the ancient duplication of EXO70A, one of which is always highly expressed in pollen, occurred independently in monocots and dicots. In summary, EXO70A2 is a crucial component of the exocyst complex in Arabidopsis pollen that is required for efficient plant sexual reproduction.

Exploring the complexity of soybean (Glycine max) transcriptional regulation using global gene co-expression networks

MAIN CONCLUSION

We report a soybean gene co-expression network built with data from 1284 RNA-Seq experiments, which was used to identify important regulators, modules and to elucidate the fates of gene duplicates. Soybean (Glycine max (L.) Merr.) is one of the most important crops worldwide, constituting a major source of protein and edible oil. Gene co-expression networks (GCN) have been extensively used to study transcriptional regulation and evolution of genes and genomes. Here, we report a soybean GCN using 1284 publicly available RNA-Seq samples from 15 distinct tissues. We found modules that are differentially regulated in specific tissues, comprising processes such as photosynthesis, gluconeogenesis, lignin metabolism, and response to biotic stress. We identified transcription factors among intramodular hubs, which probably integrate different pathways and shape the transcriptional landscape in different conditions. The top hubs for each module tend to encode proteins with critical roles, such as succinate dehydrogenase and RNA polymerase subunits. Importantly, gene essentiality was strongly correlated with degree centrality and essential hubs were enriched in genes involved in nucleic acids metabolism and regulation of cell replication. Using a guilt-by-association approach, we predicted functions for 93 of 106 hubs without functional description in soybean. Most of the duplicated genes had different transcriptional profiles, supporting their functional divergence, although paralogs originating from whole-genome duplications (WGD) are more often preserved in the same module than those from other mechanisms. Together, our results highlight the importance of GCN analysis in unraveling key functional aspects of the soybean genome, in particular those associated with hub genes and WGD events.

Genomics of Evolutionary Novelty in Hybrids and Polyploids

It has long been recognized that hybridization and polyploidy are prominent processes in plant evolution. Although classically recognized as significant in speciation and adaptation, recognition of the importance of interspecific gene flow has dramatically increased during the genomics era, concomitant with an unending flood of empirical examples, with or without genome doubling. Interspecific gene flow is thus increasingly thought to lead to evolutionary innovation and diversification, via adaptive introgression, homoploid hybrid speciation and allopolyploid speciation. Less well understood, however, are the suite of genetic and genomic mechanisms set in motion by the merger of differentiated genomes, and the temporal scale over which recombinational complexity mediated by gene flow might be expressed and exposed to natural selection. We focus on these issues here, considering the types of molecular genetic and genomic processes that might be set in motion by the saltational event of genome merger between two diverged species, either with or without genome doubling, and how these various processes can contribute to novel phenotypes. Genetic mechanisms include the infusion of new alleles and the genesis of novel structural variation including translocations and inversions, homoeologous exchanges, transposable element mobilization and novel insertional effects, presence-absence variation and copy number variation. Polyploidy generates massive transcriptomic and regulatory alteration, presumably set in motion by disrupted stoichiometries of regulatory factors, small RNAs and other genome interactions that cascade from single-gene expression change up through entire networks of transformed regulatory modules. We highlight both these novel combinatorial possibilities and the range of temporal scales over which such complexity might be generated, and thus exposed to natural selection and drift.

Multi-faceted analysis provides little evidence for recurrent whole-genome duplications during hexapod evolution

BACKGROUND

Gene duplication events play an important role in the evolution and adaptation of organisms. Duplicated genes can arise through different mechanisms, including whole-genome duplications (WGDs). Recently, WGD was suggested to be an important driver of evolution, also in hexapod animals.

RESULTS

Here, we analyzed 20 high-quality hexapod genomes using whole-paranome distributions of estimated synonymous distances (KS), patterns of within-genome co-linearity, and phylogenomic gene tree-species tree reconciliation methods. We observe an abundance of gene duplicates in the majority of these hexapod genomes, yet we find little evidence for WGD. The majority of gene duplicates seem to have originated through small-scale gene duplication processes. We did detect segmental duplications in six genomes, but these lacked the within-genome co-linearity signature typically associated with WGD, and the age of these duplications did not coincide with particular peaks in KS distributions. Furthermore, statistical gene tree-species tree reconciliation failed to support all but one of the previously hypothesized WGDs.

CONCLUSIONS

Our analyses therefore provide very limited evidence for WGD having played a significant role in the evolution of hexapods and suggest that alternative mechanisms drive gene duplication events in this group of animals. For instance, we propose that, along with small-scale gene duplication events, episodes of increased transposable element activity could have been an important source for gene duplicates in hexapods.

Dual-targeting by CRISPR/Cas9 leads to efficient point mutagenesis but only rare targeted deletions in the rice genome

The present study investigated the efficiency of CRISPR/Cas9 in creating genomic deletions as the basis of its application in removing selection marker genes or the intergenic regions. Three loci, representing a transgene and two rice genes, were targeted at two sites each, in separate experiments, and the deletion of the defined fragments was investigated by PCR and sequencing. Genomic deletions were found at a low rate among the transformed callus lines that could be isolated, cultured, and regenerated into plants harboring the deletion. However, randomly regenerated plants showed mixed genomic effects, and generally did not harbor heritable genomic deletions. To determine whether point mutations occurred at each targeted site, a total of 114 plants consisting of primary transgenic lines and their progeny were analyzed. Ninety-three plants showed targeting, 60 of which were targeted at both sites. The presence of point mutations at both sites was correlated with the guide RNA efficiency. In summary, genomic deletions through dual-targeting by the paired-guide RNAs were generally observed in callus, while de novo point mutations at one or both sites occurred at high rates in transgenic plants and their progeny, generating a variety of insertion–deletions or single-nucleotide variations. In this study, point mutations were exceedingly favored over genomic deletions; therefore, for the recovery of plant lines harboring targeted deletions, identifying early transformed clones harboring the deletions, and isolating them for plant regeneration is recommended.

Lineage-Specific Expression Divergence in Grasses Is Associated with Male Reproduction, Host-Pathogen Defense, and Domestication

Poaceae (grasses) is an agriculturally important and widely distributed family of plants with extraordinary phenotypic diversity, much of which was generated under recent lineage-specific evolution. Yet, little is known about the genes and functional modules involved in the lineage-specific divergence of grasses. Here, I address this question on a genome-wide scale by applying a novel branch-based statistic of lineage-specific expression divergence, LED, to RNA-seq data from nine tissues of the wild grass Brachypodium distachyon and its domesticated relatives Oryza sativa japonica (rice) and Sorghum bicolor (sorghum). I find that LED is generally smallest in B. distachyon and largest in O. sativa japonica, which underwent domestication earlier than S. bicolor, supporting the hypothesis that domestication may increase the rate of lineage-specific expression divergence in grasses. Moreover, in all three species, LED is positively correlated with protein-coding sequence divergence and tissue specificity, and negatively correlated with network connectivity. Further analysis reveals that genes with large LED are often primarily expressed in anther, implicating lineage-specific expression divergence in the evolution of male reproductive phenotypes. Gene ontology enrichment analysis also identifies an overrepresentation of terms related to male reproduction in the two domesticated grasses, as well as to those involved in host-pathogen defense in all three species. Last, examinations of genes with the largest LED reveal that their lineage-specific expression divergence may have contributed to antimicrobial functions in B. distachyon, to enhanced adaptation and yield during domestication in O. sativa japonica, and to defense against a widespread and devastating fungal pathogen in S. bicolor. Together, these findings suggest that lineage-specific expression divergence in grasses may increase under domestication and preferentially target rapidly evolving genes involved in male reproduction, host-pathogen defense, and the origin of domesticated phenotypes.

Transcriptome-based gene expression profiling of diploid radish (Raphanus sativus L.) and the corresponding autotetraploid

Polyploidy is an important evolutionary factor in most land plant lineages which possess more than two complete sets of chromosomes. Radish (Raphanus sativus L.) is an economically annual/biennial root vegetable crop worldwide. However, the expression patterns of duplicated homologs involved in the autopolyploidization remains unclear. In present study, the autotetraploid radish plants (2n = 4x = 36) were produced with colchicine and exhibited an increase in the size of flowers, leaves, stomata and pollen grains. The differential gene expression (DGE) profiling was performed to investigate the differences in gene expression patterns between diploid and its corresponding autotetraploid by RNA-Sequencing (RNA-Seq). Totally, 483 up-regulated differentially expressed genes (DEGs) and 408 down-regulated DEGs were detected in diploid and autotetraploid radishes, which majorly involved in the pathways of hormones, photosynthesis and stress response. Moreover, the xyloglucan endotransglucosylase/hydrolase (XTH) and pectin methylesterases (PME) family members related to cell enlargement and cell wall construction were found to be enriched in GO enrichment analysis, of which XTH family members enriched in "apoplast" and "cell wall" terms, while PME family members enriched in "cell wall" term. Reverse-transcription quantitative PCR (RT-qPCR) analysis indicated that the expression profile of DEGs were consistent with results from the RNA-Seq analysis. The DEGs involved in cell wall construction and auxin metabolism were predicted to be associated with organs size increase of autotetraploid radishes in the present study. These results could provide valuable information for elucidating the molecular mechanism underlying polyploidization and facilitating further genetic improvements of important traits in radish breeding programs.

Integrating Networks, Phylogenomics, and Population Genomics for the Study of Polyploidy

Duplication events are regarded as sources of evolutionary novelty, but our understanding of general trends for the long-term trajectory of additional genomic material is still lacking. Organisms with a history of whole genome duplication (WGD) offer a unique opportunity to study potential trends in the context of gene retention and/or loss, gene and network dosage, and changes in gene expression. In this review, we discuss the prevalence of polyploidy across the tree of life, followed by an overview of studies investigating genome evolution and gene expression. We then provide an overview of methods in network biology, phylogenomics, and population genomics that are critical for advancing our understanding of evolution post-WGD, highlighting the need for models that can accommodate polyploids. Finally, we close with a brief note on the importance of random processes in the evolution of polyploids with respect to neutral versus selective forces, ancestral polymorphisms, and the formation of autopolyploids versus allopolyploids.

Comparative analysis of selection mode reveals different evolutionary rate and expression pattern in Arachis duranensis and Arachis ipaënsis duplicated genes

Our results reveal that K_s is a determining factor affecting selective pressure and different evolution and expression patterns are detected between PSGs and NSGs in wild Arachis duplicates. Selective pressure, including purifying (negative) and positive selection, can be detected in organisms. However, studies on comparative evolutionary rates, gene expression patterns and gene features between negatively selected genes (NSGs) and positively selected genes (PSGs) are lagging in paralogs of plants. Arachis duranensis and Arachis ipaënsis are ancestors of the cultivated peanut, an important oil and protein crop. Here, we carried out a series of systematic analyses, comparing NSG and PSG in paralogs, using genome sequences and transcriptome datasets in A. duranensis and A. ipaënsis. We found that synonymous substitution rate (K_s) is a determining factor affecting selective pressure in A. duranensis and A. ipaënsis duplicated genes. Lower expression level, lower gene expression breadth, higher codon bias and shorter polypeptide length were found in PSGs and not in NSGs. The correlation analyses showed that gene expression breadth was positively correlated with polypeptide length and GC content at the first codon site (GC1) in PSGs and NSGs, respectively. There was a negative correlation between expression level and polypeptide length in PSGs. In NSGs, the K_s was positively correlated with expression level, gene expression breadth, GC1, and GC content at the third codon site (GC3), but selective pressure was negatively correlated with expression level, gene expression breadth, polypeptide length, GC1, and GC3 content. The function of most duplicated gene pairs was divergent under drought and nematode stress. Taken together, our results show that different evolution and expression patterns occur between PSGs and NSGs in paralogs of two wild Arachis species.

WRKY transcription factors in legumes

BACKGROUND

WRKY transcription factors, so named because of the WRKYGQK heptapeptide at the N-terminal end, are widely distributed in plants and play an important role in physiological changes and response to biotic and abiotic stressors. Many previous studies have focused on the evolution of WRKY transcription factors in a given plant; however, little is known about WRKY evolution in legumes. The gene expression pattern of duplicated WRKY transcription factors remains unclear.

RESULTS

We first identified the WRKY proteins in 12 legumes. We found that the WRKYGQK heptapeptide tended to mutate into WRKYGKK. The Q site in WRKYGQK preferentially mutated, while W, K, and Y were conserved. The phylogenetic tree shows that the WRKY proteins in legumes have multiple origins, especially group IIc. For example, WRKY64 from Lupinus angustifolius (LaWRKY64) contains three WRKY domains, of which the first two clustered together in the N-terminal WRKY domain of the group I WRKY protein, and the third WRKY domain grouped in the C-terminal WRKY domain of the group I WRKY protein. Orthologous WRKY genes have a faster evolutionary rate and are subject to constrained selective pressure, unlike paralogous WRKY genes. Different gene features were observed between duplicated WRKY genes and singleton WRKY genes. Duplicated Glycine max WRKY genes with similar gene features have gene expression divergence.

CONCLUSIONS

We analyzed the WRKY number and type in 12 legumes, concluding that the WRKY proteins have multiple origins. A novel WRKY protein, LaWRKY64, was found in L. angustifolius. The first two WRKY domains of LaWRKY64 have the same origin. The orthologous and paralogous WRKY proteins have different evolutionary rates. Duplicated WRKY genes have gene expression divergence under normal growth conditions in G. max. These results provide insight into understanding WRKY evolution and expression.

The Grayling Genome Reveals Selection on Gene Expression Regulation after Whole-Genome Duplication

Whole-genome duplication (WGD) has been a major evolutionary driver of increased genomic complexity in vertebrates. One such event occurred in the salmonid family ∼80 Ma (Ss4R) giving rise to a plethora of structural and regulatory duplicate-driven divergence, making salmonids an exemplary system to investigate the evolutionary consequences of WGD. Here, we present a draft genome assembly of European grayling (Thymallus thymallus) and use this in a comparative framework to study evolution of gene regulation following WGD. Among the Ss4R duplicates identified in European grayling and Atlantic salmon (Salmo salar), one-third reflect nonneutral tissue expression evolution, with strong purifying selection, maintained over ∼50 Myr. Of these, the majority reflect conserved tissue regulation under strong selective constraints related to brain and neural-related functions, as well as higher-order protein–protein interactions. A small subset of the duplicates have evolved tissue regulatory expression divergence in a common ancestor, which have been subsequently conserved in both lineages, suggestive of adaptive divergence following WGD. These candidates for adaptive tissue expression divergence have elevated rates of protein coding- and promoter-sequence evolution and are enriched for immune- and lipid metabolism ontology terms. Lastly, lineage-specific duplicate divergence points toward underlying differences in adaptive pressures on expression regulation in the nonanadromous grayling versus the anadromous Atlantic salmon. Our findings enhance our understanding of the role of WGD in genome evolution and highlight cases of regulatory divergence of Ss4R duplicates, possibly related to a niche shift in early salmonid evolution.

Genome-Wide Identification and Characterization of Aquaporins and Their Role in the Flower Opening Processes in Carnation (Dianthus caryophyllus)

Aquaporins (AQPs) are associated with the transport of water and other small solutes across biological membranes. Genome-wide identification and characterization will pave the way for further insights into the AQPs’ roles in the commercial carnation (Dianthus caryophyllus). This study focuses on the analysis of AQPs in carnation (DcaAQPs) involved in flower opening processes. Thirty DcaAQPs were identified and grouped to five subfamilies: nine PIPs, 11 TIPs, six NIPs, three SIPs, and one XIP. Subsequently, gene structure, protein motifs, and co-expression network of DcaAQPs were analyzed and substrate specificity of DcaAQPs was predicted. qRT-PCR, RNA-seq, and semi-qRTRCR were used for DcaAQP genes expression analysis. The analysis results indicated that DcaAQPs were relatively conserved in gene structure and protein motifs, that DcaAQPs had significant differences in substrate specificity among different subfamilies, and that DcaAQP genes’ expressions were significantly different in roots, stems, leaves and flowers. Five DcaAQP genes (DcaPIP1;3, DcaPIP2;2, DcaPIP2;5, DcaTIP1;4, and DcaTIP2;2) might play important roles in flower opening process. However, the roles they play are different in flower organs, namely, sepals, petals, stamens, and pistils. Overall, this study provides a theoretical basis for further functional analysis of DcaAQPs.

Loss of VPS9b enhances vps9a-2 phenotypes

Plant innate immunity enables plants to defend themselves against infectious pathogens. While membrane trafficking and release of exosomes are considered vital for correct execution of innate immunity, the mechanisms behind remain elusive. Recently, we have shown that VPS9a, the general guanine-nucleotide exchange factor activating Rab5 GTPases, is required for both pre- and post-invasive immunity against powdery mildew fungi in Arabidopsis thaliana. Yet, the Arabidopsis genome contains a close homologue of VPS9a, which potentially plays specific roles in innate immunity. Here we show that this gene, VPS9b, while weakly expressed, contributes to regulating development and disease resistance, which is predominantly regulated by VPS9a. Based on these observations, we suggest that VPS9b has no specialized functionality, but rather is becoming a non-expressed pseudogene.