Distinguishing among evolutionary models for the maintenance of gene duplicates

High-throughput amino acid-level characterization of the interactions of plasminogen activator inhibitor-1 with variably divergent proteases

While members of large paralogous protein families share structural features, their functional niches often diverge significantly. Serine protease inhibitors (SERPINs), whose members typically function as covalent inhibitors of serine proteases, are one such family. Plasminogen activator inhibitor-1 (PAI-1) is a prototypic SERPIN, which canonically inhibits tissue-and urokinase-type plasminogen activators (tPA and uPA) to regulate fibrinolysis. PAI-1 has been shown to also inhibit other serine proteases, including coagulation factor XIIa (FXIIa) and transmembrane serine protease 2 (TMPRSS2). The structural determinants of PAI-1 inhibitory function toward these non-canonical protease targets, and the biological significance of these functions, are unknown. We applied deep mutational scanning (DMS) to assess the effects of ∼80% of all possible single amino acid substitutions in PAI-1 on its ability to inhibit three putative serine protease targets (uPA, FXIIa, and TMPRSS2). Selection with each target protease generated a unique PAI-1 mutational landscape, with the determinants of protease specificity distributed throughout PAI-1's primary sequence. Next, we conducted a comparative analysis of extant orthologous sequences, demonstrating that key residues modulating PAI-1 inhibition of uPA and FXIIa, but not TMPRSS2, are maintained by purifying selection. PAI-1's activity toward FXIIa may reflect how protease evolutionary relationships predict SERPIN functional divergence, which we support via a cophylogenetic analysis of all secreted SERPINs and their cognate serine proteases. This work provides insight into the functional diversification of SERPINs and lays the framework for extending these studies to other proteases and their regulators.

Cytochrome P450 monooxygenase systems: Diversity and plasticity for adaptive stress response

Superfamily of cytochromes P450 (CYPs) is composed of heme-thiolate-containing monooxygenase enzymes, which play crucial roles in the biosynthesis, bioactivation, and detoxification of a variety of organic compounds, both endogenic and exogenic. Majority of CYP monooxygenase systems are multi-component and contain various redox partners, cofactors and auxiliary proteins, which contribute to their diversity in both prokaryotes and eukaryotes. Recent progress in bioinformatics and computational biology approaches make it possible to undertake whole-genome and phylogenetic analyses of CYPomes of a variety of organisms. Considerable variations in sequences within and between CYP families and high similarity in secondary and tertiary structures between all CYPs along with dramatic conformational changes in secondary structure elements of a substrate binding site during catalysis have been reported. This provides structural plasticity and substrate promiscuity, which underlie functional diversity of CYPs. Gene duplication and mutation events underlie CYP evolutionary diversity and emergence of novel selectable functions, which provide the involvement of CYPs in high adaptability to changing environmental conditions and dietary restrictions. In our review, we discuss the recent advancements and challenges in the elucidating the evolutionary origin and mechanisms underlying the CYP monooxygenase system diversity and plasticity. Our review is in the view of hypothesis that diversity of CYP monooxygenase systems is translated into the broad metabolic profiles, and this has been acquired during the long evolutionary time to provide structural plasticity leading to high adaptative capabilities to environmental stress conditions.

A Rigorous Framework to Classify the Postduplication Fate of Paralogous Genes

Gene duplication has a central role in evolution; still, little is known on the fates of the duplicated copies, their relative frequency, and on how environmental conditions affect them. Moreover, the lack of rigorous definitions concerning the fate of duplicated genes hinders the development of a global vision of this process. In this paper, we present a new framework aiming at characterizing and formally differentiating the fate of duplicated genes. Our framework has been tested via simulations, where the evolution of populations has been simulated using aevol, an in silico experimental evolution platform. Our results show several patterns that confirm some of the conclusions from previous studies, while also exhibiting new tendencies; this may open up new avenues to better understand the role of duplications as a driver of evolution.

The roles of gene duplications in the dynamics of evolutionary conflicts

Evolutionary conflicts occur when there is antagonistic selection between different individuals of the same or different species, life stages or between levels of biological organization. Remarkably, conflicts can occur within species or within genomes. In the dynamics of evolutionary conflicts, gene duplications can play a major role because they can bring very specific changes to the genome: changes in protein dose, the generation of novel paralogues with different functions or expression patterns or the evolution of small antisense RNAs. As we describe here, by having those effects, gene duplication might spark evolutionary conflict or fuel arms race dynamics that takes place during conflicts. Interestingly, gene duplication can also contribute to the resolution of a within-locus evolutionary conflict by partitioning the functions of the gene that is under an evolutionary trade-off. In this review, we focus on intraspecific conflicts, including sexual conflict and illustrate the various roles of gene duplications with a compilation of examples. These examples reveal the level of complexity and the differences in the patterns of gene duplications within genomes under different conflicts. These examples also reveal the gene ontologies involved in conflict and the genomic location of the elements of the conflict. The examples provide a blueprint for the direct study of these conflicts or the exploration of the presence of similar conflicts in other lineages.

Genome-wide identification and expression analysis of protein arginine methyltransferase and JmjC domain-containing family in apple

Histone methylation is an important type of histone modification that regulates gene expression in plants. In this study, we identified 14 arginine methylation-related genes (Protein Arginine Methyltransferase, MdPRMT) and 32 demethylation-related genes (JmjC Domain-Containing Family, MdJMJ) in apple. Furthermore, we investigated the phylogenetic relationship, chromosome distribution, gene structure, motif analysis, promoter sequence analysis, and expression patterns of MdPRMT and MdJMJ genes. Homology analysis showed a high degree of conservation and homology between PRMT and JMJ genes in Arabidopsis and apple. We identified the types of duplicated genes in the MdJMJ and MdPRMT gene families, found a large number of whole-genome duplicates (WGD) gene pairs and a small number of tandem duplicates (TD) pairs, transposed duplication (TRD) gene pairs as well as proximal duplicates (PD) pairs, and discussed the possible evolutionary pathways of the gene families from the perspective of duplicated genes. Homology analysis showed a high degree of conservation and homology between PRMT and JMJ genes in Arabidopsis and apple. In addition, the promoter regions of MdPRMT and MdJMJ contain numerous cis-acting elements involved in plant growth and development, hormone response, and stress responses. Based on the transcriptional profiles of MdPRMT and MdJMJ in different tissues and developmental stages, it was found that MdPRMT and MdJMJ may play multiple roles in apple growth and development, for example, MdJMJ21 may be involved in the regulation of apple endosperm formation. MdPRMT and MdJMJ exhibit different expression patterns in response to hormone signaling in apple, MdJMJ3, MdJMJ18, MdJMJ30, MdPRMT2, MdPRMT13, and MdPRMT14 may play roles in apple response to drought stress, while the expression of MdJMJ13, MdPRMT3, MdPRMT4, and MdPRMT6 is affected by cold stress. Our study provides a foundation for determining the functional roles of MdPRMT and MdJMJ genes in apple.

Adaptive evolution of chloroplast division mechanisms during plant terrestrialization

Despite extensive research, the origin and evolution of the chloroplast division machinery remain unclear. Here, we employ recently sequenced genomes and transcriptomes of Archaeplastida clades to identify the core components of chloroplast division and reconstruct their evolutionary histories, respectively. Our findings show that complete division ring structures emerged in Charophytes. We find that Glaucophytes experienced strong selection pressure, generating diverse variants adapted to the changing terrestrial environments. By integrating the functions of chloroplast division genes (CDGs) annotated in a workflow developed using large-scale multi-omics data, we further show that dispersed duplications acquire more species-specific functions under stronger selection pressures. Notably, PARC6, a dispersed duplicate CDG, regulates leaf color and plant growth in Solanum lycopersicum, demonstrating neofunctionalization. Our findings provide an integrated perspective on the functional evolution of chloroplast division machinery and highlight the potential of dispersed duplicate genes as the primary source of adaptive evolution of chloroplast division.

Models for the retention of duplicate genes and their biological underpinnings

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

Molecular cloning and tissue distribution of glucokinase and glucose-6-phosphatase catalytic subunit paralogs in largemouth bass Micropterus salmoides: Regulation by dietary starch levels and a glucose load

The dysregulation of glucose-G6P (glucose-6-phosphate) interconversion is thought to be one of the main reasons for the low glucose disposal of carnivorous fish, but is not yet well understood in largemouth bass Micropterus salmoides (LMB). In this study, the full length cDNA sequences of genes encoding glucokinase (Gck, catalyzing glucose phosphorylation) and glucose-6-phosphatase catalytic subunit (G6pc, catalyzing glucose dephosphorylation) were cloned by the RACE method from the liver of LMB. Subsequently, the distribution of g6pc and gck as well as their transcriptional regulation by dietary starch levels and a glucose load were investigated. Only one gck gene was identified, while the tandem duplication of g6pca.1 gene was named as g6pca.2 in LMB. The full cDNA sequences of g6pca.1, g6pca.2 and gck in LMB were 1585, 1813 and 2115 bp in length, encoding 478, 352 and 359 amino acids, respectively. Gck was predicted to contain two hexokinase domains, an ATP-binding domain and multiple functional sites, while G6pca.1 and G6pca.2 contained nine transmembrane helices, a PAP2 (type-2 phosphatidic acid phosphatase) domain and multiple functional amino acid sites. Both g6pca.1 and g6pca.2 were predominantly distributed in the liver and to some extent in the intraperitoneal fat, intestine and pyloric caeca, while gck was mainly transcribed in the liver and to some extent in the heart, intestine and brain. Both feeding a high starch diet and a glucose load stimulated the mRNA expression of gck in the liver of LMB. An increase of dietary starch from 9% to 14% down-regulated the transcription of g6pca.1 in the liver of LMB. However, both the mRNA levels of hepatic g6pca.1 and g6pca.2 were sharply up-regulated in LMB during 1-3 h after a glucose load. Overall, the results of this study suggested that the functions of G6pc (G6pca.1 and G6pca.2) and Gck in LMB were highly conserved in evolution. Though hepatic glucose-G6P interconversion was well regulated at the transcript level in LMB fed high starch diets, a futile cycle between glucose and G6P was induced in the liver after a glucose load.

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

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

Comparative Analysis Reveals Different Evolutionary Fates and Biological Functions in Wheat Duplicated Genes (Triticum aestivum L.)

Wheat (Triticum aestivum L.) is a staple food crop that provides 20% of total human calorie consumption. Gene duplication has been considered to play an important role in evolution by providing new genetic resources. However, the evolutionary fates and biological functions of the duplicated genes in wheat remain to be elucidated. In this study, the resulting data showed that the duplicated genes evolved faster with shorter gene lengths, higher codon usage bias, lower expression levels, and higher tissue specificity when compared to non-duplicated genes. Our analysis further revealed functions of duplicated genes in various biological processes with significant enrichment to environmental stresses. In addition, duplicated genes derived from dispersed, proximal, tandem, transposed, and whole-genome duplication differed in abundance, evolutionary rate, gene compactness, expression pattern, and genetic diversity. Tandem and proximal duplicates experienced stronger selective pressure and showed a more compact gene structure with diverse expression profiles than other duplication modes. Moreover, genes derived from different duplication modes showed an asymmetrical evolutionary pattern for wheat A, B, and D subgenomes. Several candidate duplication hotspots associated with wheat domestication or polyploidization were characterized as potential targets for wheat molecular breeding. Our comprehensive analysis revealed the evolutionary trajectory of duplicated genes and laid the foundation for future functional studies on wheat.

Regulatory and functional divergence among members of Ibβfruct2, a sweet potato vacuolar invertase gene controlling starch and glucose content

Sweet potato [Ipomoea batatas (L.) Lam.] is an important food and industrial crop. Its storage root is rich in starch, which is present in the form of granules and represents the principal storage carbohydrate in plants. Starch content is an important trait of sweet potato controlling the quality and yield of industrial products. Vacuolar invertase encoding gene Ibβfruct2 was supposed to be a key regulator of starch content in sweet potato, but its function and regulation were unclear. In this study, three Ibβfruct2 gene members were detected. Their promoters displayed differences in sequence, activity, and cis-regulatory elements and might interact with different transcription factors, indicating that the three Ibβfruct2 family members are governed by different regulatory mechanisms at the transcription level. Among them, we found that only Ibβfruct2-1 show a high expression level and promoter activity, and encodes a protein with invertase activity, and the conserved domains and three conserved motifs NDPNG, RDP, and WEC are critical to this activity. Only two and six amino acid residue variations were detected in sequences of proteins encoded by Ibβfruct2-2 and Ibβfruct2-3, respectively, compared with Ibβfruct2-1; although not within key motifs, these variations affected protein structure and affinities for the catalytic substrate, resulting in functional deficiency and low activity. Heterologous expression of Ibβfruct2-1 in Arabidopsis decreased starch content but increased glucose content in leaves, indicating Ibβfruct2-1 was a negative regulator of starch content. These findings represent an important advance in understanding the regulatory and functional divergence among duplicated genes in sweet potato, and provide critical information for functional studies and utilization of these genes in genetic improvement.

Nuclear phylogeny and insights into whole-genome duplications and reproductive development of Solanaceae plants

Solanaceae, the nightshade family, have ∼2700 species, including the important crops potato and tomato, ornamentals, and medicinal plants. Several sequenced Solanaceae genomes show evidence for whole-genome duplication (WGD), providing an excellent opportunity to investigate WGD and its impacts. Here, we generated 93 transcriptomes/genomes and combined them with 87 public datasets, for a total of 180 Solanaceae species representing all four subfamilies and 14 of 15 tribes. Nearly 1700 nuclear genes from these transcriptomic/genomic datasets were used to reconstruct a highly resolved Solanaceae phylogenetic tree with six major clades. The Solanaceae tree supports four previously recognized subfamilies (Goetzeioideae, Cestroideae, Nicotianoideae, and Solanoideae) and the designation of three other subfamilies (Schizanthoideae, Schwenckioideae, and Petunioideae), with the placement of several previously unassigned genera. We placed a Solanaceae-specific whole-genome triplication (WGT1) at ∼81 million years ago (mya), before the divergence of Schizanthoideae from other Solanaceae subfamilies at ∼73 mya. In addition, we detected two gene duplication bursts (GDBs) supporting proposed WGD events and four other GDBs. An investigation of the evolutionary histories of homologs of carpel and fruit developmental genes in 14 gene (sub)families revealed that 21 gene clades have retained gene duplicates. These were likely generated by the Solanaceae WGT1 and may have promoted fleshy fruit development. This study presents a well-resolved Solanaceae phylogeny and a new perspective on retained gene duplicates and carpel/fruit development, providing an improved understanding of Solanaceae evolution.

ADH Gene Cloning and Identification of Flooding-Responsive Genes in Taxodium distichum (L.) Rich

As a flooding-tolerant tree species, Taxodium distichum has been utilized in afforestation projects and proven to have important value in flooding areas. Alcohol dehydrogenase (ADH), which participates in ethanol fermentation, is essential for tolerance to the anaerobic conditions caused by flooding. In a comprehensive analysis of the ADH gene family in T. distichum, TdADHs were cloned on the basis of whole-genome sequencing, and then bioinformatic analysis, subcellular localization, and gene expression level analysis under flooding were conducted. The results show that the putative protein sequences of 15 cloned genes contained seven TdADHs and eight TdADH-like genes (one Class III ADH included) that were divided into five clades. All the sequences had an ADH_N domain, and except for TdADH-likeE2, all the other genes had an ADH_zinc_N domain. Moreover, the TdADHs in clades A, B, C, and D had a similar motif composition. Additionally, the number of TdADH amino acids ranged from 277 to 403, with an average of 370.13. Subcellular localization showed that, except for TdADH-likeD3, which was not expressed in the nucleus, the other genes were predominantly expressed in both the nucleus and cytosol. TdADH-likeC2 was significantly upregulated in all three organs (roots, stems, and leaves), and TdADHA3 was also highly upregulated under 24 h flooding treatment; the two genes might play key roles in ethanol fermentation and flooding tolerance. These findings offer a comprehensive understanding of TdADHs and could provide a foundation for the molecular breeding of T. distichum and current research on the molecular mechanisms driving flooding tolerance.

Carbon dioxide equivalent emissions from corn silage fermentation

The European Climate Law recently codified the goal for European climate neutrality by 2050, highlighting the need for sustainable farming practices within a robust and transparent carbon dioxide equivalent (CO₂e) accounting system. In the present study, a series of equations were proposed for the estimation of CO₂e emissions from corn silage fermentation. Systematic review of previous meta-analyses of corn silage fermentation identified the mean and standard deviation statistics for key model inputs of acetic acid, ethanol, lactic acid, ammonia, and volatile-corrected dry matter loss. Estimates of CO₂e emissions were determined for a mock dataset comprising 1,000 iterations of randomly-generated values for each metric in accordance with mean and variance statistics of the source data. Estimates for CO₂e emissions of corn silage based on meta-analysis review of laboratory experiments were 1.9 ± 5.6% (GWP₂₀) and 0.2 ± 5.5% (GWP₁₀₀) of silage dry matter. Furthermore, model results demonstrated a precedent for CO₂ recycling by silage microorganisms, which was supported by genome annotation of strains belonging to common silage species. Linear model equations for GWP₂₀ and GWP₁₀₀ with inputs and outputs in mg kg^-1 silage dry matter were developed, where inputs are acetic acid (A), ethanol (E), lactic acid (L), and volatile corrected dry matter loss (D_V). Linear equations are (for GWP₂₀; Eq. 11): GWP 20 = - 3626.1 - 0.04343 A + 0.8011 E - 0.03173 L + 1.46573 D V and for GWP₁₀₀; Eq. 12: GWP 100 = - 8526.10 - 0.22403 A - 0.11963 E - 0.03173 L + 1.46573 D V . .

Twenty-seven ZAD-ZNF genes of Drosophila melanogaster are orthologous to the embryo polarity determining mosquito gene cucoid

The C2H2 zinc finger gene cucoid establishes anterior-posterior (AP) polarity in the early embryo of culicine mosquitoes. This gene is unrelated to genes that establish embryo polarity in other fly species (Diptera), such as the homeobox gene bicoid, which serves this function in the traditional model organism Drosophila melanogaster. The cucoid gene is a conserved single copy gene across lower dipterans but nothing is known about its function in other species, and its evolution in higher dipterans, including Drosophila, is unresolved. We found that cucoid is a member of the ZAD-containing C2H2 zinc finger (ZAD-ZNF) gene family and is orthologous to 27 of the 91 members of this family in D. melanogaster, including M1BP, ranshi, ouib, nom, zaf1, odj, Nnk, trem, Zif, and eighteen uncharacterized genes. Available knowledge of the functions of cucoid orthologs in Drosophila melanogaster suggest that the progenitor of this lineage specific expansion may have played a role in regulating chromatin. We also describe many aspects of the gene duplication history of cucoid in the brachyceran lineage of D. melanogaster, thereby providing a framework for predicting potential redundancies among these genes in D. melanogaster.

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

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

Early effects of gene duplication on the robustness and phenotypic variability of gene regulatory networks

Background

Research on gene duplication is abundant and comes from a wide range of approaches, from high-throughput analyses and experimental evolution to bioinformatics and theoretical models. Notwithstanding, a consensus is still lacking regarding evolutionary mechanisms involved in evolution through gene duplication as well as the conditions that affect them. We argue that a better understanding of evolution through gene duplication requires considering explicitly that genes do not act in isolation. It demands studying how the perturbation that gene duplication implies percolates through the web of gene interactions. Due to evolution’s contingent nature, the paths that lead to the final fate of duplicates must depend strongly on the early stages of gene duplication, before gene copies have accumulated distinctive changes.

Methods

Here we use a widely-known model of gene regulatory networks to study how gene duplication affects network behavior in early stages. Such networks comprise sets of genes that cross-regulate. They organize gene activity creating the gene expression patterns that give cells their phenotypic properties. We focus on how duplication affects two evolutionarily relevant properties of gene regulatory networks: mitigation of the effect of new mutations and access to new phenotypic variants through mutation.

Results

Among other observations, we find that those networks that are better at maintaining the original phenotype after duplication are usually also better at buffering the effect of single interaction mutations and that duplication tends to enhance further this ability. Moreover, the effect of mutations after duplication depends on both the kind of mutation and genes involved in it. We also found that those phenotypes that had easier access through mutation before duplication had higher chances of remaining accessible through new mutations after duplication.

Conclusion

Our results support that gene duplication often mitigates the impact of new mutations and that this effect is not merely due to changes in the number of genes. The work that we put forward helps to identify conditions under which gene duplication may enhance evolvability and robustness to mutations.

A draft genome of the medicinal plant Cremastra appendiculata (D. Don) provides insights into the colchicine biosynthetic pathway

Cremastra appendiculata (D. Don) Makino is a rare terrestrial orchid with a high market value as an ornamental and Chinese traditional medicinal herb with a wide range of pharmacological properties. The pseudobulbs of C. appendiculata are one of the primary sources of the famous traditional Chinese medicine "Shancigu", which has been clinically used for treating many diseases, especially, as the main component to treat gout. The lack of genetic research and genome data restricts the modern development and clinical use of C. appendiculata. Here, we report a 2.3 Gb chromosome-level genome of C. appendiculata. We identify a series of candidates of 35 candidate genes responsible for colchicine biosynthesis, among which O-methyltransferase (OMT) gene exhibits an important role in colchicine biosynthesis. Co-expression analysis reveal purple and green-yellow module have close relationships with pseudobulb parts and comprise most of the colchicine pathway genes. Overall, our genome data and the candidate genes reported here set the foundation to decipher the colchicine biosynthesis pathways in medicinal plants.

Gene Recruitments and Dismissals in the Argonaut Genome Provide Insights into Pelagic Lifestyle Adaptation and Shell-like Eggcase Reacquisition

The paper nautilus or greater argonaut, Argonauta argo, is a species of octopods which is characterized by its pelagic lifestyle and by the presence of a protective spiral-shaped shell-like eggcase in females. To reveal the genomic background of how the species adapted to the pelagic lifestyle and acquired its shell-like eggcase, we sequenced the draft genome of the species. The genome size was 1.1 Gb, which is the smallest among the cephalopods known to date, with the top 215 scaffolds (average length 5,064,479 bp) covering 81% (1.09 Gb) of the total assembly. A total of 26,433 protein-coding genes were predicted from 16,802 assembled scaffolds. From these, we identified nearly intact HOX, Parahox, Wnt clusters, and some gene clusters that could probably be related to the pelagic lifestyle, such as reflectin, tyrosinase, and opsin. The gene models also revealed several homologous genes related to calcified shell formation in Conchiferan mollusks, such as Pif-like, SOD, and TRX. Interestingly, comparative genomics analysis revealed that the homologous genes for such genes were also found in the genome of the shell-less octopus, as well as Nautilus, which has a true outer shell. Therefore, the draft genome sequence of Arg. argo presented here has helped us to gain further insights into the genetic background of the dynamic recruitment and dismissal of genes to form an important, converging extended phenotypic structure such as the shell and the shell-like eggcase. Additionally, it allows us to explore the evolution of from benthic to pelagic lifestyles in cephalopods and octopods.

Functional compensation of mouse duplicates by their paralogs expressed in the same tissues

Analyses in a number of organisms have shown that duplicated genes are less likely to be essential than singletons. This implies that genes can often compensate for the loss of their paralogs. However, it is unclear why the loss of some duplicates can be compensated by their paralogs, whereas the loss of other duplicates cannot. Surprisingly, initial analyses in mice did not detect differences in the essentiality of duplicates and singletons. Only subsequent analyses, using larger gene knockout datasets and controlling for a number of confounding factors, did detect significant differences. Previous studies have not taken into account the tissues in which duplicates are expressed. We hypothesized that in complex organisms, in order for a gene's loss to be compensated by one or more of its paralogs, such paralogs need to be expressed in at least the same set of tissues as the lost gene. To test our hypothesis, we classified mouse duplicates into two categories based on the expression patterns of their paralogs: "compensable duplicates" (those with paralogs expressed in all the tissues in which the gene is expressed) and "non-compensable duplicates" (those whose paralogs are not expressed in all the tissues where the gene is expressed). In agreement with our hypothesis, the essentiality of non-compensable duplicates is similar to that of singletons, whereas compensable duplicates exhibit a substantially lower essentiality. Our results imply that duplicates can often compensate for the loss of their paralogs, but only if they are expressed in the same tissues. Indeed, the compensation ability is more dependent on expression patterns than on protein sequence similarity. The existence of these two kinds of duplicates with different essentialities, which has been overlooked by prior studies, may have hindered the detection of differences between singletons and duplicates.

A Simple Evolutionary Model of Genetic Robustness After Gene Duplication

When a dispensable gene is duplicated (referred to the ancestral dispensability denoted by O⁺), genetic buffering and duplicate compensation together maintain the duplicate redundancy, whereas duplicate compensation is the only mechanism when an essential gene is duplicated (referred to the ancestral essentiality denoted by O^-). To investigate these evolutionary scenarios of genetic robustness, I formulated a simple mixture model for analyzing duplicate pairs with one of the following states: double dispensable (DD), semi-dispensable (one dispensable one essential, DE), or double essential (EE). This model was applied to the yeast duplicate pairs from a whole-genome duplication (WGD) occurred about 100 million years ago (mya), and the mouse duplicate pairs from a WGD occurred about more than 500 mya. Both case studies revealed that the proportion of essentiality for those duplicates with ancestral essentiality [P_E(O^-)] was much higher than that for those with ancestral dispensability [P_E(O⁺)]. While it was negligible in the yeast duplicate pairs, P_E(O⁺) (about 20%) was shown statistically significant in the mouse duplicate pairs. These findings, together, support the hypothesis that both sub-functionalization and neo-functionalization may play some roles after gene duplication, though the former may be much faster than the later.

Mammalian cytochrome P450 biodiversity: Physiological importance, function, and protein and genomic structures of cytochromes P4502B in multiple species of woodrats with different dietary preferences

The vast diversity of cytochrome P450 enzymes in mammals has been proposed to result in large measure from plant-animal warfare, whereby evolution of chemical defenses such as phenolics and terpenoids in plants led to duplication and divergence of P450 genes in herbivores. Over evolutionary time, natural selection is predicted to have produced P450s with high affinity and enhanced metabolism of substrates that are ingested regularly by herbivores. Interestingly, however, almost all knowledge of the interactions of mammalian P450 enzymes with substrates stems from studies of the metabolism of drugs and model compounds rather than studies on wild mammalian herbivores and their respective PSMs. A question of particular interest centers on the role of individual P450 enzymes in the ability of certain herbivores to specialize on plants that are lethal to most other species, including those from the same genus as the specialists. We tackled this intricate problem using a tractable natural system (herbivorous woodrats, genus Neotoma) focusing on comparisons of the specialist N. stephensi, the facultative specialist N. lepida, and the generalist N. albigula, and employing a cross-disciplinary approach involving ecology, biochemistry, pharmacology, structural biology, and genomics. Based on multiple findings suggesting the importance of CYP2B enzymes for ingestion of juniper and a major constituent, α-pinene, we characterized the structure, function and activity of several CYP2B enzymes in woodrats with different dietary habits. Results to date suggest that differences in CYP2B gene copy number may contribute to differential tolerance of PSMs among woodrat species, although additional work is warranted to firmly link gene copy number to juniper tolerance.

Adaptative Mechanisms of Halophytic Eutrema salsugineum Encountering Saline Environment

Salt cress (Eutrema salsugineum), an Arabidopsis-related halophyte, can naturally adapt to various harsh climates and soil conditions; thus, it is considered a desirable model plant for deciphering mechanisms of salt and other abiotic stresses. Accumulating evidence has revealed that compared with Arabidopsis, salt cress possesses stomata that close more tightly and more succulent leaves during extreme salt stress, a noticeably higher level of proline, inositols, sugars, and organic acids, as well as stress-associated transcripts in unstressed plants, and they are induced rapidly under stress. In this review, we systematically summarize the research on the morphology, physiology, genome, gene expression and regulation, and protein and metabolite profile of salt cress under salt stress. We emphasize the latest advances in research on the genome adaptive evolution encountering saline environments, and epigenetic regulation, and discuss the mechanisms underlying salt tolerance in salt cress. Finally, we discuss the existing questions and opportunities for future research in halophytic Eutrema. Together, the review fosters a better understanding of the mechanism of plant salt tolerance and provides a reference for the research and utilization of Eutrema as a model extremophile in the future. Furthermore, the prospects for salt cress applied to explore the mechanism of salt tolerance provide a theoretical basis to develop new strategies for agricultural biotechnology.

A Population-Genetic Lens into the Process of Gene Loss Following Whole-Genome Duplication

Whole-genome duplications (WGDs) have occurred in many eukaryotic lineages. However, the underlying evolutionary forces and molecular mechanisms responsible for the long-term retention of gene duplicates created by WGDs are not well understood. We employ a population-genomic approach to understand the selective forces acting on paralogs and investigate ongoing duplicate-gene loss in multiple species of Paramecium that share an ancient WGD. We show that mutations that abolish protein function are more likely to be segregating in retained WGD paralogs than in single-copy genes, most likely because of ongoing nonfunctionalization post-WGD. This relaxation of purifying selection occurs in only one WGD paralog, accompanied by the gradual fixation of nonsynonymous mutations and reduction in levels of expression, and occurs over a long period of evolutionary time, "marking" one locus for future loss. Concordantly, the fitness effects of new nonsynonymous mutations and frameshift-causing indels are significantly more deleterious in the highly expressed copy compared with their paralogs with lower expression. Our results provide a novel mechanistic model of gene duplicate loss following WGDs, wherein selection acts on the sum of functional activity of both duplicate genes, allowing the two to wander in expression and functional space, until one duplicate locus eventually degenerates enough in functional efficiency or expression that its contribution to total activity is too insignificant to be retained by purifying selection. Retention of duplicates by such mechanisms predicts long times to duplicate-gene loss, which should not be falsely attributed to retention due to gain/change in function.

Dominance shifts increase the likelihood of soft selective sweeps

Genetic models of adaptation to a new environment have typically assumed that the alleles involved maintain a constant fitness dominance across the old and new environments. However, theories of dominance suggest that this should often not be the case. Instead, the alleles involved should frequently shift from recessive deleterious in the old environment to dominant beneficial in the new environment. Here, we study the consequences of these expected dominance shifts for the genetics of adaptation to a new environment. We find that dominance shifts increase the likelihood that adaptation occurs from standing variation, and that multiple alleles from the standing variation are involved (a soft selective sweep). Furthermore, we find that expected dominance shifts increase the haplotypic diversity of selective sweeps, rendering soft sweeps more detectable in small genomic samples. In cases where an environmental change threatens the viability of the population, we show that expected dominance shifts of newly beneficial alleles increase the likelihood of evolutionary rescue and the number of alleles involved. Finally, we apply our results to a well-studied case of adaptation to a new environment: the evolution of pesticide resistance at the Ace locus in Drosophila melanogaster. We show that, under reasonable demographic assumptions, the expected dominance shift of resistant alleles causes soft sweeps to be the most frequent outcome in this case, with the primary source of these soft sweeps being the standing variation at the onset of pesticide use, rather than recurrent mutation thereafter.

Toxicological Implications of Biological Heterogeneity

From a micro to macro scale of biological organization, macromolecular diversity and biological heterogeneity are fundamental properties of biological systems. Heterogeneity may result from genetic, epigenetic, and non-genetic characteristics (e.g., tissue microenvironment). Macromolecular diversity and biological heterogeneity are tolerated as long as the sustenance and propagation of life are not disrupted. They also provide the raw materials for microevolutionary changes that may help organisms adapt to new selection pressures arising from the environment. Sequence evolution, functional divergence, and positive selection of gene and promoter dosage play a major role in the evolution of life's diversity including complex metabolic networks, which is ultimately reflected in changes in the allele frequency over time. Robustness in evolvable biological systems is conferred by functional redundancy that is often created by macromolecular diversity and biological heterogeneity. The ability to investigate biological macromolecules at an increasingly finer level has uncovered a wealth of information in this regard. Therefore, the dynamics of biological complexity should be taken into consideration in biomedical research.

Population-genomic analysis identifies a low rate of global adaptive fixation in the proteins of the cyclical parthenogen Daphnia magna

Daphnia are well-established ecological and evolutionary models, and the interaction between D. magna and its microparasites is widely considered a paragon of the host-parasite coevolutionary process. Like other well-studied arthropods such as Drosophila melanogaster and Anopheles gambiae, D. magna is a small, widespread, and abundant species that is therefore expected to display a large long-term population size and high rates of adaptive protein evolution. However, unlike these other species, D. magna is cyclically asexual and lives in a highly structured environment (ponds and lakes) with moderate levels of dispersal, both of which are predicted to impact upon long-term effective population size and adaptive protein evolution. To investigate patterns of adaptive protein fixation, we produced the complete coding genomes of 36 D. magna clones sampled from across the European range (Western Palaearctic), along with draft sequences for the close relatives D. similis and D. lumholtzi, used as outgroups. We analyzed genome-wide patterns of adaptive fixation, with a particular focus on genes that have an a priori expectation of high rates, such as those likely to mediate immune responses, RNA interference against viruses and transposable elements, and those with a strongly male-biased expression pattern. We find that, as expected, D. magna displays high levels of diversity and that this is highly structured among populations. However, compared with Drosophila, we find that D. magna proteins appear to have a high proportion of weakly deleterious variants and do not show evidence of pervasive adaptive fixation across its entire range. This is true of the genome as a whole, and also of putative ‘arms race’ genes that often show elevated levels of adaptive substitution in other species. In addition to the likely impact of extensive, and previously documented, local adaptation, we speculate that these findings may reflect reduced efficacy of selection associated with cyclical asexual reproduction.

Gene duplication and rate variation in the evolution of plastid ACCase and Clp genes in angiosperms

While the chloroplast (plastid) is known for its role in photosynthesis, it is also involved in many other metabolic pathways essential for plant survival. As such, plastids contain an extensive suite of enzymes required for non-photosynthetic processes. The evolution of the associated genes has been especially dynamic in flowering plants (angiosperms), including examples of gene duplication and extensive rate variation. We examined the role of ongoing gene duplication in two key plastid enzymes, the acetyl-CoA carboxylase (ACCase) and the caseinolytic protease (Clp), responsible for fatty acid biosynthesis and protein turnover, respectively. In plants, there are two ACCase complexes-a homomeric version present in the cytosol and a heteromeric version present in the plastid. Duplications of the nuclear-encoded homomeric ACCase gene and retargeting of one resultant protein to the plastid have been previously reported in multiple species. We find that these retargeted homomeric ACCase proteins exhibit elevated rates of sequence evolution, consistent with neofunctionalization and/or relaxation of selection. The plastid Clp complex catalytic core is composed of nine paralogous proteins that arose via ancient gene duplication in the cyanobacterial/plastid lineage. We show that further gene duplication occurred more recently in the nuclear-encoded core subunits of this complex, yielding additional paralogs in many species of angiosperms. Moreover, in six of eight cases, subunits that have undergone recent duplication display increased rates of sequence evolution relative to those that have remained single copy. We also compared substitution patterns between pairs of Clp core paralogs to gain insight into post-duplication evolutionary routes. These results show that gene duplication and rate variation continue to shape the plastid proteome.

Species Tree Inference Methods Intended to Deal with Incomplete Lineage Sorting Are Robust to the Presence of Paralogs

Many recent phylogenetic methods have focused on accurately inferring species trees when there is gene tree discordance due to incomplete lineage sorting (ILS). For almost all of these methods, and for phylogenetic methods in general, the data for each locus are assumed to consist of orthologous, single-copy sequences. Loci that are present in more than a single copy in any of the studied genomes are excluded from the data. These steps greatly reduce the number of loci available for analysis. The question we seek to answer in this study is: what happens if one runs such species tree inference methods on data where paralogy is present, in addition to or without ILS being present? Through simulation studies and analyses of two large biological data sets, we show that running such methods on data with paralogs can still provide accurate results. We use multiple different methods, some of which are based directly on the multispecies coalescent model, and some of which have been proven to be statistically consistent under it. We also treat the paralogous loci in multiple ways: from explicitly denoting them as paralogs, to randomly selecting one copy per species. In all cases, the inferred species trees are as accurate as equivalent analyses using single-copy orthologs. Our results have significant implications for the use of ILS-aware phylogenomic analyses, demonstrating that they do not have to be restricted to single-copy loci. This will greatly increase the amount of data that can be used for phylogenetic inference.[Gene duplication and loss; incomplete lineage sorting; multispecies coalescent; orthology; paralogy.].

Retrogene survival is not impacted by linkage relationships

In retrogene evolution, the out-of-the-X pattern is the retroduplication of X-linked housekeeping genes to autosomes, hypothesized to be driven by meiotic sex chromosome inactivation during spermatogenesis. This pattern suggests that some retrogene survival is driven by selection on X-linkage. We asked if selection on linkage constitutes an important evolutionary force in retrogene survival, including for autosomal parents. Specifically, is there a correlation between retrogene survival and changes in linkage with parental gene networks? To answer this question, we compiled data on retrogenes in both Homo sapiens and Drosophila melanogaster and using Monte Carlo methods, we tested whether retrogenes exhibit significantly different linkage relationships than expected under a null assumption of uniform distribution in the genome. Overall, after excluding genes involved in the out-of-the-X pattern, no general pattern was found associating genetic linkage and retrogene survival. This demonstrates that selection on linkage may not represent an overarching force in retrogene survival. However, it remains possible that this type of selection still influences the survival of specific retrogenes.

Duplication and specialization of NUDX1 in Rosaceae led to geraniol production in rose petals

Nudix hydrolases are conserved enzymes ubiquitously present in all kingdoms of life. Recent research revealed that several Nudix hydrolases are involved in terpenoid metabolism in plants. In modern roses, RhNUDX1 is responsible for formation of geraniol, a major compound of rose scent. Nevertheless, this compound is produced by monoterpene synthases in many geraniol-producing plants. As a consequence, this raised the question about the origin of RhNUDX1 function and the NUDX1 gene evolution in Rosaceae, in wild roses or/and during the domestication process. Here, we showed that three distinct clades of NUDX1 emerged in the Rosoidae subfamily (Nudx1-1 to Nudx1-3 clades), and two subclades evolved in the Rosa genus (Nudx1-1a and Nudx1-1b subclades). We also showed that the Nudx1-1b subclade was more ancient than the Nudx1-1a subclade, and that the NUDX1-1a gene emerged by a trans-duplication of the more ancient NUDX1-1b gene. After the transposition, NUDX1-1a was cis-duplicated, leading to a gene dosage effect on the production of geraniol in different species. Furthermore, the NUDX1-1a appearance was accompanied by the evolution of its promoter, most likely from a Copia retrotransposon origin, leading to its petal-specific expression. Thus, our data strongly suggest that the unique function of NUDX1-1a in geraniol formation was evolved naturally in the genus Rosa before domestication.

The relative impact of evolving pleiotropy and mutational correlation on trait divergence

Both pleiotropic connectivity and mutational correlations can restrict the decoupling of traits under divergent selection, but it is unknown which is more important in trait evolution. To address this question, we create a model that permits within-population variation in both pleiotropic connectivity and mutational correlation, and compare their relative importance to trait evolution. Specifically, we developed an individual-based stochastic model where mutations can affect whether a locus affects a trait and the extent of mutational correlations in a population. We find that traits can decouple whether there is evolution in pleiotropic connectivity or mutational correlation, but when both can evolve, then evolution in pleiotropic connectivity is more likely to allow for decoupling to occur. The most common genotype found in this case is characterized by having one locus that maintains connectivity to all traits and another that loses connectivity to the traits under stabilizing selection (subfunctionalization). This genotype is favored because it allows the subfunctionalized locus to accumulate greater effect size alleles, contributing to increasingly divergent trait values in the traits under divergent selection without changing the trait values of the other traits (genetic modularization). These results provide evidence that partial subfunctionalization of pleiotropic loci may be a common mechanism of trait decoupling under regimes of corridor selection.

Dosage sensitivity and exon shuffling shape the landscape of polymorphic duplicates in Drosophila and humans

Despite polymorphic duplicate genes' importance for the early stages of duplicate gene evolution, they are less studied than old gene duplicates. Two essential questions thus remain poorly addressed: how does dosage sensitivity, imposed by stoichiometry in protein complexes or by X chromosome dosage compensation, affect the emergence of complete duplicate genes? Do introns facilitate intergenic and intragenic chimaerism as predicted by the theory of exon shuffling? Here, we analysed new data for Drosophila and public data for humans, to characterize polymorphic duplicate genes with respect to dosage, exon-intron structures and allele frequencies. We found that complete duplicate genes are under dosage constraint induced by protein stoichiometry but potentially tolerated by X chromosome dosage compensation. We also found that in the intron-rich human genome, gene fusions and intragenic duplications extensively use intronic breakpoints generating in-frame proteins, in accordance with the theory of exon shuffling. Finally, we found that only a small proportion of complete or partial duplicates are at high frequencies, indicating the deleterious nature of dosage or gene structural changes. Altogether, we demonstrate how mechanistic factors including dosage sensitivity and exon-intron structure shape the short-term functional consequences of gene duplication.

Defence by duplication: The relation between phenotypic glyphosate resistance and EPSPS gene copy number variation in Amaranthus palmeri

Gene copy number variation (CNV) has been increasingly associated with organismal responses to environmental stress, but we know little about the quantitative relation between CNV and phenotypic variation. In this study we quantify the relation between variation in EPSPS (5-enolpyruvylshikimate-3-phosphate synthase) copy number using digital drop PCR and variation in phenotypic glyphosate resistance in 22 populations of Amaranthus palmeri (Palmer Amaranth), a range-expanding agricultural weed. Overall, we detected a significant positive relation between population mean copy number and resistance. The majority of populations exhibited high glyphosate resistance yet maintained low-resistance individuals, resulting in bimodality in many populations. We also investigated threshold models for the relation between copy number and resistance, and found evidence for a threshold of ~15 EPSPS copies: there was a steep increase in resistance below the threshold, followed by a much shallower increase. Across 924 individuals, as copy number increased the range of variation in resistance decreased, yielding an increasing frequency of high phenotypic resistance individuals. Among populations we detected a decline in variation (s.d.) as mean phenotypic resistance increased from moderate to high, consistent with the prediction that as phenotypic resistance increases in populations, stabilizing selection decreases variation in the trait. Our study demonstrates that populations of A. palmeri can harbour wide variation in EPSPS copy number and phenotypic glyphosate resistance, reflecting the history of, and template for future, resistance evolution.

Diversity in olfactory receptor repertoires is associated with dietary specialization in a genus of frugivorous bat

Mammalian olfactory receptor genes (ORs) are a diverse family of genes encoding proteins that directly interact with environmental chemical cues. ORs evolve via gene duplication in a birth-death fashion, neofunctionalizing and pseudogenizing over time. Olfaction is a primary sense used for food detection in plant-visiting bats, but the relationship between dietary specialization and OR repertoire diversity is unclear. Within neotropical Leaf-nosed bats (Phyllostomidae), many lineages are plant specialists, and some have a distinct OR repertoire compared to insectivorous species. Yet, whether specialization on particular plant genera is associated with the evolution of specialized, less diverse OR repertoires has never been tested. Using targeted sequence capture, we sequenced the OR repertoires of three sympatric species of short-tailed fruit bats (Carollia), which vary in their degree of specialization on the fruits of Piper plants. We characterized orthologous vs duplicated receptors among Carollia species, and explored the diversity and redundancy of the receptor gene repertoire. At the species level, the most dedicated Piper specialist, Carollia castanea, had lower OR diversity compared to the two generalists (C. sowelli and C. perspicillata), but we discovered a few unique sets of ORs within C. castanea with high redundancy of similar gene duplicates. These unique receptors potentially enable C. castanea to detect Piper fruit odorants better than its two congeners. Carollia perspicillata, the species with the most generalist diet, had a higher diversity of intact receptors, suggesting the ability to detect a wider range of odorant molecules. Variation among ORs may be a factor in the coexistence of these sympatric species, facilitating the exploitation of different plant resources. Our study sheds light on how gene duplication and changes in OR diversity may play a role in dietary adaptations and underlie ecological interactions between bats and plants.

Segregational drift constrains the evolutionary rate of prokaryotic plasmids

Plasmids are extrachromosomal genetic elements in prokaryotes that have been recognized as important drivers of microbial ecology and evolution. Plasmids are found in multiple copies inside their host cell where independent emergence of mutations may lead to intracellular genetic heterogeneity. The intracellular plasmid diversity is thus subject to changes upon cell division. However, the effect of plasmid segregation on plasmid evolution remains understudied. Here, we show that genetic drift during cell division—segregational drift—leads to the rapid extinction of novel plasmid alleles. We established a novel experimental approach to control plasmid allele frequency at the levels of a single cell and the whole population. Following the dynamics of plasmid alleles in an evolution experiment, we find that the mode of plasmid inheritance—random or clustered—is an important determinant of plasmid allele dynamics. Phylogenetic reconstruction of our model plasmid in clinical isolates furthermore reveals a slow evolutionary rate of plasmid-encoded genes in comparison to chromosomal genes. Our study provides empirical evidence that genetic drift in plasmid evolution occurs at multiple levels: the host cell and the population of hosts. Segregational drift has implications for the evolutionary rate heterogeneity of extrachromosomal genetic elements.

Complete genome sequence and identification of polyunsaturated fatty acid biosynthesis genes of the myxobacterium Minicystis rosea DSM 24000T

Background

Myxobacteria harbor numerous biosynthetic gene clusters that can produce a diverse range of secondary metabolites. Minicystis rosea DSM 24000^T is a soil-dwelling myxobacterium belonging to the suborderSorangiineae and family Polyangiaceae and is known to produce various secondary metabolites as well as polyunsaturated fatty acids (PUFAs). Here, we use whole-genome sequencing to explore the diversity of biosynthetic gene clusters in M. rosea.

Results

Using PacBio sequencing technology, we assembled the 16.04 Mbp complete genome of M. rosea DSM 24000^T, the largest bacterial genome sequenced to date. About 44% of its coding potential represents paralogous genes predominantly associated with signal transduction, transcriptional regulation, and protein folding. These genes are involved in various essential functions such as cellular organization, diverse niche adaptation, and bacterial cooperation, and enable social behavior like gliding motility, sporulation, and predation, typical of myxobacteria. A profusion of eukaryotic-like kinases (353) and an elevated ratio of phosphatases (8.2/1) in M. rosea as compared to other myxobacteria suggest gene duplication as one of the primary modes of genome expansion. About 7.7% of the genes are involved in the biosynthesis of a diverse array of secondary metabolites such as polyketides, terpenes, and bacteriocins. Phylogeny of the genes involved in PUFA biosynthesis (pfa) together with the conserved synteny of the complete pfa gene cluster suggests acquisition via horizontal gene transfer from Actinobacteria.

Conclusion

Overall, this study describes the complete genome sequence of M. rosea, comparative genomic analysis to explore the putative reasons for its large genome size, and explores the secondary metabolite potential, including the biosynthesis of polyunsaturated fatty acids.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07955-x.

Functional divergence of Brassica napus BnaABI1 paralogs in the structurally conserved PP2CA gene subfamily of Brassicaceae

Group A PP2C (PP2CA) genes form a gene subfamily whose members play an important role in regulating many biological processes by dephosphorylation of target proteins. In this study we examined the effects of evolutionary changes responsible for functional divergence of BnaABI1 paralogs in Brassica napus against the background of the conserved PP2CA gene subfamily in Brassicaceae. We performed comprehensive phylogenetic analyses of 192 PP2CA genes in 15 species in combination with protein structure homology modeling. Fundamentally, the number of PP2CA genes remained relatively constant in these taxa, except in the Brassica genus and Camelina sativa. The expansion of this gene subfamily in these species has resulted from whole genome duplication. We demonstrated a high degree of structural conservation of the PP2CA genes, with a few minor variations between the different PP2CA groups. Furthermore, the pattern of conserved sequence motifs in the PP2CA proteins and their secondary and 3D structures revealed strong conservation of the key ion-binding sites. Syntenic analysis of triplicated regions including ABI1 paralogs revealed significant structural rearrangements of the Brassica genomes. The functional and syntenic data clearly show that triplication of BnaABI1 in B. napus has had an impact on its functions, as well as the positions of adjacent genes in the corresponding chromosomal regions. The expression profiling of BnaABI1 genes showed functional divergence, i.e. subfunctionalization, potentially leading to neofunctionalization. These differences in expression are likely due to changes in the promoters of the BnaABI1 paralogs. Our results highlight the complexity of PP2CA gene subfamily evolution in Brassicaceae.

Functional Diversification, Redundancy, and Epistasis among Paralogs of the Drosophila melanogaster Obp50a-d Gene Cluster

Large multigene families, such as the insect odorant-binding proteins (OBPs), are thought to arise through functional diversification after repeated gene duplications. Whereas many OBPs function in chemoreception, members of this family are also expressed in tissues outside chemosensory organs. Paralogs of the Obp50 gene cluster are expressed in metabolic and male reproductive tissues, but their functions and interrelationships remain unknown. Here, we report the genetic dissection of four members of the Obp50 cluster, which are in close physical proximity without intervening genes. We used CRISPR technology to excise the entire cluster while introducing a PhiC31 reintegration site to reinsert constructs in which different combinations of the constituent Obp genes were either intact or rendered inactive. We performed whole transcriptome sequencing and assessed sexually dimorphic changes in transcript abundances (transcriptional niches) associated with each gene-edited genotype. Using this approach, we were able to estimate redundancy, additivity, diversification, and epistasis among Obp50 paralogs. We analyzed the effects of gene editing of this cluster on organismal phenotypes and found a significant skewing of sex ratios attributable to Obp50a, and sex-specific effects on starvation stress resistance attributable to Obp50d. Thus, there is functional diversification within the Obp50 cluster with Obp50a contributing to development and Obp50d to stress resistance. The deletion-reinsertion approach we applied to the Obp50 cluster provides a general paradigm for the genetic dissection of paralogs of multigene families.

Characterization of the turbot Scophthalmus maximus (L.) myeloperoxidase. An insight into the evolution of vertebrate peroxidases

We have completed the characterization of the turbot (Scophthalmus maximus) myeloperoxidase (mpx) gene and protein, which we partially described in a previous study. The turbot mpx gene has 15 exons that encode a protein of 767 aa, with a signal peptide, propeptide and light and heavy chains, and also with haem cavities, a Ca⁺²-binding motif and several N- and O-glycosylation sites. The mature protein forms homodimers of about 150 kDa and is very abundant in turbot neutrophils. In addition to the mpx (epx2a) gene, another three peroxidase genes, named epx1, epx2b1 and epx2b2, were identified in the turbot genome. Epx1, Epx2b1 and Epx2b2 proteins also have signal peptides and many structural characteristics of mammalian MPO and eosinophil peroxidase (EPX). Mpx was strongly expressed in head kidney, while epx2b1 and epx2b2 were strongly expressed in the gills, and epx1 was not expressed in any of the tissues or organs analysed. In vitro stimulation of head kidney leucocytes with the parasite Philasterides dicentrarchi caused a decrease in mpx expression and an increase in epx2b1 expression over time. In turbot infected experimentally with P. dicentrarchi a significant increase in mpx expression in the head kidney was observed on day 7 postinfection, while the other genes were not regulated. However, mpx, epx2b1 and epx2b2 were downregulated in the gills of infected fish, and epx1 expression was not affected. These results suggest that the four genes responded differently to the same stimuli. Interestingly, BLAST analysis revealed that Epx1 and Mpx showed greater similarity to mammalian EPX than to MPO. Considering the phylogenetic and synteny data obtained, we concluded that the epx/mpx genes of Gnathostomes can be divided into three main clades: EPX1, which contains turbot epx1, EPX2, which contains turbot mpx (epx2a) and epx2b1 and epx2b2 genes, and a clade containing mammalian EPX and MPO (EPX/MPO). EPX/MPO and EPX2 clades share a common ancestor with the chondrichthyan elephant shark (Callorhinchus milii) and the coelacanth (Latimeria chalumnae) peroxidases. EPX2 was only found in fish and includes two sister groups. One of the groups includes turbot mpx and was only found in teleosts. Finally, the other group contains epx2b1 and epx2b2 genes, and epx2b1-2b2 loci share orthologous genes with other teleosts and also with holosteans, suggesting that these genes appeared earlier on than the mpx gene.

Kingdom-wide analysis of the evolution of the plant type III polyketide synthase superfamily

The emergence of type III polyketide synthases (PKSs) was a prerequisite for the conquest of land by the green lineage. Within the PKS superfamily, chalcone synthases (CHSs) provide the entry point reaction to the flavonoid pathway, while LESS ADHESIVE POLLEN 5 and 6 (LAP5/6) provide constituents of the outer exine pollen wall. To study the deep evolutionary history of this key family, we conducted phylogenomic synteny network and phylogenetic analyses of whole-genome data from 126 species spanning the green lineage including Arabidopsis thaliana, tomato (Solanum lycopersicum), and maize (Zea mays). This study thereby combined study of genomic location and context with changes in gene sequences. We found that the two major clades, CHS and LAP5/6 homologs, evolved early by a segmental duplication event prior to the divergence of Bryophytes and Tracheophytes. We propose that the macroevolution of the type III PKS superfamily is governed by whole-genome duplications and triplications. The combined phylogenetic and synteny analyses in this study provide insights into changes in the genomic location and context that are retained for a longer time scale with more recent functional divergence captured by gene sequence alterations.

Copy Number Variations of Glycoside Hydrolase 45 Genes in Bursaphelenchus xylophilus and Their Impact on the Pathogenesis of Pine Wilt Disease

The pine wood nematode Bursaphelenchus xylophilus parasitizes millions of pine trees worldwide each year, causing severe wilt and the death of host trees. Glycoside hydrolase 45 genes of B. xylophilus are reported to have been acquired by horizontal gene transfer from fungi and are responsible for cell wall degradation during nematode infection. Previous studies ignored the possibility of copy number variations of such genes. In this study, we determined that two of the glycoside hydrolase 45 genes evolved to maintain multiple copies with distinct expression levels, enabling the nematode to infect a variety of pine hosts. Additionally, tandem repeat variations within coding regions were also detected between different copies of glycoside hydrolase 45 genes that could result in changes in protein sequences and serve as an effective biological marker to detect copy number variations among different B. xylophilus populations. Consequently, we were able to further identify the copy number variations of glycoside hydrolase 45 genes among B. xylophilus strains with different virulence. Our results provide new insights into the pathogenicity of B. xylophilus, provide a practical marker to genotype copy number variations and may aid in population classification.

Copy number variations and young duplicate genes have high methylation levels in sticklebacks

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.

Gene family amplification facilitates adaptation in freshwater unionid bivalve Megalonaias nervosa

Freshwater unionid bivalves currently face severe anthropogenic challenges. Over 70% of species in the United States are threatened, endangered or extinct due to pollution, damming of waterways and overfishing. These species are notable for their unusual life history strategy, parasite-host co-evolution and biparental mitochondrial inheritance. Among this clade, the washboard mussel Megalonaias nervosa is one species that remains prevalent across the Southeastern United States, with robust population sizes. We have created a reference genome for M. nervosa to determine how genome content has evolved in the face of these widespread environmental challenges. We observe dynamic changes in genome content, with a burst of recent transposable element proliferation causing a 382 Mb expansion in genome content. Birth-death models suggest rapid expansions among gene families, with a mutation rate of 1.16 × 10^-8 duplications per gene per generation. Cytochrome P450 gene families have experienced exceptional recent amplification beyond expectations based on genome-wide birth-death processes. These genes are associated with increased rates of amino acid changes, a signature of selection driving evolution of detox genes. Fitting evolutionary models of adaptation from standing genetic variation, we can compare adaptive potential across species and mutation types. The large population size in M. nervosa suggests a 4.7-fold advantage in the ability to adapt from standing genetic variation compared with a low diversity endemic E. hopetonensis. Estimates suggest that gene family evolution may offer an exceptional substrate of genetic variation in M. nervosa, with P_sgv  = 0.185 compared with P_sgv  = 0.067 for single nucleotide changes. Hence, we suggest that gene family evolution is a source of 'hopeful monsters' within the genome that may facilitate adaptation when selective pressures shift. These results suggest that gene family expansion is a key driver of adaptive evolution in this key species of freshwater Unionidae that is currently facing widespread environmental challenges. This work has clear implications for conservation genomics on freshwater bivalves as well as evolutionary theory. This genome represents a first step to facilitate reverse ecological genomics in Unionidae and identify the genetic underpinnings of phenotypic diversity.

Signatures of selection and evolutionary relevance of cytochrome P450s in plant-insect interactions

Enzymes in the cytochrome P450 (P450) superfamily have important functions ranging from those that are essential for the physiology and development of the individual to those that mediate interactions between individuals and their biotic environment. Until recently the study of P450s had focused on single functions, substrates, or pathways. Recent advances in sequencing, genome assembly, and phylogenetic methods have returned emphasis to the adaptive value of these enzymes in the context of herbivory. Comparisons of whole repertoires of P450s across related species reveal that P450s capable of metabolizing xenobiotics have an increased rate of gains compared to losses after gene duplications. In plants, studies have focused on enzymes and end-functions that have converged to provide increased resistance to herbivory. This review summarizes the latest findings related to the ecological value of P450s in the interactions between phytophagous insects and their host plants.

Phylogeny of Slc15 family and response to Aeromonas hydrophila infection following Lactococcus lactis dietary supplementation in Cyprinus carpio

Solute carrier 15 family (Slc15) are membrane proteins that utilize the proton gradient and negative membrane protential for the transmembrane transporter of di-/tripeptide and peptide-mimetic molecules, in addition, they also play important roles in immunoreaction. In this study, 10 Slc15 genes were identified in the common carp genome database. Comparative genomics analysis showed considerable expansion of the Slc15 genes and verified the four-round whole genome duplication (WGD) event in common carp. Phylogenetic analysis revealed all Slc15 genes of common carp were clustered into orthologous groups indicating the highly conservative during evolution. Besides, the tissues and temporal expression examined by RT-PCR and qRT-PCR showed that most of the Slc15 genes had a narrow tissue distribution and exhibited tissue-specific expression patterns. Expression divergences were observed between these copies proving function divergence after the WGD. Then, we investigated the dietary supplementation effects of three Lactococcus lactis strains on the expression of Slc15 genes in common carp infected by A. hydrophila to find an effective way to treat aquatic diseases. Almost all of the Slc15 genes had an increased expression trend in the early post-challenge stage, and reached the highest expression level at 12h post-challenge. Then, the expression level showed a bluff descent at the last two stages and the expression level reached the lowest at 48 h post-challenge. Slc15 genes expression is actively up-regulated when stimulated by inflammatory factors, which can "amplify" immune signals, and improve the body's defense against foreign invasion in the early stage of the inflammatory response. So activation of the Slc15 genes may be an effective way for infectious disease treatment. As expected, three strains improved the expression of Slc15 genes variously compared with the control/infection groups. The strain 3 of L. lactis had a better induction of Slc15 genes compared with strain 1 and strain 2. It might be applied as a potential activation of Slc15 genes for disease treatment and adding befitting L. lactis may be a good way to protect aquatilia from bacillosis.

VLRs expression were significantly affected by complement C3 knockdown morphants in Lampetra morii

The complement component 3 of the lamprey, a jawless vertebrate, functions as an opsonin during the phagocytosis of rabbit red cells. Furthermore, lamprey C3 may be activated and cleaved into C3b, which is attached to the surface of target cells in the cytolytic process. However, the mechanism mediating the biological function of C3 in the lamprey is unknown. To our knowledge, this study is the first to show that variable lymphocyte receptors (VLRs) expression were significantly affected by complement C3 knockdown morphants in Lampetra morii. We identified the C3 gene in the lamprey genome based on its orthologs, conserved synteny, functional domains, phylogenetic tree, and conserved motifs. Additionally, we determined the optimal infection concentration of Aeromonas hydrophila to perform immune stimulation experiments in the lamprey larvae. The quantitative real-time polymerase chain reaction and immunofluorescence analyses revealed that the expression of Lampetra morii C3 (lmC3) was significantly upregulated in the larvae infected with 10⁷ CFU/mL of A. hydrophila. The lmC3 morphants (lmC3 MO) of lamprey larvae were generated by morpholino-mediated knockdown. The lmC3 MO larvae were highly susceptible to A. hydrophila infection, which indicated that lmC3 is critical in lamprey immune response. The expression of a selected panel of orthologous genes was comparatively analyzed in the infected wild type, infected lmC3 MO, infected control MO, uninfected wild type and uninfected lmC3 MO one-month-old ammocoete larvae. The knockdown of lmC3 strongly affected the expression of VLRA⁺/VLRB⁺/VLRC⁺-associated genes, which was also confirmed by immunohistochemical analysis. Thus, VLR expression were significantly affected by complement C3 knockdown morphants in Lampetra morii.

Genome-Wide Analyses Provide Insights into the Scavenging Lifestyle of the Striped Hyena (Hyaena hyaena)

Hyenas (family Hyaenidae) occupy a variety of different niches, of which the striped hyena (Hyaena hyaena) scavenges mainly on the carcasses of animals. We compared its genome with the genomes of nine other mammals, focusing on similarities and differences in chemoreception, detoxification, digestive, and immune systems. The results showed that the striped hyena's immune and digestive system-related gene families have significantly expanded, which was likely to be an adaptive response to its scavenging lifestyle. In addition, 88 and 26 positive selected genes (PSGs) were identified in the immune system and digestive system, respectively, which may be the molecular basis for immune defense system to effectively resist pathogen invasion. Functional enrichment analysis of PSGs revealed that most of them were involved in the immune regulation process. Among them, eight specific missense mutations were found in two PSGs (MHC class II antigen DOA and MHC class II antigen DOB), suggesting important reorganization of the immune system in the striped hyena. Moreover, we identified one cathelicidin gene and four defensin genes in the striped hyenas by genome mining, which have high-efficiency and broad-spectrum antimicrobial activity. Of particular interest, a striped hyena-specific missense mutation was found in the cathelicidin gene. PolyPhen-2 classified the missense mutation as a harmful mutation, which may have aided in immune adaptation to carrion feeding. Our genomic analyses on the striped hyena provided insights into its success in the adaptation to the scavenging lifestyle.

Dynamic and Differential Expression of Duplicated Cxcr4/Cxcl12 Genes Facilitates Antiviral Response in Hexaploid Gibel Carp

Chemokine receptor cxcr4 and its ligand cxcl12 have evolved two paralogs in the teleost lineage. In this study, we have identified four duplicated cxcr4 and cxcl12 genes from hexaploid gibel carp, Carassius gibelio, respectively. Cgcxcr4bs and Cgcxcl12as were dynamically and differentially expressed in immune-related tissues, and significantly up-regulated in head kidney and spleen after crucian carp herpesvirus (CaHV) infection. Blocking Cxcr4/Cxcl12 axis by injecting AMD3100 brought more severe bleeding symptom and lower survival rate in CaHV-infected fish. AMD3100 treatment also suppressed the up-regulation of key antiviral genes in head kidney and spleen, and resulted in more acute replication of CaHV in vivo. Consistently, the similar suppression of up-regulated expression of key antiviral genes were also observed in CAB cells treated by AMD3100 after poly(I:C) stimulation. Finally, MAPK3 and JAK/STAT were identified as the possible pathways that CgCxcr4s and CgCxcl12s participate in to promote the antiviral response in vitro.