Modeling gene and genome duplications in eukaryotes

Inverse and Proportional Trans Modulation of Gene Expression in Human Aneuploidies

Genomic imbalance in aneuploidy is often detrimental to organisms. To gain insight into the molecular basis of aneuploidies in humans, we analyzed transcriptome data from several autosomal and sex chromosome aneuploidies. The results showed that in human aneuploid cells, genes located on unvaried chromosomes are inversely or proportionally trans-modulated, while a subset of genes on the varied chromosomes are compensated. Less genome-wide modulation is found for sex chromosome aneuploidy compared with autosomal aneuploidy due to X inactivation and the retention of dosage sensitive regulators on both sex chromosomes to limit the effective dosage change. We also found that lncRNA and mRNA can have different responses to aneuploidy. Furthermore, we analyzed the relationship between dosage-sensitive transcription factors and their targets, which illustrated the modulations and indicates genomic imbalance is related to stoichiometric changes in components of gene regulatory complexes.In summary, this study demonstrates the existence of trans-acting effects and compensation mechanisms in human aneuploidies and contributes to our understanding of gene expression regulation in unbalanced genomes and disease states.

wgd v2: a suite of tools to uncover and date ancient polyploidy and whole-genome duplication

MOTIVATION

Major improvements in sequencing technologies and genome sequence assembly have led to a huge increase in the number of available genome sequences. In turn, these genome sequences form an invaluable source for evolutionary, ecological, and comparative studies. One kind of analysis that has become routine is the search for traces of ancient polyploidy, particularly for plant genomes, where whole-genome duplication (WGD) is rampant.

RESULTS

Here, we present a major update of a previously developed tool wgd, namely wgd v2, to look for remnants of ancient polyploidy, or WGD. We implemented novel and improved previously developed tools to (a) construct KS age distributions for the whole-paranome (collection of all duplicated genes in a genome), (b) unravel intragenomic and intergenomic collinearity resulting from WGDs, (c) fit mixture models to age distributions of gene duplicates, (d) correct substitution rate variation for phylogenetic placement of WGDs, and (e) date ancient WGDs via phylogenetic dating of WGD-retained gene duplicates. The applicability and feasibility of wgd v2 for the identification and the relative and absolute dating of ancient WGDs is demonstrated using different plant genomes.

AVAILABILITY AND IMPLEMENTATION

wgd v2 is open source and available at https://github.com/heche-psb/wgd.

Genome organization and botanical diversity

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

Halophytes and heavy metals: A multi-omics approach to understand the role of gene and genome duplication in the abiotic stress tolerance of Cakile maritima

PREMISE

The origin of diversity is a fundamental biological question. Gene duplications are one mechanism that provides raw material for the emergence of novel traits, but evolutionary outcomes depend on which genes are retained and how they become functionalized. Yet, following different duplication types (polyploidy and tandem duplication), the events driving gene retention and functionalization remain poorly understood. Here we used Cakile maritima, a species that is tolerant to salt and heavy metals and shares an ancient whole-genome triplication with closely related salt-sensitive mustard crops (Brassica), as a model to explore the evolution of abiotic stress tolerance following polyploidy.

METHODS

Using a combination of ionomics, free amino acid profiling, and comparative genomics, we characterize aspects of salt stress response in C. maritima and identify retained duplicate genes that have likely enabled adaptation to salt and mild levels of cadmium.

RESULTS

Cakile maritima is tolerant to both cadmium and salt treatments through uptake of cadmium in the roots. Proline constitutes greater than 30% of the free amino acid pool in C. maritima and likely contributes to abiotic stress tolerance. We find duplicated gene families are enriched in metabolic and transport processes and identify key transport genes that may be involved in C. maritima abiotic stress tolerance.

CONCLUSIONS

These findings identify pathways and genes that could be used to enhance plant resilience and provide a putative understanding of the roles of duplication types and retention on the evolution of abiotic stress response.

Origin, evolution, and diversification of inositol 1,4,5-trisphosphate 3-kinases in plants and animals

BACKGROUND

In Eukaryotes, inositol polyphosphates (InsPs) represent a large family of secondary messengers and play crucial roes in various cellular processes. InsPs are synthesized through a series of pohophorylation reactions catalyzed by various InsP kinases in a sequential manner. Inositol 1,4,5-trisphosphate 3-kinase (IP3 3-kinase/IP3K), one member of InsP kinase, plays important regulation roles in InsPs metabolism by specifically phosphorylating inositol 1,4,5-trisphosphate (IP3) to inositol 1,3,4,5-tetrakisphosphate (IP4) in animal cells. IP3Ks were widespread in fungi, plants and animals. However, its evolutionary history and patterns have not been examined systematically.

RESULTS

A total of 104 and 31 IP3K orthologues were identified across 57 plant genomes and 13 animal genomes, respectively. Phylogenetic analyses indicate that IP3K originated in the common ancestor before the divergence of fungi, plants and animals. In most plants and animals, IP3K maintained low-copy numbers suggesting functional conservation during plant and animal evolution. In Brassicaceae and vertebrate, IP3K underwent one and two duplication events, respectively, resulting in multiple gene copies. Whole-genome duplication (WGD) was the main mechanism for IP3K duplications, and the IP3K duplicates have experienced functional divergence. Finally, a hypothetical evolutionary model for the IP3K proteins is proposed based on phylogenetic theory.

CONCLUSION

Our study reveals the evolutionary history of IP3K proteins and guides the future functions of animal, plant, and fungal IP3K proteins.

Ancient Duplication and Lineage-Specific Transposition Determine Evolutionary Trajectory of ERF Subfamily across Angiosperms

AP2/ERF transcription factor family plays an important role in plant development and stress responses. Previous studies have shed light on the evolutionary trajectory of the AP2 and DREB subfamilies. However, knowledge about the evolutionary history of the ERF subfamily in angiosperms still remains limited. In this study, we performed a comprehensive analysis of the ERF subfamily from 107 representative angiosperm species by combining phylogenomic and synteny network approaches. We observed that the expansion of the ERF subfamily was driven not only by whole-genome duplication (WGD) but also by tandem duplication (TD) and transposition duplication events. We also found multiple transposition events in Poaceae, Brassicaceae, Poales, Brassicales, and Commelinids. These events may have had notable impacts on copy number variation and subsequent functional divergence of the ERF subfamily. Moreover, we observed a number of ancient tandem duplications occurred in the ERF subfamily across angiosperms, e.g., in Subgroup IX, IXb originated from ancient tandem duplication events within IXa. These findings together provide novel insights into the evolution of this important transcription factor family.

Genome-wide analysis of PvMADS in common bean and functional characterization of PvMADS31 in Arabidopsis thaliana as a player in abiotic stress responses

Changing climatic conditions with rising temperatures and altered precipitation patterns pose significant challenges to agricultural productivity, particularly for common bean crops. Transcription factors (TFs) are crucial regulators that can mitigate the impact of biotic and abiotic stresses on crop production. The MADS-box TFs family has been implicated in various plant physiological processes, including stress-responsive mechanisms. However, their role in common bean and their response to stressful conditions remain poorly understood. Here, we identified 35 MADS-box gene family members in common bean, with conserved MADS-box domains and other functional domains. Gene duplication events were observed, suggesting the significance of duplication in the evolutionary development of gene families. The analysis of promoter regions revealed diverse elements, including stress-responsive elements, indicating their potential involvement in stress responses. Notably, PvMADS31, a member of the PvMADS-box gene family, demonstrated rapid upregulation under various abiotic stress conditions, including NaCl, polyethylene glycol, drought, and abscisic acid (ABA) treatments. Transgenic plants overexpressing PvMADS31 displayed enhanced lateral root development, root elongation, and seed germination under stress conditions. Furthermore, PvMADS31 overexpression in Arabidopsis resulted in improved drought tolerance, likely attributed to the enhanced scavenging of ROS and increased proline accumulation. These findings suggest that PvMADS31 might play a crucial role in modulating seed germination, root development, and stress responses, potentially through its involvement in auxin and ABA signaling pathways. Overall, this study provides valuable insights into the potential roles of PvMADS-box genes in abiotic stress responses in common bean, offering prospects for crop improvement strategies to enhance resilience under changing environmental conditions.

Enzymes Involved in Antioxidant and Detoxification Processes Present Changes in the Expression Levels of Their Coding Genes under the Stress Caused by the Presence of Antimony in Tomato

Currently, there is an increasing presence of heavy metals and metalloids in soils and water due to anthropogenic activities. However, the biggest problem caused by this increase is the difficulty in recycling these elements and their high permanence in soils. There are plants with great capacity to assimilate these elements or make them less accessible to other organisms. We analyzed the behavior of Solanum lycopersicum L., a crop with great agronomic interest, under the stress caused by antimony (Sb). We evaluated the antioxidant response throughout different exposure times to the metalloid. Our results showed that the enzymes involved in the AsA-GSH cycle show changes in their expression level under the stress caused by Sb but could not find a relationship between the NITROSOGLUTATHIONE REDUCTASE (GSNOR) expression data and nitric oxide (NO) content in tomato roots exposed to Sb. We hypothesize that a better understanding of how these enzymes work could be key to develop more tolerant varieties to this kind of abiotic stress and could explain a greater or lesser phytoremediation capacity. Moreover, we deepened our knowledge about Glutathione S-transferase (GST) and Glutathione Reductase (GR) due to their involvement in the elimination of the xenobiotic component.

Models for the retention of duplicate genes and their biological underpinnings

Gene content in genomes changes through several different processes, with gene duplication being an important contributor to such changes. Gene duplication occurs over a range of scales from individual genes to whole genomes, and the dynamics of this process can be context dependent. Still, there are rules by which genes are retained or lost from genomes after duplication, and probabilistic modeling has enabled characterization of these rules, including their context-dependence. Here, we describe the biology and corresponding mathematical models that are used to understand duplicate gene retention and its contribution to the set of biochemical functions encoded in a genome.

PbRbohH/J mediates ROS generation to regulate the growth of pollen tube in pear

Respiratory burst oxidase homolog (Rboh) family genes play crucial functions in development and growth. However, comprehensive and systematic investigation of Rboh family members in Rosaceae and their specific functions during pear pollen development are still limited. In the study, 63 Rboh genes were identified from eight Rosaceae genomes (Malus domestica, Pyrus bretschneideri, Pyrus communis, Prunus persica, Rubus occidentalis, Fragaria vesca, Prunus mume and Prunus avium) and divided into seven main subfamilies (I-VII) according to phylogenetic and structural features. Different modes of gene duplication led to the expansion of Rboh family, with purifying selection playing a vital role in the evolution of Rboh genes. In addition, RNA sequencing and qRT-PCR results indicated that PbRbohH and PbRbohJ were specifically high-expressed in pear pollen. Subsequently, subcellular localization revealed that PbRbohH/J distributed at the plasma membrane. Furthermore, by pharmacological analysis and antisense oligodeoxynucleotide assay, PbRbohH/J were demonstrated to mediate the formation of reactive oxygen species (ROS) to manage pollen tube growth. In conclusion, our results provide useful insights into the functions, expression patterns, evolutionary history of the Rboh genes in pear and other Rosaceae species.

Chromosome-level assembly and analysis of Camelina neglecta: a novel diploid model for Camelina biotechnology research

Camelina neglecta is a new diploid Brassicaceae species, which has great research value because of its close relationship with the hexaploid oilseed crop Camelina sativa. Here, we report a chromosome-level assembly of C. neglecta with a total length of 210 Mb. By adopting PacBio sequencing and Hi-C technology, the C. neglecta genome was assembled into 6 chromosomes with scaffold N50 of 29.62 Mb. C. neglecta has undergone the whole-genome triplication (γ) shared among eudicots and two whole-genome duplications (α and β) shared by crucifers, but it has not undergone a specific whole-genome duplication event. By synteny analysis between C. neglecta and C. sativa, we successfully used the method of calculating Ks to distinguish the three subgenomes of C. sativa and determined that C. neglecta was closest to the first subgenome (SG1) of C. sativa. Further, transcriptomic analysis revealed the key genes associated with seed oil biosynthesis and its transcriptional regulation, including SAD, FAD2, FAD3, FAE1, ABI3, WRI1 and FUS3 displaying high expression levels in C. neglecta seeds. The high representability of C. neglecta as a model species for Camelina-based biotechnology research has been demonstrated for the first time. In particular, floral Agrobacterium tumefaciens infiltration-based transformation of C. neglecta, leading to overexpression of CvLPAT2, CpDGAT1 and CvFatB1 transgenes, was demonstrated for medium-chain fatty acid accumulation in C. neglecta seed oil. This study provides an important genomic resource and establishes C. neglecta as a new model for oilseed biotechnology research.

Comprehensive Identification of the β-Amylase (BAM) Gene Family in Response to Cold Stress in White Clover

White clover (Trifolium repens L.) is an allopolyploid plant and an excellent perennial legume forage. However, white clover is subjected to various stresses during its growth, with cold stress being one of the major limiting factors affecting its growth and development. Beta-amylase (BAM) is an important starch-hydrolyzing enzyme that plays a significant role in starch degradation and responses to environmental stress. In this study, 21 members of the BAM gene family were identified in the white clover genome. A phylogenetic analysis using BAMs from Arabidopsis divided TrBAMs into four groups based on sequence similarity. Through analysis of conserved motifs, gene duplication, synteny analysis, and cis-acting elements, a deeper understanding of the structure and evolution of TrBAMs in white clover was gained. Additionally, a gene regulatory network (GRN) containing TrBAMs was constructed; gene ontology (GO) annotation analysis revealed close interactions between TrBAMs and AMY (α-amylase) and DPE (4-alpha-glucanotransferase). To determine the function of TrBAMs under various tissues and stresses, RNA-seq datasets were analyzed, showing that most TrBAMs were significantly upregulated in response to biotic and abiotic stresses and the highest expression in leaves. These results were validated through qRT-PCR experiments, indicating their involvement in multiple gene regulatory pathways responding to cold stress. This study provides new insights into the structure, evolution, and function of the white clover BAM gene family, laying the foundation for further exploration of the functional mechanisms through which TrBAMs respond to cold stress.

The genome of the toxic invasive species Heracleum sosnowskyi carries an increased number of genes despite absence of recent whole-genome duplications

Heracleum sosnowskyi, belonging to a group of giant hogweeds, is a plant with large effects on ecosystems and human health. It is an invasive species that contributes to the deterioration of grassland ecosystems. The ability of H. sosnowskyi to produce linear furanocoumarins (FCs), photosensitizing compounds, makes it very dangerous. At the same time, linear FCs are compounds with high pharmaceutical value used in skin disease therapies. Despite this high importance, it has not been the focus of genetic and genomic studies. Here, we report a chromosome-scale assembly of Sosnowsky's hogweed genome. Genomic analysis revealed an unusually high number of genes (55106) in the hogweed genome, in contrast to the 25-35 thousand found in most plants. However, we did not find any traces of recent whole-genome duplications not shared with its confamiliar, Daucus carota (carrot), which has approximately thirty thousand genes. The analysis of the genomic proximity of duplicated genes indicates on tandem duplications as a main reason for this increase. We performed a genome-wide search of the genes of the FC biosynthesis pathway and surveyed their expression in aboveground plant parts. Using a combination of expression data and phylogenetic analysis, we found candidate genes for psoralen synthase and experimentally showed the activity of one of them using a heterologous yeast expression system. These findings expand our knowledge on the evolution of gene space in plants and lay a foundation for further analysis of hogweed as an invasive plant and as a source of FCs.

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.

Expectations of duplicate gene retention under the gene duplicability hypothesis

BACKGROUND

Gene duplication is an important process in evolution. What causes some genes to be retained after duplication and others to be lost is a process not well understood. The most prevalent theory is the gene duplicability hypothesis, that something about the function and number of interacting partners (number of subunits of protein complex, etc.), determines whether copies have more opportunity to be retained for long evolutionary periods. Some genes are also more susceptible to dosage balance effects following WGD events, making them more likely to be retained for longer periods of time. One would expect these processes that affect the retention of duplicate copies to affect the conditional probability ratio after consecutive whole genome duplication events. The probability that a gene will be retained after a second whole genome duplication event (WGD2), given that it was retained after the first whole genome duplication event (WGD1) versus the probability a gene will be retained after WGD2, given it was lost after WGD1 defines the probability ratio that is calculated.

RESULTS

Since duplicate gene retention is a time heterogeneous process, the time between the events (t1) and the time since the most recent event (t2) are relevant factors in calculating the expectation for observation in any genome. Here, we use a survival analysis framework to predict the probability ratio for genomes with different values of t1 and t2 under the gene duplicability hypothesis, that some genes are more susceptible to selectable functional shifts, some more susceptible to dosage compensation, and others only drifting. We also predict the probability ratio with different values of t1 and t2 under the mutational opportunity hypothesis, that probability of retention for certain genes changes in subsequent events depending upon how they were previously retained. These models are nested such that the mutational opportunity model encompasses the gene duplicability model with shifting duplicability over time. Here we present a formalization of the gene duplicability and mutational opportunity hypotheses to characterize evolutionary dynamics and explanatory power in a recently developed statistical framework.

CONCLUSIONS

This work presents expectations of the gene duplicability and mutational opportunity hypotheses over time under different sets of assumptions. This expectation will enable formal testing of processes leading to duplicate gene retention.

Identification and analysis of BAHD superfamily related to malonyl ginsenoside biosynthesis in Panax ginseng

Introduction

The BAHD (benzylalcohol O-acetyl transferase, anthocyanin O-hydroxycinnamoyl transferase, N-hydroxycinnamoyl anthranilate benzoyl transferase and deacetylvindoline 4-O-acetyltransferase), has various biological functions in plants, including catalyzing the biosynthesis of terpenes, phenolics and esters, participating in plant stress response, affecting cell stability, and regulating fruit quality.

Methods

Bioinformatics methods, real-time fluorescence quantitative PCR technology, and ultra-high-performance liquid chromatography combined with an Orbitrap mass spectrometer were used to explore the relationship between the BAHD gene family and malonyl ginsenosides in Panax ginseng.

Results

In this study, 103 BAHD genes were identified in P. ginseng, mainly distributed in three major clades. Most PgBAHDs contain cis-acting elements associated with abiotic stress response and plant hormone response. Among the 103 genes, 68 PgBAHDs are WGD (whole-genome duplication) genes. The significance of malonylation in biosynthesis has garnered considerable attention in the study of malonyltransferases. The phylogenetic tree results showed 34 PgBAHDs were clustered with genes that have malonyl characterization. Among them, seven PgBAHDs (PgBAHD4, 45, 65, 74, 90, 97, and 99) showed correlations > 0.9 with crucial enzyme genes involved in ginsenoside biosynthesis and > 0.8 with malonyl ginsenosides. These seven genes were considered potential candidates involved in the biosynthesis of malonyl ginsenosides.

Discussion

These results help elucidate the structure, evolution, and functions of the P. ginseng BAHD gene family, and establish the foundation for further research on the mechanism of BAHD genes in ginsenoside biosynthesis.

Comparative transcriptomics reveals divergence in pathogen response gene families amongst 20 forest tree species

Forest trees provide critical ecosystem services for humanity that are under threat due to ongoing global change. Measuring and characterizing genetic diversity are key to understanding adaptive potential and developing strategies to mitigate negative consequences arising from climate change. In the area of forest genetic diversity, genetic divergence caused by large-scale changes at the chromosomal level has been largely understudied. In this study, we used the RNA-seq data of 20 co-occurring forest trees species from genera including Acer, Alnus, Amelanchier, Betula, Cornus, Corylus, Dirca, Fraxinus, Ostrya, Populus, Prunus, Quercus, Ribes, Tilia, and Ulmus sampled from Upper Peninsula of Michigan. These data were used to infer the origin and maintenance of gene family variation, species divergence time, as well as gene family expansion and contraction. We identified a signal of common whole genome duplication events shared by core eudicots. We also found rapid evolution, namely fast expansion or fast contraction of gene families, in plant-pathogen interaction genes amongst the studied diploid species. Finally, the results lay the foundation for further research on the genetic diversity and adaptive capacity of forest trees, which will inform forest management and conservation policies.

Genome-Wide Identification and Characterization of Gibberellic Acid-Stimulated Arabidopsis Gene Family in Pineapple (Ananas comosus)

The gibberellic acid-stimulated Arabidopsis (GASA) gene family plays a crucial role in growth, development, and stress response, and it is specific to plants. This gene family has been extensively studied in various plant species, and its functional role in pineapple has yet to be characterized. In this study, 15 AcGASA genes were identified in pineapple through a genome-wide scan and categorized into three major branches based on a phylogenetic tree. All AcGASA proteins share a common structural domain with 12 cysteine residues, but they exhibit slight variations in their physicochemical properties and motif composition. Predictions regarding subcellular localization suggest that AcGASA proteins are present in the cell membrane, Golgi apparatus, nucleus, and cell wall. An analysis of gene synteny indicated that both tandem and segmental repeats have a significant impact on the expansion of the AcGASA gene family. Our findings demonstrate the differing regulatory effects of these hormones (GA, NAA, IAA, MeJA, and ABA) on the AcGASA genes. We analyzed the expression profiles of GASA genes in different pineapple tissue parts, and the results indicated that AcGASA genes exhibit diverse expression patterns during the development of different plant tissues, particularly in the regulation of floral organ development. This study provides a comprehensive understanding of GASA family genes in pineapple. It serves as a valuable reference for future studies on the functional characterization of GASA genes in other perennial herbaceous plants.

Gene Duplication, Shifting Selection, and Dosage Balance of Silicon Transporter Proteins in Marine and Freshwater Diatoms

Numerous factors shape the evolution of protein-coding genes, including shifts in the strength or type of selection following gene duplications or changes in the environment. Diatoms and other silicifying organisms use a family of silicon transporters (SITs) to import dissolved silicon from the environment. Freshwaters contain higher silicon levels than oceans, and marine diatoms have more efficient uptake kinetics and less silicon in their cell walls, making them better competitors for a scarce resource. We compiled SITs from 37 diatom genomes to characterize shifts in selection following gene duplications and marine-freshwater transitions. A deep gene duplication, which coincided with a whole-genome duplication, gave rise to two gene lineages. One of them (SIT1-2) is present in multiple copies in most species and is known to actively import silicon. These SITs have evolved under strong purifying selection that was relaxed in freshwater taxa. Episodic diversifying selection was detected but not associated with gene duplications or habitat shifts. In contrast, genes in the second SIT lineage (SIT3) were present in just half the species, the result of multiple losses. Despite conservation of SIT3 in some lineages for the past 90-100 million years, repeated losses, relaxed selection, and low expression highlighted the dispensability of SIT3, consistent with a model of deterioration and eventual loss due to relaxed selection on SIT3 expression. The extensive but relatively balanced history of duplications and losses, together with paralog-specific expression patterns, suggest diatoms continuously balance gene dosage and expression dynamics to optimize silicon transport across major environmental gradients.

Tackling redundancy: genetic mechanisms underlying paralog compensation in plants

Gene duplication is a powerful source of biological innovation giving rise to paralogous genes that undergo diverse fates. Redundancy between paralogous genes is an intriguing outcome of duplicate gene evolution, and its maintenance over evolutionary time has long been considered a paradox. Redundancy can also be dubbed 'a geneticist's nightmare': It hinders the predictability of genome editing outcomes and limits our ability to link genotypes to phenotypes. Genetic studies in yeast and plants have suggested that the ability of ancient redundant duplicates to compensate for dosage perturbations resulting from a loss of function depends on the reprogramming of gene expression, a phenomenon known as active compensation. Starting from considerations on the stoichiometric constraints that drive the evolutionary stability of redundancy, this review aims to provide insights into the mechanisms of active compensation between duplicates that could be targeted for breaking paralog dependencies - the next frontier in plant functional studies.

Identification and characterization of the cupin_1 domain-containing proteins in ma bamboo (Dendrocalamus latiflorus) and their potential role in rhizome sprouting

Cupin_1 domain-containing protein (CDP) family, which is a member of the cupin superfamily with the most diverse functions in plants, has been found to be involved in hormone pathways that are closely related to rhizome sprouting (RS), a vital form of asexual reproduction in plants. Ma bamboo is a typical clumping bamboo, which mainly reproduces by RS. In this study, we identified and characterized 53 Dendrocalamus latiflorus CDP genes and divided them into seven subfamilies. Comparing the genetic structures among subfamilies showed a relatively conserved gene structure within each subfamily, and the number of cupin_1 domains affected the conservation among D. latiflorus CDP genes. Gene collinearity results showed that segmental duplication and tandem duplication both contributed to the expansion of D. latiflorus CDP genes, and lineage-specific gene duplication was an important factor influencing the evolution of CDP genes. Expression patterns showed that CDP genes generally had higher expression levels in germinating underground buds, indicating that they might play important roles in promoting shoot sprouting. Transcription factor binding site prediction and co-expression network analysis indicated that D. latiflorus CDPs were regulated by a large number of transcription factors, and collectively participated in rhizome buds and shoot development. This study significantly provided new insights into the evolutionary patterns and molecular functions of CDP genes, and laid a foundation for further studying the regulatory mechanisms of plant rhizome sprouting.

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

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

Ohnologs and SSD Paralogs Differ in Genomic and Expression Features Related to Dosage Constraints

Gene duplication is recognized as a critical process in genome evolution; however, many questions about this process remain unanswered. Although gene duplicability has been observed to differ by duplication mechanism and evolutionary rate, there is so far no broad characterization of its determinants. Many features correlate with this difference in duplicability; however, our ability to exploit these observations to advance our understanding of the role of duplication in evolution is hampered by limitations within existing work. In particular, the existence of methodological differences across studies impedes meaningful comparison. Here, we use consistent definitions of duplicability in the human lineage to explore these associations, allow resolution of the impact of confounding factors, and define the overall relevance of individual features. Using a classifier approach and controlling for the confounding effect of duplicate longevity, we find a subset of gene features important in differentiating genes duplicable by small-scale duplication from those duplicable by whole-genome duplication, revealing critical roles for gene dosage and expression costs in duplicability. We further delve into patterns of functional enrichment and find a lack of constraint on duplicate retention in any context for genes duplicable by small-scale duplication.

Evolutionary Consequences of Functional and Regulatory Divergence of HD-Zip I Transcription Factors as a Source of Diversity in Protein Interaction Networks in Plants

The HD superfamily has been studied in detail for several decades. The plant-specific HD-Zip I subfamily attracts the most attention because of its involvement in plant development and stress responses. In this review, we provide a comprehensive insight into the evolutionary events responsible for the functional redundancy and diversification of the HD-Zip I genes in regulating various biological processes. We summarized the evolutionary history of the HD-Zip family, highlighting the important role of WGDs in its expansion and divergence of retained duplicates in the genome. To determine the relationship between the evolutionary origin and functional conservation of HD-Zip I in different species, we performed a phylogenetic analysis, compared their expression profiles in different tissues and under stress and traced the role of orthologs and paralogs in regulating developmental processes. We found that HD-Zip I from different species have similar gene structures with a highly conserved HD and Zip, bind to the same DNA sequences and are involved in similar biological processes. However, they exhibit a functional diversity, which is manifested in altered expression patterns. Some of them are involved in the regulation of species-specific leaf morphology and phenotypes. Here, we discuss the role of changes in functional domains involved in DNA binding and protein interaction of HD-Zip I and in cis-regulated regions of its target genes in promoting adaptive innovations through the formation of de novo regulatory systems. Understanding the role of the HD-Zip I subfamily in organism-environment interactions remains a challenge for evolutionary developmental biology (evo-devo).

Genome-wide identification and expression analysis of the trihelix transcription factor family in sesame (Sesamum indicum L.) under abiotic stress

BACKGROUND

The plant trihelix gene family is among the earliest discovered transcription factor families, and it is vital in modulating light, plant growth, and stress responses.

METHODS

The identification and characterization of trihelix family members in the sesame genome were analyzed by bioinformatics methods, and the expression patterns of sesame trihelix genes were assessed by quantitative real-time PCR.

RESULTS

There were 34 trihelix genes discovered in the genome of sesame, which were irregularly distributed among 10 linkage groups. Also, the genome contained 5 duplicate gene pairs. The 34 trihelix genes were divided into six sub-families through a phylogenetic study. A tissue-specific expression revealed that SiTH genes exhibited spatial expression patterns distinct from other trihelix genes in the same subfamily. The cis-element showed that the SiTHs gene promoter contained various elements associated with responses to hormones and multiple abiotic stresses. Additionally, the expression patterns of 8 SiTH genes in leaves under abiotic stresses demonstrated that all selected genes were significantly upregulated or downregulated at least once in the stress period. Furthermore, the SiTH4 gene was significantly induced in response to drought and salt stress, showing that SiTH genes may be engaged in the stress response mechanisms of sesame.

CONCLUSION

These findings establish a foundation for further investigation of the trihelix gene-mediated response to abiotic stress in sesame.

Genomic and Transcriptomic Insights into the Evolution and Divergence of MIKC-Type MADS-Box Genes in Carica papaya

MIKC-type MADS-box genes, also known as type II genes, play a crucial role in regulating the formation of floral organs and reproductive development in plants. However, the genome-wide identification and characterization of type II genes as well as a transcriptomic survey of their potential roles in Carica papaya remain unresolved. Here, we identified and characterized 24 type II genes in the C. papaya genome, and investigated their evolutional scenario and potential roles with a widespread expression profile. The type II genes were divided into thirteen subclades, and gene loss events likely occurred in papaya, as evidenced by the contracted member size of most subclades. Gene duplication mainly contributed to MIKC-type gene formation in papaya, and the duplicated gene pairs displayed prevalent expression divergence, implying the evolutionary significance of gene duplication in shaping the diversity of type II genes in papaya. A large-scale transcriptome analysis of 152 samples indicated that different subclasses of these genes showed distinct expression patterns in various tissues, biotic stress response, and abiotic stress response, reflecting their divergent functions. The hub-network of male and female flowers and qRT-PCR suggested that TT16-3 and AGL8 participated in male flower development and seed germination. Overall, this study provides valuable insights into the evolution and functions of MIKC-type genes in C. papaya.

Expression analysis of alpha keratins and corneous beta-protein genes during embryonic development of Gekko japonicus

Epidermal appendages of birds and reptiles, including claws, feathers, scales, and setae, are primarily composed of alpha keratins (KRTs) and corneous beta-proteins (CBPs). A comprehensive and systematic knowledge of KRTs and CBPs in Schlegel's Japanese gecko (Gekko japonicus) is still lacking. In this study, 22 candidate Gecko japonicus keratin (GjKRT) family genes (12 type I genes, 10 type II genes) were identified in the G. japonicus genome. The majority of GjKRT genes across various subgroups had undergone a prolonged and highly conservative evolutionary process. Through a combination of morphological observation, RNA-seq analysis, and qRT-PCR assay, it was possible to discern the dynamic alterations in the expression of GjKRTs and Gecko japonicus corneous beta-proteins genes (GjCBPs). These findings strongly indicate that GjKRTs gradually accumulate to constitute an α-layer, which is subsequently succeeded by the formation of the corneous beta layer containing GjCBPs at late stages (40-42) of embryonic development. The epidermal appendages in G. japonicus may result from the joint accumulation of KRTs and CBPs, with stages 40-42 being critical for their development. These findings provide novel insights into KRTs and CBPs of G. japonicus and offer a foundation for investigating the functions of GjKRT and GjCBP gene families. Furthermore, this knowledge contributes to unraveling the molecular mechanisms underlying the formation of epidermal appendages in G. japonicus.

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

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

Identification and evolution analysis of YUCCA genes of Medicago sativa and Medicago truncatula and their expression profiles under abiotic stress

The YUCCAs (YUC) are functionally identified flavin-containing monooxidases (FMOs) in plants that act as an important rate-limiting enzyme functioning in the auxin synthesis IPA (indole-3-pyruvic acid) pathway. In this study, 12 MsYUCs and 15 MtYUCs containing characteristic conserved motifs were identified in M. sativa (Medicago sativa L.) and M. truncatula (Medicago truncatula Gaertn.), respectively. Phylogenetic analysis revealed that YUC proteins underwent an evolutionary divergence. Both tandem and segmental duplication events were presented in MsYUC and MtYUC genes. Comparative syntenic maps of M. sativa with M. truncatula, Arabidopsis (Arabidopsis thaliana), or rice (Oryza sativa L.) were constructed to illustrate the evolution relationship of the YUC gene family. A large number of cis-acting elements related to stress response and hormone regulation were revealed in the promoter sequences of MsYUCs. Expression analysis showed that MsYUCs had a tissue-specific, genotype-differential expression and a differential abiotic stress response pattern based on transcriptome data analysis of M. sativa online. In addition, RT-qPCR confirmed that salt stress significantly induced the expression of MsYUC1/MsYUC10 but significantly inhibited MsYUC2/MsYUC3 expression and the expression of MsYUC10/MsYUC11/MsYUC12 was significantly induced by cold treatment. These results could provide valuable information for functional analysis of YUC genes via gene engineering of the auxin synthetic IPA pathway in Medicago.

Analysis of gene duplication within the Arabidopsis NUCLEAR FACTOR Y, subunit B (NF-YB) protein family reveals domains under both purifying and diversifying selection

Gene duplication is an evolutionary mechanism that provides new genetic material. Since gene duplication is a major driver for molecular evolution, examining the fate of duplicated genes is an area of active research. The fate of duplicated genes can include loss, subfunctionalization, and neofunctionalization. In this manuscript, we chose to experimentally study the fate of duplicated genes using the Arabidopsis NUCLEAR FACTOR Y (NF-Y) transcription factor family. NF-Y transcription factors are heterotrimeric complexes, composed of NF-YA, NF-YB, and NF-YC. NF-YA subunits are responsible for nucleotide-specific binding to a CCAAT cis-regulatory element. NF-YB and NF-YC subunits make less specific, but essential complex-stabilizing contacts with the DNA flanking the core CCAAT pentamer. While ubiquitous in eukaryotes, each NF-Y family has expanded by duplication in the plant lineage. For example, the model plant Arabidopsis contains 10 each of the NF-Y subunits. Here we examine the fate of duplicated NF-YB proteins in Arabidopsis, which are composed of central histone fold domains (HFD) and less conserved flanking regions (N- and C-termini). Specifically, the principal question we wished to address in this manuscript was to what extent can the 10 Arabidopsis NF-YB paralogs functionally substitute the genes NF-YB2 and NF-YB3 in the promotion of photoperiodic flowering? Our results demonstrate that the conserved histone fold domains (HFD) may be under pressure for purifying (negative) selection, while the non-conserved N- and C-termini may be under pressure for diversifying (positive) selection, which explained each paralog's ability to substitute. In conclusion, our data demonstrate that the N- and C-termini may have allowed the duplicated genes to undergo functional diversification, allowing the retention of the duplicated genes.

Whole-genome Duplications and the Long-term Evolution of Gene Regulatory Networks in Angiosperms

Angiosperms have a complex history of whole-genome duplications (WGDs), with varying numbers and ages of WGD events across clades. These WGDs have greatly affected the composition of plant genomes due to the biased retention of genes belonging to certain functional categories following their duplication. In particular, regulatory genes and genes encoding proteins that act in multiprotein complexes have been retained in excess following WGD. Here, we inferred protein-protein interaction (PPI) networks and gene regulatory networks (GRNs) for seven well-characterized angiosperm species and explored the impact of both WGD and small-scale duplications (SSDs) in network topology by analyzing changes in frequency of network motifs. We found that PPI networks are enriched in WGD-derived genes associated with dosage-sensitive intricate systems, and strong selection pressures constrain the divergence of WGD-derived genes at the sequence and PPI levels. WGD-derived genes in network motifs are mostly associated with dosage-sensitive processes, such as regulation of transcription and cell cycle, translation, photosynthesis, and carbon metabolism, whereas SSD-derived genes in motifs are associated with response to biotic and abiotic stress. Recent polyploids have higher motif frequencies than ancient polyploids, whereas WGD-derived network motifs tend to be disrupted on the longer term. Our findings demonstrate that both WGD and SSD have contributed to the evolution of angiosperm GRNs, but in different ways, with WGD events likely having a more significant impact on the short-term evolution of polyploids.

Whole-genome identification and expression profiling of growth-regulating factor (GRF) and GRF-interacting factor (GIF) gene families in Panax ginseng

BACKGROUND

Panax ginseng is a perennial herb and one of the most widely used traditional medicines in China. During its long growth period, it is affected by various environmental factors. Past studies have shown that growth-regulating factors (GRFs) and GRF-interacting factors (GIFs) are involved in regulating plant growth and development, responding to environmental stress, and responding to the induction of exogenous hormones. However, GRF and GIF transcription factors in ginseng have not been reported.

RESULTS

In this study, 20 GRF gene members of ginseng were systematically identified and found to be distributed on 13 chromosomes. The ginseng GIF gene family has only ten members, which are distributed on ten chromosomes. Phylogenetic analysis divided these PgGRFs into six clades and PgGIFs into two clades. In total, 18 of the 20 PgGRFs and eight of the ten PgGIFs are segmental duplications. Most PgGRF and PgGIF gene promoters contain some hormone- and stress- related cis-regulatory elements. Based on the available public RNA-Seq data, the expression patterns of PgGRF and PgGIF genes were analysed from 14 different tissues. The responses of the PgGRF gene to different hormones (6-BA, ABA, GA3, IAA) and abiotic stresses (cold, heat, drought, and salt) were studied. The expression of the PgGRF gene was significantly upregulated under GA3 induction and three weeks of heat treatment. The expression level of the PgGIF gene changed only slightly after one week of heat treatment.

CONCLUSIONS

The results of this study may be helpful for further study of the function of PgGRF and PgGIF genes and lay a foundation for further study of their role in the growth and development of Panax ginseng.

Genome-Wide Identification of PEBP Gene Family in Solanum lycopersicum

The PEBP gene family is crucial for the growth and development of plants, the transition between vegetative and reproductive growth, the response to light, the production of florigen, and the reaction to several abiotic stressors. The PEBP gene family has been found in numerous species, but the SLPEBP gene family has not yet received a thorough bioinformatics investigation, and the members of this gene family are currently unknown. In this study, bioinformatics was used to identify 12 members of the SLPEBP gene family in tomato and localize them on the chromosomes. The physicochemical characteristics of the proteins encoded by members of the SLPEBP gene family were also examined, along with their intraspecific collinearity, gene structure, conserved motifs, and cis-acting elements. In parallel, a phylogenetic tree was built and the collinear relationships of the PEBP gene family among tomato, potato, pepper, and Arabidopsis were examined. The expression of 12 genes in different tissues and organs of tomato was analyzed using transcriptomic data. It was also hypothesized that SLPEBP3, SLPEBP5, SLPEBP6, SLPEBP8, SLPEBP9, and SLPEBP10 might be related to tomato flowering and that SLPEBP2, SLPEBP3, SLPEBP7, and SLPEBP11 might be related to ovary development based on the tissue-specific expression analysis of SLPEBP gene family members at five different stages during flower bud formation to fruit set. This article's goal is to offer suggestions and research directions for further study of tomato PEBP gene family members.

Functional and evolutionary study of MLO gene family in the regulation of Sclerotinia stem rot resistance in Brassica napus L

BACKGROUND

Oilseed rape (Brassica napus L.) is known as one of the most important oilseed crops cultivated around the world. However, its production continuously faces a huge challenge of Sclerotinia stem rot (SSR), a destructive disease caused by the fungus Sclerotinia sclerotiorum, resulting in huge yield loss annually. The SSR resistance in B. napus is quantitative and controlled by a set of minor genes. Identification of these genes and pyramiding them into a variety are a major strategy for SSR resistance breeding in B. napus.

RESULTS

Here, we performed a genome-wide association study (GWAS) using a natural population of B. napus consisting of 222 accessions to identify BnaA08g25340D (BnMLO2_2) as a candidate gene that regulates the SSR resistance. BnMLO2_2 was a member of seven homolog genes of Arabidopsis Mildew Locus O 2 (MLO2) and the significantly SNPs were mainly distributed in the promoter of BnMLO2_2, suggesting a role of BnMLO2_2 expression level in the regulation of SSR resistance. We expressed BnMLO2_2 in Arabidopsis and the transgenic plants displayed an enhanced SSR resistance. Transcriptome profiling of different tissues of B. napus revealed that BnMLO2_2 had the most expression level in leaf and silique tissues among all the 7 BnMLO2 members and also expressed higher in the SSR resistant accession than in the susceptible accession. In Arabidopsis, mlo2 plants displayed reduced resistance to SSR, whereas overexpression of MLO2 conferred plants an enhanced SSR resistance. Moreover, a higher expression level of MLO2 showed a stronger SSR resistance in the transgenic plants. The regulation of MLO2 in SSR resistance may be associated with the cell death. Collinearity and phylogenetic analysis revealed a large expansion of MLO family in Brassica crops.

CONCLUSION

Our study revealed an important role of BnMLO2 in the regulation of SSR resistance and provided a new gene candidate for future improvement of SSR resistance in B. napus and also new insights into understanding of MLO family evolution in Brassica crops.

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

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

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Genome-Wide Identification and Expression Analysis of the SWEET Gene Family in Annual Alfalfa (Medicago polymorpha)

SWEET (Sugars will eventually be exported transporter) proteins are a group of sugar transporters that are involved in sugar efflux, phloem loading, reproductive development, plant senescence, and stress responses. In this study, 23 SWEET transporter members were identified in the Medicago polymorpha genome, heterogeneously distributed on seven chromosomes. These MpSWEET genes were divided into four subfamilies, which showed similar gene structure and motif composition within the same subfamily. Seventeen MpSWEET genes encode seven transmembrane helices (TMHs), and all MpSWEET proteins possess conserved membrane domains and putative serine phosphorylation sites. Four and three pairs of MpSWEET genes were predicted to be segmentally and tandemly duplicated, respectively, which may have contributed to their evolution of M. polymorpha. The results of microarray and RNA-Seq data showed that some MpSWEET genes were specifically expressed in disparate developmental stages (including seedling stage, early flowering stage, and late flowering stage) or tissues such as flower and large pod. Based on protein network interaction and expression patterns of MpSWEET genes, six MpSWEET genes were selected for further quantitative real-time PCR validation in different stress treatments. qRT-PCR results showed that MpSWEET05, MpSWEET07, MpSWEET12, MpSWEET15, and MpSWEET21 were significantly upregulated for at least two of the three abiotic stress treatments. These findings provide new insights into the complex transcriptional regulation of MpSWEET genes, which facilitates future research to elucidate the function of MpSWEET genes in M. polymorpha and other legume crops.

Genome-Wide Characterization and Expression of the Hsf Gene Family in Salvia miltiorrhiza (Danshen) and the Potential Thermotolerance of SmHsf1 and SmHsf7 in Yeast

Salvia miltiorrhiza Bunge (Danshen) is a traditional Chinese herb with significant medicinal value. The yield and quality of Danshen are greatly affected by climatic conditions, in particular high temperatures. Heat shock factors (Hsfs) play important regulatory roles in plant response to heat and other environmental stresses. However, little is currently known about the role played by the Hsf gene family in S. miltiorrhiza. Here, we identified 35 SmHsf genes and classified them into three major groups: SmHsfA (n = 22), SmHsfB (n = 11), and SmHsfC (n = 2) using phylogenetic analysis. The gene structure and protein motifs were relatively conserved within subgroups but diverged among the different groups. The expansion of the SmHsf gene family was mainly driven by whole-genome/segmental and dispersed gene duplications. The expression profile of SmHsfs in four distinct organs revealed its members (23/35) are predominantly expressed in the root. The expression of a large number of SmHsfs was regulated by drought, ultraviolet, heat and exogenous hormones. Notably, the SmHsf1 and SmHsf7 genes in SmHsfB2 were the most responsive to heat and are conserved between dicots and monocots. Finally, heterologous expression analysis showed that SmHsf1 and SmHsf7 enhance thermotolerance in yeast. Our results provide a solid foundation for further functional investigation of SmHsfs in Danshen plants as a response to abiotic stresses.

Adaptive Selection in the Evolution of Aquaglyceroporins in Mammals

Aquaporins (AQPs) are integral membrane proteins responsible for water transport across cellular membranes in both prokaryotes and eukaryotes. A subfamily of AQPs, known as aquaglyceroporins (AQGPs), facilitate the transport of small solutes such as glycerol, water, and other solutes across cellular membranes. These proteins are involved in a variety of physiological processes, such as organogenesis, wound healing, and hydration. Although AQPs have been studied extensively in different species, their conservation patterns, phylogenetic relationships, and evolution in mammals remain unexplored. In the present study, 119 AQGP coding sequences from 31 mammalian species were analysed to identify conserved residues, gene organisation, and most importantly, the nature of AQGP gene selection. Repertoire analysis revealed the absence of AQP7, 9, and 10 genes in certain species of Primates, Rodentia, and Diprotodontia, although not all three genes were absent in a single species. Two Asparagine-Proline-Alanine (NPA) motifs located at the N- and C-terminal ends, aspartic acid (D) residues, and the ar/R region were conserved in AQP3, 9, and 10. Six exons encoding the functional MIP domain of AQGP genes were found to be conserved across mammalian species. Evolutionary analysis indicated signatures of positive selection in AQP7, 9, and 10 amongst different mammalian lineages. Furthermore, substitutions of certain amino acids located close to critical residues may alter AQGP functionality, which is crucial for substrate selectivity, pore formation, and transport efficiency required for the maintenance of homeostasis in different mammalian species.

Dosage balance acts as a time-dependent selective barrier to subfunctionalization

BACKGROUND

Gene duplication is an important process for genome expansion, sometimes allowing for new gene functions to develop. Duplicate genes can be retained through multiple processes, either for intermediate periods of time through processes such as dosage balance, or over extended periods of time through processes such as subfunctionalization and neofunctionalization.

RESULTS

Here, we built upon an existing subfunctionalization Markov model by incorporating dosage balance to describe the interplay between subfunctionalization and dosage balance to explore selective pressures on duplicate copies. Our model incorporates dosage balance using a biophysical framework that penalizes the fitness of genetic states with stoichiometrically imbalanced proteins. These imbalanced states cause increased concentrations of exposed hydrophobic surface areas, which cause deleterious mis-interactions. We draw comparison between our Subfunctionalization + Dosage-Balance Model (Sub + Dos) and the previous Subfunctionalization-Only (Sub-Only) Model. This comparison includes how the retention probabilities change over time, dependent upon the effective population size and the selective cost associated with spurious interaction of dosage-imbalanced partners. We show comparison between Sub-Only and Sub + Dos models for both whole-genome duplication and small-scale duplication events.

CONCLUSION

These comparisons show that following whole-genome duplication, dosage balance serves as a time-dependent selective barrier to the subfunctionalization process, by causing an overall delay but ultimately leading to a larger portion of the genome retained through subfunctionalization. This higher percentage of the genome that is ultimately retained is caused by the alternative competing process, nonfunctionalization, being selectively blocked to a greater extent. In small-scale duplication, the reverse pattern is seen, where dosage balance drives faster rates of subfunctionalization, but ultimately leads to a smaller portion of the genome retained as duplicates. This faster rate of subfunctionalization is because the dosage balance of interacting gene products is negatively affected immediately after duplication and the loss of a duplicate restores the stoichiometric balance. Our findings provide support that the subfunctionalization of genes that are susceptible to dosage balance effects, such as proteins involved in complexes, is not a purely neutral process. With stronger selection against stoichiometrically imbalanced gene partners, the rates of subfunctionalization and nonfunctionalization slow; however, this ultimately leads to a greater proportion of subfunctionalized gene pairs.

A chromosome-level genome assembly of radish (Raphanus sativus L.) reveals insights into genome adaptation and differential bolting regulation

High-quality radish (Raphanus sativus) genome represents a valuable resource for agronomical trait improvements and understanding genome evolution among Brassicaceae species. However, existing radish genome assembly remains fragmentary, which greatly hampered functional genomics research and genome-assisted breeding. Here, using a NAU-LB radish inbred line, we generated a reference genome of 476.32 Mb with a scaffold N50 of 56.88 Mb by incorporating Illumina, PacBio and BioNano optical mapping techniques. Utilizing Hi-C data, 448.12 Mb (94.08%) of the assembled sequences were anchored to nine radish chromosomes with 40 306 protein-coding genes annotated. In total, 249.14 Mb (52.31%) comprised the repetitive sequences, among which long terminal repeats (LTRs, 30.31%) were the most abundant class. Beyond confirming the whole-genome triplication (WGT) event in R. sativus lineage, we found several tandem arrayed genes were involved in stress response process, which may account for the distinctive phenotype of high disease resistance in R. sativus. By comparing against the existing Xin-li-mei radish genome, a total of 2 108 573 SNPs, 7740 large insertions, 7757 deletions and 84 inversions were identified. Interestingly, a 647-bp insertion in the promoter of RsVRN1 gene can be directly bound by the DOF transcription repressor RsCDF3, resulting into its low promoter activity and late-bolting phenotype of NAU-LB cultivar. Importantly, introgression of this 647-bp insertion allele, RsVRN1^In-536 , into early-bolting genotype could contribute to delayed bolting time, indicating that it is a potential genetic resource for radish late-bolting breeding. Together, this genome resource provides valuable information to facilitate comparative genomic analysis and accelerate genome-guided breeding and improvement in radish.

The Gynandropsis gynandra genome provides insights into whole-genome duplications and the evolution of C4 photosynthesis in Cleomaceae

Gynandropsis gynandra (Cleomaceae) is a cosmopolitan leafy vegetable and medicinal plant, which has also been used as a model to study C4 photosynthesis due to its evolutionary proximity to C3 Arabidopsis (Arabidopsis thaliana). Here, we present the genome sequence of G. gynandra, anchored onto 17 main pseudomolecules with a total length of 740 Mb, an N50 of 42 Mb and 30,933 well-supported gene models. The G. gynandra genome and previously released genomes of C3 relatives in the Cleomaceae and Brassicaceae make an excellent model for studying the role of genome evolution in the transition from C3 to C4 photosynthesis. Our analyses revealed that G. gynandra and its C3 relative Tarenaya hassleriana shared a whole-genome duplication event (Gg-α), then an addition of a third genome (Th-α, +1×) took place in T. hassleriana but not in G. gynandra. Analysis of syntenic copy number of C4 photosynthesis-related gene families indicates that G. gynandra generally retained more duplicated copies of these genes than C3T. hassleriana, and also that the G. gynandra C4 genes might have been under positive selection pressure. Both whole-genome and single-gene duplication were found to contribute to the expansion of the aforementioned gene families in G. gynandra. Collectively, this study enhances our understanding of the polyploidy history, gene duplication and retention, as well as their impact on the evolution of C4 photosynthesis in Cleomaceae.

DNA methylation signatures of duplicate gene evolution in angiosperms

Gene duplication is a source of evolutionary novelty. DNA methylation may play a role in the evolution of duplicate genes (paralogs) through its association with gene expression. While this relationship has been examined to varying extents in a few individual species, the generalizability of these results at either a broad phylogenetic scale with species of differing duplication histories or across a population remains unknown. We applied a comparative epigenomics approach to 43 angiosperm species across the phylogeny and a population of 928 Arabidopsis (Arabidopsis thaliana) accessions, examining the association of DNA methylation with paralog evolution. Genic DNA methylation was differentially associated with duplication type, the age of duplication, sequence evolution, and gene expression. Whole genome duplicates were typically enriched for CG-only gene-body methylated or unmethylated genes, while single-gene duplications were typically enriched for non-CG methylated or unmethylated genes. Non-CG methylation, in particular, was characteristic of more recent single-gene duplicates. Core angiosperm gene families differentiated into those which preferentially retain paralogs and 'duplication-resistant' families, which convergently reverted to singletons following duplication. Duplication-resistant families that still have paralogous copies were, uncharacteristically for core angiosperm genes, enriched for non-CG methylation. Non-CG methylated paralogs had higher rates of sequence evolution, higher frequency of presence-absence variation, and more limited expression. This suggests that silencing by non-CG methylation may be important to maintaining dosage following duplication and be a precursor to fractionation. Our results indicate that genic methylation marks differing evolutionary trajectories and fates between paralogous genes and have a role in maintaining dosage following duplication.

A2TEA: Identifying trait-specific evolutionary adaptations

Background: Plants differ in their ability to cope with external stresses (e.g., drought tolerance). Genome duplications are an important mechanism to enable plant adaptation. This leads to characteristic footprints in the genome, such as protein family expansion. We explore genetic diversity and uncover evolutionary adaptation to stresses by exploiting genome comparisons between stress tolerant and sensitive species and RNA-Seq data sets from stress experiments. Expanded gene families that are stress-responsive based on differential expression analysis could hint at species or clade-specific adaptation, making these gene families exciting candidates for follow-up tolerance studies and crop improvement. Software: Integration of such cross-species omics data is a challenging task, requiring various steps of transformation and filtering. Ultimately, visualization is crucial for quality control and interpretation. To address this, we developed A2TEA: Automated Assessment of Trait-specific Evolutionary Adaptations, a Snakemake workflow for detecting adaptation footprints in silico. It functions as a one-stop processing pipeline, integrating protein family, phylogeny, expression, and protein function analyses. The pipeline is accompanied by an R Shiny web application that allows exploring, highlighting, and exporting the results interactively. This allows the user to formulate hypotheses regarding the genomic adaptations of one or a subset of the investigated species to a given stress. Conclusions: While our research focus is on crops, the pipeline is entirely independent of the underlying species and can be used with any set of species. We demonstrate pipeline efficiency on real-world datasets and discuss the implementation and limits of our analysis workflow as well as planned extensions to its current state. The A2TEA workflow and web application are publicly available at: https://github.com/tgstoecker/A2TEA.Workflow and https://github.com/tgstoecker/A2TEA.WebApp, respectively.

Genomic imbalance modulates transposable element expression in maize

Genomic imbalance refers to the more severe phenotypic consequences of changing part of a chromosome compared with the whole genome set. Previous genome imbalance studies in maize have identified prevalent inverse modulation of genes on the unvaried chromosomes (trans) with both the addition or subtraction of chromosome arms. Transposable elements (TEs) comprise a substantial fraction of the genome, and their reaction to genomic imbalance is therefore of interest. Here, we analyzed TE expression using RNA-seq data of aneuploidy and ploidy series and found that most aneuploidies showed an inverse modulation of TEs, but reductions in monosomy and increases in disomy and trisomy were also common. By contrast, the ploidy series showed little TE modulation. The modulation of TEs and genes in the same experimental group were compared, and TEs showed greater modulation than genes, especially in disomy. Class I and II TEs were differentially modulated in most aneuploidies, and some superfamilies in each TE class also showed differential modulation. Finally, the significantly upregulated TEs in three disomies (TB-7Lb, TB9Lc, and TB-10L19) did not increase the proportion of adjacent gene expression when compared with non-differentially expressed TEs, indicating that modulations of TEs do not compound the effect on genes. These results suggest that the prevalent inverse TE modulation in aneuploidy results from stoichiometric upset of the regulatory machinery used by TEs, similar to the response of core genes to genomic imbalance.

Identification of Camellia oleifera WRKY transcription factor genes and functional characterization of CoWRKY78

Camellia oleifera Abel is a highly valued woody edible oil tree, which is endemic to China. It has great economic value because C. oleifera seed oil contains a high proportion of polyunsaturated fatty acids. C. oleifera anthracnose caused by Colletotrichum fructicola, poses a serious threat to C. oleifera growth and yield and causes the benefit of the C. oleifera industry to suffer directly. The WRKY transcription factor family members have been widely characterized as vital regulators in plant response to pathogen infection. Until now, the number, type and biological function of C. oleifera WRKY genes are remains unknown. Here, we identified 90 C. oleifera WRKY members, which were distributed across 15 chromosomes. C. oleifera WRKY gene expansion was mainly attributed to segmental duplication. We performed transcriptomic analyses to verify the expression patterns of CoWRKYs between anthracnose-resistant and -susceptible cultivars of C. oleifera. These results demonstrated that multiple candidate CoWRKYs can be induced by anthracnose and provide useful clues for their functional studies. CoWRKY78, an anthracnose-induced WRKY gene, was isolated from C. oleifera. It was significantly down-regulated in anthracnose-resistant cultivars. Overexpression of CoWRKY78 in tobacco markedly reduced resistance to anthracnose than WT plants, as evidenced by more cell death, higher malonaldehyde content and reactive oxygen species (ROS), but lower activities of superoxide dismutase (SOD), peroxidase (POD), as well as phenylalanine ammonia-lyase (PAL). Furthermore, the expression of multiple stress-related genes, which are associated with ROS-homeostasis (NtSOD and NtPOD), pathogen challenge (NtPAL), and pathogen defense (NtPR1, NtNPR1, and NtPDF1.2) were altered in the CoWRKY78-overexpressing plants. These findings increase our understanding of the CoWRKY genes and lay the foundation for the exploration of anthracnose resistance mechanisms and expedite the breeding of anthracnose-resistant C. oleifera cultivars.

Genome-wide analysis of R2R3-MYB genes in cultivated peanut (Arachis hypogaea L.): Gene duplications, functional conservation, and diversification

The cultivated Peanut (Arachis hypogaea L.), an important oilseed and edible legume, are widely grown worldwide. The R2R3-MYB transcription factor, one of the largest gene families in plants, is involved in various plant developmental processes and responds to multiple stresses. In this study we identified 196 typical R2R3-MYB genes in the genome of cultivated peanut. Comparative phylogenetic analysis with Arabidopsis divided them into 48 subgroups. The motif composition and gene structure independently supported the subgroup delineation. Collinearity analysis indicated polyploidization, tandem, and segmental duplication were the main driver of the R2R3-MYB gene amplification in peanut. Homologous gene pairs between the two subgroups showed tissue specific biased expression. In addition, a total of 90 R2R3-MYB genes showed significant differential expression levels in response to waterlogging stress. Furthermore, we identified an SNP located in the third exon region of AdMYB03-18 (AhMYB033) by association analysis, and the three haplotypes of the SNP were significantly correlated with total branch number (TBN), pod length (PL) and root-shoot ratio (RS ratio), respectively, revealing the potential function of AdMYB03-18 (AhMYB033) in improving peanut yield. Together, these studies provide evidence for functional diversity in the R2R3-MYB genes and will contribute to understanding the function of R2R3-MYB genes in peanut.

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.

Inference of Ancient Polyploidy from Genomic Data

Whole-genome sequence data have revealed that numerous eukaryotic organisms derive from distant polyploid ancestors, even when these same organisms are genetically and karyotypically diploid. Such ancient whole-genome duplications (WGDs) have been important for long-term genome evolution and are often speculatively associated with important evolutionary events such as key innovations, adaptive radiations, or survival after mass extinctions. Clearly, reliable methods for unveiling ancient WGDs are key toward furthering understanding of the long-term evolutionary significance of polyploidy. In this chapter, we describe a set of basic established comparative genomics approaches for the inference of ancient WGDs from genomic data based on empirical age distributions and collinearity analyses, explain the principles on which they are based, and illustrate a basic workflow using the software "wgd," geared toward a typical exploratory analysis of a newly obtained genome sequence.

MicroRNAs play regulatory roles in genomic balance

Classic genetics studies found that genomic imbalance caused by changing the dosage of part of the genome (aneuploidy) has more detrimental effects than altering the dosage of the whole genome (ploidy). Previous analysis revealed global modulation of gene expression triggered by aneuploidy across various species, including maize (Zea mays), Arabidopsis, yeast, mammals, etc. Plant microRNAs (miRNAs) are a class of 20- to 24-nt endogenous small noncoding RNAs that carry out post-transcriptional gene expression regulation. That miRNAs and their putative targets are preferentially retained as duplicates after whole-genome duplication, as are many transcription factors and signaling components, indicates miRNAs are likely to be dosage-sensitive and potentially involved in genomic balance networks. This review addresses the following questions regarding the role of miRNAs in genomic imbalance. (1) How do aneuploidy and polyploidy impact the expression of miRNAs? (2) Do miRNAs play a regulatory role in modulating the expression of their targets under genomic imbalance?