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

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

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

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

Duplicate Gene Expression and Possible Mechanisms of Paralog Retention During Bacterial Genome Expansion.

Gene duplication contributes to the evolution of expression and the origin of new genes, but the relative importance of different patterns of duplicate gene expression and mechanisms of retention remains debated and particularly poorly understood in bacteria. Here, we investigated gene expression patterns for two lab strains of the cyanobacterium Acaryochloris marina with expanding genomes that contain about 10-fold more gene duplicates compared with most bacteria. Strikingly, we observed a generally stoichiometric pattern of greater combined duplicate transcript dosage with increased gene copy number, in contrast to the prevalence of expression reduction reported for many eukaryotes. We conclude that increased transcript dosage is likely an important mechanism of initial duplicate retention in these bacteria and may persist over long periods of evolutionary time. However, we also observed that paralog expression can diverge rapidly, including possible functional partitioning, for which different copies were respectively more highly expressed in at least one condition. Divergence may be promoted by the physical separation of most Acaryochloris duplicates on different genetic elements. In addition, expression pattern for ancestrally shared duplicates could differ between strains, emphasizing that duplicate expression fate need not be deterministic. We further observed evidence for context-dependent transcript dosage, where the aggregate expression of duplicates was either greater or lower than their single-copy homolog depending on physiological state. Finally, we illustrate how these different expression patterns of duplicated genes impact Acaryochloris biology for the innovation of a novel light-harvesting apparatus and for the regulation of recA paralogs in response to environmental change.

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.

Genome-Wide Investigation of the PLD Gene Family in Tomato: Identification, Analysis, and Expression.

Phospholipase Ds (PLDs) are important phospholipid hydrolases in plants that play crucial roles in the regulation of plant growth, development, and stress tolerance. In this study, 14 PLD genes were identified in the tomato genome and were localized on eight chromosomes, and one tandem-duplicated gene pair was identified. According to a phylogenetic analysis, the genes were categorized into four subtypes: SlPLDα, ß, and δ belonged to the C2-PLD subfamily, while SlPLDζ belonged to the PXPH-PLD subfamily. The gene structure and protein physicochemical properties were highly conserved within the same subtype. The promoter of all the SlPLD genes contained hormone-, light-, and stress-responsive cis-acting regulatory elements, but no significant correlation between the number, distribution, and type of cis-acting elements was observed among the members of the same subtype. Transcriptome data showed that the expression of the SlPLD genes was different in multiple tissues. A quantitative RT-PCR analysis revealed that the SlPLD genes responded positively to cold, salt, drought, and abscisic acid treatments, particularly to salt stress. Different expression patterns were observed for different genes under the same stress, and for the same gene under different stresses. The results provide important insights into the functions of SlPLD genes and lay a foundation for further studies of the response of SlPLD genes to abiotic stresses.

Young duplicate genes show developmental stage- and cell type-specific expression and function in Caenorhabditis elegans.

Gene duplication produces the material that fuels evolutionary innovation. The "out-of-testis" hypothesis suggests that sperm competition creates selective pressure encouraging the emergence of new genes in male germline, but the somatic expression and function of the newly evolved genes are not well understood. We systematically mapped the expression of young duplicate genes throughout development in Caenorhabditis elegans using both whole-organism and single-cell transcriptomic data. Based on the expression dynamics across developmental stages, young duplicate genes fall into three clusters that are preferentially expressed in early embryos, mid-stage embryos, and late-stage larvae. Early embryonic genes are involved in protein degradation and develop essentiality comparable to the genomic average. In mid-to-late embryos and L4-stage larvae, young genes are enriched in intestine, epidermal cells, coelomocytes, and amphid chemosensory neurons. Their molecular functions and inducible expression indicate potential roles in innate immune response and chemosensory perceptions, which may contribute to adaptation outside of the sperm.

Characterization of the pyruvate kinase gene family in soybean and identification of a putative salt responsive gene GmPK21.

BACKGROUND: As a key regulatory enzyme in the glycolysis pathway, pyruvate kinase (PK) plays crucial roles in multiple physiological processes during plant growth and is also involved in the abiotic stress response. However, little information is known about PKs in soybean. RESULTS: In this study, we identified 27 PK family genes against the genome of soybean cultivar Zhonghuang13. They were classified into 2 subfamilies including PKc and PKp. 22 segmental duplicated gene pairs and 1 tandem duplicated gene pair were identified and all of them experienced a strong purifying selective pressure during evolution. Furthermore, the abiotic stresses (especially salt stress) and hormone responsive cis-elements were present in the promoters of GmPK genes, suggesting their potential roles in abiotic stress tolerance. By performing the qRT-PCR, 6 GmPK genes that continuously respond to both NaCl and ABA were identified. Subsequently, GmPK21, which represented the most significant change under NaCl treatment was chosen for further study. Its encoded protein GmPK21 was localized in the cytoplasm and plasma membrane. The transgenic Arabidopsis overexpressing GmPK21 exhibited weakened salinity tolerance. CONCLUSIONS: This study provides genomic information of soybean PK genes and a molecular basis for mining salt tolerance function of PKs in the future.

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.

Why old duplicated genes are not thrown away in (paleo)polyploids? Example from the petC gene in Brassica napus.

Duplication of genes at different time period, through recurrent and frequent polyploidization events, have played a major role in plant evolution, adaptation and diversification. Interestingly, some of the ancestral duplicated genes (referred as paleologs), have been maintained for millions of years, and there is still a poor knowledge of the reasons of their retention, especially when testing the phenotypic effect of individual copies by using functional genetic approaches. To fill this gap, we performed functional genetic (CRISPR-Cas9), physiological, transcriptomic and evolutionary studies to finely investigate this open question, taking the example of the petC gene (involved in cytochrome b6/f and thus impacting photosynthesis) that is present in four paleologous copies in the oilseed crop Brassica napus. RNA-Seq and selective pressure analyses suggested that all paleologous copies conserved the same function and that they were all highly transcribed. Thereafter, the Knock Out (K.O.) of one, several or all petC copies highlighted that all paleologous copies have to be K.O. to suppress the gene function. In addition, we could determine that phenotypic effects in single and double mutants could only be deciphered in high light conditions. Interestingly, we did not detect any significant differences between single mutants K.O. for either the A03 or A09 copy (despite being differentially transcribed), or even between mutants for a single or two petC copies. Altogether, this work revealed that petC paleologs have retained their ancestral function and that the retention of these copies is explained by their compensatory role, especially in optimal environmental conditions.

Expression profiles in knock-down transgenic plants of high and low diversified duplicate genes in Arabidopsis thaliana.

Duplicated genes show various degrees of functional diversification in plants. We previously identified 1,052 pairs of high diversified duplicates (HDDs) and 600 pairs of low diversified duplicates (LDDs) in Arabidopsis thaliana. Single knock-down of HDDs induced abnormal phenotypic changes because the other gene copy could not compensate for the knock-down effect, while single knock-down of LDDs did not induce abnormal phenotypic changes because of functional compensation by the copy gene. Here, focusing on one pair each of HDDs and LDDs, we performed transcriptome analyses in single-knock-down transgenic plants. The numbers of differentially expressed genes in single-knock-down transgenic plants were not different between HDDs and LDDs. Thus, functional compensation inferred by transcriptomics was similar between HDDs and LDDs. However, the trend of differentially expressed genes was similar in the pair of LDDs, while expression profiles were dissimilar in the pair of HDDs. This result indicates that a pair of LDDs tends to share similar functions but a pair of HDDs tends to have undergone functional divergence. Taking these findings together, as the reason for no phenotypic changes in single knock-down of LDDs but phenotypic changes in double knock-down of LDDs, we concluded that phenotypic changes of LDDs were induced by decreasing gene dosage.

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

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

Interlocus Gene Conversion, Natural Selection, and Paralog Homogenization.

Following a duplication, the resulting paralogs tend to diverge. While mutation and natural selection can accelerate this process, they can also slow it. Here, we quantify the paralog homogenization that is caused by point mutations and interlocus gene conversion (IGC). Among 164 duplicated teleost genes, the median percentage of postduplication codon substitutions that arise from IGC rather than point mutation is estimated to be between 7% and 8%. By differentiating between the nonsynonymous codon substitutions that homogenize the protein sequences of paralogs and the nonhomogenizing nonsynonymous substitutions, we estimate the homogenizing nonsynonymous rates to be higher for 163 of the 164 teleost data sets as well as for all 14 data sets of duplicated yeast ribosomal protein-coding genes that we consider. For all 14 yeast data sets, the estimated homogenizing nonsynonymous rates exceed the synonymous rates.

Genome-Wide Identification and Expression Patterns of Cucumber Invertases and Their Inhibitor Genes.

Invertases and their inhibitors play important roles in sucrose metabolism, growth and development, signal transduction, and biotic and abiotic stress tolerance in many plant species. However, in cucumber, both the gene members and functions of invertase and its inhibitor families remain largely unclear. In this study, in comparison with the orthologues of Citrullus lanatus (watermelon), Cucumis melo (melon), and Arabidopsis thaliana (Arabidopsis), 12 invertase genes and 12 invertase inhibitor genes were identified from the genome of Cucumis sativus (cucumber). Among them, the 12 invertase genes were classified as 4 cell wall invertases, 6 cytoplasmic invertases, and 2 vacuolar invertases. Most invertase genes were conserved in cucumber, melon, and watermelon, with several duplicate genes in melon and watermelon. Transcriptome analysis distinguished these genes into various expression patterns, which included genes CsaV3_2G025540 and CsaV3_2G007220, which were significantly expressed in different tissues, organs, and development stages, and genes CsaV3_7G034730 and CsaV3_5G005910, which might be involved in biotic and abiotic stress. Six genes were further validated in cucumber based on quantitative real-time PCR (qRT-PCR), and three of them showed consistent expression patterns as revealed in the transcriptome. These results provide important information for further studies on the physiological functions of cucumber invertases (CSINVs) and their inhibitors (CSINHs).

Novel genes and alleles of the BTB/POZ protein family in Oryza rufipogon.

The BTB/POZ family of proteins is widespread in plants and animals, playing important roles in development, growth, metabolism, and environmental responses. Although members of the expanded BTB/POZ gene family (OsBTB) have been identified in cultivated rice (Oryza sativa), their conservation, novelty, and potential applications for allele mining in O. rufipogon, the direct progenitor of O. sativa ssp. japonica and potential wide-introgression donor, are yet to be explored. This study describes an analysis of 110 BTB/POZ encoding gene loci (OrBTB) across the genome of O. rufipogon as outcomes of tandem duplication events. Phylogenetic grouping of duplicated OrBTB genes was supported by the analysis of gene sequences and protein domain architecture, shedding some light on their evolution and functional divergence. The O. rufipogon genome encodes nine novel BTB/POZ genes with orthologs in its distant cousins in the family Poaceae (Sorghum bicolor, Brachypodium distachyon), but such orthologs appeared to have been lost in its domesticated descendant, O. sativa ssp. japonica. Comparative sequence analysis and structure comparisons of novel OrBTB genes revealed that diverged upstream regulatory sequences and regulon restructuring are the key features of the evolution of this large gene family. Novel genes from the wild progenitor serve as a reservoir of potential new alleles that can bring novel functions to cultivars when introgressed by wide hybridization. This study establishes a foundation for hypothesis-driven functional genomic studies and their applications for widening the genetic base of rice cultivars through the introgression of novel genes or alleles from the exotic gene pool.

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.

Spatial Colinear but Broken Temporal Expression of Duplicated ParaHox Genes in Asexually Reproducing Annelids, Nais communis and Pristina longiseta.

ParaHox genes are key developmental regulators involved in the patterning of the digestive tract along the anteroposterior axis and the development of the nervous system. Most studies have focused on the function of these genes in embryogenesis, while their expression patterns in postembryonic development often remain unknown. In this study, we identified for the first time all ParaHox orthologs in two naidid oligochaetes, N. communis and P. longiseta, and described their expression patterns during normal growth and fission in these animals. We showed that Gsx and Cdx are presented by two paralogs, while Xlox is a single copy gene in both species. Using whole-mount in situ hybridization, we also found that orthologs, except for the Xlox gene, have similar activity patterns with minor differences in details, while the expression patterns of paralogs can differ significantly. However, all these genes are involved in axial patterning and/or in tissue remodeling during growth and asexual reproduction in naidids. Moreover, during paratomic fission, these genes are expressed with spatial colinearity but temporal colinearity is broken. The results of this study may be evidence of the functional diversification of duplicated genes and suggest involvement of the ParaHox genes in whole-body patterning during growth and asexual reproduction in annelids.

Genome-Wide Analysis of the MADS-box Gene Family and Expression Analysis during Anther Development in Salvia miltiorrhiza.

MADS-box genes constitute a large family of transcription factors that play important roles in plant growth and development. However, our understanding of MADS-box genes involved in anther development and male sterility in Salvia miltiorrhiza is still limited. In this study, 63 MADS-box genes were identified from the genome of the male sterility ecotype Sichuan S. miltiorrhiza (S. miltiorrhiza_SC) unevenly distributed among eight chromosomes. Phylogenetic analysis classified them into two types and 17 subfamilies. They contained 1 to 12 exons and 10 conserved motifs. Evolution analysis showed that segmental duplication was the main force for the expansion of the SmMADS gene family, and duplication gene pairs were under purifying selection. Cis-acting elements analysis demonstrated that the promoter of SmMADS genes contain numerous elements associated with plant growth and development, plant hormones, and stress response. RNA-seq showed that the expression levels of B-class and C-class SmMADS genes were highly expressed during anther development, with SmMADS11 likely playing an important role in regulating anther development and male fertility in S. miltiorrhiza_SC. Overall, this study provides a comprehensive analysis of the MADS-box gene family in S. miltiorrhiza, shedding light on their potential role in anther development and male sterility.

Genomic and Transcriptional Profiles of Kelch-like (klhl) Gene Family in Polyploid Carassius Complex.

Genome duplication supplies raw genetic materials and has been thought to be essential for evolutionary innovation and ecological adaptation. Here, we select Kelch-like (klhl) genes to study the evolution of the duplicated genes in the polyploid Carassius complex, including amphidiploid C. auratus and amphitriploid C. gibelio. Phylogenetic, chromosomal location and read coverage analyses indicate that most of Carassius klhl genes exhibit a 2:1 relationship with zebrafish orthologs and confirm two rounds of polyploidy, an allotetraploidy followed by an autotriploidy, occurred during Carassius evolution. The lineage-specific expansion and biased retention/loss of klhl genes are also found in Carassius. Transcriptome analyses across eight adult tissues and seven embryogenesis stages reveal varied expression dominance and divergence between the two species. The expression of klhls in response to Carassius herpesvirus 2 infection shows different expression changes corresponding to distinct herpesvirus resistances in three C. gibelio gynogenetic clones. Finally, we find that most C. gibelio klhl genes possess three alleles except eight genes that have lost one or two alleles due to genome rearrangement. The allele expression bias is prosperous for Cgklhl genes and varies during embryogenesis owning to the sequential expression manner of the alleles. The current study provides global insights into the genomic and transcriptional evolution of duplicated genes in a given superfamily resulting from multiple rounds of polyploidization.

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.