Genomic evidence for adaptation by gene duplication

In Silico Analysis of Glutamate Receptors in Capsicum chinense: Structure, Evolution, and Molecular Interactions

Plant glutamate receptors (GLRs) are integral membrane proteins that function as non-selective cation channels, involved in the regulation of developmental events crucial in plants. Knowledge of these proteins is restricted to a few species and their true agonists are still unknown in plants. Using tomato SlGLRs, a search was performed in the pepper database to identify GLR sequences in habanero pepper (Capsicum chinense Jacq.). Structural, phylogenetic, and orthology analysis of the CcGLRs, as well as molecular docking and protein interaction networks, were conducted. Seventeen CcGLRs were identified, which contained the characteristic domains of GLR. The variation of conserved residues in the M2 transmembrane domain between members suggests a difference in ion selectivity and/or conduction. Also, new conserved motifs in the ligand-binding regions are reported. Duplication events seem to drive the expansion of the species, and these were located in the evolution by using orthologs. Molecular docking analysis allowed us to identify differences in the agonist binding pocket between CcGLRs, which suggest the existence of different affinities for amino acids. The possible interaction of some CcGLRs with proteins leads to suggesting specific functions for them within the plant. These results offer important functional clues for CcGLR, probably extrapolated to other Solanaceae.

Aberrant light sensing and motility in the green alga Chlamydomonas priscuii from the ice-covered Antarctic Lake Bonney

The Antarctic green alga Chlamydomonas priscuii is an obligate psychrophile and an emerging model for photosynthetic adaptation to extreme conditions. Endemic to the ice-covered Lake Bonney, this alga thrives at highly unusual light conditions characterized by very low light irradiance (<15 μmol m-2 s-1), a narrow wavelength spectrum enriched in blue light, and an extreme photoperiod. Genome sequencing of C. priscuii exposed an unusually large genome, with hundreds of highly similar gene duplicates and expanded gene families, some of which could be aiding its survival in extreme conditions. In contrast to the described expansion in the genetic repertoire in C. priscuii, here we suggest that the gene family encoding for photoreceptors is reduced when compared to related green algae. This alga also possesses a very small eyespot and exhibits an aberrant phototactic response, compared to the model Chlamydomonas reinhardtii. We also investigated the genome and behavior of the closely related psychrophilic alga Chlamydomonas sp. ICE-MDV, that is found throughout the photic zone of Lake Bonney and is naturally exposed to higher light levels. Our analyses revealed a photoreceptor gene family and a robust phototactic response similar to those in the model Chlamydomonas reinhardtii. These results suggest that the aberrant phototactic response in C. priscuii is a result of life under extreme shading rather than a common feature of all psychrophilic algae. We discuss the implications of these results on the evolution and survival of shade adapted polar algae.

The structured coalescent in the context of gene copy number variation

The Structured Coalescent was introduced to describe the coalescent process in spatially subdivided populations with migration. Here, we re-interpret migration routes of individuals in the original model as "migration routes" of single genes in tandemly arranged gene arrays. A gene copy may change its position within the array via unequal recombination. Hence, in a coalescent framework, two copies sampled from two chromosomes may coalesce only if they are at exactly homologous positions. Otherwise, one or multiple recombination events have to occur before they can coalesce, thereby increasing mean coalescence time and expected genetic diversity among the copies in a gene array. We explicitly calculate the transition probabilities on these routes backward in time. We simulate the structured coalescent with migration and coalescence rates informed by the unequal recombination process of gene copies. With this novel interpretation of population structure models we determine coalescence times and expected genetic diversity in samples of orthologous and paralogous copies from a gene family. As a case study, we discuss the site frequency spectrum of a small gene family in the two scenarios of high and of no gene copy number variation among individuals. These examples underline the significance of our model, since standard test-statistics may lead to misinterpretations when analyzing sequence data of multi-copy genes due to their different expected genetic diversity.

Venomous Noodles: The Evolution of Toxins in Nemertea through Positive Selection and Gene Duplication

Some, probably most and perhaps all, members of the phylum Nemertea are poisonous, documented so far from marine and benthic specimens. Although the toxicity of these animals has been long known, systematic studies on the characterization of toxins, mechanisms of toxicity, and toxin evolution for this group are scarce. Here, we present the first investigation of the molecular evolution of toxins in Nemertea. Using a proteo-transcriptomic approach, we described toxins in the body and poisonous mucus of the pilidiophoran Lineus sanguineus and the hoplonemertean Nemertopsis pamelaroeae. Using these new and publicly available transcriptomes, we investigated the molecular evolution of six selected toxin gene families. In addition, we also characterized in silico the toxin genes found in the interstitial hoplonemertean, Ototyphlonemertes erneba, a meiofaunal taxa. We successfully identified over 200 toxin transcripts in each of these species. Evidence of positive selection and gene duplication was observed in all investigated toxin genes. We hypothesized that the increased rates of gene duplications observed for Pilidiophora could be involved with the expansion of toxin genes. Studies concerning the natural history of Nemertea are still needed to understand the evolution of their toxins. Nevertheless, our results show evolutionary mechanisms similar to other venomous groups.

Gene duplication and functional divergence of new genes contributed to the polar acclimation of Antarctic green algae

Psychrophilic microalgae successfully survive in the extreme and highly variable polar ecosystems, which represent the energy base of most food webs and play a fundamental role in nutrient cycling. The success of microalgae is rooted in their adaptive evolution. Revealing how they have evolved to thrive in extreme polar environments will help us better understand the origin of life in polar ecosystems. We isolated a psychrophilic unicellular green alga, Microglena sp. YARC, from Antarctic sea ice which has a huge genome. Therefore, we predicted that gene replication may play an important role in its polar adaptive evolution. We found that its protein-coding gene number significantly increased and the duplication time was dated between 37 and 48 million years ago, which is consistent with the formation of the circumpolar Southern Ocean. Most duplicated paralogous genes were enriched in pathways related to photosynthesis, DNA repair, and fatty acid metabolism. Moreover, there were a total of 657 Microglena-specific families, including collagen-like proteins. The divergence in the expression patterns of the duplicated and species-specific genes reflects sub- and neo-functionalization during stress acclimation. Overall, key findings from this study provide new information on how gene duplication and their functional novelty contributed to polar algae adaptation to the highly variable polar environmental conditions.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00203-z.

The genome of Anoplarchus purpurescens (Stichaeidae) reflects its carnivorous diet

Digestion is driven by digestive enzymes and digestive enzyme gene copy number can provide insights on the genomic underpinnings of dietary specialization. The "Adaptive Modulation Hypothesis" (AMH) proposes that digestive enzyme activity, which increases with increased gene copy number, should correlate with substrate quantity in the diet. To test the AMH and reveal some of the genetics of herbivory vs carnivory, we sequenced, assembled, and annotated the genome of Anoplarchus purpurescens, a carnivorous prickleback fish in the family Stichaeidae, and compared the gene copy number for key digestive enzymes to that of Cebidichthys violaceus, a herbivorous fish from the same family. A highly contiguous genome assembly of high quality (N50 = 10.6 Mb) was produced for A. purpurescens, using combined long-read and short-read technology, with an estimated 33,842 protein-coding genes. The digestive enzymes that we examined include pancreatic α-amylase, carboxyl ester lipase, alanyl aminopeptidase, trypsin, and chymotrypsin. Anoplarchus purpurescens had fewer copies of pancreatic α-amylase (carbohydrate digestion) than C. violaceus (1 vs. 3 copies). Moreover, A. purpurescens had one fewer copy of carboxyl ester lipase (plant lipid digestion) than C. violaceus (4 vs. 5). We observed an expansion in copy number for several protein digestion genes in A. purpurescens compared to C. violaceus, including trypsin (5 vs. 3) and total aminopeptidases (6 vs. 5). Collectively, these genomic differences coincide with measured digestive enzyme activities (phenotypes) in the two species and they support the AMH. Moreover, this genomic resource is now available to better understand fish biology and dietary specialization.

Chromosome-Level Analysis of the Pelochelys cantorii Genome Provides Insights to Its Immunity, Growth and Longevity

The Asian giant soft-shelled turtle, Pelochelys cantorii (Trionychidae), is one of the largest aquatic turtles in China and was designated as a First-Grade Protected Animal in China in 1989. Previous investigation based on a combination of Illumina short-read, PacBio long-read and Hi-C scaffolding technologies acquired a high-quality chromosome-level genome of Pc. cantorii. In this study, comparative genomic analysis between Pc. cantorii and 16 other vertebrate genomes indicated that turtles separated from the ancestor of archosaurians approximately 256.6 (95% highest posterior density interval, 263.6-251.9) million years ago (Mya) (Upper Permian to Triassic) and that Pc. cantorii separated from the ancestor of Pd. sinensis and R. swinhoei approximately 59.3 (95% highest posterior density interval, 64.3-54.3) Mya. Moreover, several candidate genes, such as VWA5A, ABCG2, A2M and IGSF1, associated with tumor suppression, growth and age were expanded, implicating their potential roles in the exceptional longevity of turtles. This new chromosome-level assembly has important scientific value in the study of conservation of Pc. cantorii and also enriches the evolutionary investigation of turtle species.

Parallel expansion and divergence of an adhesin family in pathogenic yeasts

Opportunistic yeast pathogens arose multiple times in the Saccharomycetes class, including the recently emerged, multidrug-resistant (MDR) Candida auris. We show that homologs of a known yeast adhesin family in Candida albicans, the Hyr/Iff-like (Hil) family, are enriched in distinct clades of Candida species as a result of multiple, independent expansions. Following gene duplication, the tandem repeat-rich region in these proteins diverged extremely rapidly and generated large variations in length and β-aggregation potential, both of which are known to directly affect adhesion. The conserved N-terminal effector domain was predicted to adopt a β-helical fold followed by an α-crystallin domain, making it structurally similar to a group of unrelated bacterial adhesins. Evolutionary analyses of the effector domain in C. auris revealed relaxed selective constraint combined with signatures of positive selection, suggesting functional diversification after gene duplication. Lastly, we found the Hil family genes to be enriched at chromosomal ends, which likely contributed to their expansion via ectopic recombination and break-induced replication. Combined, these results suggest that the expansion and diversification of adhesin families generate variation in adhesion and virulence within and between species and are a key step toward the emergence of fungal pathogens.

Evolutionary conservation of sequence motifs at sites of protein modification

Gene duplications are common in biology and are likely to be an important source of functional diversification and specialization. The yeast Saccharomyces cerevisiae underwent a whole genome duplication event early in evolution, and a substantial number of duplicated genes have been retained. We identified more than 3,500 instances where only one of two paralogous proteins undergoes post-translational modification despite having retained the same amino acid residue in both. We also developed a web-based search algorithm (CoSMoS.c.) that scores conservation of amino acid sequences based on 1011 wild and domesticated yeast isolates and used it to compare differentially-modified pairs of paralogous proteins. We found that the most common modifications - phosphorylation, ubiquitylation and acylation but not N-glycosylation - occur in regions of high sequence conservation. Such conservation is evident even for ubiquitylation and succinylation, where there is no established 'consensus site' for modification. Differences in phosphorylation were not associated with predicted secondary structure or solvent accessibility, but did mirror known differences in kinase-substrate interactions. Thus, differences in post-translational modification likely result from differences in adjoining amino acids and their interactions with modifying enzymes. By integrating data from large scale proteomics and genomics analysis, in a system with such substantial genetic diversity, we obtained a more comprehensive understanding of the functional basis for genetic redundancies that have persisted for 100 million years.

Protein S-Acyl Transferase GhPAT27 Was Associated with Verticillium wilt Resistance in Cotton

Protein palmitoylation is an ability of the frame of the cell marker protein is one of the most notable reversible changes after translation. However, studies on protein palmitoylation in cotton have not yet been performed. In our current research, the PAT gene family was systematically identified and bioinformatically analyzed in G. arboreum, G. raimondii, G. barbadense and G. hirsutum, and 211 PAT genes were authenticated and classified into six subfamilies. Sixty-nine PAT genes were identified in upland cotton, mainly at the ends of its the 26 chromosomes of upland cotton. The majority of these genes are located in the nucleus of the plant. Gene structure analysis revealed that each member encodes a protein that which contains at least one DHHC structural domain. Cis-acting element analysis indicated that GhPATs genes are mainly involved in hormone production, light response and stress response. Gene expression pattern analysis indicated that most GhPATs genes were differentially expressed upon induction by pathogenic bacteria, drought, salt, hot and cold stresses, and some GhPATs could be induced by multiple abiotic stresses simultaneously. GhPATs genes showed different expression patterns in tissue-specific assays and were found to be preferentially expressed in roots, followed by expression in stems and leaves. Virus-induced gene silencing (VIGS) experiments showed that cotton was significantly less resistant to Verticillium dahliae when GhPAT27 was silenced. We conclude that the GhPAT27 gene, which mediates S-palmitoylation acetylation, may be involved in the regulation of upland cotton resistance to Verticillium wilt (VW). Overall, this work has provided a fundamental framework for understanding the latent capabilities of GhPATs and a solid foundation for molecular breeding and plant pathogen resistance in cotton.

An overview of online resources for intra-species detection of gene duplications

Gene duplication plays an important role in evolutionary mechanism, which can act as a new source of genetic material in genome evolution. However, detecting duplicate genes from genomic data can be challenging. Various bioinformatics resources have been developed to identify duplicate genes from single and/or multiple species. Here, we summarize the metrics used to measure sequence identity among gene duplicates within species, compare several computational approaches that have been used to predict gene duplicates, and review recent advancements of a Basic Local Alignment Search Tool (BLAST)-based web tool and database, allowing future researchers to easily identify intra-species gene duplications. This article is a quick reference guide for research tools used for detecting gene duplicates.

Genome-guided prediction of acid resistance mechanisms in acidophilic methanotrophs of phylogenetically deep-rooted Verrucomicrobia isolated from geothermal environments

Verrucomicrobia are a group of microorganisms that have been proposed to be deeply rooted in the Tree of Life. Some are methanotrophs that oxidize the potent greenhouse gas methane and are thus important in decreasing atmospheric concentrations of the gas, potentially ameliorating climate change. They are widespread in various environments including soil and fresh or marine waters. Recently, a clade of extremely acidophilic Verrucomicrobia, flourishing at pH < 3, were described from high-temperature geothermal ecosystems. This novel group could be of interest for studies about the emergence of life on Earth and to astrobiologists as homologs for possible extraterrestrial life. In this paper, we describe predicted mechanisms for survival of this clade at low pH and suggest its possible evolutionary trajectory from an inferred neutrophilic ancestor. Extreme acidophiles are defined as organisms that thrive in extremely low pH environments (≤ pH 3). Many are polyextremophiles facing high temperatures and high salt as well as low pH. They are important to study for both providing fundamental insights into biological mechanisms of survival and evolution in such extreme environments and for understanding their roles in biotechnological applications such as industrial mineral recovery (bioleaching) and mitigation of acid mine drainage. They are also, potentially, a rich source of novel genes and pathways for the genetic engineering of microbial strains. Acidophiles of the Verrucomicrobia phylum are unique as they are the only known aerobic methanotrophs that can grow optimally under acidic (pH 2–3) and moderately thermophilic conditions (50–60°C). Three moderately thermophilic genera, namely Methylacidiphilum, Methylacidimicrobium, and Ca. Methylacidithermus, have been described in geothermal environments. Most of the investigations of these organisms have focused on their methane oxidizing capabilities (methanotrophy) and use of lanthanides as a protein cofactor, with no extensive study that sheds light on the mechanisms that they use to flourish at extremely low pH. In this paper, we extend the phylogenetic description of this group of acidophiles using whole genome information and we identify several mechanisms, potentially involved in acid resistance, including “first line of defense” mechanisms that impede the entry of protons into the cell. These include the presence of membrane-associated hopanoids, multiple copies of the outer membrane protein (Slp), and inner membrane potassium channels (kup, kdp) that generate a reversed membrane potential repelling the intrusion of protons. Acidophilic Verrucomicrobia also display a wide array of proteins potentially involved in the “second line of defense” where protons that evaded the first line of defense and entered the cell are expelled or neutralized, such as the glutamate decarboxylation (gadAB) and phosphate-uptake systems. An exclusive N-type ATPase F₀-F₁ was identified only in acidophiles of Verrucomicrobia and is predicted to be a specific adaptation in these organisms. Phylogenetic analyses suggest that many predicted mechanisms are evolutionarily conserved and most likely entered the acidophilic lineage of Verrucomicrobia by vertical descent from a common ancestor. However, it is likely that some defense mechanisms such as gadA and kup entered the acidophilic Verrucomicrobia lineage by horizontal gene transfer.

Bifurcation drives the evolution of assembly-line biosynthesis

Reprogramming biosynthetic assembly-lines is a topic of intense interest. This is unsurprising as the scaffolds of most antibiotics in current clinical use are produced by such pathways. The modular nature of assembly-lines provides a direct relationship between the sequence of enzymatic domains and the chemical structure of the product, but rational reprogramming efforts have been met with limited success. To gain greater insight into the design process, we wanted to examine how Nature creates assembly-lines and searched for biosynthetic pathways that might represent evolutionary transitions. By examining the biosynthesis of the anti-tubercular wollamides, we uncover how whole gene duplication and neofunctionalization can result in pathway bifurcation. We show that, in the case of the wollamide biosynthesis, neofunctionalization is initiated by intragenomic recombination. This pathway bifurcation leads to redundancy, providing the genetic robustness required to enable large structural changes during the evolution of antibiotic structures. Should the new product be non-functional, gene loss can restore the original genotype. However, if the new product confers an advantage, depreciation and eventual loss of the original gene creates a new linear pathway. This provides the blind watchmaker equivalent to the design, build, test cycle of synthetic biology.

Functional Characterization of Soybean Diacylglycerol Acyltransferase 3 in Yeast and Soybean

Diacylglycerol acyltransferases (DGAT) function as the key rate-limiting enzymes in de novo biosynthesis of triacylglycerol (TAG) by transferring an acyl group from acyl-CoA to sn-3 of diacylglycerol (DAG) to form TAG. Here, two members of the type 3 DGAT gene family, GmDGAT3-1 and GmDGAT3-2, were identified from the soybean (Glycine max) genome. Both of them were predicted to encode soluble cytosolic proteins containing the typical thioredoxin-like ferredoxin domain. Quantitative PCR analysis revealed that GmDGAT3-2 expression was much higher than GmDGAT3-1's in various soybean tissues such as leaves, flowers, and seeds. Functional complementation assay using TAG-deficient yeast (Saccharomyces cerevisiae) mutant H1246 demonstrated that GmDGAT3-2 fully restored TAG biosynthesis in the yeast and preferentially incorporated monounsaturated fatty acids (MUFAs), especially oleic acid (C18:1) into TAGs. This substrate specificity was further verified by fatty-acid feeding assays and in vitro enzyme activity characterization. Notably, transgenic tobacco (Nicotiana benthamiana) data showed that heterogeneous expression of GmDGAT3-2 resulted in a significant increase in seed oil and C18:1 levels but little change in contents of protein and starch compared to the EV-transformed tobacco plants. Taken together, GmDGAT3-2 displayed a strong enzymatic activity to catalyze TAG assembly with high substrate specificity for MUFAs, particularly C18:1, playing an important role in the cytosolic pathway of TAG synthesis in soybean. The present findings provide a scientific reference for improving oil yield and FA composition in soybean through gene modification, further expanding our knowledge of TAG biosynthesis and its regulatory mechanism in oilseeds.

Convergent Loss of Chemoreceptors across Independent Origins of Slave-Making in Ants

The evolution of an obligate parasitic lifestyle often leads to the reduction of morphological and physiological traits, which may be accompanied by loss of genes and functions. Slave-making ants are social parasites that exploit the work force of closely related ant species for social behaviors such as brood care and foraging. Recent divergence between these social parasites and their hosts enables comparative studies of gene family evolution. We sequenced the genomes of eight ant species, representing three independent origins of ant slavery. During the evolution of eusociality, chemoreceptor genes multiplied due to the importance of chemical communication in insect societies. We investigated the evolutionary fate of these chemoreceptors and found that slave-making ant genomes harbored only half as many gustatory receptors as their hosts', potentially mirroring the outsourcing of foraging tasks to host workers. In addition, parasites had fewer odorant receptors and their loss shows striking patterns of convergence across independent origins of parasitism, in particular in orthologs often implicated in sociality like the 9-exon odorant receptors. These convergent losses represent a rare case of convergent molecular evolution at the level of individual genes. Thus, evolution can operate in a way that is both repeatable and reversible when independent ant lineages lose important social traits during the transition to a parasitic lifestyle.

A constitutive stress response is a result of low temperature growth in the Antarctic green alga Chlamydomonas sp. UWO241

The Antarctic green alga Chlamydomonas sp. UWO241 is an obligate psychrophile that thrives in the cold (4-6°C) but is unable to survive at temperatures ≥18°C. Little is known how exposure to heat affects its physiology or whether it mounts a heat stress response in a manner comparable to mesophiles. Here, we dissect the responses of UWO241 to temperature stress by examining its growth, primary metabolome and transcriptome under steady-state low temperature and heat stress conditions. In comparison with Chlamydomonas reinhardtii, UWO241 constitutively accumulates metabolites and proteins commonly considered as stress markers, including soluble sugars, antioxidants, polyamines, and heat shock proteins to ensure efficient protein folding at low temperatures. We propose that this results from life at extreme conditions. A shift from 4°C to a non-permissive temperature of 24°C alters the UWO241 primary metabolome and transcriptome, but growth of UWO241 at higher permissive temperatures (10 and 15°C) does not provide enhanced heat protection. UWO241 also fails to induce the accumulation of HSPs when exposed to heat, suggesting that it has lost the ability to fine-tune its heat stress response. Our work adds to the growing body of research on temperature stress in psychrophiles, many of which are threatened by climate change.

HSDFinder: A BLAST-Based Strategy for Identifying Highly Similar Duplicated Genes in Eukaryotic Genomes

Gene duplication is an important evolutionary mechanism capable of providing new genetic material for adaptive and nonadaptive evolution. However, bioinformatics tools for identifying duplicate genes are often limited to the detection of paralogs in multiple species or to specific types of gene duplicates, such as retrocopies. Here, we present a user-friendly, BLAST-based web tool, called HSDFinder, which can identify, annotate, categorize, and visualize highly similar duplicate genes (HSDs) in eukaryotic nuclear genomes. HSDFinder includes an online heatmap plotting option, allowing users to compare HSDs among different species and visualize the results in different Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway functional categories. The external software requirements are BLAST, InterProScan, and KEGG. The utility of HSDFinder was tested on various model eukaryotic species, including Chlamydomonas reinhardtii, Arabidopsis thaliana, Oryza sativa, and Zea mays as well as the psychrophilic green alga Chlamydomonas sp. UWO241, and was proven to be a practical and accurate tool for gene duplication analyses. The web tool is free to use at http://hsdfinder.com. Documentation and tutorials can be found via the GitHub: https://github.com/zx0223winner/HSDFinder.

Helicobacter pylori virulence factors: relationship between genetic variability and phylogeographic origin

Background

Helicobacter pylori is a pathogenic bacteria that colonize the gastrointestinal tract from human stomachs and causes diseases including gastritis, peptic ulcers, gastric lymphoma (MALT), and gastric cancer, with a higher prevalence in developing countries. Its high genetic diversity among strains is caused by a high mutation rate, observing virulence factors (VFs) variations in different geographic lineages. This study aimed to postulate the genetic variability associated with virulence factors present in the Helicobacter pylori strains, to identify the relationship of these genes with their phylogeographic origin.

Methods

The complete genomes of 135 strains available in NCBI, from different population origins, were analyzed using bioinformatics tools, identifying a high rate; as well as reorganization events in 87 virulence factor genes, divided into seven functional groups, to determine changes in position, number of copies, nucleotide identity and size, contrasting them with their geographical lineage and pathogenic phenotype.

Results

Bioinformatics analyses show a high rate of gene annotation errors in VF. Analysis of genetic variability of VFs shown that there is not a direct relationship between the reorganization and geographic lineage. However, regarding the pathogenic phenotype demonstrated in the analysis of many copies, size, and similarity when dividing the strains that possess and not the cag pathogenicity island (cagPAI), having a higher risk of developing gastritis and peptic ulcer was evidenced. Our data has shown that the analysis of the overall genetic variability of all VFs present in each strain of H. pylori is key information in understanding its pathogenic behavior.

The Evolution of Hemocyanin Genes in Caenogastropoda: Gene Duplications and Intron Accumulation in Highly Diverse Gastropods

Hemocyanin is the oxygen transport protein of most molluscs and represents an important physiological factor that has to be well-adapted to their environments because of the strong influences of abiotic factors on its oxygen affinity. Multiple independent gene duplications and intron gains have been reported for hemocyanin genes of Tectipleura (Heterobranchia) and the caenogastropod species Pomacea canaliculata, which contrast with the uniform gene architectures of hemocyanins in Vetigastropoda. The goal of this study was to analyze hemocyanin gene evolution within the diverse group of Caenogastropoda in more detail. Our findings reveal multiple gene duplications and intron gains and imply that these represent general features of Apogastropoda hemocyanins. Whereas hemocyanin exon–intron structures are identical within different Tectipleura lineages, they differ strongly within Caenogastropoda among phylogenetic groups as well as between paralogous hemocyanin genes of the same species. Thus, intron accumulation took place more gradually within Caenogastropoda but finally led to a similar consequence, namely, a multitude of introns. Since both phenomena occurred independently within Heterobranchia and Caenogastropoda, the results support the hypothesis that introns may contribute to adaptive radiation by offering new opportunities for genetic variability (multiple paralogs that may evolve differently) and regulation (multiple introns). Our study indicates that adaptation of hemocyanin genes may be one of several factors that contributed to the evolution of the large diversity of Apogastropoda. While questions remain, this hypothesis is presented as a starting point for the further study of hemocyanin genes and possible correlations between hemocyanin diversity and adaptive radiation.

Investigation of RNA metabolism through large-scale genetic interaction profiling in yeast

Gene deletion and gene expression alteration can lead to growth defects that are amplified or reduced when a second mutation is present in the same cells. We performed 154 genetic interaction mapping (GIM) screens with query mutants related with RNA metabolism and estimated the growth rates of about 700 000 double mutant Saccharomyces cerevisiae strains. The tested targets included the gene deletion collection and 900 strains in which essential genes were affected by mRNA destabilization (DAmP). To analyze the results, we developed RECAP, a strategy that validates genetic interaction profiles by comparison with gene co-citation frequency, and identified links between 1471 genes and 117 biological processes. In addition to these large-scale results, we validated both enhancement and suppression of slow growth measured for specific RNA-related pathways. Thus, negative genetic interactions identified a role for the OCA inositol polyphosphate hydrolase complex in mRNA translation initiation. By analysis of suppressors, we found that Puf4, a Pumilio family RNA binding protein, inhibits ribosomal protein Rpl9 function, by acting on a conserved UGUAcauUA motif located downstream the stop codon of the RPL9B mRNA. Altogether, the results and their analysis should represent a useful resource for discovery of gene function in yeast.

Expression attenuation as a mechanism of robustness against gene duplication

Many studies have focused on the mechanisms of long-term retention of gene duplicates, such as the gain of functions or reciprocal losses. However, such changes are more likely to occur if the duplicates are maintained for a long period. This time span will be short if duplication is immediately deleterious. We measured the distribution of fitness effects of gene duplication for 899 genes in budding yeast. We find that gene duplication is more likely to be deleterious than beneficial. However, contrary to previous models, in general, gene duplication does not affect fitness by altering the organization of protein complexes. We show that expression attenuation may protect complexes from the effects of gene duplication.

Gene duplication is ubiquitous and a major driver of phenotypic diversity across the tree of life, but its immediate consequences are not fully understood. Deleterious effects would decrease the probability of retention of duplicates and prevent their contribution to long-term evolution. One possible detrimental effect of duplication is the perturbation of the stoichiometry of protein complexes. Here, we measured the fitness effects of the duplication of 899 essential genes in the budding yeast using high-resolution competition assays. At least 10% of genes caused a fitness disadvantage when duplicated. Intriguingly, the duplication of most protein complex subunits had small to nondetectable effects on fitness, with few exceptions. We selected four complexes with subunits that had an impact on fitness when duplicated and measured the impact of individual gene duplications on their protein–protein interactions. We found that very few duplications affect both fitness and interactions. Furthermore, large complexes such as the 26S proteasome are protected from gene duplication by attenuation of protein abundance. Regulatory mechanisms that maintain the stoichiometric balance of protein complexes may protect from the immediate effects of gene duplication. Our results show that a better understanding of protein regulation and assembly in complexes is required for the refinement of current models of gene duplication.

Diversification of DNA-Binding Specificity by Permissive and Specificity-Switching Mutations in the ParB/Noc Protein Family

Specific interactions between proteins and DNA are essential to many biological processes. Yet, it remains unclear how the diversification in DNA-binding specificity was brought about, and the mutational paths that led to changes in specificity are unknown. Using a pair of evolutionarily related DNA-binding proteins, each with a different DNA preference (ParB [Partitioning Protein B] and Noc [Nucleoid Occlusion Factor], which both play roles in bacterial chromosome maintenance), we show that specificity is encoded by a set of four residues at the protein-DNA interface. Combining X-ray crystallography and deep mutational scanning of the interface, we suggest that permissive mutations must be introduced before specificity-switching mutations to reprogram specificity and that mutational paths to new specificity do not necessarily involve dual-specificity intermediates. Overall, our results provide insight into the possible evolutionary history of ParB and Noc and, in a broader context, might be useful for understanding the evolution of other classes of DNA-binding proteins.

•

DNA-binding specificity for parS and NBS is conserved within ParB and Noc family

•

Specificity is encoded by a set of four residues at the protein-DNA interface

•

Mutations must be introduced in a defined order to reprogram specificity

The evolution of opsin genes in five species of mirid bugs: duplication of long-wavelength opsins and loss of blue-sensitive opsins

BACKGROUND

Color vision and phototactic behavior based on opsins are important for the fitness of insects because of their roles in foraging and mate choice. Related topics, including the duplication and loss of opsin genes, have been well investigated in insect orders such as Coleoptera, Lepidoptera, Hymenoptera, Odonata and Orthoptera, and the findings have been used to develop pest management strategies involving light trapping. Mirid bugs of Hemiptera, which are pests that cause heavy economic losses, show capacity for color discrimination and phototaxis. However, the opsins in mirid bugs remain uncharacterized. Herein, we examined five species to investigate the evolution of opsins in the family Miridae.

RESULTS

Using RNA-seq, we identified several contigs showing high identity with opsins, including four contigs in Apolygus lucorum and three contigs each in Adelphocoris suturalis, Adelphocoris fasciaticollis, Adelphocoris lineolatus and Nesidiocoris tenuis. Phylogenetic analyses indicated that one of these genes clustered with ultraviolet-sensitive (UV) opsins and that the others clustered with long-wavelength (LW) opsins, suggesting that duplication of LW opsins and loss of blue light-sensitive (B) opsins occurred in mirid bugs. The existence of introns in the LW opsins of mirid bugs suggested that the duplication events were DNA based. Both LW1 and LW2 opsins of mirid bugs were found to be under strong purifying selection. The LW1 opsins were significantly more highly expressed than the LW2 and UV opsins.

CONCLUSIONS

We identified the opsins of mirid bugs using five selected mirid species as a representative sample. Phylogenetic analyses clustered one of the genes with UV opsins and the others with LW opsins, suggesting the occurrence of LW opsin duplication and B opsin loss during the evolution of mirid bugs. Intron detection suggested that the identified duplication event was DNA based. The evidence of strong purifying selection and the relatively high expression levels suggested that these opsins exhibit fundamental functions in mirid bugs.

Reduced Translational Efficiency of Eukaryotic Genes after Duplication Events

Control of gene expression has been found to be predominantly determined at the level of protein translation. However, to date, reduced expression from duplicated genes in eukaryotes for dosage maintenance has only been linked to transcriptional control involving epigenetic mechanisms. Here, we hypothesize that dosage maintenance following gene duplication also involves regulation at the protein level. To test this hypothesis, we compared transcriptome and proteome data of yeast models, Saccharomyces cerevisiae and Schizosaccharomyces pombe, and worm models, Caenorhabditis elegans and Caenorhabditis briggsae, to investigate lineage-specifically duplicated genes. Duplicated genes in both eukaryotic models exhibited a reduced protein-to-mRNA abundance ratio. Moreover, dosage sensitive genes, represented by genes encoding protein complex subunits, reduced their protein-to-mRNA abundance ratios more significantly than the other genes after duplication events. An analysis of ribosome profiling (Ribo-Seq) data further showed that reduced translational efficiency was more prominent for dosage sensitive genes than for the other genes. Meanwhile, no difference in protein degradation rate was associated with duplication events. Translationally repressed duplicated genes were also more likely to be inhibited at the level of transcription. Taken together, these results suggest that translation-mediated dosage control is partially contributed by natural selection and it enhances transcriptional control in maintaining gene dosage after gene duplication events during eukaryotic genome evolution.

Global network analysis in Schizosaccharomyces pombe reveals three distinct consequences of the common 1-kb deletion causing juvenile CLN3 disease

Juvenile CLN3 disease is a recessively inherited paediatric neurodegenerative disorder, with most patients homozygous for a 1-kb intragenic deletion in CLN3. The btn1 gene is the Schizosaccharomyces pombe orthologue of CLN3. Here, we have extended the use of synthetic genetic array (SGA) analyses to delineate functional signatures for two different disease-causing mutations in addition to complete deletion of btn1. We show that genetic-interaction signatures can differ for mutations in the same gene, which helps to dissect their distinct functional effects. The mutation equivalent to the minor transcript arising from the 1-kb deletion (btn1^102–208del) shows a distinct interaction pattern. Taken together, our results imply that the minor 1-kb deletion transcript has three consequences for CLN3: to both lose and retain some inherent functions and to acquire abnormal characteristics. This has particular implications for the therapeutic development of juvenile CLN3 disease. In addition, this proof of concept could be applied to conserved genes for other mendelian disorders or any gene of interest, aiding in the dissection of their functional domains, unpacking the global consequences of disease pathogenesis, and clarifying genotype–phenotype correlations. In doing so, this detail will enhance the goals of personalised medicine to improve treatment outcomes and reduce adverse events.

Thermal adaptation rather than demographic history drives genetic structure inferred by copy number variants in a marine fish

Increasing evidence shows that structural variants represent an overlooked aspect of genetic variation with consequential evolutionary roles. Among those, copy number variants (CNVs), including duplicated genomic regions and transposable elements (TEs), may contribute to local adaptation and/or reproductive isolation among divergent populations. Those mechanisms suppose that CNVs could be used to infer neutral and/or adaptive population genetic structure, whose study has been restricted to microsatellites, mitochondrial DNA and Amplified fragment length polymorphism markers in the past and more recently the use of single nucleotide polymorphisms (SNPs). Taking advantage of recent developments allowing CNV analysis from RAD-seq data, we investigated how variation in fitness-related traits, local environmental conditions and demographic history are associated with CNVs, and how subsequent copy number variation drives population genetic structure in a marine fish, the capelin (Mallotus villosus). We collected 1538 DNA samples from 35 sampling sites in the north Atlantic Ocean and identified 6620 putative CNVs. We found associations between CNVs and the gonadosomatic index, suggesting that six duplicated regions could affect female fitness by modulating oocyte production. We also detected 105 CNV candidates associated with water temperature, among which 20% corresponded to genomic regions located within the sequence of protein-coding genes, suggesting local adaptation to cold water by means of gene sequence amplification. We also identified 175 CNVs associated with the divergence of three previously defined parapatric glacial lineages, of which 24% were located within protein-coding genes, making those loci potential candidates for reproductive isolation. Lastly, our analyses unveiled a hierarchical, complex CNV population structure determined by temperature and local geography, which was in stark contrast to that inferred based on SNPs in a previous study. Our findings underline the complementarity of those two types of genomic variation in population genomics studies.

Draft genome sequence of the Antarctic green alga Chlamydomonas sp. UWO241

Antarctica is home to an assortment of psychrophilic algae, which have evolved various survival strategies for coping with their frigid environments. Here, we explore Antarctic psychrophily by examining the ∼212 Mb draft nuclear genome of the green alga Chlamydomonas sp. UWO241, which resides within the water column of a perennially ice-covered, hypersaline lake. Like certain other Antarctic algae, UWO241 encodes a large number (≥37) of ice-binding proteins, putatively originating from horizontal gene transfer. Even more striking, UWO241 harbors hundreds of highly similar duplicated genes involved in diverse cellular processes, some of which we argue are aiding its survival in the Antarctic via gene dosage. Gene and partial gene duplication appear to be an ongoing phenomenon within UWO241, one which might be mediated by retrotransposons. Ultimately, we consider how such a process could be associated with adaptation to extreme environments but explore potential non-adaptive hypotheses as well.

•

Chlamydomonas sp. UWO241 is a green alga originating from Lake Bonney, Antarctica

•

We present a draft nuclear genome sequence of UWO241 (∼212 Mb).

•

The UWO genome contains hundreds of highly similar duplicated genes

•

These duplicates, we argue, might be involved in cold adaptation

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.

Pervasive duplication of tumor suppressors in Afrotherians during the evolution of large bodies and reduced cancer risk

The risk of developing cancer is correlated with body size and lifespan within species. Between species, however, there is no correlation between cancer and either body size or lifespan, indicating that large, long-lived species have evolved enhanced cancer protection mechanisms. Elephants and their relatives (Proboscideans) are a particularly interesting lineage for the exploration of mechanisms underlying the evolution of augmented cancer resistance because they evolved large bodies recently within a clade of smaller-bodied species (Afrotherians). Here, we explore the contribution of gene duplication to body size and cancer risk in Afrotherians. Unexpectedly, we found that tumor suppressor duplication was pervasive in Afrotherian genomes, rather than restricted to Proboscideans. Proboscideans, however, have duplicates in unique pathways that may underlie some aspects of their remarkable anti-cancer cell biology. These data suggest that duplication of tumor suppressor genes facilitated the evolution of increased body size by compensating for decreasing intrinsic cancer risk.

A network-based comparative framework to study conservation and divergence of proteomes in plant phylogenies

Comparative functional genomics offers a powerful approach to study species evolution. To date, the majority of these studies have focused on the transcriptome in mammalian and yeast phylogenies. Here, we present a novel multi-species proteomic dataset and a computational pipeline to systematically compare the protein levels across multiple plant species. Globally we find that protein levels diverge according to phylogenetic distance but is more constrained than the mRNA level. Module-level comparative analysis of groups of proteins shows that proteins that are more highly expressed tend to be more conserved. To interpret the evolutionary patterns of conservation and divergence, we develop a novel network-based integrative analysis pipeline that combines publicly available transcriptomic datasets to define co-expression modules. Our analysis pipeline can be used to relate the changes in protein levels to different species-specific phenotypic traits. We present a case study with the rhizobia-legume symbiosis process that supports the role of autophagy in this symbiotic association.

Genomic signatures of globally enhanced gene duplicate accumulation in the megadiverse higher Diptera fueling intralocus sexual conflict resolution

Gene duplication is an important source of evolutionary innovation. To explore the relative impact of gene duplication during the diversification of major insect model system lineages, we performed a comparative analysis of lineage-specific gene duplications in the fruit fly Drosophila melanogaster (Diptera: Brachycera), the mosquito Anopheles gambiae (Diptera: Culicomorpha), the red flour beetle Tribolium castaneum (Coleoptera), and the honeybee Apis mellifera (Hymenoptera). Focusing on close to 6,000 insect core gene families containing maximally six paralogs, we detected a conspicuously higher number of lineage-specific duplications in Drosophila (689) compared to Anopheles (315), Tribolium (386), and Apis (223). Based on analyses of sequence divergence, phylogenetic distribution, and gene ontology information, we present evidence that an increased background rate of gene duplicate accumulation played an exceptional role during the diversification of the higher Diptera (Brachycera), in part by providing enriched opportunities for intralocus sexual conflict resolution, which may have boosted speciation rates during the early radiation of the megadiverse brachyceran subclade Schizophora.

Autosomal sdY Pseudogenes Explain Discordances Between Phenotypic Sex and DNA Marker for Sex Identification in Atlantic Salmon

Despite the key role that sex-determination plays in evolutionary processes, it is still poorly understood in many species. In salmonids, which are among the best studied fishes, the master sex-determining gene sexually dimorphic on the Y-chromosome (sdY) has been identified. However, sdY displays unexplained discordance to the phenotypic sex, with a variable frequency of phenotypic females being reported as genetic males. Multiple sex determining loci in Atlantic salmon have also been reported, possibly as a result of recent transposition events in this species. We hypothesized the existence of an autosomal copy of sdY, causing apparent discordance between phenotypic and genetic sex, that is transmitted in accordance with autosomal inheritance. To test this, we developed a qPCR methodology to detect the total number of sdY copies present in the genome. Based on the observed phenotype/genotype frequencies and linkage analysis among 2,025 offspring from 64 pedigree-controlled families of accurately phenotyped Atlantic salmon, we identified both males and females carrying one or two autosomal copies of sdY in addition to the Y-specific copy present in males. Patterns across families were highly consistent with autosomal inheritance. These autosomal sdY copies appear to have lost the ability to function as a sex determining gene and were only occasionally assigned to the actual sex chromosome in any of the affected families.

Molecular Evolution of the Glutathione S-Transferase Family in the Bemisia tabaci Species Complex

The glutathione S-transferase (GST) family plays an important role in the adaptation of herbivorous insects to new host plants and other environmental constrains. The family codes for enzymes that neutralize reactive oxygen species and phytotoxins through the conjugation of reduced glutathione. Here, we studied the molecular evolution of the GST family in Bemisia tabaci, a complex of >35 sibling species, differing in their geographic and host ranges. We tested if some enzymes evolved different functionality, by comparing their sequences in six species, representing five of the six major genetic clades in the complex. Comparisons of the nonsynonymous to synonymous substitution ratios detected positive selection events in 11 codons of 5 cytosolic GSTs. Ten of them are located in the periphery of the GST dimer, suggesting a putative involvement in interactions with other proteins. Modeling the tertiary structure of orthologous enzymes, identified additional 19 mutations in 9 GSTs, likely affecting the enzymes' functionality. Most of the mutation events were found in the environmentally responsive classes Delta and Sigma, indicating a slightly different delta/sigma tool box in each species. At a broader genomic perspective, our analyses indicated a significant expansion of the Delta GST class in B. tabaci and a general association between the diet breadth of hemipteran species and their total number of GST genes. We raise the possibility that at least some of the identified changes improve the fitness of the B. tabaci species carrying them, leading to their better adaptation to specific environments.

Structural Insights into the Molecular Evolution of the Archaeal Exo-β-d-Glucosaminidase

The archaeal exo-β-d-glucosaminidase (GlmA), a thermostable enzyme belonging to the glycosidase hydrolase (GH) 35 family, hydrolyzes chitosan oligosaccharides into monomer glucosamines. GlmA is a novel enzyme in terms of its primary structure, as it is homologous to both GH35 and GH42 β-galactosidases. The catalytic mechanism of GlmA is not known. Here, we summarize the recent reports on the crystallographic analysis of GlmA. GlmA is a homodimer, with each subunit comprising three distinct domains: a catalytic TIM-barrel domain, an α/β domain, and a β1 domain. Surprisingly, the structure of GlmA presents features common to GH35 and GH42 β-galactosidases, with the domain organization resembling that of GH42 β-galactosidases and the active-site architecture resembling that of GH35 β-galactosidases. Additionally, the GlmA structure also provides critical information about its catalytic mechanism, in particular, on how the enzyme can recognize glucosamine. Finally, we postulate an evolutionary pathway based on the structure of an ancestor GlmA to extant GH35 and GH42 β-galactosidases.

Genome-wide identification and comparative evolutionary analysis of the Dof transcription factor family in physic nut and castor bean

DNA-binding with one finger (Dof) proteins comprise a plant-specific transcription factor family involved in plant growth, development and stress responses. This study presents a genome-wide comparison of Dof family genes in physic nut (Jatropha curcas) and castor bean (Ricinus communis), two Euphorbiaceae plants that have not experienced any recent whole-genome duplication. A total of 25 or 24 Dof genes were identified from physic nut and castor genomes, respectively, where JcDof genes are distributed across nine out of 11 chromosomes. Phylogenetic analysis assigned these genes into nine groups representing four subfamilies, and 24 orthologous groups were also proposed based on comparison of physic nut, castor, Arabidopsis and rice Dofs. Conserved microsynteny was observed between physic nut and castor Dof-coding scaffolds, which allowed anchoring of 23 RcDof genes to nine physic nut chromosomes. In contrast to how no recent duplicate was present in castor, two tandem duplications and one gene loss were found in the Dof gene family of physic nut. Global transcriptome profiling revealed diverse patterns of Jc/RcDof genes over various tissues, and key Dof genes involved in flower development and stress response were also identified in physic nut. These findings provide valuable information for further studies of Dof genes in physic nut and castor.

Duplication and expression of horizontally transferred polygalacturonase genes is associated with host range expansion of mirid bugs

BACKGROUD

Horizontal gene transfer and gene duplication are two major mechanisms contributing to the evolutionary adaptation of organisms. Previously, polygalacturonase genes (PGs) were independently horizontally transferred and underwent multiple duplications in insects (e.g., mirid bugs and beetles). Here, we chose three phytozoophagous mirid bugs (Adelphocoris suturalis, A. fasciaticollis, A. lineolatus) and one zoophytophagous mirid bug (Nesidiocoris tenuis) to detect whether the duplication, molecular evolution, and expression levels of PGs were related to host range expansion in mirid bugs.

RESULTS

By RNA-seq, we reported 30, 20, 19 and 8 PGs in A. suturalis, A. fasciaticollis, A. lineolatus and N. tenuis, respectively. Interestingly, the number of PGs was significantly positive correlation to the number of host plants (P = 0.0339) in mirid bugs. Most PGs (> 17) were highly expressed in the three phytozoophagous mirid bugs, while only one PG was relatively highly expressed in the zoophytophagous mirid bug. Natural selection analysis clearly showed that a significant relaxation of selection pressure acted on the PGs in zoophytophagous mirid bugs (K = 0.546, P = 0.0158) rather than in phytozoophagous mirid bugs (K = 1, P = 0.92), suggesting a function constraint of PGs in phytozoophagous mirid bugs.

CONCLUSION

Taken together with gene duplication, molecular evolution, and expression levels, our results suggest that PGs are more strictly required by phytozoophagous than by zoophytophagous mirid bugs and that the duplication of PGs is associated with the expansion of host plant ranges in mirid bugs.

A two-level model for the role of complex and young genes in the formation of organism complexity and new insights into the relationship between evolution and development

Background

How genome complexity affects organismal phenotypic complexity is a fundamental question in evolutionary developmental biology. Previous studies proposed various contributing factors of genome complexity and tried to find the connection between genomic complexity and organism complexity. However, a general model to answer this question is lacking. Here, we introduce a ‘two-level’ model for the realization of genome complexity at phenotypic level.

Results

Five representative species across Protostomia and Deuterostomia were involved in this study. The intrinsic gene properties contributing to genome complexity were classified into two generalized groups: the complexity and age degree of both protein-coding and noncoding genes. We found that young genes tend to be simpler; however, the mid-age genes, rather than the oldest genes, show the highest proportion of high complexity. Complex genes tend to be utilized preferentially in each stage of embryonic development, with maximum representation during the late stage of organogenesis. This trend is mainly attributed to mid-age complex genes. In contrast, young genes tend to be expressed in specific spatiotemporal states. An obvious correlation between the time point of the change in over- and under-representation and the order of gene age was observed, which supports the funnel-like model of the conservation pattern of development. In addition, we found some probable causes for the seemingly contradictory ‘funnel-like’ or ‘hourglass’ model.

Conclusions

These results indicate that complex and young genes contribute to organismal complexity at two different levels: Complex genes contribute to the complexity of individual proteomes in certain states, whereas young genes contribute to the diversity of proteomes in different spatiotemporal states. This conclusion is valid across the five species investigated, indicating it is a conserved model across Protostomia and Deuterostomia. The results in this study also support ‘funnel-like model’ from a new viewpoint and explain why there are different evo–devo relation models.

Electronic supplementary material

The online version of this article (10.1186/s13227-018-0111-4) contains supplementary material, which is available to authorized users.

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

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

Differential Substrate Recognition by Maltose Binding Proteins Influenced by Structure and Dynamics

The genome of the hyperthermophile Thermotoga maritima contains three isoforms of maltose binding protein (MBP) that are high-affinity receptors for di-, tri-, and tetrasaccharides. Two of these proteins (tmMBP1 and tmMBP2) share significant sequence identity, approximately 90%, while the third (tmMBP3) shares less than 40% identity. MBP from Escherichia coli (ecMBP) shares 35% sequence identity with the tmMBPs. This subset of MBP isoforms offers an interesting opportunity to investigate the mechanisms underlying the evolution of substrate specificity and affinity profiles in a genome where redundant MBP genes are present. In this study, the X-ray crystal structures of tmMBP1, tmMBP2, and tmMBP3 are reported in the absence and presence of oligosaccharides. tmMBP1 and tmMBP2 have binding pockets that are larger than that of tmMBP3, enabling them to bind to larger substrates, while tmMBP1 and tmMBP2 also undergo substrate-induced hinge bending motions (∼52°) that are larger than that of tmMBP3 (∼35°). Small-angle X-ray scattering was used to compare protein behavior in solution, and computer simulations provided insights into dynamics of these proteins. Comparing quantitative protein-substrate interactions and dynamical properties of tmMBPs with those of the promiscuous ecMBP and disaccharide selective Thermococcus litoralis MBP provides insights into the features that enable selective binding. Collectively, the results provide insights into how the structure and dynamics of tmMBP homologues enable them to differentiate between a myriad of chemical entities while maintaining their common fold.

Reduced intrinsic DNA curvature leads to increased mutation rate

BACKGROUND

Mutation rates vary across the genome. Many trans factors that influence mutation rates have been identified, as have specific sequence motifs at the 1-7-bp scale, but cis elements remain poorly characterized. The lack of understanding regarding why different sequences have different mutation rates hampers our ability to identify positive selection in evolution and to identify driver mutations in tumorigenesis.

RESULTS

Here, we use a combination of synthetic genes and sequences of thousands of isolated yeast colonies to show that intrinsic DNA curvature is a major cis determinant of mutation rate. Mutation rate negatively correlates with DNA curvature within genes, and a 10% decrease in curvature results in a 70% increase in mutation rate. Consistently, both yeast and humans accumulate mutations in regions with small curvature. We further show that this effect is due to differences in the intrinsic mutation rate, likely due to differences in mutagen sensitivity and not due to differences in the local activity of DNA repair.

CONCLUSIONS

Our study establishes a framework for understanding the cis properties of DNA sequence in modulating the local mutation rate and identifies a novel causal source of non-uniform mutation rates across the genome.

Genetic Interaction Network as an Important Determinant of Gene Order in Genome Evolution

Although it is generally accepted that eukaryotic gene order is not random, the basic principles of gene arrangement on a chromosome remain poorly understood. Here, we extended existing population genetics theories that were based on two-locus models and proposed a hypothesis that genetic interaction networks drive the evolution of eukaryotic gene order. We predicted that genes with positive epistasis would move toward each other in evolution, during which a negative correlation between epistasis and gene distance formed. We tested and confirmed our prediction with computational simulations and empirical data analyses. Importantly, we demonstrated that gene order in the budding yeast could be successfully predicted from the genetic interaction network. Taken together, our study reveals the role of the genetic interaction network in the evolution of gene order, extends our understanding of the encoding principles in genomes, and potentially offers new strategies to improve synthetic biology.

Comparative genomics and the nature of placozoan species

Placozoans are a phylum of nonbilaterian marine animals currently represented by a single described species, Trichoplax adhaerens, Schulze 1883. Placozoans arguably show the simplest animal morphology, which is identical among isolates collected worldwide, despite an apparently sizeable genetic diversity within the phylum. Here, we use a comparative genomics approach for a deeper appreciation of the structure and causes of the deeply diverging lineages in the Placozoa. We generated a high-quality draft genome of the genetic lineage H13 isolated from Hong Kong and compared it to the distantly related T. adhaerens. We uncovered substantial structural differences between the two genomes that point to a deep genomic separation and provide support that adaptation by gene duplication is likely a crucial mechanism in placozoan speciation. We further provide genetic evidence for reproductively isolated species and suggest a genus-level difference of H13 to T. adhaerens, justifying the designation of H13 as a new species, Hoilungia hongkongensis nov. gen., nov. spec., now the second described placozoan species and the first in a new genus. Our multilevel comparative genomics approach is, therefore, likely to prove valuable for species distinctions in other cryptic microscopic animal groups that lack diagnostic morphological characters, such as some nematodes, copepods, rotifers, or mites.

Analyses of the genetic diversity and protein expression variation of the acyl: CoA medium-chain ligases, ACSM2A and ACSM2B

Benzoate (found in milk and widely used as preservative), salicylate (present in fruits and the active component of aspirin), dietary polyphenols produced by gut microbiota, metabolites from organic acidemias, and medium-chain fatty acids (MCFAs) are all metabolised/detoxified by the glycine conjugation pathway. Xenobiotics are first activated to an acyl-CoA by the mitochondrial xenobiotic/medium-chain fatty acid: CoA ligases (ACSMs) and subsequently conjugated to glycine by glycine N-acyltransferase (GLYAT). The MCFAs are activated to acyl-CoA by the ACSMs before entering mitochondrial β-oxidation. This two-step enzymatic pathway has, however, not been thoroughly investigated and the biggest gap in the literature remains the fact that studies continuously characterise the pathway as a one-step reaction. There are no studies available on the interaction/competition of the various substrates involved in the pathway, whilst very little research has been done on the ACSM ligases. To identify variants/haplotypes that should be characterised in future detoxification association studies, this study assessed the naturally observed sequence diversity and protein expression variation of ACSM2A and ACSM2B. The allelic variation, haplotype diversity, Tajima's D values, and phylogenetic analyses indicated that ACSM2A and ACSM2B are highly conserved. This confirmed an earlier hypothesis that the glycine conjugation pathway is highly conserved and essential for life as it maintains the CoA and glycine homeostasis in the liver mitochondria. The protein expression analyses showed that ACSM2A is the predominant transcript in liver. Future studies should investigate the effect of the variants identified in this study on the substrate specificity of these proteins.

MicroRNA duplication accelerates the recruitment of new targets during vertebrate evolution

The repertoire of miRNAs has considerably expanded during metazoan evolution, and duplication is an important mechanism for generating new functional miRNAs. However, relatively little is known about the functional divergence between paralogous miRNAs and the possible coevolution between duplicated miRNAs and the genomic contexts. By systematically examining small RNA expression profiles across various human tissues and interrogating the publicly available miRNA:mRNA pairing chimeras, we found that changes in expression patterns and targeting preferences are widespread for duplicated miRNAs in vertebrates. Both the empirical interactions and target predictions suggest that evolutionarily conserved homo-seed duplicated miRNAs pair with significantly higher numbers of target sites compared to the single-copy miRNAs. Our birth-and-death evolutionary analysis revealed that the new target sites of miRNAs experienced frequent gains and losses during function development. Our results suggest that a newly emerged target site has a higher probability to be functional and maintained by natural selection if it is paired to a seed shared by multiple paralogous miRNAs rather than being paired to a single-copy miRNA. We experimentally verified the divergence in target repression between two paralogous miRNAs by transfecting let-7a and let-7b mimics into kidney-derived cell lines of four mammalian species and measuring the resulting transcriptome alterations by extensive high-throughput sequencing. Our results also suggest that the gains and losses of let-7 target sites might be associated with the evolution of repressiveness of let-7 across mammalian species.

Selective Constraints on Coding Sequences of Nervous System Genes Are a Major Determinant of Duplicate Gene Retention in Vertebrates

The evolutionary history of vertebrates is marked by three ancient whole-genome duplications: two successive rounds in the ancestor of vertebrates, and a third one specific to teleost fishes. Biased loss of most duplicates enriched the genome for specific genes, such as slow evolving genes, but this selective retention process is not well understood. To understand what drives the long-term preservation of duplicate genes, we characterized duplicated genes in terms of their expression patterns. We used a new method of expression enrichment analysis, TopAnat, applied to in situ hybridization data from thousands of genes from zebrafish and mouse. We showed that the presence of expression in the nervous system is a good predictor of a higher rate of retention of duplicate genes after whole-genome duplication. Further analyses suggest that purifying selection against the toxic effects of misfolded or misinteracting proteins, which is particularly strong in nonrenewing neural tissues, likely constrains the evolution of coding sequences of nervous system genes, leading indirectly to the preservation of duplicate genes after whole-genome duplication. Whole-genome duplications thus greatly contributed to the expansion of the toolkit of genes available for the evolution of profound novelties of the nervous system at the base of the vertebrate radiation.

CNVs are associated with genomic architecture in a songbird

Background

Understanding variation in genome structure is essential to understand phenotypic differences within populations and the evolutionary history of species. A promising form of this structural variation is copy number variation (CNV). CNVs can be generated by different recombination mechanisms, such as non-allelic homologous recombination, that rely on specific characteristics of the genome architecture. These structural variants can therefore be more abundant at particular genes ultimately leading to variation in phenotypes under selection. Detailed characterization of CNVs therefore can reveal evolutionary footprints of selection and provide insight in their contribution to phenotypic variation in wild populations.

Results

Here we use genotypic data from a long-term population of great tits (Parus major), a widely studied passerine bird in ecology and evolution, to detect CNVs and identify genomic features prevailing within these regions. We used allele intensities and frequencies from high-density SNP array data from 2,175 birds. We detected 41,029 CNVs concatenated into 8,008 distinct CNV regions (CNVRs). We successfully validated 93.75% of the CNVs tested by qPCR, which were sampled at different frequencies and sizes. A mother-daughter family structure allowed for the evaluation of the inheritance of a number of these CNVs. Thereby, only CNVs with 40 probes or more display segregation in accordance with Mendelian inheritance, suggesting a high rate of false negative calls for smaller CNVs. As CNVRs are a coarse-grained map of CNV loci, we also inferred the frequency of coincident CNV start and end breakpoints. We observed frequency-dependent enrichment of these breakpoints at homologous regions, CpG sites and AT-rich intervals. A gene ontology enrichment analyses showed that CNVs are enriched in genes underpinning neural, cardiac and ion transport pathways.

Conclusion

Great tit CNVs are present in almost half of the genes and prominent at repetitive-homologous and regulatory regions. Although overlapping genes under selection, the high number of false negatives make neutrality or association tests on CNVs detected here difficult. Therefore, CNVs should be further addressed in the light of their false negative rate and architecture to improve the comprehension of their association with phenotypes and evolutionary history.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4577-1) contains supplementary material, which is available to authorized users.

Interactome evolution: insights from genome-wide analyses of protein-protein interactions

We highlight new evolutionary insights enabled by recent genome-wide studies on protein-protein interaction (PPI) networks ('interactomes'). While most PPIs are mediated by a single sequence region promoting or inhibiting interactions, many PPIs are mediated by multiple sequence regions acting cooperatively. Most PPIs perform important functions maintained by negative selection: we estimate that less than ∼10% of the human interactome is effectively neutral upon perturbation (i.e. 'junk' PPIs), and the rest are deleterious upon perturbation; interfacial sites evolve more slowly than other sites; many conserved PPIs show signatures of co-evolution at the interface; PPIs evolve more slowly than protein sequence. At the same time, many PPIs undergo rewiring during evolution for lineage-specific adaptation. Finally, chaperone-protein and host-pathogen interactomes are governed by distinct evolutionary principles.

Detection of Pathways Affected by Positive Selection in Primate Lineages Ancestral to Humans

Gene set enrichment approaches have been increasingly successful in finding signals of recent polygenic selection in the human genome. In this study, we aim at detecting biological pathways affected by positive selection in more ancient human evolutionary history. Focusing on four branches of the primate tree that lead to modern humans, we tested all available protein coding gene trees of the Primates clade for signals of adaptation in these branches, using the likelihood-based branch site test of positive selection. The results of these locus-specific tests were then used as input for a gene set enrichment test, where whole pathways are globally scored for a signal of positive selection, instead of focusing only on outlier "significant" genes. We identified signals of positive selection in several pathways that are mainly involved in immune response, sensory perception, metabolism, and energy production. These pathway-level results are highly significant, even though there is no functional enrichment when only focusing on top scoring genes. Interestingly, several gene sets are found significant at multiple levels in the phylogeny, but different genes are responsible for the selection signal in the different branches. This suggests that the same function has been optimized in different ways at different times in primate evolution.