Accelerated evolution after gene duplication: a time-dependent process affecting just one copy

Genome-wide analysis reveals the contributors to fast molecular evolution of the Chinese hook snout carp (Opsariichthys bidens)

Variations in molecular evolutionary rate have been widely investigated among lineages and genes. However, it remains an open question whether fast rate of molecular evolution is driven by natural selection or random drift, and how the fast rate is linked to metabolic rate. Additionally, previous studies on fast molecular evolution have been largely restricted to concatenated matrix of genes or a few specifically selected genes, but less is known for individual genes at the genome-wide level. Here we addressed these questions using more than 5000 single-copy orthologous (SCO) genes through comparative genomic and phylogenetic analyses among fishes, with a special focus on a newly-sequenced clupeocephalan fish the Chinese hook snout carp Opsariichthys bidens. We showed O. bidens displays significantly higher mean substitution rate and more fast-evolving SCO genes (2172 genes) than most fishes studied here. The rapidly evolving genes are enriched in highly conserved and very basic functions such as translation and ribosome that are critical for biological fitness. We further revealed that ∼25 % of these fast-evolving genes exhibit a constant increase of substitution rate from the common ancestor down to the present, suggesting a neglected but important contribution from ancestral states. Model fitting showed that ∼85 % of fast-evolving genes exclusive to O. bidens and related species follow the adaptive evolutionary model rather than random-drift model, and 7.6 % of fast-evolving genes identified in O. bidens have experienced positive selection, both indicating the reflection of adaptive selection. Finally, metabolic rate was observed to be linked with substitution rate in a gene-specific manner. Overall, our findings reveal fast molecular evolution of SCO genes at genome-wide level in O. bidens, and uncover the evolutionary and ecological contributors to it.

Expansion of MHC-IIB Has Constrained the Evolution of MHC-IIA in Passerines

The major histocompatibility complex (MHC) is central in adaptive immunity, with the highly polymorphic MHC genes encoding antigen-presenting molecules. Two MHC class II (MHC-II) loci, DA1 and DA2, predate the radiation of extant birds and persist throughout much of the avian phylogeny. Within each locus, the MHC-II molecules are encoded by A-genes (DAA) and B-genes (DAB), which are arranged in A-B dyads. However, in passerines (order Passeriformes), the DA2 locus has been lost, and the ancestral A-B dyad at the DA1 locus has been replaced by a putatively single A-gene (DAA1) and an array of highly polymorphic B-genes (DAB1). In this study, we genotyped the DAA1 gene of 15 passerine species and confirmed that passerines possess just one copy of DAA1. We then compared selection patterns in DAA1 between passerines and nonpasserines and found that exon 2, which encodes the antigen-presenting domain, has been subject to weaker positive selection and stronger negative selection in passerines compared with nonpasserines. Additional comparisons showed that the patterns of selection in the passerine DAA1 gene are unlikely to be related to the loss of the DA2 locus. Instead, our findings suggest that the expansion of DAB1 (MHC-IIB) has imposed an evolutionary constraint on the passerine DAA1 (MHC-IIA) gene. We speculate that this constraint may be the result of each DAA1 chain forming heterodimers with many different DAB1 chains.

A telomere-to-telomere Eucalyptus regnans genome: unveiling haplotype variance in structure and genes within one of the world's tallest trees

BACKGROUND

Eucalyptus regnans (Mountain Ash) is an Australian native giant tree species which form forests that are among the highest known carbon-dense biomasses in the world. To enhance genomic studies in this ecologically important species, we assembled a high-quality, mostly telomere-to-telomere complete, chromosome-level, haplotype-resolved reference genome. We sampled a single tree, the Centurion, which is currently a contender for the world's tallest flowering plant.

RESULTS

Using long-read sequencing data (PacBio HiFi, Oxford Nanopore ultra-long reads) and chromosome conformation capture data (Hi-C), we assembled the most contiguous and complete Eucalyptus reference genome to date. For each haplotype, we observed contig N50s exceeding 36 Mbp, scaffold N50s exceeding 43 Mbp, and genome BUSCO completeness exceeding 99%. The assembled genome revealed extensive structural variations between the two haplotypes, consisting mostly of insertions, deletions, duplications and translocations. Analysis of gene content revealed haplotype-specific genes, which were enriched in functional categories related to transcription, energy production and conservation. Additionally, many genes reside within structurally rearranged regions, particularly duplications, suggesting that haplotype-specific variation may contribute to environmental adaptation in the species.

CONCLUSIONS

Our study provides a foundation for future research into E. regnans environmental adaptation, and the high-quality genome will be a powerful resource for conservation of carbon-dense giant tree forests.

Conservation and selective pressures shaping baleen whale olfactory receptor genes supports their use of olfaction in the marine environment

The relative importance of various sensory modalities can shift in response to evolutionary transitions, resulting in changes to underlying gene families encoding their reception systems. The rapid birth-and-death process underlying the evolution of the large olfactory receptor (OR) gene family has accelerated genomic-level change for the sense of smell in particular. The transition from the land to sea in marine mammals is an attractive model for understanding the influence of habitat shifts on sensory systems, with the retained OR repertoire of baleen whales contrasting with its loss in toothed whales. In this study, we examine to what extent the transition from a terrestrial to a marine environment has influenced the evolution of baleen whale OR repertoires. We developed Gene Mining Pipeline (GMPipe) (https://github.com/AprilJauhal/GMPipe), which can accurately identify large numbers of candidate OR genes. GMPipe identified 707 OR sequences from eight baleen whale species. These repertoires exhibited distinct family count distributions compared to terrestrial mammals, including signs of relative expansion in families OR10, OR11 and OR13. While many receptors have been lost or show signs of random drift in baleen whales, others exhibit signs of evolving under purifying or positive selection. Over 85% of OR genes could be sorted into orthologous groups of sequences containing at least four homologous sequences. Many of these groups, particularly from family OR10, presented signs of relative expansion and purifying selective pressure. Overall, our results suggest that the relatively small size of baleen whale OR repertoires result from specialisation to novel olfactory landscapes, as opposed to random drift.

Dynamic Expansions and Retinal Expression of Spectrally Distinct Short-Wavelength Opsin Genes in Sea Snakes

The photopigment-encoding visual opsin genes that mediate color perception show great variation in copy number and adaptive function across vertebrates. An open question is how this variation has been shaped by the interaction of lineage-specific structural genomic architecture and ecological selection pressures. We contribute to this issue by investigating the expansion dynamics and expression of the duplicated Short-Wavelength-Sensitive-1 opsin (SWS1) in sea snakes (Elapidae). We generated one new genome, 45 resequencing datasets, 10 retinal transcriptomes, and 81 SWS1 exon sequences for sea snakes, and analyzed these alongside 16 existing genomes for sea snakes and their terrestrial relatives. Our analyses revealed multiple independent transitions in SWS1 copy number in the marine Hydrophis clade, with at least three lineages having multiple intact SWS1 genes: the previously studied Hydrophis cyanocinctus and at least two close relatives of this species; Hydrophis atriceps and Hydrophis fasciatus; and an individual Hydrophis curtus. In each lineage, gene copy divergence at a key spectral tuning site resulted in distinct UV and Violet/Blue-sensitive SWS1 subtypes. Both spectral variants were simultaneously expressed in the retinae of H. cyanocinctus and H. atriceps, providing the first evidence that these SWS1 expansions confer novel phenotypes. Finally, chromosome annotation for nine species revealed shared structural features in proximity to SWS1 regardless of copy number. If these features are associated with SWS1 duplication, expanded opsin complements could be more common in snakes than is currently recognized. Alternatively, selection pressures specific to aquatic environments could favor improved chromatic distinction in just some lineages.

A comprehensive evolutionary scenario for the origin and neofunctionalization of the Drosophila speciation gene Odysseus (OdsH)

Odysseus (OdsH) was the first speciation gene described in Drosophila related to hybrid sterility in offspring of mating between Drosophila mauritiana and Drosophila simulans. Its origin is attributed to the duplication of the gene unc-4 in the subgenus Sophophora. By using a much larger sample of Drosophilidae species, we showed that contrary to what has been previously proposed, OdsH origin occurred 62 MYA. Evolutionary rates, expression, and transcription factor-binding sites of OdsH evidence that it may have rapidly experienced neofunctionalization in male sexual functions. Furthermore, the analysis of the OdsH peptide allowed the identification of mutations of D. mauritiana that could result in incompatibility in hybrids. In order to find if OdsH could be related to hybrid sterility, beyond Sophophora, we explored the expression of OdsH in Drosophila arizonae and Drosophila mojavensis, a pair of sister species with incomplete reproductive isolation. Our data indicated that OdsH expression is not atypical in their male-sterile hybrids. In conclusion, we have proposed that the origin of OdsH occurred earlier than previously proposed, followed by neofunctionalization. Our results also suggested that its role as a speciation gene might be restricted to D. mauritiana and D. simulans.

Degeneration of the Olfactory System in a Murid Rodent that Evolved Diurnalism

In mammalian research, it has been debated what can initiate an evolutionary tradeoff between different senses, and the phenomenon of sensory tradeoff in rodents, the most abundant mammalian clade, is not evident. The Nile rat (Arvicanthis niloticus), a murid rodent, recently adapted to a diurnal niche through an evolutionary acquisition of daylight vision with enhanced visual acuity. As such, this model provides an opportunity for a cross-species investigation where comparative morphological and multi-omic analyses of the Nile rat are made with its closely related nocturnal species, e.g. the mouse (Mus musculus) and the rat (Rattus norvegicus). Thus, morphological examinations were performed, and evolutionary reductions in relative sizes of turbinal bone surfaces, the cribriform plate, and the olfactory bulb were discovered in Nile rats. Subsequently, we compared multiple murid genomes, and profiled olfactory epithelium transcriptomes of mice and Nile rats at various ages with RNA sequencing. The results further demonstrate that, in comparison with mouse olfactory receptor (OR) genes, Nile rat OR genes have experienced less frequent gain, more frequent loss, and more frequent expression reduction during their evolution. Furthermore, functional degeneration of coding sequences in the Nile rat lineage was found in OR genes, yet not in other genes. Taken together, these results suggest that acquisition of improved vision in the Nile rat has been accompanied by degeneration of both olfaction-related anatomical structures and OR gene repertoires, consistent with the hypothesis of an olfaction-vision tradeoff initiated by the switch from a nocturnal to a diurnal lifestyle in mammals.

Functional Divergence and Origin of the Vertebrate Praja Family

The Praja family is an E3 ubiquitin ligase, promoting polyubiquitination and subsequent degradation of substrates. It comprises two paralogs, praja1 and praja2. Prior research suggests these paralogs have undergone functional divergence, with examples, such as their distinct roles in neurite outgrowth. However, the specific evolutionary trajectories of each paralog remain largely unexplored preventing mechanistic understanding of functional differences between paralogs. Here, we investigated the phylogeny and divergence of the vertebrate Praja family through molecular evolutionary analysis. Phylogenetic examination of the vertebrate praja revealed that praja1 and praja2 originated from the common ancestor of placentals via gene duplication, with praja1 evolving at twice the rate of praja2 shortly after the duplication. Moreover, a unique evolutionary trajectory for praja1 relative to other vertebrate Praja was indicated, as evidenced by principal component analysis on GC content, codon usage frequency, and amino acid composition. Subsequent motif/domain comparison revealed conserved N terminus and C terminus in praja1 and praja2, together with praja1-specific motifs, including nuclear localization signal and Ala-Gly-Ser repeats. The nuclear localization signal was demonstrated to be functional in human neuroblastoma SH-SY5Y cells using deletion mutant, while praja2 was exclusively expressed in the nucleus. These discoveries contribute to a more comprehensive understanding of the Praja family's phylogeny and suggest a functional divergence between praja1 and praja2. Specifically, the shift of praja1 into the nucleus implies the degradation of novel substrates located in the nucleus as an evolutionary consequence.

Leveraging genomic redundancy to improve inference and alignment of orthologous proteins

Identifying protein sequences with common ancestry is a core task in bioinformatics and evolutionary biology. However, methods for inferring and aligning such sequences in annotated genomes have not kept pace with the increasing scale and complexity of the available data. Thus, in this work, we implemented several improvements to the traditional methodology that more fully leverage the redundancy of closely related genomes and the organization of their annotations. Two highlights include the application of the more flexible k-clique percolation algorithm for identifying clusters of orthologous proteins and the development of a novel technique for removing poorly supported regions of alignments with a phylogenetic hidden Markov model (phylo-HMM). In making the latter, we wrote a fully documented Python package Homomorph that implements standard HMM algorithms and created a set of tutorials to promote its use by a wide audience. We applied the resulting pipeline to a set of 33 annotated Drosophila genomes, generating 22,813 orthologous groups and 8,566 high-quality alignments.

The evolution of bitter taste receptor gene in primates: Gene duplication and selection

Bitter taste perception plays an important role in preventing animals from digesting poisonous and harmful substances. In primates, especially the Cercopithecidae species, most species feed on plants; thus, it is reasonable to speculate that most of the bitter taste receptor genes (T2Rs) of primates are under purifying selection to maintain the functional stability of bitter taste perception. Gene duplication has happened in T2Rs frequently, and what will be the fate of T2Rs copies is another question we are concerned about. To answer these questions, we selected the T2Rs of primates reported in another study and conducted corresponding selective pressure analyses to determine what kind of selective pressure was acting on them. Further, we carried out selective pressure analyses on gene copies and their corresponding ancestors by considering several possible situations. The results showed that among the 25 gene groups examined here, 15 groups are subject to purifying selection and others are under relaxed selection, with many positively selected sites detected. Gene copies existed in several groups, but only some groups (clade1_a1-b2, clade1_c-c2, clade1_d1-d3, clade1_f1-f2, T2R10, T2R13, and T2R42) have positively selected sites, inferring that they may have some relation to functional divergence. Taken together, T2Rs in primates are under diverse selective pressures, and most gene copies are subject to the same selective pressures. In such cases, the copies may be just to keep the function conservative, and more copies can increase the quantity of the bitter taste receptor, raise the efficiency of bitter substance recognition, and finally enhance the fitness of feeding during the evolutionary course of primates. This study can improve our understanding of T2Rs evolution in primates.

Population History Shapes Responses to Different Temperature Regimes in Drosophila subobscura

Drosophila subobscura is considered a good model species for investigation of a population's ability to adapt and cope with climate changes. Decade long research has shown that inversion frequencies change in response to environmental factors indicating their role in adaptation to novel environments. The mechanisms behind organisms' responses to temperature are complex, involving changes in physiology, behavior, gene expression and regulation. On the other hand, a population's ability to respond to suboptimal conditions depends on standing genetic variation and population history. In order to elucidate the role of local adaptation in population response to the changing temperature, we investigated the response to temperature in D. subobscura individuals originating from two different altitudes by combining traditional cytogenetic techniques with assessing the levels of Hsp70 protein expression. Inversion polymorphism was assessed in the flies sampled from natural populations and in flies reared in laboratory conditions at three different temperatures after five and sixteen generations and Hsp70 protein expression profile in 12th generation flies at the basal level and after heat shock induction. Our results indicate that local adaptation and population history influence population response to the changing temperature.

Interspecies genome divergence is predominantly due to frequent small scale rearrangements in Eucalyptus

Synteny, the ordering of sequences within homologous chromosomes, must be maintained within the genomes of sexually reproducing species for the sharing of alleles and production of viable, reproducing offspring. However, when the genomes of closely related species are compared, a loss of synteny is often observed. Unequal homologous recombination is the primary mechanism behind synteny loss, occurring more often in transposon rich regions, and resulting in the formation of chromosomal rearrangements. To examine patterns of synteny among three closely related, interbreeding, and wild Eucalyptus species, we assembled their genomes using long-read DNA sequencing and de novo assembly. We identify syntenic and rearranged regions between these genomes and estimate that ~48% of our genomes remain syntenic while ~36% is rearranged. We observed that rearrangements highly fragment microsynteny. Our results suggest that synteny between these species is primarily lost through small-scale rearrangements, not through sequence loss, gain, or sequence divergence. Further examination of identified rearrangements suggests that rearrangements may be altering the phenotypes of Eucalyptus species. Our study also underscores that the use of single reference genomes in genomic variation studies could lead to reference bias, especially given the scale at which we show potentially adaptive loci have highly diverged, deleted, duplicated and/or rearranged. This study provides an unbiased framework to look at potential speciation and adaptive loci among a rapidly radiating foundation species of woodland trees that are free from selective breeding seen in most crop species.

Overexpression of TgERF1, a Transcription Factor from Tectona grandis, Increases Tolerance to Drought and Salt Stress in Tobacco

Teak (Tectona grandis) is one of the most important wood sources, and it is cultivated in tropical regions with a significant market around the world. Abiotic stresses are an increasingly common and worrying environmental phenomenon because it causes production losses in both agriculture and forestry. Plants adapt to these stress conditions by activation or repression of specific genes, and they synthesize numerous stress proteins to maintain their cellular function. For example, APETALA2/ethylene response factor (AP2/ERF) was found to be involved in stress signal transduction. A search in the teak transcriptome database identified an AP2/ERF gene named TgERF1 with a key AP2/ERF domain. We then verified that the TgERF1 expression is rapidly induced by Polyethylene Glycol (PEG), NaCl, and exogenous phytohormone treatments, suggesting a potential role in drought and salt stress tolerance in teak. The full-length coding sequence of TgERF1 gene was isolated from teak young stems, characterized, cloned, and constitutively overexpressed in tobacco plants. In transgenic tobacco plants, the overexpressed TgERF1 protein was localized exclusively in the cell nucleus, as expected for a transcription factor. Furthermore, functional characterization of TgERF1 provided evidence that TgERF1 is a promising candidate gene to be used as selective marker on plant breeding intending to improve plant stress tolerance.

The evolutionary patterns, expression profiles, and genetic diversity of expanded genes in barley

Gene duplication resulting from whole-genome duplication (WGD), small-scale duplication (SSD), or unequal hybridization plays an important role in the expansion of gene families. Gene family expansion can also mediate species formation and adaptive evolution. Barley (Hordeum vulgare) is the world's fourth largest cereal crop, and it contains valuable genetic resources due to its ability to tolerate various types of environmental stress. In this study, 27,438 orthogroups in the genomes of seven Poaceae were identified, and 214 of them were significantly expanded in barley. The evolutionary rates, gene properties, expression profiles, and nucleotide diversity between expanded and non-expanded genes were compared. Expanded genes evolved more rapidly and experienced lower negative selection. Expanded genes, including their exons and introns, were shorter, they had fewer exons, their GC content was lower, and their first exons were longer compared with non-expanded genes. Codon usage bias was also lower for expanded genes than for non-expanded genes; the expression levels of expanded genes were lower than those of non-expanded genes, and the expression of expanded genes showed higher tissue specificity than that of non-expanded genes. Several stress-response-related genes/gene families were identified, and these genes could be used to breed barley plants with greater resistance to environmental stress. Overall, our analysis revealed evolutionary, structural, and functional differences between expanded and non-expanded genes in barley. Additional studies are needed to clarify the functions of the candidate genes identified in our study and evaluate their utility for breeding barley plants with greater stress resistance.

Genome-Wide Identification and Expression Analysis of Kinesin Family in Barley (Hordeum vulgare)

Kinesin, as a member of the molecular motor protein superfamily, plays an essential function in various plants' developmental processes. Especially at the early stages of plant growth, including influences on plants' growth rate, yield, and quality. In this study, we did a genome-wide identification and expression profile analysis of the kinesin family in barley. Forty-two HvKINs were identified and screened from the barley genome, and a generated phylogenetic tree was used to compare the evolutionary relationships between Rice and Arabidopsis. The protein structure prediction, physicochemical properties, and bioinformatics of the HvKINs were also dissected. Our results reveal the important regulatory roles of HvKIN genes in barley growth. We found many cis- elements related to GA3 and ABA in homeopathic elements of the HvKIN gene and verified them by QRT-PCR, indicating their potential role in the barley kinesin family. The current study revealed the biological functions of barley kinesin genes in barley and will aid in further investigating the kinesin in other plant species.

The final piece of the Triangle of U: Evolution of the tetraploid Brassica carinata genome

Ethiopian mustard (Brassica carinata) is an ancient crop with remarkable stress resilience and a desirable seed fatty acid profile for biofuel uses. Brassica carinata is one of six Brassica species that share three major genomes from three diploid species (AA, BB, and CC) that spontaneously hybridized in a pairwise manner to form three allotetraploid species (AABB, AACC, and BBCC). Of the genomes of these species, that of B. carinata is the least understood. Here, we report a chromosome scale 1.31-Gbp genome assembly with 156.9-fold sequencing coverage for B. carinata, completing the reference genomes comprising the classic Triangle of U, a classical theory of the evolutionary relationships among these six species. Our assembly provides insights into the hybridization event that led to the current B. carinata genome and the genomic features that gave rise to the superior agronomic traits of B. carinata. Notably, we identified an expansion of transcription factor networks and agronomically important gene families. Completion of the Triangle of U comparative genomics platform has allowed us to examine the dynamics of polyploid evolution and the role of subgenome dominance in the domestication and continuing agronomic improvement of B. carinata and other Brassica species.

Nramp: Deprive and conquer?

Solute carriers 11 (Slc11) evolved from bacterial permease (MntH) to eukaryotic antibacterial defense (Nramp) while continuously mediating proton (H⁺)-dependent manganese (Mn²⁺) import. Also, Nramp horizontal gene transfer (HGT) toward bacteria led to mntH polyphyly. Prior demonstration that evolutionary rate-shifts distinguishing Slc11 from outgroup carriers dictate catalytic specificity suggested that resolving Slc11 family tree may provide a function-aware phylogenetic framework. Hence, MntH C (MC) subgroups resulted from HGTs of prototype Nramp (pNs) parologs while archetype Nramp (aNs) correlated with phagocytosis. PHI-Blast based taxonomic profiling confirmed MntH B phylogroup is confined to anaerobic bacteria vs. MntH A (MA)’s broad distribution; suggested niche-related spread of MC subgroups; established that MA-variant MH, which carries ‘eukaryotic signature’ marks, predominates in archaea. Slc11 phylogeny shows MH is sister to Nramp. Site-specific analysis of Slc11 charge network known to interact with the protonmotive force demonstrates sequential rate-shifts that recapitulate Slc11 evolution. 3D mapping of similarly coevolved sites across Slc11 hydrophobic core revealed successive targeting of discrete areas. The data imply that pN HGT could advantage recipient bacteria for H⁺-dependent Mn²⁺ acquisition and Alphafold 3D models suggest conformational divergence among MC subgroups. It is proposed that Slc11 originated as a bacterial stress resistance function allowing Mn²⁺-dependent persistence in conditions adverse for growth, and that archaeal MH could contribute to eukaryogenesis as a Mn²⁺ sequestering defense perhaps favoring intracellular growth-competent bacteria.

Multiple Pristionchus pacificus genomes reveal distinct evolutionary dynamics between de novo candidates and duplicated genes

The birth of new genes is a major molecular innovation driving phenotypic diversity across all domains of life. Although repurposing of existing protein-coding material by duplication is considered the main process of new gene formation, recent studies have discovered thousands of transcriptionally active sequences as a rich source of new genes. However, differential loss rates have to be assumed to reconcile the high birth rates of these incipient de novo genes with the dominance of ancient gene families in individual genomes. Here, we test this rapid turnover hypothesis in the context of the nematode model organism Pristionchus pacificus We extended the existing species-level phylogenomic framework by sequencing the genomes of six divergent P. pacificus strains. We used these data to study the evolutionary dynamics of different age classes and categories of origin at a population level. Contrasting de novo candidates with new families that arose by duplication and divergence from known genes, we find that de novo candidates are typically shorter, show less expression, and are overrepresented on the sex chromosome. Although the contribution of de novo candidates increases toward young age classes, multiple comparisons within the same age class showed significantly higher attrition in de novo candidates than in known genes. Similarly, young genes remain under weak evolutionary constraints with de novo candidates representing the fastest evolving subcategory. Altogether, this study provides empirical evidence for the rapid turnover hypothesis and highlights the importance of the evolutionary timescale when quantifying the contribution of different mechanisms toward new gene formation.

Evolution of binding preferences among whole-genome duplicated transcription factors

Throughout evolution, new transcription factors (TFs) emerge by gene duplication, promoting growth and rewiring of transcriptional networks. How TF duplicates diverge was studied in a few cases only. To provide a genome-scale view, we considered the set of budding yeast TFs classified as whole-genome duplication (WGD)-retained paralogs (~35% of all specific TFs). Using high-resolution profiling, we find that ~60% of paralogs evolved differential binding preferences. We show that this divergence results primarily from variations outside the DNA-binding domains (DBDs), while DBD preferences remain largely conserved. Analysis of non-WGD orthologs revealed uneven splitting of ancestral preferences between duplicates, and the preferential acquiring of new targets by the least conserved paralog (biased neo/sub-functionalization). Interactions between paralogs were rare, and, when present, occurred through weak competition for DNA-binding or dependency between dimer-forming paralogs. We discuss the implications of our findings for the evolutionary design of transcriptional networks.

Systematic identification and expression profiles of the BAHD superfamily acyltransferases in barley (Hordeum vulgare)

BAHD superfamily acyltransferases play an important role in catalyzing and regulating secondary metabolism in plants. Despite this, there is relatively little information regarding the BAHD superfamily in barley. In this study, we identified 116 HvBAHD acyltransferases from the barley genome. Based on phylogenetic analysis and classification in model monocotyledonous and dicotyledonous plants, we divided the genes into eight groups, I-a, I-b, II, III-a, III-b, IV, V-a and V-b. The Clade IV genes, including Agmatine Coumarol Transferase (ACT) that is associated with resistance of plants to Gibberella fungi, were absent in Arabidopsis. Cis-regulatory element analysis of the HvBAHDs showed that the genes respond positively to GA3 treatment. In-silico expression and qPCR analysis showed the HvBAHD genes are expressed in a range of tissues and developmental stages, and highly enriched in the seedling stage, consistent with diverse roles. Single nucleotide polymorphism (SNP) scanning analysis revealed that the natural variation in the coding regions of the HvBAHDs is low and the sequences have been conserved during barley domestication. Our results reveal the complexity of the HvBAHDs and will help facilitate their analysis in further studies.

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.

Evidence from Drosophila Supports Higher Duplicability of Faster Evolving Genes

The faster rate of evolution of duplicated genes relative to singletons has been well documented in multiple lineages. This observation has generally been attributed to a presumed release from constraint following creation of a redundant, duplicate copy. However, it is not obvious that the relationship operates in this direction. An alternative possibility—that the faster rate of evolution predates the duplication event and the observed differences result from a higher propensity to duplicate in fast-evolving genes—has been tested in primates and in insects. However, these studies arrived at different conclusions and clarity is needed on whether these contrasting results relate to differences in methodology or legitimate biological differences between the lineages selected. Here, we test whether duplicable genes are faster evolving independent of duplication in the Drosophila lineage and find that our results support the conclusion that faster evolving genes are more likely to duplicate, in agreement with previous work in primates. Our findings indicate that this characteristic of gene duplication is not restricted to a single lineage and has broad implications for the interpretation of the impact of gene duplication. We identify a subset of “singletons” which defy the general trends and appear to be faster evolving. Further investigation implicates homology detection failure and suggests that these may be duplicable genes with unidentifiable paralogs.

Comparable Number of Genes Having Experienced Positive Selection among Great Ape Species

Simple Summary

It is of great interest to quantify adaptive evolution in human lineage by studying genes under positive selection, since these genes could reveal insights into our own adaptive evolutionary history compared to our closely related species and often these genes are functionally important. We used the great apes as the subjects to detect gene-level adaptive evolution signals in all the great ape lineages and investigated the evolutionary patterns and functional relevance of these adaptive evolution signals. Even the differences in population size among these closely related great apes have resulted in differences in their ability to remove deleterious alleles and to adapt to changing environments, we found that they experienced comparable numbers of positive selection. Notably, we identified several genes that offer insights into great ape and human evolution. For example, SOD1, a gene associated with aging in humans, experienced positive selection in the common ancestor of the great ape and this positive selection may contribute to the aging evolution in great apes. Overall, an updated list of positively selected genes reported by this study not only informs us of adaptive evolution during great ape evolution, but is also helpful to the further study of non-human primate models for disease and other fields.

Abstract

Alleles that cause advantageous phenotypes with positive selection contribute to adaptive evolution. Investigations of positive selection in protein-coding genes rely on the accuracy of orthology, models, the quality of assemblies, and alignment. Here, based on the latest genome assemblies and gene annotations, we present a comparative analysis on positive selection in four great ape species and identify 211 high-confidence positively selected genes (PSGs). Even the differences in population size among these closely related great apes have resulted in differences in their ability to remove deleterious alleles and to adapt to changing environments, we found that they experienced comparable numbers of positive selection. We also uncovered that more than half of multigene families exhibited signals of positive selection, suggesting that imbalanced positive selection resulted in the functional divergence of duplicates. Moreover, at the expression level, although positive selection led to a more non-uniform pattern across tissues, the correlation between positive selection and expression patterns is diverse. Overall, this updated list of PSGs is of great significance for the further study of the phenotypic evolution in great apes.

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.

Simultaneous expression of UV and violet SWS1 opsins expands the visual palette in a group of freshwater snakes

Snakes are known to express a rod visual opsin and two cone opsins, only (SWS1, LWS), a reduced palette resulting from their supposedly fossorial origins. Dipsadid snakes in the genus Helicops are highly visual predators that successfully invaded freshwater habitats from ancestral terrestrial-only habitats. Here, we report the first case of multiple SWS1 visual pigments in a vertebrate, simultaneously expressed in different photoreceptors and conferring both UV and violet sensitivity to Helicops snakes. Molecular analysis and in vitro expression confirmed the presence of two functional SWS1 opsins, likely the result of recent gene duplication. Evolutionary analyses indicate that each sws1 variant has undergone different evolutionary paths with strong purifying selection acting on the UV-sensitive copy and d_N/d_S ∼1 on the violet-sensitive copy. Site-directed mutagenesis points to the functional role of a single amino acid substitution, Phe86Val, in the large spectral shift between UV and violet opsins. In addition, higher densities of photoreceptors and SWS1 cones in the ventral retina suggest improved acuity in the upper visual field possibly correlated with visually guided behaviors. The expanded visual opsin repertoire and specialized retinal architecture are likely to improve photon uptake in underwater and terrestrial environments, and provide the neural substrate for a gain in chromatic discrimination, potentially conferring unique color vision in the UV–violet range. Our findings highlight the innovative solutions undertaken by a highly specialized lineage to tackle the challenges imposed by the invasion of novel photic environments and the extraordinary diversity of evolutionary trajectories taken by visual opsin-based perception in vertebrates.

Olfactory receptor gene evolution is unusually rapid across Tetrapoda and outpaces chemosensory phenotypic change

Chemosensation is the most ubiquitous sense in animals, enacted by the products of complex gene families that detect environmental chemical cues and larger-scale sensory structures that process these cues. While there is a general conception that olfactory receptor (OR) genes evolve rapidly, the universality of this phenomenon across vertebrates, and its magnitude, are unclear. The supposed correlation between molecular rates of chemosensory evolution and phenotypic diversity of chemosensory systems is largely untested. We combine comparative genomics and sensory morphology to test whether OR genes and olfactory phenotypic traits evolve at faster rates than other genes or traits. Using published genomes, we identified ORs in 21 tetrapods, including amphibians, reptiles, birds, and mammals and compared their rates of evolution to those of orthologous non-OR protein-coding genes. We found that, for all clades investigated, most OR genes evolve nearly an order of magnitude faster than other protein-coding genes, with many OR genes showing signatures of diversifying selection across nearly all taxa in this study. This rapid rate of evolution suggests that chemoreceptor genes are in "evolutionary overdrive," perhaps evolving in response to the ever-changing chemical space of the environment. To obtain complementary morphological data, we stained whole fixed specimens with iodine, µCT-scanned the specimens, and digitally segmented chemosensory and nonchemosensory brain regions. We then estimated phenotypic variation within traits and among tetrapods. While we found considerable variation in chemosensory structures, they were no more diverse than nonchemosensory regions. We suggest chemoreceptor genes evolve quickly in reflection of an ever-changing chemical space, whereas chemosensory phenotypes and processing regions are more conserved because they use a standardized or constrained architecture to receive and process a range of chemical cues.

Special Care Is Needed in Applying Phylogenetic Comparative Methods to Gene Trees with Speciation and Duplication Nodes

How gene function evolves is a central question of evolutionary biology. It can be investigated by comparing functional genomics results between species and between genes. Most comparative studies of functional genomics have used pairwise comparisons. Yet it has been shown that this can provide biased results, as genes, like species, are phylogenetically related. Phylogenetic comparative methods should be used to correct for this, but they depend on strong assumptions, including unbiased tree estimates relative to the hypothesis being tested. Such methods have recently been used to test the “ortholog conjecture,” the hypothesis that functional evolution is faster in paralogs than in orthologs. Although pairwise comparisons of tissue specificity (τ) provided support for the ortholog conjecture, phylogenetic independent contrasts did not. Our reanalysis on the same gene trees identified problems with the time calibration of duplication nodes. We find that the gene trees used suffer from important biases, due to the inclusion of trees with no duplication nodes, to the relative age of speciations and duplications, to systematic differences in branch lengths, and to non-Brownian motion of tissue specificity on many trees. We find that incorrect implementation of phylogenetic method in empirical gene trees with duplications can be problematic. Controlling for biases allows successful use of phylogenetic methods to study the evolution of gene function and provides some support for the ortholog conjecture using three different phylogenetic approaches.

Rapid evolution of mammalian APLP1 as a synaptic adhesion molecule

Amyloid precursor protein (APP) family members are involved in essential neuronal development including neurite outgrowth, neuronal migration and maturation of synapse and neuromuscular junction. Among the APP gene family members, amyloid precursor-like protein 1 (APLP1) is selectively expressed in neurons and has specialized functions during synaptogenesis. Although a potential role for APLP1 in neuronal evolution has been indicated, its precise evolutionary and functional contributions are unknown. This study shows the molecular evolution of the vertebrate APP family based on phylogenetic analysis, while contrasting the evolutionary differences within the APP family. Phylogenetic analysis showed 15 times higher substitution rate that is driven by positive selection at the stem branch of the mammalian APLP1, resulting in dissimilar protein sequences compared to APP/APLP2. Docking simulation identified one positively selected site in APLP1 that alters the heparin-binding site, which could affect its function, and dimerization rate. Furthermore, the evolutionary rate covariation between the mammalian APP family and synaptic adhesion molecules (SAMs) was confirmed, indicating that only APLP1 has evolved to gain synaptic adhesion property. Overall, our results suggest that the enhanced synaptogenesis property of APLP1 as one of the SAMs may have played a role in mammalian brain evolution.

Selection following Gene Duplication Shapes Recent Genome Evolution in the Pea Aphid Acyrthosiphon pisum

Ecology of insects is as wide as their diversity, which reflects their high capacity of adaptation in most of the environments of our planet. Aphids, with over 4,000 species, have developed a series of adaptations including a high phenotypic plasticity and the ability to feed on the phloem sap of plants, which is enriched in sugars derived from photosynthesis. Recent analyses of aphid genomes have indicated a high level of shared ancestral gene duplications that might represent a basis for genetic innovation and broad adaptations. In addition, there are a large number of recent, species-specific gene duplications whose role in adaptation remains poorly understood. Here, we tested whether duplicates specific to the pea aphid Acyrthosiphon pisum are related to genomic innovation by combining comparative genomics, transcriptomics, and chromatin accessibility analyses. Consistent with large levels of neofunctionalization, we found that most of the recent pairs of gene duplicates evolved asymmetrically, showing divergent patterns of positive selection and gene expression. Genes under selection involved a plethora of biological functions, suggesting that neofunctionalization and tissue specificity, among other evolutionary mechanisms, have orchestrated the evolution of recent paralogs in the pea aphid and may have facilitated host-symbiont cooperation. Our comprehensive phylogenomics analysis allowed us to tackle the history of duplicated genes to pave the road toward understanding the role of gene duplication in ecological adaptation.

Divergence in Coding Sequence and Expression of Different Functional Categories of Immune Genes between Two Wild Rodent Species

Differences in immune function between species could be a result of interspecific divergence in coding sequence and/or expression of immune genes. Here, we investigate how the degree of divergence in coding sequence and expression differs between functional categories of immune genes, and if differences between categories occur independently of other factors (expression level, pleiotropy). To this end, we compared spleen transcriptomes of wild-caught yellow-necked mice and bank voles. Immune genes expressed in the spleen were divided into four categories depending on the function of the encoded protein: pattern recognition receptors (PRR); signal transduction proteins; transcription factors; and cyto- and chemokines and their receptors. Genes encoding PRR and cyto-/chemokines had higher sequence divergence than genes encoding signal transduction proteins and transcription factors, even when controlling for potentially confounding factors. Genes encoding PRR also had higher expression divergence than genes encoding signal transduction proteins and transcription factors. There was a positive correlation between expression divergence and coding sequence divergence, in particular for PRR genes. We propose that this is a result of that divergence in PRR coding sequence leads to divergence in PRR expression through positive feedback of PRR ligand binding on PRR expression. When controlling for sequence divergence, expression divergence of PRR genes did not differ from other categories. Taken together, the results indicate that coding sequence divergence of PRR genes is a major cause of differences in immune function between species.

Diverse Evolution in 111 Plant Genomes Reveals Purifying and Dosage Balancing Selection Models for F-Box Genes

The F-box proteins function as substrate receptors to determine the specificity of Skp1-Cul1-F-box ubiquitin ligases. Genomic studies revealed large and diverse sizes of the F-box gene superfamily across plant species. Our previous studies suggested that the plant F-box gene superfamily is under genomic drift evolution promoted by epigenomic programming. However, how the size of the superfamily drifts across plant genomes is currently unknown. Through a large-scale genomic and phylogenetic comparison of the F-box gene superfamily covering 110 green plants and one red algal species, I discovered four distinct groups of plant F-box genes with diverse evolutionary processes. While the members in Clusters 1 and 2 are species/lineage-specific, those in Clusters 3 and 4 are present in over 46 plant genomes. Statistical modeling suggests that F-box genes from the former two groups are skewed toward fewer species and more paralogs compared to those of the latter two groups whose presence frequency and sizes in plant genomes follow a random statistical model. The enrichment of known Arabidopsis F-box genes in Clusters 3 and 4, along with comprehensive biochemical evidence showing that Arabidopsis members in Cluster 4 interact with the Arabidopsis Skp1-like 1 (ASK1), demonstrates over-representation of active F-box genes in these two groups. Collectively, I propose purifying and dosage balancing selection models to explain the lineage/species-specific duplications and expansions of F-box genes in plant genomes. The purifying selection model suggests that most, if not all, lineage/species-specific F-box genes are detrimental and are thus kept at low frequencies in plant genomes.

Dynamic Evolution of the Cthrc1 Genes, a Newly Defined Collagen-Like Family

Collagen triple helix repeat containing protein 1 (Cthrc1) is a secreted glycoprotein reported to regulate collagen deposition and to be linked to the Transforming growth factor β/Bone morphogenetic protein and the Wnt/planar cell polarity pathways. It was first identified as being induced upon injury to rat arteries and was found to be highly expressed in multiple human cancer types. Here, we explore the phylogenetic and evolutionary trends of this metazoan gene family, previously studied only in vertebrates. We identify Cthrc1 orthologs in two distant cnidarian species, the sea anemone Nematostella vectensis and the hydrozoan Clytia hemisphaerica, both of which harbor multiple copies of this gene. We find that Cthrc1 clade-specific diversification occurred multiple times in cnidarians as well as in most metazoan clades where we detected this gene. Many other groups, such as arthropods and nematodes, have entirely lost this gene family. Most vertebrates display a single highly conserved gene, and we show that the sequence evolutionary rate of Cthrc1 drastically decreased within the gnathostome lineage. Interestingly, this reduction coincided with the origin of its conserved upstream neighboring gene, Frizzled 6 (FZD6), which in mice has been shown to functionally interact with Cthrc1. Structural modeling methods further reveal that the yet uncharacterized C-terminal domain of Cthrc1 is similar in structure to the globular C1q superfamily domain, also found in the C-termini of collagens VIII and X. Thus, our studies show that the Cthrc1 genes are a collagen-like family with a variable short collagen triple helix domain and a highly conserved C-terminal domain structure resembling the C1q family.

Richard Dickerson, Molecular Clocks, and Rates of Protein Evolution

Proteins approximately behave as molecular clocks, accumulating amino acid replacements at a more or less constant rate. Nonetheless, each protein displays a characteristic rate of evolution: whereas some proteins remain largely unaltered over large periods of time, others can rapidly accumulate amino acid replacements. An article by Richard Dickerson, published in the first issue of the Journal of Molecular Evolution (J Mol Evol 1:26-45, 1971), described the first analysis in which the rates of evolution of many proteins were compared, and the differences were interpreted in the light of their function. When comparing the sequences of fibrinopeptides, hemoglobin, and cytochrome c of different species, he observed a linear relationship between the number of amino acid replacements and divergence time. Remarkably, fibrinopeptides had evolved fast, cytochrome c had evolved slowly, and hemoglobin exhibited an intermediate rate of evolution. As the Journal of Molecular Evolution celebrates its 50th anniversary, I highlight this landmark article and reflect on its impact on the field of Molecular Evolution.

Dense sampling of bird diversity increases power of comparative genomics

Whole-genome sequencing projects are increasingly populating the tree of life and characterizing biodiversity1-4. Sparse taxon sampling has previously been proposed to confound phylogenetic inference5, and captures only a fraction of the genomic diversity. Here we report a substantial step towards the dense representation of avian phylogenetic and molecular diversity, by analysing 363 genomes from 92.4% of bird families-including 267 newly sequenced genomes produced for phase II of the Bird 10,000 Genomes (B10K) Project. We use this comparative genome dataset in combination with a pipeline that leverages a reference-free whole-genome alignment to identify orthologous regions in greater numbers than has previously been possible and to recognize genomic novelties in particular bird lineages. The densely sampled alignment provides a single-base-pair map of selection, has more than doubled the fraction of bases that are confidently predicted to be under conservation and reveals extensive patterns of weak selection in predominantly non-coding DNA. Our results demonstrate that increasing the diversity of genomes used in comparative studies can reveal more shared and lineage-specific variation, and improve the investigation of genomic characteristics. We anticipate that this genomic resource will offer new perspectives on evolutionary processes in cross-species comparative analyses and assist in efforts to conserve species.

Genetic, transcriptional, and regulatory landscape of monolignol biosynthesis pathway in Miscanthus × giganteus

BACKGROUND

Miscanthus × giganteus is widely recognized as a promising lignocellulosic biomass crop due to its advantages of high biomass production, low environmental impacts, and the potential to be cultivated on marginal land. However, the high costs of bioethanol production still limit the current commercialization of lignocellulosic bioethanol. The lignin in the cell wall and its by-products released in the pretreatment step is the main component inhibiting the enzymatic reactions in the saccharification and fermentation processes. Hence, genetic modification of the genes involved in lignin biosynthesis could be a feasible strategy to overcome this barrier by manipulating the lignin content and composition of M. × giganteus. For this purpose, the essential knowledge of these genes and understanding the underlying regulatory mechanisms in M. × giganteus is required.

RESULTS

In this study, MgPAL1, MgPAL5, Mg4CL1, Mg4CL3, MgHCT1, MgHCT2, MgC3'H1, MgCCoAOMT1, MgCCoAOMT3, MgCCR1, MgCCR2, MgF5H, MgCOMT, and MgCAD were identified as the major monolignol biosynthetic genes in M. × giganteus based on genetic and transcriptional evidence. Among them, 12 genes were cloned and sequenced. By combining transcription factor binding site prediction and expression correlation analysis, MYB46, MYB61, MYB63, WRKY24, WRKY35, WRKY12, ERF021, ERF058, and ERF017 were inferred to regulate the expression of these genes directly. On the basis of these results, an integrated model was summarized to depict the monolignol biosynthesis pathway and the underlying regulatory mechanism in M. × giganteus.

CONCLUSIONS

This study provides a list of potential gene targets for genetic improvement of lignocellulosic biomass quality of M. × giganteus, and reveals the genetic, transcriptional, and regulatory landscape of the monolignol biosynthesis pathway in M. × giganteus.

Recurrent sequence evolution after independent gene duplication

Background

Convergent and parallel evolution provide unique insights into the mechanisms of natural selection. Some of the most striking convergent and parallel (collectively recurrent) amino acid substitutions in proteins are adaptive, but there are also many that are selectively neutral. Accordingly, genome-wide assessment has shown that recurrent sequence evolution in orthologs is chiefly explained by nearly neutral evolution. For paralogs, more frequent functional change is expected because additional copies are generally not retained if they do not acquire their own niche. Yet, it is unknown to what extent recurrent sequence differentiation is discernible after independent gene duplications in different eukaryotic taxa.

Results

We develop a framework that detects patterns of recurrent sequence evolution in duplicated genes. This is used to analyze the genomes of 90 diverse eukaryotes. We find a remarkable number of families with a potentially predictable functional differentiation following gene duplication. In some protein families, more than ten independent duplications show a similar sequence-level differentiation between paralogs. Based on further analysis, the sequence divergence is found to be generally asymmetric. Moreover, about 6% of the recurrent sequence evolution between paralog pairs can be attributed to recurrent differentiation of subcellular localization. Finally, we reveal the specific recurrent patterns for the gene families Hint1/Hint2, Sco1/Sco2 and vma11/vma3.

Conclusions

The presented methodology provides a means to study the biochemical underpinning of functional differentiation between paralogs. For instance, two abundantly repeated substitutions are identified between independently derived Sco1 and Sco2 paralogs. Such identified substitutions allow direct experimental testing of the biological role of these residues for the repeated functional differentiation. We also uncover a diverse set of families with recurrent sequence evolution and reveal trends in the functional and evolutionary trajectories of this hitherto understudied phenomenon.

Divergent Traits and Ligand-Binding Properties of the Cytomegalovirus CD48 Gene Family

The genesis of gene families by the capture of host genes and their subsequent duplication is a crucial process in the evolution of large DNA viruses. CD48 is a cell surface molecule that interacts via its N-terminal immunoglobulin (Ig) domain with the cell surface receptor 2B4 (CD244), regulating leukocyte cytotoxicity. We previously reported the presence of five CD48 homologs (vCD48s) in two related cytomegaloviruses, and demonstrated that one of them, A43, binds 2B4 and acts as a soluble CD48 decoy receptor impairing NK cell function. Here, we have characterized the rest of these vCD48s. We show that they are highly glycosylated proteins that display remarkably distinct features: divergent biochemical properties, cellular locations, and temporal expression kinetics. In contrast to A43, none of them interacts with 2B4. Consistent with this, molecular modeling of the N-terminal Ig domains of these vCD48s evidences notable changes as compared to CD48, suggesting that they interact with alternative targets. Accordingly, we demonstrate that one of them, S30, tightly binds CD2, a crucial T- and NK-cell adhesion and costimulatory molecule. Thus, our findings show how a key host immune receptor gene captured by a virus can be subsequently remodeled to evolve new immunoevasins with altered binding properties.

Duplication and functional divergence of a calcium sensor in the Brassicaceae

The presence of varied numbers of CALCINEURIN B-LIKE10 (CBL10) calcium sensor genes in species across the Brassicaceae and the demonstrated role of CBL10 in salt tolerance in Arabidopsis thaliana and Eutrema salsugineum provided a unique opportunity to determine if CBL10 function is modified in different species and linked to salt tolerance. Salinity effects on species growth and cross-species complementation were used to determine the extent of conservation and divergence of CBL10 function in four species representing major lineages within the core Brassicaceae (A. thaliana, E. salsugineum, Schrenkiella parvula, and Sisymbrium irio) as well as the first diverging lineage (Aethionema arabicum). Evolutionary and functional analyses indicate that CBL10 duplicated within expanded lineage II of the Brassicaceae and that, while portions of CBL10 function are conserved across the family, there are species-specific variations in CBL10 function. Paralogous CBL10 genes within a species diverged in expression and function probably contributing to the maintenance of the duplicated gene pairs. Orthologous CBL10 genes diverged in function in a species-specific manner, suggesting that functions arose post-speciation. Multiple CBL10 genes and their functional divergence may have expanded calcium-mediated signaling responses and contributed to the ability of certain members of the Brassicaceae to maintain growth in salt-affected soils.

Molecular Chaperones Accelerate the Evolution of Their Protein Clients in Yeast

Protein stability is a major constraint on protein evolution. Molecular chaperones, also known as heat-shock proteins, can relax this constraint and promote protein evolution by diminishing the deleterious effect of mutations on protein stability and folding. This effect, however, has only been stablished for a few chaperones. Here, we use a comprehensive chaperone–protein interaction network to study the effect of all yeast chaperones on the evolution of their protein substrates, that is, their clients. In particular, we analyze how yeast chaperones affect the evolutionary rates of their clients at two very different evolutionary time scales. We first study the effect of chaperone-mediated folding on protein evolution over the evolutionary divergence of Saccharomyces cerevisiae and S. paradoxus. We then test whether yeast chaperones have left a similar signature on the patterns of standing genetic variation found in modern wild and domesticated strains of S. cerevisiae. We find that genes encoding chaperone clients have diverged faster than genes encoding non-client proteins when controlling for their number of protein–protein interactions. We also find that genes encoding client proteins have accumulated more intraspecific genetic diversity than those encoding non-client proteins. In a number of multivariate analyses, controlling by other well-known factors that affect protein evolution, we find that chaperone dependence explains the largest fraction of the observed variance in the rate of evolution at both evolutionary time scales. Chaperones affecting rates of protein evolution mostly belong to two major chaperone families: Hsp70s and Hsp90s. Our analyses show that protein chaperones, by virtue of their ability to buffer destabilizing mutations and their role in modulating protein genotype–phenotype maps, have a considerable accelerating effect on protein evolution.

Comparative Genomic Analysis of the Pheromone Receptor Class 1 Family (V1R) Reveals Extreme Complexity in Mouse Lemurs (Genus, Microcebus) and a Chromosomal Hotspot across Mammals

Sensory gene families are of special interest for both what they can tell us about molecular evolution and what they imply as mediators of social communication. The vomeronasal type-1 receptors (V1Rs) have often been hypothesized as playing a fundamental role in driving or maintaining species boundaries given their likely function as mediators of intraspecific mate choice, particularly in nocturnal mammals. Here, we employ a comparative genomic approach for revealing patterns of V1R evolution within primates, with a special focus on the small-bodied nocturnal mouse and dwarf lemurs of Madagascar (genera Microcebus and Cheirogaleus, respectively). By doubling the existing genomic resources for strepsirrhine primates (i.e. the lemurs and lorises), we find that the highly speciose and morphologically cryptic mouse lemurs have experienced an elaborate proliferation of V1Rs that we argue is functionally related to their capacity for rapid lineage diversification. Contrary to a previous study that found equivalent degrees of V1R diversity in diurnal and nocturnal lemurs, our study finds a strong correlation between nocturnality and V1R elaboration, with nocturnal lemurs showing elaborate V1R repertoires and diurnal lemurs showing less diverse repertoires. Recognized subfamilies among V1Rs show unique signatures of diversifying positive selection, as might be expected if they have each evolved to respond to specific stimuli. Furthermore, a detailed syntenic comparison of mouse lemurs with mouse (genus Mus) and other mammalian outgroups shows that orthologous mammalian subfamilies, predicted to be of ancient origin, tend to cluster in a densely populated region across syntenic chromosomes that we refer to as a V1R "hotspot."

Omics Technologies for Profiling Toxin Diversity and Evolution in Snake Venom: Impacts on the Discovery of Therapeutic and Diagnostic Agents

Snake venoms are primarily composed of proteins and peptides, and these toxins have developed high selectivity to their biological targets. This makes venoms interesting for exploration into protein evolution and structure-function relationships. A single venom protein superfamily can exhibit a variety of pharmacological effects; these variations in activity originate from differences in functional sites, domains, posttranslational modifications, and the formations of toxin complexes. In this review, we discuss examples of how the major venom protein superfamilies have diversified, as well as how newer technologies in the omics fields, such as genomics, transcriptomics, and proteomics, can be used to characterize both known and unknown toxins.Because toxins are bioactive molecules with a rich diversity of activities, they can be useful as therapeutic and diagnostic agents, and successful examples of toxin applications in these areas are also reviewed. With the current rapid pace of technology, snake venom research and its applications will only continue to expand.

Greater genetic and regulatory plasticity of retained duplicates in Epichloë endophytic fungi

Gene duplicates can act as a source of genetic material from which new functions arise. Most duplicated genes revert to single copy genes and only a small proportion are retained. However, it remains unclear why some duplicate genes persist in the genome for an extended time. We investigate this question by analysing retained gene duplicates in the fungal genus Epichloë, ascomycete fungi that form close endophytic symbioses with their host grasses. Retained duplicates within this genus have two independent origins, but both long pre-date the origin and diversification of the genus Epichloë. We find that loss of retained duplicates within the genus is frequent and often associated with speciation. Retained duplicates have faster evolutionary rates (Ka) and show relaxed selection (Ka/Ks) compared to single copy genes. Both features are time-dependent. Through comparison of conspecific strains, we find greater evolutionary rates in coding regions and sequence divergence in regulatory regions of retained duplicates than single copy genes, with this pattern more pronounced for strains adapted to different grass host species. Consistent with this sequence divergence in regulatory regions, transcriptome analyses show greater expression variation of retained duplicates than single copy genes. This suggest that cis-regulatory changes make important contributions to the expression patterns of retained duplicates. Coupled with supporting observations from the model yeast Saccharomyces cerevisiae, these data suggest that genetic robustness and regulatory plasticity are common drivers behind the retention of duplicated genes in fungi.

Unraveling the evolutionary origin of ELR motif using fish CXC chemokine CXCL8

Chemokines are chemotactic proteins involved in host defense through the migration of immune-regulatory cells to the site of infection. Interleukin-8 (CXCL8/IL8) is the most studied "ELR-CXC chemokine/neutrophil activating chemokine (NAC) that regulate neutrophil trafficking during infections and inflammation by binding to its cognate G-protein coupled receptors CXCR1/CXCR2. The "ELR" motif of NAC chemokines is essential for the CXCR1/CXCR2 receptor activation. In order to understand the evolutionary origin of "ELR" motif in the CXC chemokines, a thorough evolutionary study of CXCL8 gene from various fishes and primates was performed. Phylogenetic analysis revealed that the CXCL8 gene can be classified into four distinct lineages (CXCL8-L1a, CXCL8-L1b, CXCL8-L2, and CXCL8-L3), where CXCL8-L1a is the fastest evolving lineage and CXCL8-L3 is the slowest. Selection analysis suggested that The "ELR/DLR" motif containing branches (gadoid and coelacanth) are positively selected. The probable evolutionary trend of "ELR" motif suggested that this motif in ancestor CXCL8 is evolved from the GGR of Lamprey (Agnatha), followed by duplication giving rise to two main motifs in CXCL8 "NXH" in L3 lineage and "ELR/DLR" in L1a/L1b lineages. Although, structural analysis suggested that the overall topology of the CXCL8 proteins is similar, differences do exist at the individual structural elements among the members of different lineages. Functional distance analysis suggested that the CXCL8-L3 lineage is more distant compared to the CXCL8-L1a and L1b lineages from the inferred ancestor. Functional divergence analysis between different lineages suggested that most of the selected residues are important for receptor or glycosaminoglycan binding. Such a functional diversification can be attributed to the novel set of functions adopted by CXCL8 in various species.

Chromatin Evolution-Key Innovations Underpinning Morphological Complexity

The history of life consists of a series of major evolutionary transitions, including emergence and radiation of complex multicellular eukaryotes from unicellular ancestors. The cells of multicellular organisms, with few exceptions, contain the same genome, however, their organs are composed of a variety of cell types that differ in both structure and function. This variation is largely due to the transcriptional activity of different sets of genes in different cell types. This indicates that complex transcriptional regulation played a key role in the evolution of complexity in eukaryotes. In this review, we summarize how gene duplication and subsequent evolutionary innovations, including the structural evolution of nucleosomes and chromatin-related factors, contributed to the complexity of the transcriptional system and provided a basis for morphological diversity.

Role for Cela1 in Postnatal Lung Remodeling and Alpha-1 Antitrypsin-Deficient Emphysema

Alpha-1 antitrypsin (AAT) deficiency-related emphysema is the fourth leading indication for lung transplant. Chymotrypsin-like elastase 1 (Cela1) is a digestive protease that is expressed during lung development in association with regions of elastin remodeling, exhibits stretch-dependent expression during lung regeneration, and binds lung elastin in a stretch-dependent manner. AAT covalently neutralizes Cela1 in vitro. We sought to determine the role of Cela1 in postnatal lung physiology, whether it interacted with AAT in vivo, and to detect any effects it may have in the context of AAT deficiency. The lungs of Cela1^-/- mice had aberrant lung elastin structure and higher elastance as assessed with the flexiVent system. On the basis of in situ zymography with ex vivo lung stretch, Cela1 was solely responsible for stretch-inducible lung elastase activity. By mass spectrometry, Cela1 degraded mature elastin similarly to pancreatic elastase. Cela1 promoter and protein sequences were phylogenetically distinct in the placental mammal lineage, suggesting an adaptive role for lung-expressed Cela1 in this clade. A 6-week antisense oligonucleotide mouse model of AAT deficiency resulted in emphysema with increased Cela1 mRNA and reduction of approximately 70 kD Cela1, consistent with covalent binding of Cela1 by AAT. Cela1^-/- mice were completely protected against emphysema in this model. Cela1 was increased in human AAT-deficient emphysema. Cela1 is important in physiologic and pathologic stretch-dependent remodeling processes in the postnatal lung. AAT is an important regulator of this process. Our findings provide proof of concept for the development of anti-Cela1 therapies to prevent and/or treat AAT-deficient emphysema.

Faster Evolving Primate Genes Are More Likely to Duplicate

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

Transcriptomic analyses reveal groups of co-expressed, syntenic lncRNAs in four species of the genus Caenorhabditis

Long non-coding RNAs (lncRNAs) are a heterogeneous class of genes that do not code for proteins. Since lncRNAs (or a fraction thereof) are expected to be functional, many efforts have been dedicated to catalog lncRNAs in numerous organisms, but our knowledge of lncRNAs in non vertebrate species remains very limited. Here, we annotated lncRNAs using transcriptomic data from the same larval stage of four Caenorhabditis species. The number of annotated lncRNAs in self-fertile nematodes was lower than in out-crossing species. We used a combination of approaches to identify putatively homologous lncRNAs: synteny, sequence conservation, and structural conservation. We classified a total of 1,532 out of 7,635 genes from the four species into families of lncRNAs with conserved synteny and expression at the larval stage, suggesting that a large fraction of the predicted lncRNAs may be species specific. Despite both sequence and local secondary structure seem to be poorly conserved, sequences within families frequently shared BLASTn hits and short sequence motifs, which were more likely to be unpaired in the predicted structures. We provide the first multi-species catalog of lncRNAs in nematodes and identify groups of lncRNAs with conserved synteny and expression, that share exposed motifs.

Angiosperm-Wide and Family-Level Analyses of AP2/ERF Genes Reveal Differential Retention and Sequence Divergence After Whole-Genome Duplication

Plants are immobile and often face stressful environmental conditions, prompting the evolution of genes regulating environmental responses. Such evolution is achieved largely through gene duplication and subsequent divergence. One of the most important gene families involved in regulating plant environmental responses and development is the AP2/ERF superfamily; however, the evolutionary history of these genes is unclear across angiosperms and in major angiosperm families adapted to various ecological niches. Specifically, the impact on gene copy number of whole-genome duplication events occurring around the time of the origins of several plant families is unknown. Here, we present the first angiosperm-wide comparative study of AP2/ERF genes, identifying 75 Angiosperm OrthoGroups (AOGs), each derived from an ancestral angiosperm gene copy. Among these AOGs, 21 retain duplicates with increased copy number in many angiosperm lineages, while the remaining 54 AOGs tend to maintain low copy number. Further analyses of multiple species in the Brassicaceae family indicated that family-specific duplicates experienced differential selective pressures in coding regions, with some paralogs showing signs of positive selection. Further, cis regulatory elements also exhibit extensive divergence between duplicates in Arabidopsis. Moreover, comparison of expression levels suggested that AP2/ERF genes with frequently retained duplicates are enriched for broad expression patterns, offering increased opportunities for functional diversification via changes in expression patterns, and providing a mechanism for repeated duplicate retention in some AOGs. Our results represent the most comprehensive evolutionary history of the AP2/ERF gene family, and support the hypothesis that AP2/ERF genes with broader expression patterns are more likely to be retained as duplicates than those with narrower expression profiles, which could lead to a higher chance of duplicate gene subfunctionalization. The greater tendency of some AOGs to retain duplicates, allowing expression and functional divergence, may facilitate the evolution of complex signaling networks in response to new environmental conditions.