Search

Duplicated genes are thought to follow one of three evolutionary trajectories that resolve their redundancy: neofunctionalization, subfunctionalization, or pseudogenization. Differences in expression patterns have been documented for many duplicated gene pairs and interpreted as evidence of subfunctionalization and a loss of redundancy. However, little is known about the functional impact of such differences and about their molecular basis. Here, we investigate the genetic and molecular basis for the partial loss of redundancy between the two BLADE-ON-PETIOLE genes BOP1 and BOP2 in red shepherd's purse (Capsella rubella) compared to Arabidopsis (Arabidopsis thaliana). While both genes remain almost fully redundant in A. thaliana, BOP1 in C. rubella can no longer ensure wild-type floral organ numbers and suppress bract formation, due to an altered expression pattern in the region of the cryptic bract primordium. We use two complementary approaches, transgenic rescue of A. thaliana atbop1 atbop2 double mutants and deletions in the endogenous AtBOP1 promoter, to demonstrate that several BOP1 promoter regions containing conserved noncoding sequences interact in a nonadditive manner to control BOP1 expression in the bract primordium and that changes in these interactions underlie the evolutionary divergence between C. rubella and A. thaliana BOP1 expression and activity. Similarly, altered interactions between cis-regulatory regions underlie the divergence in functional promoter architecture related to the control of floral organ abscission by BOP1. These findings highlight the complexity of promoter architecture in plants and suggest that changes in the interactions between cis-regulatory elements are key drivers for evolutionary divergence in gene expression and the loss of redundancy., Competing Interests: Conflict of interest statement. None declared., (© The Author(s) 2024. Published by Oxford University Press on behalf of American Society of Plant Biologists.)

Polyploidy, a prevalent event in plant evolution, drives phenotypic diversification and speciation. While transcriptional changes and regulation in polyploids have been extensively studied, the translational level impact remains largely unexplored. To address this gap, we conducted a comparative transcriptomic and translatomic analysis of cotton leaves from allopolyploid species G. hirsutum (AD 1 ) and G. barbadense (AD 2 ) relative to their model A-genome and D-genome diploid progenitors. Our data revealed that while allopolyploidization significantly affects the transcriptional landscape, its impact on translation was relatively modest, evidenced by a narrower expression range and fewer expression changes in ribosome-protected fragments than in mRNA levels. Allopolyploid-specific changes commonly identified in both AD 1 and AD 2 were observed in 7393 genes at either transcriptional or translational levels. Interestingly, the majority of translational changes exhibited concordant down-regulation in both ribosome-protected fragments and mRNA, particularly associated with terpenoid synthesis and metabolism (352 genes). Regarding translational efficiency (TE), at least one-fifth of cotton genes exhibit translational level regulation, with a general trend of more down-regulation (13.9-15.1%) than up-regulation (7.3-11.2%) of TE. The magnitude of translational regulation was slightly reduced in allopolyploids compared with diploids, and allopolyploidy tends to have a more profound impact on genes and functional associations with ultra-low TE. Moreover, we demonstrated a reduced extent of homeolog expression biases during translation compared with transcription. Our study provides insights into the regulatory consequences of allopolyploidy post-transcription, contributing to a comprehensive understanding of regulatory mechanisms of duplicated gene expression evolution.

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., (© 2024. The Author(s).)

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.

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., Competing Interests: No competing interests were disclosed., (Copyright: © 2024 Assis R et al.)

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., (© 2023. The Author(s).)

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.

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., (© 2023. The Author(s), under exclusive licence to Springer Nature B.V.)

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.

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 epigenomic 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 a characteristic of more recent single-gene duplicates. Core angiosperm gene families were 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., Competing Interests: Conflict of interest statement. None declared., (© The Author(s) 2023. Published by Oxford University Press on behalf of American Society of Plant Biologists.)

Tandem duplication, one of the major types of duplication, provides the raw material for the evolution of divergent functions. In this study, we identified 1 pair of tandem duplicate genes (AT5G12950 and AT5G12960) in Arabidopsis (Arabidopsis thaliana) that originated within the last 16 million years after the split of Arabidopsis from the Capsella-Boechera ancestor. We systematically used bioinformatic tools to redefine their putative biochemical function as β-L-arabinofuranosidases that release L-Arabinose from the β-L-Araf-containing molecules in Arabidopsis. Comprehensive transcriptomic and proteomic analyses using various datasets showed divergent expression patterns among tissues between the 2 duplicate genes. We further collected phenotypic data from 2 types of measurements to indicate that AT5G12950 and AT5G12960 have different roles resulting in divergent phenotypic effects. Overall, AT5G12950 and AT5G12960 represent putative β-L-arabinofuranosidase encoding genes in Arabidopsis. After duplication, 1 duplicate copy developed diverged biological functions and contributed to a different phenotypic evolution in Arabidopsis., Competing Interests: Conflict of interest statement. All authors declared there are no conflicts of interest., (© American Society of Plant Biologists 2023. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

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.

Many tools have been developed to measure the degree of similarity between gene duplicates within and between species. Here, we present HSDecipher, a bioinformatics pipeline to assist users in the analysis and visualization of highly similar duplicate genes (HSDs). We describe the steps for analysis of HSDs statistics, expanding HSD gene sets, and visualizing the results of comparative genomic analyses. HSDecipher represents a useful tool for researchers exploring the evolution of duplicate genes in select eukaryotic species. For complete details on the use and execution of this protocol, please refer to Zhang et al. (2021) 1 and Zhang et al. (2022). 2 ., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Gene duplications have long been recognized as a contributor to the evolution of genes with new functions. Multiple copies of genes can result from tandem duplication, from transposition to new chromosomes, or from whole-genome duplication (polyploidy). The most common fate is that one member of the pair is deleted to return the gene to the singleton state. Other paths involve the reduced expression of both copies (hypofunctionalization) that are held in duplicate to maintain sufficient quantity of function. The two copies can split functions (subfunctionalization) or can diverge to generate a new function (neofunctionalization). Retention of duplicates resulting from doubling of the whole genome occurs for genes involved with multicomponent interactions such as transcription factors and signal transduction components. In contrast, these classes of genes are underrepresented in small segmental duplications. This complementary pattern suggests that the balance of interactors affects the fate of the duplicate pair. We discuss the different mechanisms that maintain duplicated genes, which may change over time and intersect., (� The Author(s) 2022. Published by Oxford University Press on behalf of American Society of Plant Biologists.)

Background: Gene conversion has an important effect on duplicate genes produced by polyploidization. Poplar (Populus trichocarpa) and willow (Salix brachista) are leading models and excellent green plants in the Salicaceae. Although much attention has been paid to the evolution of duplicated genes in poplar and willow, the role of conversion between duplicates generated from polyploidization remains poorly understood., Results: Here, through genomic synteny analyses, we identified duplicate genes generated by the Salicaceae common tetraploidization (SCT) in the poplar and willow genomes. We estimated that at least 0.58% and 0.25% of poplar and willow duplicates were affected by whole-gene conversion after the poplar-willow divergence, with more (5.73% and 2.66%) affected by partial-gene conversion. Moreover, we found that the converted duplicated genes were unevenly distributed on each chromosome in the two genomes, and the well-preserved homoeologous chromosome regions may facilitate the conversion of duplicates. Notably, we found that conversion maintained the similarity of duplicates, likely contributing to the conservation of certain sequences, but is essentially accelerated the rate of evolution and increased species divergence. In addition, we found that converted duplicates tended to have more similar expression patterns than nonconverted duplicates. We found that genes associated with multigene families were preferentially converted. We also found that the genes encoding conserved structural domains associated with specific traits exhibited a high frequency of conversion., Conclusions: Extensive conversion between duplicate genes generated from the SCT contributes to the diversification of the family Salicaceae and has had long-lasting effects on those genes with important biological functions., (© 2022. The Author(s).)

Bactrocera dorsalis is an invasive polyphagous pest causing considerable ecological and economic damage worldwide. We report a high-quality chromosome-level genome assembly and combine various transcriptome data to explore the molecular mechanisms of its rapid adaptation to new environments. The expansions of the DDE transposase superfamily and key gene families related to environmental adaptation and enrichment of the expanded and unique gene families in metabolism and defence response pathways explain its environmental adaptability. The relatively high but not significantly different expression of heat-shock proteins, regardless of the environmental conditions, suggests an intrinsic mechanism underlying its adaptation to high temperatures. The mitogen-activated protein kinase pathway plays a key role in adaptation to new environments. The prevalence of duplicated genes in its genome explains the diversity in the B. dorsalis complex. These findings provide insights into the genetic basis of the invasiveness and diversity of B. dorsalis, explaining its rapid adaptation and expansion., (© 2022. The Author(s).)

Gene duplication is a prevalent phenomenon across the tree of life. The processes that lead to the retention of duplicated genes are not well understood. Functional genomics approaches in model organisms, such as yeast, provide useful tools to test the mechanisms underlying retention with functional redundancy and divergence of duplicated genes, including fates associated with neofunctionalization, subfunctionalization, back-up compensation, and dosage amplification. Duplicated genes may also be retained as a consequence of structural and functional entanglement. Advances in human gene editing have enabled the interrogation of duplicated genes in the human genome, providing new tools to evaluate the relative contributions of each of these factors to duplicate gene retention and the evolution of genome structure., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

The 14-3-3 family genes are highly conserved regulatory factors in eukaryotes with involvement in multiple important cellular processes. However, detailed investigations of this family in fishes are very limited. Here, a comparative genomic and transcriptomic survey were performed to investigate the 14-3-3 family in fishes. We confirmed that the numbers of 14-3-3 genes ranged from 5 to 7 in non-teleost fishes, as well as additional 14-3-3 genes (9 to 11) in teleost fishes. In addition, some special teleost fishes possess 17 to 25 14-3-3s, which undergone the fourth whole-genome duplication (WGD). We also found that six pairs of fish 14-3-3 genes were clustered with mammalian ε, γ, ς, η, τand β isotypes, respectively, while σ was absent with a potential specificity within mammals, on the basis of their phylogenetic and synteny analyses. According to our results, we inferred that the diversity of 14-3-3 genes in fishes seems to be generated from a combination of WGD and gene loss. Comparative transcriptomic analysis revealed that there are differences in tissue distribution, and we speculated that 14-3-3 genes may contribute to terrestrial adaptations in mudskippers. In addition, protein sequence alignments of 14-3-3s supported their differential roles in fishes. In summary, our present comparative genomic and transcriptomic survey will benefit for further functional investigations of these fish genes., (Copyright © 2021. Published by Elsevier B.V.)

The tea plant (Camellia sinensis) is a thermophilic cash crop and contains a highly duplicated and repeat-rich genome. It is still unclear how DNA methylation regulates the evolution of duplicated genes and chilling stress in tea plants. We therefore generated a single-base-resolution DNA methylation map of tea plants under chilling stress. We found that, compared with other plants, the tea plant genome is highly methylated in all three sequence contexts, including CG, CHG and CHH (where H = A, T, or C), which is further proven to be correlated with its repeat content and genome size. We show that DNA methylation in the gene body negatively regulates the gene expression of tea plants, whereas non-CG methylation in the flanking region enables a positive regulation of gene expression. We demonstrate that transposable element-mediated methylation dynamics significantly drives the expression divergence of duplicated genes in tea plants. The DNA methylation and expression divergence of duplicated genes in the tea plant increases with evolutionary age and selective pressure. Moreover, we detect thousands of differentially methylated genes, some of which are functionally associated with chilling stress. We also experimentally reveal that DNA methyltransferase genes of tea plants are significantly downregulated, whereas demethylase genes are upregulated at the initial stage of chilling stress, which is in line with the significant loss of DNA methylation of three well-known cold-responsive genes at their promoter and gene body regions. Overall, our findings underscore the importance of DNA methylation regulation and offer new insights into duplicated gene evolution and chilling tolerance in tea plants., (© 2021 Society for Experimental Biology and John Wiley & Sons Ltd.)

Gene duplication is an important driver of genomic diversity that can promote adaptive evolution. However, like most mutations, a newly duplicated gene is often deleterious and removed from the genome by drift or natural selection. The early molecular changes that occur soon after duplication therefore may influence the long-term survival of gene duplicates, but relatively little empirical data exist on the events near the onset of duplication before mutations have time to accumulate. In this study, we contrast gene expression and DNA methylation levels of duplicate genes in the threespine stickleback, Gasterosteus aculeatus, including recently emerged duplications that segregate as copy number variations (CNVs). We find that younger duplicate genes have higher levels of promoter methylation than older genes, and that gene CNVs have higher promoter methylation than non-CNVs. These results suggest preferential duplication of highly methylated genes or rapid methylation changes soon after duplication. We also find a negative association between methylation and expression, providing a putative role for methylation in suppressing transcription that compensates for increases in gene copy numbers and promoting paralog retention. We propose that methylation contributes to the longevity of young duplicate genes, extending the window of opportunity for functional divergence via mutation., (© 2021 The Authors. Evolution © 2021 The Society for the Study of Evolution.)

Vector population control using insecticides is a key element of current strategies to prevent malaria transmission in Africa. The introduction of effective insecticides, such as the organophosphate pirimiphos-methyl, is essential to overcome the recurrent emergence of resistance driven by the highly diverse Anopheles genomes. Here, we use a population genomic approach to investigate the basis of pirimiphos-methyl resistance in the major malaria vectors Anopheles gambiae and A. coluzzii. A combination of copy number variation and a single non-synonymous substitution in the acetylcholinesterase gene, Ace1, provides the key resistance diagnostic in an A. coluzzii population from Côte d'Ivoire that we used for sequence-based association mapping, with replication in other West African populations. The Ace1 substitution and duplications occur on a unique resistance haplotype that evolved in A. gambiae and introgressed into A. coluzzii, and is now common in West Africa primarily due to selection imposed by other organophosphate or carbamate insecticides. Our findings highlight the predictive value of this complex resistance haplotype for phenotypic resistance and clarify its evolutionary history, providing tools to for molecular surveillance of the current and future effectiveness of pirimiphos-methyl based interventions., Competing Interests: The authors have declared that no competing interests exist.

Background: Alu elements are most abundant retrotransposons with > 1.2 million copies in the primate genome. AluYb8 subfamily was diverged from AluY lineage, and has accumulated eight diagnostic mutations and 7-bp duplication during primate evolution. A total of 1851 AluYb copies are present in the human genome, and most of them are human-specific. On the other hand, only a few AluYb8 copies were identified in the chimpanzee genome by previous studies on AluYb8. The significantly different number of species-specific AluYb8 elements between human and chimpanzee might result from the incompletion of chimpanzee reference genome sequences at the time of the previous study., Objective: AluYb8 elements could generate genomic structural variations in the chimpanzee genome. This study aimed to identify and characterize chimpanzee-specific AluYb elements using the most updated chimpanzee reference genome sequences (Jan. 2018, panTro6)., Methods: To identify chimpanzee-specific AluYb8, we carried out genomic comparison with non-chimpanzee primate genome using the UCSC table browser. In addition, chimpanzee-specific AluYb8 candidates were manually inspected and experimentally verified using PCR and Sanger sequencing., Results: Among a total of 231 chimpanzee-specific AluYb8 candidates, 11 of the candidates are chimpanzee-specific AluYb8, and 29 elements are shared between the chimpanzee and non-chimpanzee primate genomes. Through the sequence analysis of AluYb8 and other Alu subfamilies, we were able to observe various diagnostic mutations and variable length duplications in 7-bp duplication region of AluYb8 element. In addition, we further validated two of the chimpanzee-specific AluYb8 elements (CS8 and CS20) that were not previously discovered by display PCR and Sanger sequencing. Interestingly, we identified a AluYb8 insertion-mediated deletion (CS8 locus) in the chimpanzee genome., Conclusion: Our study found that AluYb8 elements are much more abundant in the human genome than chimpanzee genome, and that it could be due to the absence of hyperactive "master" AluYb8 elements in the chimpanzee genome.

Achieving complete and precise genome duplication requires that each genomic segment be replicated only once per cell division cycle. Protecting large eukaryotic genomes from re-replication requires an overlapping set of molecular mechanisms that prevent the first DNA replication step, the DNA loading of MCM helicase complexes to license replication origins, after S phase begins. Previous reports have defined many such origin licensing inhibition mechanisms, but the temporal relationships among them are not clear, particularly with respect to preventing re-replication in G2 and M phases. Using a combination of mutagenesis, biochemistry, and single cell analyses in human cells, we define a new mechanism that prevents re-replication through hyperphosphorylation of the essential MCM loading protein, Cdt1. We demonstrate that Cyclin A/CDK1 can hyperphosphorylate Cdt1 to inhibit MCM re-loading in G2 phase. The mechanism of inhibition is to block Cdt1 binding to MCM independently of other known Cdt1 inactivation mechanisms such as Cdt1 degradation during S phase or Geminin binding. Moreover, our findings suggest that Cdt1 dephosphorylation at the mitosis-to-G1 phase transition re-activates Cdt1. We propose that multiple distinct, non-redundant licensing inhibition mechanisms act in a series of sequential relays through each cell cycle phase to ensure precise genome duplication., Competing Interests: The authors have declared that no competing interests exist.

The Cpi-17 ( ppp1r14 ) gene family is an evolutionarily conserved, vertebrate specific group of protein phosphatase 1 (PP1) inhibitors. When phosphorylated, Cpi-17 is a potent inhibitor of myosin phosphatase (MP), a holoenzyme complex of the regulatory subunit Mypt1 and the catalytic subunit PP1. Myosin phosphatase dephosphorylates the regulatory myosin light chain (Mlc2) and promotes actomyosin relaxation, which in turn, regulates numerous cellular processes including smooth muscle contraction, cytokinesis, cell motility, and tumor cell invasion. We analyzed zebrafish homologs of the Cpi-17 family, to better understand the mechanisms of myosin phosphatase regulation. We found single homologs of both Kepi ( ppp1r14c ) and Gbpi ( ppp1r14d ) in silico, but we detected no expression of these genes during early embryonic development. Cpi-17 ( ppp1r14a ) and Phi-1 ( ppp1r14b ) each had two duplicate paralogs, ( ppp1r14aa and ppp1r14ab ) and ( ppp1r14ba and ppp1r14bb ), which were each expressed during early development. The spatial expression pattern of these genes has diverged, with ppp1r14aa and ppp1r14bb expressed primarily in smooth muscle and skeletal muscle, respectively, while ppp1r14ab and ppp1r14ba are primarily expressed in neural tissue. We observed that, in in vitro and heterologous cellular systems, the Cpi-17 paralogs both acted as potent myosin phosphatase inhibitors, and were indistinguishable from one another. In contrast, the two Phi-1 paralogs displayed weak myosin phosphatase inhibitory activity in vitro, and did not alter myosin phosphorylation in cells. Through deletion and chimeric analysis, we identified that the difference in specificity for myosin phosphatase between Cpi-17 and Phi-1 was encoded by the highly conserved PHIN (phosphatase holoenzyme inhibitory) domain, and not the more divergent N- and C- termini. We also showed that either Cpi-17 paralog can rescue the knockdown phenotype, but neither Phi-1 paralog could do so. Thus, we provide new evidence about the biochemical and developmental distinctions of the zebrafish Cpi-17 protein family.

Internal tandem duplication (ITD) mutations of FMS-like tyrosine kinase-3 (FLT3) are the most frequent genetic alterations in acute myeloid leukemia (AML) and predict a poor prognosis. FLT3 tyrosine kinase inhibitors (TKIs) provide short-term clinical responses, but the long-term prognosis of FLT3/ITD + AML patients remains poor. Notch signaling is important in numerous types of tumors. However, the role of Notch signaling in FLT3/ITD + AML remains to be elucidated. In the current study, we found that Notch signaling was activated upon FLT3-TKI treatment in FLT3/ITD + cell lines and primary cells. As Notch signaling can be blocked by γ-secretase inhibitors (GSIs), we examined the combinatorial antitumor efficacy of FLT3-TKIs and GSIs against FLT3/ITD + AML and explored the underlying molecular mechanisms. As a result, we observed synergistic cytotoxic effects, and the treatment preferentially reduced cell proliferation and induced apoptosis in FLT3/ITD + AML cell lines and in primary AML cells. Furthermore, the combination of FLT3-TKI and GSI eradicated leukemic cells and prolonged survival in an FLT3/ITD + patient-derived xenograft AML model. Mechanistically, differential expression analysis suggested that CXCR3 may be partially responsible for the observed synergy, possibly through ERK signaling. Our findings suggest that combined therapies of FLT3-TKIs with GSI may be exploited as a potential therapeutic strategy to treat FLT3/ITD + AML.

We describe a 5-year-old male with developmental delay, behavioral problems, and dysmorphic features who was found by microarray to have a 93-kb duplication of uncertain significance that fully encompasses the third exon of CTNND2 (delta catenin). Mate-pair sequencing was used to determine that the duplication is tandem and is predicted to lead to CTNND2 haploinsufficiency. Haploinsufficiency for CTNND2 has been shown to result in developmental delay and intellectual disability, providing a unifying diagnosis for this patient. His features overlap those associated with the larger cri-du-chat deletion of this region, expanding the clinical phenotype of isolated CTNND2 variants. The use of mate-pair sequencing to determine the orientation of the small duplication was essential to the diagnosis and avoided the use of exome sequencing, which would not have defined the arrangement of the duplication. This is only the second reported patient, to our knowledge, with a single exon duplication of CTNND2., (© 2019 Wiley Periodicals, Inc.)

Mitochondrial function is critical in eukaryotes. To maintain an adequate supply of energy, precise interactions must be maintained between nuclear- and mitochondrial-encoded gene products. Such interactions are paramount in chimeric enzymes such as the oxidative phosphorylation (OXPHOS) complexes. Mutualistic coevolution between the two genomes has therefore been suggested to be a critical, ubiquitous feature of eukaryotes that acts to maintain cellular function. However, mitochondrial genomes can also act selfishly and increase their own transmission at the expense of organismal function. For example, male-harming mutations are predisposed to accumulate in mitochondrial genomes due to their maternal inheritance ("mother's curse"). Here, we investigate sexually antagonistic mitonuclear coevolution in nuclear-encoded OXPHOS paralogs from mammals and Drosophila. These duplicate genes are highly divergent but must interact with the same set of mitochondrial-encoded genes. Many such paralogs show testis-specific expression, prompting previous hypotheses suggesting they may have evolved under selection to counteract male-harming mitochondrial mutations. We found increased rates of evolution in OXPHOS paralogs with testis-specific expression in mammals and Drosophila, supporting this hypothesis. However, further analyses suggested such patterns may be due to relaxed, not positive selection, especially in Drosophila. Structural data also suggest that mitonuclear interactions do not play a major role in the evolution of many OXPHOS paralogs in a consistent way. In conclusion, no single OXPHOS paralog met all our criteria for being under selection to counteract male-harming mitochondrial mutations. We discuss alternative explanations for the drastic patterns of evolution in these genes, including mutualistic mitonuclear coevolution, adaptive subfunctionalization after gene duplication, and relaxed selection on OXPHOS in male tissues., (© The Author(s) 2019. Published by Oxford University Press on behalf of the Society for Integrative and Comparative Biology. All rights reserved. For permissions please email: journals.permissions@oup.com.)

Old and young duplicate genes have been reported in some organisms. However, little is known about the properties of old and young duplicate genes in Arachis. Here, we have identified old and young duplicate genes in Arachis duranensis, and analyzed the evolution, gene complexity, gene expression pattern, and functional divergence between old and young duplicate genes. Our results showed different evolutionary, gene complexity and gene expression patterns, as well as differing correlations between old and young duplicate genes. Gene ontology results showed that old duplicate genes play a crucial role in lipid and amino acid biosynthesis and the oxidation-reduction process and that young duplicate genes are preferentially involved in photosynthesis and response to biotic stimulus. Transcriptome data sets revealed that most old and young duplicate genes had asymmetric function, and only a few duplicate genes exhibited symmetric function under drought and nematode stress. We found that old duplicate genes are preferentially involved in lipid and amino acid metabolism and response to abiotic stress, while young duplicate genes are likely to participate in photosynthesis and response to biotic stress. This work provides a better understanding of the evolution and functional divergence of old and young duplicate genes in A. duranensis.

The uniquely human α7-nAChR gene (CHRFAM7A) is evolved from the fusion of two partially duplicated genes, FAM7 and α7-nAChR gene (CHRNA7), and is inserted on same chromosome 15, 5' end of the CHRNA7 gene. Transcription of CHRFAM7A gene produces a 1256-bp open reading frame encoding dup-α7-nAChR, where a 27-aminoacid residues from FAM7 replaced the 146-aminoacid residues of the N-terminal extracellular ligand binding domain of α7-nAChR. In vitro, dup-α7-nAChR has been shown to form hetero-pentamer with α7-nAChR and dominant-negatively regulates the channel functions of α7-nAChR. However, the contribution of CHRFAM7A gene to the biology of α7-nAChR in the brain in vivo remains largely a matter of conjecture. CHRFAM7A transgenic mouse was created and differentially expressed proteins were profiled from the whole brain using iTRAQ-2D-LC-MS/MS proteomic technology. Proteins with a fold change of ≥1.2 or ≤0.83 and p < 0.05 were considered to be significant. Bioinformatics analysis showed that over-expression of the CHRFAM7A gene significantly modulated the proteins commonly involved in the signaling pathways of α7-nAChR-mediated neuropsychiatric disorders including Parkinson's disease, Alzheimer's disease, Huntington's disease, and alcoholism, suggesting that the CHRFAM7A gene contributes to the pathogenesis of neuropsychiatric disorders mostly likely through fine-tuning the functions of α7-nAChR in the brain., (Copyright © 2019 Elsevier B.V. All rights reserved.)

Pilocytic astrocytoma (PA) is emerging as a tumor entity with dysregulated RAS/RAF/MEK/ERK signaling. In this study, we report the identification of a novel recurrent BRAF insertion (p.V504&#95;R506dup) in five PA cases harboring exclusively this somatic tandem duplication. This recurrent alteration leads to an addition of three amino acids in the kinase domain of BRAF and has functional impact on activating MAPK phosphorylation. Importantly, we show that this mutation confers resistance to RAF inhibitors without changing effectiveness while downstream MEK inhibitors remain effective. Our results further emphasize the importance of BRAF alterations in PA and the need to characterize them in a given tumor as this can affect therapeutic strategies and their potential use as tumor marker in molecular diagnostics.

During the evolution of multicellular eukaryotes, gene duplication occurs frequently to generate new genes and/or functions. A duplicated gene may have a similar function to its ancestral gene. Therefore, it may be expected that duplicated genes are less likely to be critical for the survival of an organism, since there are multiple copies of the gene rendering each individual copy redundant. In this study, we explored the developmental expression patterns of duplicate gene pairs and the relationship between development co-expression and phenotypes resulting from the knockout of duplicate genes in the mouse. We define genes that generate lethal phenotypes in single gene knockout experiments as essential genes. We found that duplicate gene pairs comprised of two essential genes tend to be expressed at different stages of development, compared to duplicate gene pairs with at least one non-essential member, showing that the timing of developmental expression affects the ability of one paralogue to compensate for the loss of the other. Gene essentiality, developmental expression and gene duplication are thus closely linked.

Gene duplication is an important process in the evolution of gene content in eukaryotic genomes. Understanding when gene duplicates contribute new molecular functions to genomes through molecular adaptation is one important goal in comparative genomics. In large gene families, however, characterizing adaptation and neofunctionalization across species is challenging, as models have traditionally quantified the timing of duplications without considering underlying gene trees. This protocol combines multiple approaches to detect adaptation in protein duplicates at a phylogenetic scale. We include a description of models for gene tree-species tree reconciliation that enable different types of inference, as well as a practical guide to their use. Although simulation-based approaches successfully detect shifts in the rate of duplication/retention, the conflation between the duplication and retention processes, the distinct trajectories of duplicates under non-, sub-, and neofunctionalization, as well as dosage effects offer hitherto unexplored analytical avenues. We introduce mathematical descriptions of these probabilities and offer a road map to computational implementation whose starting point is parsimony reconciliation. Sequence evolution information based on the ratio of nonsynonymous to synonymous nucleotide substitution rates (dN/dS) can be combined with duplicate survival probabilities to better predict the emergence of new molecular functions in retained duplicates. Together, these methods enable characterization of potentially adaptive candidate duplicates whose neofunctionalization may contribute to phenotypic divergence across species.

Whole genome duplication, a prevalent force of evolution in plants, results in massive genome restructuring in different organisms. Roles of the resultant duplicated genes are poorly understood, both functionally and evolutionarily. In the present study, differentially expressed ethylene responsive factors (GhERF1s), anchored on Chr-A07 and Chr-D07, were isolated from a high-yielding cotton hybrid (XZM2) and its parents. The GhERF1 was located in the B3 subgroup of the ethylene responsive factors subfamily involved in conferring tolerance to abiotic stress. Nucleotide sequence analysis of 524 diverse accessions, together with quantitative real-time polymerase chain reaction analysis, elucidated that de-functionalization of GhERF1-7A occurred due to one base insertion following formation of the allotetraploid cotton. Our quantitative trait loci and association mapping analyses highlighted a role for GhERF1-7A in conferring high boll number per plant in modern cotton cultivars. Overexpression of GhERF1-7A in transgenic Arabidopsis resulted in a substantial increase in the number of siliques and total seed yield. Neo-functionalization of GhERF1-7A was also observed in modern cultivars rather than in races and/or landraces, further supporting its role in the development of high-yielding cotton cultivars. Both de- and neo-functionalization occurred in one of the duplicate genes, thus providing new genomic insight into the evolution of allotetraploid cotton species., (© 2018 Institute of Botany, Chinese Academy of Sciences.)

FACT (facilitates chromatin transcription) is a histone chaperone complex important in genomic processes including transcription, DNA replication, and DNA repair. FACT is composed of two proteins, SSRP1 and SPT16, which are highly conserved across eukaryotes. While the mechanisms for FACT in nucleosome reorganization and its relationship to DNA processes is well established, how these roles impact coordination in multicellular animal development are less well understood. Here we characterize the genes encoding FACT complex proteins in the nematode C. elegans. We show that whereas C. elegans includes one SPT16 gene (spt-16), two genes (hmg-3 and hmg-4) encode SSRP1 proteins. Depletion of FACT complex genes interferes with embryonic cell division and cell cycle timing generally, with anterior pharynx development especially sensitive to these defects. hmg-3 and hmg-4 exhibit redundancy for these maternally-provided embryonic functions, but are each uniquely required zygotically for normal germline development. This work provides a framework to study FACT gene function in developmental processes, and identifies that distinct functional requirements for gene duplicates can be manifest within a single tissue., (Copyright © 2018 Elsevier Inc. All rights reserved.)

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

In salmonids, the majority of circulating insulin-like growth factor-I (IGF-I) is bound to IGF binding proteins (IGFBP), with IGFBP-2b being the most abundant in circulation. We used CRISPR/Cas9 methodology to disrupt expression of a functional IGFBP-2b protein by co-targeting for gene editing IGFBP-2b1 and IGFBP-2b2 subtypes, which represent salmonid-specific gene duplicates. Twenty-four rainbow trout were produced with mutations in the IGFBP-2b1 and IGFBP-2b2 genes. Mutant fish exhibited between 8-100% and 2-83% gene disruption for IGFBP-2b1 and IGFBP-2b2, respectively, with a positive correlation (P < 0.001) in gene mutation rate between individual fish. Analysis of IGFBP-2b protein indicated reductions in plasma IGFBP-2b abundance to between 0.04-0.96-fold of control levels. Plasma IGF-I, body weight, and fork length were reduced in mutants at 8 and 10 months post-hatch, which supports that IGFBP-2b is significant for carrying IGF-I. Despite reduced plasma IGF-I and IGFBP-2b in mutants, growth retardation in mutants was less severe between 10 and 12 months post-hatch (P < 0.05), suggesting a compensatory growth response occurred. These findings indicate that gene editing using CRISPR/Cas9 and ligand blotting is a feasible approach for characterizing protein-level functions of duplicated IGFBP genes in salmonids and is useful to unravel IGF-related endocrine mechanisms.

Variation in genome content is a potent mechanism of microbial adaptation. The genomes of members of the cyanobacterial genus Acaryochloris vary greatly in gene content as a consequence of the idiosyncratic retention of both recent gene duplicates and plasmid-encoded genes acquired by horizontal transfer. For example, the genome of Acaryochloris strain MBIC11017, which was isolated from an iron-limited environment, is enriched in duplicated and novel genes involved in iron assimilation. Here, we took an integrative approach to characterize the adaptation of Acaryochloris MBIC11017 to low environmental iron availability and the relative contributions of the expression of duplicated versus novel genes. We observed that Acaryochloris MBIC11017 grew faster and to a higher yield in the presence of nanomolar concentrations of iron than did a closely related strain. These differences were associated with both a higher rate of iron assimilation and a greater abundance of iron assimilation transcripts. However, recently duplicated genes contributed little to increased transcript dosage; rather, the maintenance of these duplicates in the MBIC11017 genome is likely due to the sharing of ancestral dosage by expression reduction. Instead, novel, horizontally transferred genes are responsible for the differences in transcript abundance. The study provides insights on the mechanisms of adaptive genome evolution and gene expression in Acaryochloris.

Males and females have different fitness optima but share the vast majority of their genomes, causing an inherent genetic conflict between the two sexes that must be resolved to achieve maximal population fitness. We show that two tandem duplicate genes found specifically in Drosophila melanogaster are sexually antagonistic, but rapidly evolved sex-specific functions and expression patterns that mitigate their antagonistic effects. We use copy-specific knockouts and rescue experiments to show that Apollo (Apl) is essential for male fertility but detrimental to female fertility, in addition to its important role in development, while Artemis (Arts) is essential for female fertility but detrimental to male fertility. Further analyses show that Apl and Arts have essential roles in spermatogenesis and oogenesis. These duplicates formed ~200,000 years ago, underwent a strong selective sweep and lost most expression in the antagonized sex. These data provide direct evidence that gene duplication allowed rapid mitigation of sexual conflict by allowing Apl and Arts to evolve essential sex-specific reproductive functions and complementary expression in male and female gonads.

Key Message: Using a gRNA and Agrobacterium-mediated transformation, we performed simultaneous site-directed mutagenesis of two GmPPD loci in soybean. Mutations in GmPPD loci were confirmed in at least 33% of T 2 seeds. The clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated endonuclease 9 (Cas9) system is a powerful tool for site-directed mutagenesis in crops. Using a single guide RNA (gRNA) and Agrobacterium-mediated transformation, we performed simultaneous site-directed mutagenesis of two homoeologous loci in soybean (Glycine max), GmPPD1 and GmPPD2, which encode the orthologs of Arabidopsis thaliana PEAPOD (PPD). Most of the T 1 plants had heterozygous and/or chimeric mutations for the targeted loci. The sequencing analysis of T 1 and T 2 generations indicates that putative mutation induced in the T 0 plant is transmitted to the T 1 generation. The inheritable mutation induced in the T 1 plant was also detected. This result indicates that continuous induction of mutations during T 1 plant development increases the occurrence of mutations in germ cells, which ensures the transmission of mutations to the next generation. Simultaneous site-directed mutagenesis in both GmPPD loci was confirmed in at least 33% of T 2 seeds examined. Approximately 19% of double mutants did not contain the Cas9/gRNA expression construct. Double mutants with frameshift mutations in both GmPPD1 and GmPPD2 had dome-shaped trifoliate leaves, extremely twisted pods, and produced few seeds. Taken together, our data indicate that continuous induction of mutations in the whole plant and advancing generations of transgenic plants enable efficient simultaneous site-directed mutagenesis in duplicated loci in soybean.

An attractive and long-standing hypothesis regarding the evolution of genes after duplication posits that the duplication event creates new evolutionary possibilities by releasing a copy of the gene from constraint. Apparent support was found in numerous analyses, particularly, the observation of higher rates of evolution in duplicated as compared with singleton genes. Could it, instead, be that more duplicable genes (owing to mutation, fixation, or retention biases) are intrinsically faster evolving? To uncouple the measurement of rates of evolution from the determination of duplicate or singleton status, we measure the rates of evolution in singleton genes in outgroup primate lineages but classify these genes as to whether they have duplicated or not in a crown group of great apes. We find that rates of evolution are higher in duplicable genes prior to the duplication event. In part this is owing to a negative correlation between coding sequence length and rate of evolution, coupled with a bias toward smaller genes being more duplicable. The effect is masked by difference in expression rate between duplicable genes and singletons. Additionally, in contradiction to the classical assumption, we find no convincing evidence for an increase in dN/dS after duplication, nor for rate asymmetry between duplicates. We conclude that high rates of evolution of duplicated genes are not solely a consequence of the duplication event, but are rather a predictor of duplicability. These results are consistent with a model in which successful gene duplication events in mammals are skewed toward events of minimal phenotypic impact., (© The Author 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Salmonella enterica serovar Dublin is adapted to cattle but is able to infect humans with high invasiveness. An acute inflammatory response at the intestine helps to prevent Salmonella dissemination to systemic sites. Flagella contribute to this response by providing motility and FliC-mediated signaling through pattern recognition receptors. In a previous work, we reported a high frequency (11 out of 25) of S Dublin isolates lacking flagella in a collection obtained from humans and cattle. The aflagellate strains were impaired in their proinflammatory properties in vitro and in vivo The aim of this work was to elucidate the underlying cause of the absence of flagella in S Dublin isolates. We report here that class 3 flagellar genes are repressed in the human aflagellate isolates, due to impaired secretion of FliA anti-sigma factor FlgM. This phenotype is due to an in-frame 42-nucleotide deletion in the fliE gene, which codes for a protein located in the flagellar basal body. The deletion is predicted to produce a protein lacking amino acids 18 to 31. The aflagellate phenotype was highly stable; revertants were obtained only when fliA was artificially overexpressed combined with several successive passages in motility agar. DNA sequence analysis revealed that motile revertants resulted from duplications of DNA sequences in fliE adjacent to the deleted region. These duplications produced a FliE protein of similar length to the wild type and demonstrate that amino acids 18 to 31 of FliE are not essential. The same deletion was detected in S Dublin isolates obtained from cattle, indicating that this mutation circulates in nature., (Copyright © 2017 American Society for Microbiology.)

Pairs of duplicated genes generally display a combination of conserved expression patterns inherited from their unduplicated ancestor and newly acquired domains. However, how the cis-regulatory architecture of duplicated loci evolves to produce these expression patterns is poorly understood. We have directly examined the gene-regulatory evolution of two tandem duplicates, the Drosophila Ly6 genes CG9336 and CG9338, which arose at the base of the drosophilids between 40 and 60 Ma. Comparing the expression patterns of the two paralogs in four Drosophila species with that of the unduplicated ortholog in the tephritid Ceratitis capitata, we show that they diverged from each other as well as from the unduplicated ortholog. Moreover, the expression divergence appears to have occurred close to the duplication event and also more recently in a lineage-specific manner. The comparison of the tissue-specific cis-regulatory modules (CRMs) controlling the paralog expression in the four Drosophila species indicates that diverse cis-regulatory mechanisms, including the novel tissue-specific enhancers, differential inactivation, and enhancer sharing, contributed to the expression evolution. Our analysis also reveals a surprisingly variable cis-regulatory architecture, in which the CRMs driving conserved expression domains change in number, location, and specificity. Altogether, this study provides a detailed historical account that uncovers a highly dynamic picture of how the paralog expression patterns and their underlying cis-regulatory landscape evolve. We argue that our findings will encourage studying cis-regulatory evolution at the whole-locus level to understand how interactions between enhancers and other regulatory levels shape the evolution of gene expression., (© The Author 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Gene duplication is thought to play a major role in phenotypic evolution. Yet the forces involved in the functional divergence of young duplicate genes remain unclear. Here, we use population-genetic inference to elucidate the role of natural selection in the functional evolution of young duplicate genes in Drosophila melanogaster. We find that negative selection acts on young duplicates with ancestral functions, and positive selection on those with novel functions, suggesting that natural selection may determine whether and how young duplicate genes are retained. Moreover, evidence of natural selection is strongest in protein-coding regions and 3' UTRs of young duplicates, indicating that selection may primarily target encoded proteins and regulatory sequences specific to 3' UTRs. Further analysis reveals that natural selection acts immediately after duplication and weakens over time, possibly explaining the observed bias toward the acquisition of new functions by young, rather than old, duplicate gene copies. Last, we find an enrichment of testis-related functions in young duplicates that underwent recent positive selection, but not in young duplicates that did not undergo recent positive selection, or in old duplicates that either did or did not undergo recent positive selection. Thus, our findings reveal that natural selection is a key player in the functional evolution of young duplicate genes, acts rapidly and in a region-specific manner, and may underlie the origin of novel testis-specific phenotypes in Drosophila., (© The Author 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

In several organisms, particular functional categories of genes, such as regulatory and complex-forming genes, are preferentially retained after whole-genome multiplications but rarely duplicate through small-scale duplication, a pattern referred to as reciprocal retention. This peculiar duplication behavior is hypothesized to stem from constraints on the dosage balance between the genes concerned and their interaction context. However, the evidence for a relationship between reciprocal retention and dosage balance sensitivity remains fragmentary. Here, we identified which gene families are most strongly reciprocally retained in the angiosperm lineage and studied their functional and evolutionary characteristics. Reciprocally retained gene families exhibit stronger sequence divergence constraints and lower rates of functional and expression divergence than other gene families, suggesting that dosage balance sensitivity is a general characteristic of reciprocally retained genes. Gene families functioning in regulatory and signaling processes are much more strongly represented at the top of the reciprocal retention ranking than those functioning in multiprotein complexes, suggesting that regulatory imbalances may lead to stronger fitness effects than classical stoichiometric protein complex imbalances. Finally, reciprocally retained duplicates are often subject to dosage balance constraints for prolonged evolutionary times, which may have repercussions for the ease with which genome multiplications can engender evolutionary innovation., (© 2017 American Society of Plant Biologists. All rights reserved.)

Gene and genome duplications have been rampant during the evolution of flowering plants. Unlike small-scale gene duplications, whole-genome duplications (WGDs) copy entire pathways or networks, and as such create the unique situation in which such duplicated pathways or networks could evolve novel functionality through the coordinated sub- or neofunctionalization of its constituent genes. Here, we describe a remarkable case of coordinated gene expression divergence following WGDs in Arabidopsis thaliana We identified a set of 92 homoeologous gene pairs that all show a similar pattern of tissue-specific gene expression divergence following WGD, with one homoeolog showing predominant expression in aerial tissues and the other homoeolog showing biased expression in tip-growth tissues. We provide evidence that this pattern of gene expression divergence seems to involve genes with a role in cell polarity and that likely function in the maintenance of cell wall integrity. Following WGD, many of these duplicated genes evolved separate functions through subfunctionalization in growth/development and stress response. Uncoupling these processes through genome duplications likely provided important adaptations with respect to growth and morphogenesis and defense against biotic and abiotic stress., (© 2017 American Society of Plant Biologists. All rights reserved.)

Single-gene knockout experiments can fail to reveal function in the context of redundancy, which is frequently observed among duplicated genes (paralogs) with overlapping functions. We discuss the complexity associated with studying paralogs and outline how recent advances in CRISPR will help address the "phenotype gap" and impact biomedical research., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Increasing organismal complexity during the evolution of life has been attributed to the duplication of genes and entire genomes. More recently, theoretical models have been proposed that postulate the fate of duplicated genes, among them the duplication-degeneration-complementation (DDC) model. In the DDC model, the common fate of a duplicated gene is lost from the genome owing to nonfunctionalization. Duplicated genes are retained in the genome either by subfunctionalization, where the functions of the ancestral gene are sub-divided between the sister duplicate genes, or by neofunctionalization, where one of the duplicate genes acquires a new function. Both processes occur either by loss or gain of regulatory elements in the promoters of duplicated genes. Here, we review the genomic organization, evolution, and transcriptional regulation of the multigene family of intracellular lipid-binding protein (iLBP) genes from teleost fishes. Teleost fishes possess many copies of iLBP genes owing to a whole genome duplication (WGD) early in the teleost fish radiation. Moreover, the retention of duplicated iLBP genes is substantially higher than the retention of all other genes duplicated in the teleost genome. The fatty acid-binding protein genes, a subfamily of the iLBP multigene family in zebrafish, are differentially regulated by peroxisome proliferator-activated receptor (PPAR) isoforms, which may account for the retention of iLBP genes in the zebrafish genome by the process of subfunctionalization of cis-acting regulatory elements in iLBP gene promoters.

Evolution of the phenolic metabolism was critical for the transition of plants from water to land. A cytochrome P450, CYP73, with cinnamate 4-hydroxylase (C4H) activity, catalyzes the first plant-specific and rate-limiting step in this pathway. The CYP73 gene is absent from green algae, and first detected in bryophytes. A CYP73 duplication occurred in the ancestor of seed plants and was retained in Taxaceae and most angiosperms. In spite of a clear divergence in primary sequence, both paralogs can fulfill comparable cinnamate hydroxylase roles both in vitro and in vivo. One of them seems dedicated to the biosynthesis of lignin precursors. Its N-terminus forms a single membrane spanning helix and its properties and length are highly constrained. The second is characterized by an elongated and variable N-terminus, reminiscent of ancestral CYP73s. Using as proxies the Brachypodium distachyon proteins, we show that the elongation of the N-terminus does not result in an altered subcellular localization, but in a distinct membrane topology. Insertion in the membrane of endoplasmic reticulum via a double-spanning open hairpin structure allows reorientation to the lumen of the catalytic domain of the protein. In agreement with participation to a different functional unit and supramolecular organization, the protein displays modified heme proximal surface. These data suggest the evolution of divergent C4H enzymes feeding different branches of the phenolic network in seed plants. It shows that specialization required for retention of gene duplicates may result from altered protein topology rather than change in enzyme activity., (© The Author 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

It is now recognized that several rounds of whole genome duplication (WGD) have occurred during the evolution of vertebrates, but the link between WGDs and phenotypic diversification remains unsolved. We have investigated in this study the impact of the teleost-specific WGD on the evolution of the sox gene family in teleostean fishes. The sox gene family, which encodes for transcription factors, has essential role in morphology, physiology and behavior of vertebrates and teleosts, the current largest group of vertebrates. We have first redrawn the evolution of all sox genes identified in eleven teleost genomes using a comparative genomic approach including phylogenetic and synteny analyses. We noticed, compared to tetrapods, an important expansion of the sox family: 58% (11/19) of sox genes are duplicated in teleost genomes. Furthermore, all duplicated sox genes, except sox17 paralogs, are derived from the teleost-specific WGD. Then, focusing on five sox genes, analyzing the evolution of coding and non-coding sequences, as well as the expression patterns in fish embryos and adult tissues, we demonstrated that these paralogs followed lineage-specific evolutionary trajectories in teleost genomes. This work, based on whole genome data from multiple teleostean species, supports the contribution of WGDs to the expansion of gene families, as well as to the emergence of genomic differences between lineages that might promote genetic and phenotypic diversity in teleosts.

Common models for the evolution of duplicated genes after genome duplication are subfunctionalization, neofunctionalization, and pseudogenization. Although the crucial roles of cis-regulatory mutations in subfunctionalization are well-documented, their involvement in pseudogenization and/or neofunctionalization remains unclear. We addressed this issue by investigating the evolution of duplicated homeobox genes, six6.L and six6.S, in the allotetraploid frog Xenopus laevis. Based on a comparative expression analysis, we observed similar eye-specific expression patterns for the two loci and their single ortholog in the ancestral-type diploid species Xenopus tropicalis. However, we detected lower levels of six6.S expression than six6.L expression. The six6.S enhancer sequence was more highly diverged from the orthologous enhancer of X. tropicalis than the six6.L enhancer, and showed weaker activity in a transgenic reporter assay. Based on a phylogenetic analysis of the protein sequences, we observed greater divergence between X. tropicalis Six6 and Six6.S than between X. tropicalis Six6 and Six6.L, and the observed mutations were reminiscent of a microphthalmia mutation in human SIX6. Misexpression experiments showed that six6.S has weaker eye-enlarging activity than six6.L, and targeted disruption of six6.L reduced the eye size more significantly than that of six6.S. These results suggest that enhancer attenuation stimulates the accumulation of hypomorphic coding mutations, or vice versa, in one duplicated gene copy and facilitates pseudogenization. We also underscore the value of the allotetraploid genome of X. laevis as a resource for studying latent pathogenic mutations., (Copyright © 2017 Elsevier Inc. All rights reserved.)