Weak Purifying Selection Research Articles

Reversions mask the contribution of adaptive evolution in microbiomes.

When examining bacterial genomes for evidence of past selection, the results obtained depend heavily on the mutational distance between chosen genomes. Even within a bacterial species, genomes separated by larger mutational distances exhibit stronger evidence of purifying selection as assessed by , the normalized ratio of nonsynonymous to synonymous mutations. Here, we show that the classical interpretation of this scale-dependence, weak purifying selection, leads to problematic mutation accumulation when applied to available gut microbiome data. We propose an alternative, adaptive reversion model with exactly opposite implications for dynamical intuition and applications of . Reversions that occur and sweep within-host populations are nearly guaranteed in microbiomes due to large population sizes, short generation times, and variable environments. Using analytical and simulation approaches, we show that adaptive reversion can explain the decay given only dozens of locally-fluctuating selective pressures, which is realistic in the context of Bacteroides genomes. The success of the adaptive reversion model argues for interpreting low values of obtained from long-time scales with caution, as they may emerge even when adaptive sweeps are frequent. Our work thus inverts the interpretation of an old observation in bacterial evolution, illustrates the potential of mutational reversions to shape genomic landscapes over time, and highlights the importance of studying bacterial genomic evolution on short time scales.

Weak purifying selection in allelic diversity of the ADSL gene in indigenous and local chicken breeds and red junglefowl in Thailand

High levels of purine and uric acid, which are associated with health issues such as gout and cardiovascular disease, are found in the meat of fast-growing broiler chickens, which raises concerns about the quality of chicken meat and the health of the consumers who consume it. High genetic homogeneity and uniformity, particularly in genes involved in the synthesis of inosine monophosphate (IMP) and subsequent process of purine synthesis, which are associated with the meat quality, are exhibited in commercial broiler chickens owing to intensive inbreeding programs. Adenosine succinate lyase (ADSL) is a key enzyme involved in de novo purine biosynthetic pathway and its genetic polymorphisms affect IMP metabolism and purine content. In this study, we investigated the polymorphism of the ADSL gene in indigenous and local chicken breeds and red junglefowl in Thailand, using metabarcoding and genetic diversity analyses. Five alleles with 73 single nucleotide polymorphisms in exon 2, including missense and silent mutations, which may act on the synthesis efficiency of IMP and purine. Their protein structures revealed changes in amino acid composition that may affect ADSL enzyme activity. Weak purifying selection in these ADSL alleles was observed in the chicken population studied, implying that the variants have minor fitness impacts and a greater probability of fixation of beneficial mutations than strong purifying selection. A potential selective sweep was observed in Mae Hong Son chickens, whose purine content was lower than that in other breeds. This suggests a potential correlation between variations of the ADSL gene and reduced purine content and an impact of ADSL expression on the quality of chicken meat. However, further studies are required to validate its potential availability as a genetic marker for selecting useful traits that are beneficial to human health and well-being.

Pan-genome insights into adaptive evolution of bacterial symbionts in mixed host-microbe symbioses represented by human gut microbiota Bacteroides cellulosilyticus

Animal hosts harbor diverse assemblages of microbial symbionts that play crucial roles in the host's lifestyle. The link between microbial symbiosis and host development remains poorly understood. In particular, little is known about the adaptive evolution of gut bacteria in host–microbe symbioses. Recently, symbiotic relationships have been categorized as open, closed, or mixed, reflecting their modes of inter-host transmission and resulting in distinct genomic features. Members of the genus Bacteroides are the most abundant human gut microbiota and possess both probiotic and pathogenic potential, providing an excellent model for studying pan-genome evolution in symbiotic systems. Here, we determined the complete genome of an novel clinical strain PL2022, which was isolated from a blood sample and performed pan-genome analyses on a representative set of Bacteroides cellulosilyticus strains to quantify the influence of the symbiotic relationship on the evolutionary dynamics. B. cellulosilyticus exhibited correlated genomic features with both open and closed symbioses, suggesting a mixed symbiosis. An open pan-genome is characterized by abundant accessory gene families, potential horizontal gene transfer (HGT), and diverse mobile genetic elements (MGEs), indicating an innovative gene pool, mainly associated with genomic islands and plasmids. However, massive parallel gene loss, weak purifying selection, and accumulation of positively selected mutations were the main drivers of genome reduction in B. cellulosilyticus. Metagenomic read recruitment analyses showed that B. cellulosilyticus members are globally distributed and active in human gut habitats, in line with predominant vertical transmission in the human gut. However, existence and/or high abundance were also detected in non-intestinal tissues, other animal hosts, and non-host environments, indicating occasional horizontal transmission to new niches, thereby creating arenas for the acquisition of novel genes. This case study of adaptive evolution under a mixed host–microbe symbiosis advances our understanding of symbiotic pan-genome evolution. Our results highlight the complexity of genetic evolution in this unusual intestinal symbiont.

Wolbachia-mediated resistance to Zika virus infection in Aedes aegypti is dominated by diverse transcriptional regulation and weak evolutionary pressures.

A promising candidate for arbovirus control and prevention relies on replacing arbovirus-susceptible Aedes aegypti populations with mosquitoes that have been colonized by the intracellular bacterium Wolbachia and thus have a reduced capacity to transmit arboviruses. This reduced capacity to transmit arboviruses is mediated through a phenomenon referred to as pathogen blocking. Pathogen blocking has primarily been proposed as a tool to control dengue virus (DENV) transmission, however it works against a range of viruses, including Zika virus (ZIKV). Despite years of research, the molecular mechanisms underlying pathogen blocking still need to be better understood. Here, we used RNA-seq to characterize mosquito gene transcription dynamics in Ae. aegypti infected with the wMel strain of Wolbachia that are being released by the World Mosquito Program in Medellín, Colombia. Comparative analyses using ZIKV-infected, uninfected tissues, and mosquitoes without Wolbachia revealed that the influence of wMel on mosquito gene transcription is multifactorial. Importantly, because Wolbachia limits, but does not completely prevent, replication of ZIKV and other viruses in coinfected mosquitoes, there is a possibility that these viruses could evolve resistance to pathogen blocking. Therefore, to understand the influence of Wolbachia on within-host ZIKV evolution, we characterized the genetic diversity of molecularly barcoded ZIKV virus populations replicating in Wolbachia-infected mosquitoes and found that within-host ZIKV evolution was subject to weak purifying selection and, unexpectedly, loose anatomical bottlenecks in the presence and absence of Wolbachia. Together, these findings suggest that there is no clear transcriptional profile associated with Wolbachia-mediated ZIKV restriction, and that there is no evidence for ZIKV escape from this restriction in our system.

Genetic Drift Shapes the Evolution of a Highly Dynamic Metapopulation.

The dynamics of extinction and (re)colonization in habitat patches are characterizing features of dynamic metapopulations, causing them to evolve differently than large, stable populations. The propagule model, which assumes genetic bottlenecks during colonization, posits that newly founded subpopulations have low genetic diversity and are genetically highly differentiated from each other. Immigration may then increase diversity and decrease differentiation between subpopulations. Thus, older and/or less isolated subpopulations are expected to have higher genetic diversity and less genetic differentiation. We tested this theory using whole-genome pool-sequencing to characterize nucleotide diversity and differentiation in 60 subpopulations of a natural metapopulation of the cyclical parthenogen Daphnia magna. For comparison, we characterized diversity in a single, large, and stable D. magna population. We found reduced (synonymous) genomic diversity, a proxy for effective population size, weak purifying selection, and low rates of adaptive evolution in the metapopulation compared with the large, stable population. These differences suggest that genetic bottlenecks during colonization reduce effective population sizes, which leads to strong genetic drift and reduced selection efficacy in the metapopulation. Consistent with the propagule model, we found lower diversity and increased differentiation in younger and also in more isolated subpopulations. Our study sheds light on the genomic consequences of extinction-(re)colonization dynamics to an unprecedented degree, giving strong support for the propagule model. We demonstrate that the metapopulation evolves differently from a large, stable population and that evolution is largely driven by genetic drift.

The western redcedar genome reveals low genetic diversity in a self-compatible conifer.

We assembled the 9.8-Gbp genome of western redcedar (WRC; Thuja plicata), an ecologically and economically important conifer species of the Cupressaceae. The genome assembly, derived from a uniquely inbred tree produced through five generations of self-fertilization (selfing), was determined to be 86% complete by BUSCO analysis, one of the most complete genome assemblies for a conifer. Population genomic analysis revealed WRC to be one of the most genetically depauperate wild plant species, with an effective population size of approximately 300 and no significant genetic differentiation across its geographic range. Nucleotide diversity, π, is low for a continuous tree species, with many loci showing zero diversity, and the ratio of π at zero- to fourfold degenerate sites is relatively high (approximately 0.33), suggestive of weak purifying selection. Using an array of genetic lines derived from up to five generations of selfing, we explored the relationship between genetic diversity and mating system. Although overall heterozygosity was found to decline faster than expected during selfing, heterozygosity persisted at many loci, and nearly 100 loci were found to deviate from expectations of genetic drift, suggestive of associative overdominance. Nonreference alleles at such loci often harbor deleterious mutations and are rare in natural populations, implying that balanced polymorphisms are maintained by linkage to dominant beneficial alleles. This may account for how WRC remains responsive to natural and artificial selection, despite low genetic diversity.

Evolution and functionalization of vitellogenin genes in the termite Reticulitermes speratus.

Eusociality has been commonly observed in distinct animal lineages. The reproductive division of labor is a particular feature, achieved by the coordination between fertile and sterile castes within the same nest. The sociogenomic approach in social hymenopteran insects indicates that vitellogenin (Vg) has undergone neo-functionalization in sterile castes. Here, to know whether Vgs have distinct roles in nonreproductive castes in termites, we investigated the unique characteristics of Vgs in the rhinotermitid termite Reticulitermes speratus. The four Vgs were identified from R. speratus (RsVg1-4), and RsVg3 sequences were newly identified using the RACE method. Molecular phylogenetic analysis supported the monophyly of the four termite Vgs. Moreover, the termites Vg1-3 and Vg4 were positioned in two different clades. The dN/dSratios indicated that the branch leading to the common ancestor of termite Vg4 was under weak purifying selection. Expression analyses among castes (reproductives, workers, and soldiers) and females (nymphs, winged alates, and queens) showed that RsVg1-3 was highly expressed in fertile queens. In contrast, RsVg4 was highly expressed in workers and female nonreproductives (nymphs and winged adults). Localization of RsVg4messenger RNAwas confirmed in the fat body of worker heads and abdomens. These results suggest that Vg genes are functionalized after gene duplication during termite eusocial transition and that Vg4 is involved in nonreproductive roles in termites.

Excessive replacement changes drive evolution of global sheep prion protein (PRNP) sequences.

Sheep prion protein (PRNP) is the major host genetic factor responsible for susceptibility to scrapie. We aimed to understand the evolutionary history of sheep PRNP, and primarily focused on breeds from Turkey and Ethiopia, representing genome-wise ancient sheep populations. Population molecular genetic analyses are extended to European, South Asian, and East Asian populations, and for the first time to scrapie associated haplotypes. 1178 PRNP coding region nucleotide sequences were analyzed. High levels of nucleotide diversity driven by extensive low-frequency replacement changes are observed in all populations. Interspecific analyses were conducted using mouflon and domestic goat as outgroup species. Despite an abundance of silent and replacement changes, lack of silent or replacement fixations was observed. All scrapie-associated haplotype analyses from all populations also showed extensive low-frequency replacement changes. Neutrality tests did not indicate positive (directional), balancing or strong negative selection or population contraction for any of the haplotypes in any population. A simple negative selection history driven by prion disease susceptibility is not supported by the population and haplotype based analyses. Molecular function, biological process enrichment, and protein-protein interaction analyses suggested functioning of PRNP protein in multiple pathways, and possible other functional constraint selections. In conclusion, a complex selection history favoring excessive replacement changes together with weak purifying selection possibly driven by frequency-dependent selection is driving PRNP sequence evolution. Our results is not unique only to the Turkish and Ethiopian samples, but can be generalized to global sheep populations.

Continuous Variation in an Aposematic Pattern Affects Background Contrast, but Is Not Associated With Differences in Microhabitat Use

Variation in aposematic signals was once predicted to be rare, yet in recent years it has become increasingly well documented. Despite increases in the frequency with which polytypism and polymorphism have been suggested to occur, population-wide variance is rarely quantified. We comprehensively sampled a subpopulation of the poison frog Oophaga sylvatica, a species which is polytypic across its distribution and also shows considerable within-population polymorphism. On one hand, color pattern polymorphism could be the result of multifarious selection acting to balance different signaling functions and leading to the evolution of discrete sub-morphs which occupy different fitness peaks. Alternatively, variance could simply be due to relaxed selection, where variation would be predicted to be continuous. We used visual modeling of conspecific and heterospecific observers to quantify the extent of within population phenotypic variation and assess whether this variation produced distinct signals. We found that, despite considerable color pattern variation, variance could not be partitioned into distinct groups, but rather all viewers would be likely to perceive variation as continuous. Similarly, we found no evidence that frog color pattern contrast was either enhanced or diminished in the frogs’ chosen microhabitats compared to alternative patches in which conspecifics were observed. Within population phenotypic variance therefore does not seem to be indicative of strong selection toward multiple signaling strategies, but rather pattern divergence has likely arisen due to weak purifying selection, or neutral processes, on a signal that is highly salient to both conspecifics and predators.

Characterization of the Complete Mitochondrial Genome of the Bromeliad Crab Metopaulias depressus (Rathbun, 1896) (Crustacea: Decapoda: Brachyura: Sesarmidae).

Metopaulias depressus is a non-marine crab endemic to Jamaica that dwells in rainforest bromeliads and exhibits elaborate active parental care behavior. Current genomic resources on M. depressus are rare, limiting the understanding of its adaptation to terrestrial life in species that evolved from marine ancestors. This study reports the complete mitochondrial genome of M. depressus assembled using Sanger sequencing. The AT-rich mitochondrial genome of M. depressus is 15,765 bp in length and comprises 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22 transfer RNA genes. A single 691 bp-long intergenic space is assumed to be the control region (CR) or D-loop. A set of selective pressure analyses indicate that the entirety of the PCGs experience purifying selection. Cox1, cox2, nad5, cox3, and atp6 experience strong purifying selection, and atp8 experiences weak purifying selection compared to the rest of the PCGs. The secondary structures of most tRNA genes exhibit a standard ‘cloverleaf’ structure, with the exception of trnS1, which lacks the dihydroxyuridine (DHU) arm but not the loop, the trnH gene, which lacks the thymine pseudouracil cytosine (T) loop but not the arm, and trnM, which exhibits an overly developed T loop. A maximum likelihood phylogenetic analysis based on all PCGs indicated that M. depressus is more closely related to the genera Clistocoeloma, Nanosesarma, and Parasesarma than to Chiromantes, Geosesarma, and Orisarma. This study contributes to deciphering the phylogenetic relationships within the family Sesarmidae and represents a new genomic resource for this iconic crab species.

Transmission of SARS-CoV-2 in domestic cats imposes a narrow bottleneck.

The evolutionary mechanisms by which SARS-CoV-2 viruses adapt to mammalian hosts and, potentially, undergo antigenic evolution depend on the ways genetic variation is generated and selected within and between individual hosts. Using domestic cats as a model, we show that SARS-CoV-2 consensus sequences remain largely unchanged over time within hosts, while dynamic sub-consensus diversity reveals processes of genetic drift and weak purifying selection. We further identify a notable variant at amino acid position 655 in Spike (H655Y), which was previously shown to confer escape from human monoclonal antibodies. This variant arises rapidly and persists at intermediate frequencies in index cats. It also becomes fixed following transmission in two of three pairs. These dynamics suggest this site may be under positive selection in this system and illustrate how a variant can quickly arise and become fixed in parallel across multiple transmission pairs. Transmission of SARS-CoV-2 in cats involved a narrow bottleneck, with new infections founded by fewer than ten viruses. In RNA virus evolution, stochastic processes like narrow transmission bottlenecks and genetic drift typically act to constrain the overall pace of adaptive evolution. Our data suggest that here, positive selection in index cats followed by a narrow transmission bottleneck may have instead accelerated the fixation of S H655Y, a potentially beneficial SARS-CoV-2 variant. Overall, our study suggests species- and context-specific adaptations are likely to continue to emerge. This underscores the importance of continued genomic surveillance for new SARS-CoV-2 variants as well as heightened scrutiny for signatures of SARS-CoV-2 positive selection in humans and mammalian model systems.

Pervasive Strong Selection at the Level of Codon Usage Bias in Drosophila melanogaster.

Codon usage bias (CUB), where certain codons are used more frequently than expected by chance, is a ubiquitous phenomenon and occurs across the tree of life. The dominant paradigm is that the proportion of preferred codons is set by weak selection. While experimental changes in codon usage have at times shown large phenotypic effects in contrast to this paradigm, genome-wide population genetic estimates have supported the weak selection model. Here we use deep genomic population sequencing of two Drosophila melanogaster populations to measure selection on synonymous sites in a way that allowed us to estimate the prevalence of both weak and strong purifying selection. We find that selection in favor of preferred codons ranges from weak (|Nes| ∼ 1) to strong (|Nes| > 10), with strong selection acting on 10–20% of synonymous sites in preferred codons. While previous studies indicated that selection at synonymous sites could be strong, this is the first study to detect and quantify strong selection specifically at the level of CUB. Further, we find that CUB-associated polymorphism accounts for the majority of strong selection on synonymous sites, with secondary contributions of splicing (selection on alternatively spliced genes, splice junctions, and spliceosome-bound sites) and transcription factor binding. Our findings support a new model of CUB and indicate that the functional importance of CUB, as well as synonymous sites in general, have been underestimated.

Population genomics of the facultatively asexual duckweed Spirodela polyrhiza.

Clonal propagation allows some plant species to achieve massive population sizes quickly but also reduces the evolutionary independence of different sites in the genome. We examine genome-wide genetic diversity in Spirodela polyrhiza, a duckweed that reproduces primarily asexually. We find that this geographically widespread and numerically abundant species has very low levels of genetic diversity. Diversity at nonsynonymous sites relative to synonymous sites is high, suggesting that purifying selection is weak. A potential explanation for this observation is that a very low frequency of sex renders selection ineffective. However, there is a pronounced decay in linkage disequilibrium over 40kb, suggesting that though sex may be rare at the individual level it is not too infrequent at the population level. In addition, neutral diversity is affected by the physical proximity of selected sites, which would be unexpected if sex was exceedingly rare at the population level. The amount of genetic mixing as assessed by the decay in linkage disequilibrium is not dissimilar from selfing species such as Arabidopsis thaliana, yet selection appears to be much less effective in duckweed. We discuss alternative explanations for the signature of weak purifying selection.

OR23 Whole gene analysis of HLA-A, -B and -C alleles observed in 100,000 individuals demonstrate distinct evolutionary and selection patterns in each Class I gene

Aim Classical HLA Class I genes (HLA-A, -B and -C) is thought to be evolved from ancestral gene(s) via duplication events and considered to be homologous. Exon 2 and 3 of those genes encodes Antigen Recognition Sites (ARS). These exons are highly polymorphic and it was believed due to positive selection. We have analyzed over 100,000 human samples that were whole gene sequenced for HLA-A, -B and -C and contained 1188, 1237 and 1185 unique patterns for each gene respectively. Methods We performed whole gene HLA Class I gene sequencing on PacificBio platform as described earlier (Cereb et al. Hum Immunol. 2015 Dec;76(12):963–74). Mean values of nonsynonymous and synonymous divergence (dN and dS, respectively) for all codons for all pairwise comparisons among alleles were determined using SNPGenie, an implementation of the Nei-Gojobori method (Nei and Gojobori, Mol. Biol. Evol. 3:418–426.) for use with next-generation sequencing data (Mol. Biol. Evol. 3:418–426.). Results In summary regarding positive selection, the following exons give signatures of positive selection: (1) HLA-A exons 2, 3, and 5 (weak); (2) HLA-B exons 2 and 3 (weak); and (3) HLA-C exons 1, 6, and 7. Regarding HLA-C, this seems quite revolutionary, and may help to explain previous studies that failed to identify evidence of positive selection specifically favoring amino acid charge differences in the peptide binding residues (PBR) of HLA-C only. HLA-C has the most extreme signatures of positive selection heretofore observed, but in non-PBR exons: dN/dS for exon 1 is 17.94, and dN/dS for exon 7 is 19.02 — both very high values. However, it is noteworthy that these huge ratios are not due to a great increase in dN, but rather to a great decrease in d. By contrast, PBR exons 2 and 3 have overall signatures of weak purifying selection at the HLA-C locus.Overall values for dN at HLA-A, -B, and -C are 0.0305, 0.0306, and 0.0240, respectively. Overall values for dS are 0.0395, 0.0344, and 0.0300, respectively. Therefore, in terms of selectively neutral diversity, diversity is greatest in HLA-A, less in -B, and least in -C. Conclusions Whole gene sequencing of HLA Class I will provide us the foundation to understand the evolutionary forces shaping the regions and the functions of various regions of those genes. N. Cereb: 6. Stock Shareholder; Company/Organization; PacBio and Hitogenetics. C.W. Nelson: 2. Consultant; Company/Organization; Histogenetics. S. Yang: 6. Stock Shareholder; Company/Organization; PacBio and Histogenetics.

Genomic analysis reveals major determinants of cis-regulatory variation in Capsella grandiflora

Understanding the causes of cis-regulatory variation is a long-standing aim in evolutionary biology. Although cis-regulatory variation has long been considered important for adaptation, we still have a limited understanding of the selective importance and genomic determinants of standing cis-regulatory variation. To address these questions, we studied the prevalence, genomic determinants, and selective forces shaping cis-regulatory variation in the outcrossing plant Capsella grandiflora We first identified a set of 1,010 genes with common cis-regulatory variation using analyses of allele-specific expression (ASE). Population genomic analyses of whole-genome sequences from 32 individuals showed that genes with common cis-regulatory variation (i) are under weaker purifying selection and (ii) undergo less frequent positive selection than other genes. We further identified genomic determinants of cis-regulatory variation. Gene body methylation (gbM) was a major factor constraining cis-regulatory variation, whereas presence of nearby transposable elements (TEs) and tissue specificity of expression increased the odds of ASE. Our results suggest that most common cis-regulatory variation in C. grandiflora is under weak purifying selection, and that gene-specific functional constraints are more important for the maintenance of cis-regulatory variation than genome-scale variation in the intensity of selection. Our results agree with previous findings that suggest TE silencing affects nearby gene expression, and provide evidence for a link between gbM and cis-regulatory constraint, possibly reflecting greater dosage sensitivity of body-methylated genes. Given the extensive conservation of gbM in flowering plants, this suggests that gbM could be an important predictor of cis-regulatory variation in a wide range of plant species.

The Impact of Selection, Gene Conversion, and Biased Sampling on the Assessment of Microbial Demography.

Recent studies have linked demographic changes and epidemiological patterns in bacterial populations using coalescent-based approaches. We identified 26 studies using skyline plots and found that 21 inferred overall population expansion. This surprising result led us to analyze the impact of natural selection, recombination (gene conversion), and sampling biases on demographic inference using skyline plots and site frequency spectra (SFS). Forward simulations based on biologically relevant parameters from Escherichia coli populations showed that theoretical arguments on the detrimental impact of recombination and especially natural selection on the reconstructed genealogies cannot be ignored in practice. In fact, both processes systematically lead to spurious interpretations of population expansion in skyline plots (and in SFS for selection). Weak purifying selection, and especially positive selection, had important effects on skyline plots, showing patterns akin to those of population expansions. State-of-the-art techniques to remove recombination further amplified these biases. We simulated three common sampling biases in microbiological research: uniform, clustered, and mixed sampling. Alone, or together with recombination and selection, they further mislead demographic inferences producing almost any possible skyline shape or SFS. Interestingly, sampling sub-populations also affected skyline plots and SFS, because the coalescent rates of populations and their sub-populations had different distributions. This study suggests that extreme caution is needed to infer demographic changes solely based on reconstructed genealogies. We suggest that the development of novel sampling strategies and the joint analyzes of diverse population genetic methods are strictly necessary to estimate demographic changes in populations where selection, recombination, and biased sampling are present.

Small Things Considered

By now, many of us are aware that ∼8% of the human genome is retroviral, compared to only ∼1.5% that codes for human proteins. These sequences are the result of ancient infections of germ cells by retroviruses, which insert into DNA as part of their life cycle. Upon insertion into the genome, such endogenous retroviruses (ERVs) evolve neutrally and are eventually inactivated by mutations and by host repression of their expression. Therefore, these viruses are no longer viruses but a collection of broken parts. Considering that we carry out ∼10,000 billion cell divisions in a lifetime, the maintenance of this viral junk heap requires considerable resources. So why do we keep it around? Probably because of weak purifying selection combined with the hunch that “this stuff might come in handy someday.” In fact, it appears that an important evolutionary innovation—the placenta—has depended on parts pulled from this “genomic junk drawer.”

Tissue-Specific Evolution of Protein Coding Genes in Human and Mouse.

Protein-coding genes evolve at different rates, and the influence of different parameters, from gene size to expression level, has been extensively studied. While in yeast gene expression level is the major causal factor of gene evolutionary rate, the situation is more complex in animals. Here we investigate these relations further, especially taking in account gene expression in different organs as well as indirect correlations between parameters. We used RNA-seq data from two large datasets, covering 22 mouse tissues and 27 human tissues. Over all tissues, evolutionary rate only correlates weakly with levels and breadth of expression. The strongest explanatory factors of purifying selection are GC content, expression in many developmental stages, and expression in brain tissues. While the main component of evolutionary rate is purifying selection, we also find tissue-specific patterns for sites under neutral evolution and for positive selection. We observe fast evolution of genes expressed in testis, but also in other tissues, notably liver, which are explained by weak purifying selection rather than by positive selection.

Identification and expression dynamics of three WUSCHEL related homeobox 13 (WOX13) genes in peanut.

WUSCHEL-related homeobox (WOX) genes play key roles in plant stem cell maintenance and development. WOX genes showed specific expression patterns which are important for their functions. WOX13 subfamily genes as the ancestor genes of this family were less studied in the past. In this study, we cloned three Arachis hypogaea (peanut) WOX13 (AhWOX13) subfamily genes from peanut: WOX13A and WOX13B1, 2. WOX13B1 encoded a same protein as WOX13B2, and there were only two-base difference between these two genes. Differential expression patterns were observed for these three AhWOX13 subfamily genes in different tissues and developmental stages. Phylogenic trees analysis showed that these AhWOX13 subfamily genes were the most conserved WOX genes and belonged to the ancient clade of WOX family. This was also supported by the conserved motif analysis. Selective pressure analysis showed that the WOX family genes mainly underwent weak purifying selection (ω = 0.58097), while many positive mutations accumulated during the evolution history. Under the purifying selection, gene duplication event and loss of duplicated gene play important roles in the expansion and evolution of WOX family.

Molecular Polymorphism in the MTA and MTB Mating Type Genes of Tetrahymena thermophila and Related Asexual Species.

Each of the seven mating types of Tetrahymena thermophila is determined by a pair of large genes, MTA and MTB, whose expression peaks at early conjugation. Each protein consists of a mating-type specific domain and a common transmembrane domain. To assess variation in natural populations, regions of both domains from wild isolates expressing mating types V and VII were analyzed. Corresponding regions of amicronucleates incapable of mating also were examined. MTA and MTB showed high haplotype diversity, with greater sequence variation in MTB. Mating type VII was less variable than mating type V, suggesting more recent origin. No polymorphism distinguished between mat1- and mat2-like alleles encoding different arrays of mating types, nor did polymorphisms give evidence of population structure. MTA and MTB variants have different phylogenies, suggesting independent rather than concerted evolution, and are under weak purifying selection. Codon usage is less biased than for housekeeping genes, and reassigned glutamine encoding stop codons are preferentially used. Amicronucleate T. thermophila and closely related nsp15 and nsp25 have higher levels of nucleotide and amino acid substitution, consistent with cox1 distances. The results suggest that complete sequencing of mating type genes of wild isolates coupled with functional analysis will be informative.