ALBERT

Beschreibung: Deciphering the evolutionary processes driving nucleotide variation in multiallelic genes is limited by the number of genetic systems in which such genes occur. The complementary sex determiner ( csd ) gene in the honey bee Apis mellifera is an informative example for studying allelic diversity and the underlying evolutionary forces in a well-described model of balancing selection. Acting as the primary signal of sex determination, diploid individuals heterozygous for csd develop into females, whereas csd homozygotes are diploid males that have zero fitness. Examining 77 of the functional heterozygous csd allele pairs, we established a combinatorical criteria that provide insights into the minimum number of amino acid differences among those pairs. Given a data set of 244 csd sequences, we show that the total number of csd alleles found in A. mellifera ranges from 53 (locally) to 87 (worldwide), which is much higher than was previously reported (20). Using a coupon-collector model, we extrapolate the presence of in total 116–145 csd alleles worldwide. The hypervariable region (HVR) is of particular importance in determining csd allele specificity, and we provide for this region evidence of high evolutionary rate for length differences exceeding those of microsatellites. The proportion of amino acids driven by positive selection and the rate of nonsynonymous substitutions in the HVR-flanking regions reach values close to 1 but differ with respect to the HVR length. Using a model of csd coalescence, we identified the high originating rate of csd specificities as a major evolutionary force, leading to an origin of a novel csd allele every 400,000 years. The csd polymorphism frequencies in natural populations indicate an excess of new mutations, whereas signs of ancestral transspecies polymorphism can still be detected. This study provides a comprehensive view of the enormous diversity and the evolutionary forces shaping a multiallelic gene.

Beschreibung: In sexually reproducing animals, primordial germ cells (PGCs) are often set aside early in embryogenesis, a strategy that minimizes the risk of genomic damage associated with replication and mitosis during the cell cycle. Here, we have used germ line markers ( piwi , vasa , and nanos ) and microinjected cell lineage tracers to show that PGC specification in the leech genus Helobdella follows a different scenario: in this hermaphrodite, the male and female PGCs segregate from somatic lineages only after more than 20 rounds of zygotic mitosis; the male and female PGCs share the same (mesodermal) cell lineage for 19 rounds of zygotic mitosis. Moreover, while all three markers are expressed in both male and female reproductive tissues of the adult, they are expressed differentially between the male and female PGCs of the developing embryo: piwi and vasa are expressed preferentially in female PGCs at a time when nanos is expressed preferentially in male PGCs. A priori, the delayed segregation of male and female PGCs from somatic tissues and from one another increases the probability of mutations affecting both male and female PGCs of a given individual. We speculate that this suite of features, combined with a capacity for self-fertilization, may contribute to the dramatically rearranged genome of Helobdella robusta relative to other animals.

Beschreibung: Experimental evolution in combination with whole-genome sequencing (evolve and resequence [E&R]) is a promising approach to define the genotype–phenotype map and to understand adaptation in evolving populations. Many previous studies have identified a large number of putative selected sites (i.e., candidate loci), but it remains unclear to what extent these loci are genuine targets of selection or experimental noise. To address this question, we exposed the same founder population to two different selection regimes—a hot environment and a cold environment—and quantified the genomic response in each. We detected large numbers of putative selected loci in both environments, albeit with little overlap between the two sets of candidates, indicating that most resulted from habitat-specific selection. By quantifying changes across multiple independent biological replicates, we demonstrate that most of the candidate SNPs were false positives that were linked to selected sites over distances much larger than the typical linkage disequilibrium range of Drosophila melanogaster . We show that many of these mid- to long-range associations were attributable to large segregating inversions and confirm by computer simulations that such patterns could be readily replicated when strong selection acts on rare haplotypes. In light of our findings, we outline recommendations to improve the performance of future Drosophila E&R studies which include using species with negligible inversion loads, such as D. mauritiana and D. simulans , instead of D. melanogaster .

Beschreibung: Immunization with the vOka vaccine prevents varicella (chickenpox) in children and susceptible adults. The vOka vaccine strain comprises a mixture of genotypes and, despite attenuation, causes rashes in small numbers of recipients. Like wild-type virus, the vaccine establishes latency in neuronal tissue and can later reactivate to cause Herpes zoster (shingles). Using hybridization-based methodologies, we have purified and sequenced vOka directly from skin lesions. We show that alleles present in the vaccine can be recovered from the lesions and demonstrate the presence of a severe bottleneck between inoculation and lesion formation. Genotypes in any one lesion appear to be descended from one to three vaccine-genotypes with a low frequency of novel mutations. No single vOka haplotype and no novel mutations are consistently present in rashes, indicating that neither new mutations nor recombination with wild type are critical to the evolution of vOka rashes. Instead, alleles arising from attenuation (i.e., not derived from free-living virus) are present at lower frequencies in rash genotypes. We identify 11 loci at which the ancestral allele is selected for in vOka rash formation and show genotypes in rashes that have reactivated from latency cannot be distinguished from rashes occurring immediately after inoculation. We conclude that the vOka vaccine, although heterogeneous, has not evolved to form rashes through positive selection in the mode of a quasispecies, but rather alleles that were essentially neutral during the vaccine production have been selected against in the human subjects, allowing us to identify key loci for rash formation.

Beschreibung: Bitter taste perception influences human nutrition and health, and the genetic variation underlying this trait may play a role in disease susceptibility. To better understand the genetic architecture and patterns of phenotypic variability of bitter taste perception, we sequenced a 996 bp region, encompassing the coding exon of TAS2R16 , a bitter taste receptor gene, in 595 individuals from 74 African populations and in 94 non-Africans from 11 populations. We also performed genotype–phenotype association analyses of threshold levels of sensitivity to salicin, a bitter anti-inflammatory compound, in 296 individuals from Central and East Africa. In addition, we characterized TAS2R16 mutants in vitro to investigate the effects of polymorphic loci identified at this locus on receptor function. Here, we report striking signatures of positive selection, including significant Fay and Wu’s H statistics predominantly in East Africa, indicating strong local adaptation and greater genetic structure among African populations than expected under neutrality. Furthermore, we observed a "star-like" phylogeny for haplotypes with the derived allele at polymorphic site 516 associated with increased bitter taste perception that is consistent with a model of selection for "high-sensitivity" variation. In contrast, haplotypes carrying the "low-sensitivity" ancestral allele at site 516 showed evidence of strong purifying selection. We also demonstrated, for the first time, the functional effect of nonsynonymous variation at site 516 on salicin phenotypic variance in vivo in diverse Africans and showed that most other nonsynonymous substitutions have weak or no effect on cell surface expression in vitro, suggesting that one main polymorphism at TAS2R16 influences salicin recognition. Additionally, we detected geographic differences in levels of bitter taste perception in Africa not previously reported and infer an East African origin for high salicin sensitivity in human populations.

Beschreibung: Gene conversion is the nonreciprocal exchange of genetic material between homologous chromosomes. Multiple lines of evidence from a variety of taxa strongly suggest that gene conversion events are biased toward GC-bearing alleles. However, in Drosophila , the data have largely been indirect and unclear, with some studies supporting the predictions of a GC-biased gene conversion model and other data showing contradictory findings. Here, we test whether gene conversion events are GC-biased in Drosophila melanogaster using whole-genome polymorphism and divergence data. Our results provide no support for GC-biased gene conversion and thus suggest that this process is unlikely to significantly contribute to patterns of polymorphism and divergence in this system.

Beschreibung: Substitution rates vary between species, and many explanations regarding the causes of this variation have been proposed. Here we consider how new genomic data on the per-generation mutation rate impinge on proposed hypotheses for substitution rate variation in primates. We propose that the generation-time effect as it is usually understood cannot explain the observed rate variation, but instead that selection for decreased somatic mutation rates can. By considering the disparate causes underlying mutation rate changes in recent human history, we also show that the per-generation mutation rate is increasing even as the per-cell-division rate is decreasing.

Beschreibung: Gene looping, defined as the physical interaction between the promoter and terminator regions of a RNA polymerase II-transcribed gene, is widespread in yeast and mammalian cells. Gene looping has been shown to play important roles in transcription. Gene-loop formation is dependent on regulatory proteins localized at the 5' and 3' ends of genes, such as TFIIB. However, whether other factors contribute to gene looping remains to be elucidated. Here, we investigated the contribution of intrinsic DNA and chromatin structures to gene looping. We found that Saccharomyces cerevisiae looped genes show high DNA bendability around middle and 3/4 regions in open reading frames (ORFs). This bendability pattern is conserved between yeast species, whereas the position of bendability peak varies substantially among species. Looped genes in human cells also show high DNA bendability. Nucleosome positioning around looped ORF middle regions is unstable. We also present evidence indicating that this unstable nucleosome positioning is involved in gene looping. These results suggest a mechanism by which DNA bendability and unstable nucleosome positioning could assist in the formation of gene loops.

Beschreibung: Although efficient influenza vaccines are designed on a regular basis, the only protection of human populations against an unforeseen virus such as during the H1N1 pandemic in 2009 might be antiviral drugs. Adamantanes and neuraminidase inhibitors (Oseltamivir) represent two classes of such drugs that target the viral matrix protein 2 and neuraminidase, respectively. Although the emergence of resistance to both drugs has been described, the timing and spread of the acquisition of either single or dual resistances by different hosts is still unclear. Using a multilayered phylogenetic approach based on relaxed molecular clocks and large-scale maximum likelihood approaches, we show that Adamantane resistance evolved multiple times in various subtypes and hosts, possibly in breeding contexts (swine); and Oseltamivir resistance was also found in different subtypes and hosts, but its transmission is only sustained in humans. Furthermore, the dynamics of the emergence of antiviral resistance were examined for each drug. This showed that although the first mutations conferring resistance to Adamantanes precede US Food and Drug Administration (FDA) approval, general resistance emerged 15–38 years post-drug approval. This is in contrast to Oseltamivir resistance mutations that emerged at most 7 years after FDA approval of the drug. This study demonstrates the power of large-scale analyses to uncover and monitor the emergence dynamics of drug resistance.

Beschreibung: Horizontal gene transfer (HGT) has been described as a common mechanism of transferring genetic material between prokaryotes, whereas genetic transfers from eukaryotes to prokaryotes have been rarely documented. Here we report a rare case of HGT in which plant expansin genes that code for plant cell-wall loosening proteins were transferred from plants to bacteria, fungi, and amoebozoa. In several cases, the species in which the expansin gene was found is either in intimate association with plants or is a known plant pathogen. Our analyses suggest that at least two independent genetic transfers occurred from plants to bacteria and fungi. These events were followed by multiple HGT events within bacteria and fungi. We have also observed that in bacteria expansin genes have been independently fused to DNA fragments that code for an endoglucanase domain or for a carbohydrate binding module, pointing to functional convergence at the molecular level. Furthermore, the functional similarities between microbial expansins and their plant xenologs suggest that these proteins mediate microbial–plant interactions by altering the plant cell wall and therefore may provide adaptive advantages to these species. The evolution of these nonplant expansins represents a unique case in which bacteria and fungi have found innovative and adaptive ways to interact with and infect plants by acquiring genes from their host. This evolutionary paradigm suggests that despite their low frequency such HGT events may have significantly contributed to the evolution of prokaryotic and eukaryotic species.

Beschreibung: Vertebrate Tas2r taste receptors bind to bitter compounds, which are typically poisonous, to elicit bitter sensation to prevent the ingestion of toxins. Previous studies noted a marked variation in the number of Tas2r genes among species, but the underlying cause is unclear. To address this question, we compile the Tas2r gene repertoires from 41 mammals, 4 birds, 2 reptiles, 1 amphibian, and 6 fishes. The number of intact Tas2r genes varies from 0 in the bottlenose dolphin to 51 in the Western clawed frog, with numerous expansions and contractions of the gene family throughout vertebrates, especially among tetrapods. The Tas2r gene number in a species correlates with the fraction of plants in its diet. Because plant tissues contain more toxic compounds than animal tissues do, our observation supports the hypothesis that dietary toxins are a major selective force shaping the diversity of the Tas2r repertoire.

Beschreibung: The human APOBEC3 gene cluster locus encodes polynucleotide cytidine deaminases. Although many act as viral restriction factors through mutation of single-stranded DNA, recent reports have shown that human APOBEC3A was capable of efficiently hypermutating nuclear DNA and inducing DNA breaks in genomic DNA. In addition, the enzyme was unique in efficiently deaminating 5-methylcytidine in single-stranded DNA. To appreciate the evolutionary relevance of these activities, we analyzed A3A-related enzymes from the rhesus and tamarin monkey, horse, sheep, dog, and panda. All proved to be orthologous to the human enzyme in all these activities revealing strong conservation more than 148 My. Hence, their singular role in DNA catabolism is a well-established mechanism probably outweighing any deleterious or pathological roles such as genomic instability and cancer formation.

Beschreibung: DNA methylation at CpG dinucleotides can significantly increase the rate of cytosine-to-thymine mutations and the level of sequence divergence. Although the correlations between DNA methylation and genomic sequence evolution have been widely studied, an unaddressed yet fundamental question is how DNA methylation is associated with the conservation of individual nucleotides in different sequence contexts. Here, we demonstrate that in mammalian exons, the correlations between DNA methylation and the conservation of individual nucleotides are dependent on the type of exonic sequence (coding or untranslated), the degeneracy of coding nucleotides, background selection pressure, and the relative position (first or nonfirst exon in the transcript) where the nucleotides are located. For untranslated and nonzero-fold degenerate nucleotides, methylated sites are less conserved than unmethylated sites regardless of background selection pressure and the relative position of the exon. For zero-fold degenerate (or nondegenerate) nucleotides, however, the reverse trend is observed in nonfirst coding exons and first coding exons that are under stringent background selection pressure. Furthermore, cytosine-to-thymine mutations at methylated zero-fold degenerate nucleotides are predicted to be more detrimental than those that occur at unmethylated nucleotides. As zero-fold and nonzero-fold degenerate nucleotides are very close to each other, our results suggest that the "functional resolution" of DNA methylation may be finer than previously recognized. In addition, the positive correlation between CpG methylation and the level of conservation at zero-fold degenerate nucleotides implies that CpG methylation may serve as an "indicator" of functional importance of these nucleotides.

Beschreibung: Eukaryotic genomes are mosaics of genes acquired from their prokaryotic ancestors, the eubacterial endosymbiont that gave rise to the mitochondrion and its archaebacterial host. Genomic footprints of the prokaryotic merger at the origin of eukaryotes are still discernable in eukaryotic genomes, where gene expression and function correlate with their prokaryotic ancestry. Molecular chaperones are essential in all domains of life as they assist the functional folding of their substrate proteins and protect the cell against the cytotoxic effects of protein misfolding. Eubacteria and archaebacteria code for slightly different chaperones, comprising distinct protein folding pathways. Here we study the evolution of the eukaryotic protein folding pathways following the endosymbiosis event. A phylogenetic analysis of all 64 chaperones encoded in the Saccharomyces cerevisiae genome revealed 25 chaperones of eubacterial ancestry, 11 of archaebacterial ancestry, 10 of ambiguous prokaryotic ancestry, and 18 that may represent eukaryotic innovations. Several chaperone families (e.g., Hsp90 and Prefoldin) trace their ancestry to only one prokaryote group, while others, such as Hsp40 and Hsp70, are of mixed ancestry, with members contributed from both prokaryotic ancestors. Analysis of the yeast chaperone–substrate interaction network revealed no preference for interaction between chaperones and substrates of the same origin. Our results suggest that the archaebacterial and eubacterial protein folding pathways have been reorganized and integrated into the present eukaryotic pathway. The highly integrated chaperone system of yeast is a manifestation of the central role of chaperone-mediated folding in maintaining cellular fitness. Most likely, both archaebacterial and eubacterial chaperone systems were essential at the very early stages of eukaryogenesis, and the retention of both may have offered new opportunities for expanding the scope of chaperone-mediated folding.

Beschreibung: Coalescent methods are widely used to infer the demographic history of populations from gene genealogies. These approaches—often referred to as phylodynamic methods—have proven especially useful for reconstructing the dynamics of rapidly evolving viral pathogens. Yet, population dynamics inferred from viral genealogies often differ widely from those observed from other sources of epidemiological data, such as hospitalization records. We demonstrate how a modeling framework that allows for the direct fitting of mechanistic epidemiological models to genealogies can be used to test different hypotheses about what ecological factors cause phylodynamic inferences to differ from observed dynamics. We use this framework to test different hypotheses about why dengue serotype 1 (DENV-1) population dynamics in southern Vietnam inferred using existing phylodynamic methods differ from hospitalization data. Specifically, we consider how factors such as seasonality, vector dynamics, and spatial structure can affect inferences drawn from genealogies. The coalescent models we derive to take into account vector dynamics and spatial structure reveal that these ecological complexities can substantially affect coalescent rates among lineages. We show that incorporating these additional ecological complexities into coalescent models can also greatly improve estimates of historical population dynamics and lead to new insights into the factors shaping viral genealogies.

Beschreibung: Complex metazoan bodies require cell-to-cell communication for development, a process often mediated by signaling molecules binding to specific receptors. Relatively few signaling pathways have been recruited during evolution to build multicellular animals from unicellular zygotes. Of these few signaling pathways, one of particular importance is the receptor tyrosine kinase (RTK) pathway. In metazoans, fibroblast growth factors (FGFs) bind to receptors in the RTK family, but the origin of the FGF gene family has so far remained a mystery. Here we show that extant bona fide FGFs most likely originated from proteins bearing an FGF-like domain that arose in a choanoflagellate/metazoan ancestor. We found orthologous genes closely related to FGF in choanoflagellates as well as in many metazoans such as sponges, acoels, protostomes, or nonvertebrate deuterostomes. We also show that these genes have a common evolutionary history with Retinitis Pigmentosa 1 ( RP1) . Even if some metazoan signaling pathways emerged long before multicellularity, we show that FGFs, like their receptors, originated in a eumetazoan ancestor.

Beschreibung: Whole genome duplications (WGDs) occurred in the distant evolutionary history of many lineages and are particularly frequent in the flowering plant lineages. Following paleopolyploidization in plants, most duplicated genes are deleted by intrachromosomal recombination, a process referred to as fractionation. In the examples studied so far, genes are disproportionately lost from one of the parental subgenomes (biased fractionation) and the subgenome having lost the lowest number of genes is more expressed (genome dominance). In the present study, we analyzed the pattern of gene deletion and gene expression following the most recent WGD in banana (alpha event) and extended our analyses to seven other sequenced plant genomes: poplar, soybean, medicago, arabidopsis, sorghum, brassica, and maize. We propose a new class of ancient WGD, with Musa (alpha), poplar, and soybean as members, where genes are both deleted and expressed to an equal extent (unbiased fractionation and genome equivalence). We suggest that WGDs with genome dominance and biased fractionation (Class I) may result from ancient allotetraploidies, while WGDs without genome dominance or biased fractionation (Class II) may result from ancient autotetraploidies.

Beschreibung: Ribosomal and small nuclear RNAs (snRNAs) comprise numerous modified nucleotides. The modification patterns are retained during evolution, making it even possible to project them from yeast onto human. The stringent conservation of modification sites and the slow evolution of rRNAs and snRNAs contradicts the rapid evolution of small nucleolar RNA (snoRNA) sequences. To explain this discrepancy, we investigated the coevolution of snoRNAs and their targeted sites throughout vertebrates. To measure and evaluate the conservation of RNA-RNA interactions, we defined the interaction conservation index (ICI). It combines the quality of individual interaction with the scope of its conservation in a set of species and serves as an efficient measure to evaluate the conservation of the interaction of snoRNA and target. We show that functions of homologous snoRNAs are evolutionarily stable, thus, members of the same snoRNA family guide equivalent modifications. The conservation of snoRNA sequences is high at target binding regions while the remaining sequence varies significantly. In addition to elucidating principles of correlated evolution, we were able, with the help of the ICI measure, to assign functions to previously orphan snoRNAs and to associate snoRNAs as partners to known chemical modifications unassigned to a given snoRNA. Furthermore, we used predictions of snoRNA functions in conjunction with sequence conservation to identify distant homologies. Because of the high overall entropy of snoRNA sequences, such relationships are hard to detect by means of sequence homology search methods alone.

Beschreibung: The branch-site model is a widely popular approach that accommodates for the lineage- and the site-specific heterogeneity of natural selection regimes among coding sequences. This model relies on prior knowledge of the (foreground) lineage(s) evolving under positive selection at some sites. Unfortunately, such prior information is not always available in practice. A more recent technique (Guindon S, Rodrigo A, Dyer K, Huelsenbeck J. 2004. Modeling the site-specific variation of selection patterns along lineages. Proc Natl Acad Sci USA 101:12957–12962) alleviates this issue by explicitly modeling the variability of selection patterns using a stochastic process. However, the performance of this approach for deciding whether a set of homologous sequences evolved under positive selection at some point has not been assessed yet. This study compares the sensitivity and specificity of tests for positive selection derived from both the standard and the stochastic approaches using extensive simulations. We show that the two methods have low proportions of type I errors, that is, they tend to be conservative when testing the null hypothesis of no positive selection if sequences truly evolve under neutral or negative selection regimes. Also, the standard approach is more powerful than the stochastic one when the prior knowledge on foreground lineages is correct. When this prior is incorrect, however, the stochastic approach outperforms the standard model in a broad range of conditions. Additional comparisons also suggest that the stochastic branch-site method compares favorably with the recently proposed mixed-effects model of evolution of Murrell et al. (Murrell B, Wertheim JO, Moola S, Weighill T, Scheffler K, Pond SLK. 2012. Detecting individual sites subject to episodic diversifying selection. PLoS Genet . 8:e1002764). Altogether, our results show that the standard branch-site model is well suited to confirmatory analyses, whereas the stochastic approach should be preferred over the standard or the mixed-effects ones for exploratory studies.

Beschreibung: The cation chloride cotransporter (CCCs) family comprises of four subfamilies—K + -Cl – cotransporters (KCCs), Na + -K + -2Cl – cotransporters (NKCCs), and Na + -Cl – cotransporters (NCCs)—and possibly two additional members—CCC interacting protein (CIP1) and polyamine transporters (CCC9)—as well. Altogether, CCCs can play essential physiological roles in transepithelial ion reabsorption and secretion, cell volume regulation, and inhibitory neurotransmission and so are present across all domains of life. To gain insight into the evolution of this family, we performed a comprehensive phylogenetic analysis using publically available genomic information. Our results clearly support CIP1 as being a true CCC based on shared evolutionary history. By contrast, the status of CCC9 in this regard remains equivocal. We also reveal the existence of a single ancestral CCC gene present in Archaea, from which numerous duplication events at the base of archaeans and eukaryotes lead to the divergence and subsequent neofunctionalization of the paralogous CCC subfamilies. A diversity of ensuing gene-loss events resulted in the complex distribution of CCCs present across the different taxa. Importantly, the occurrence of KCCs in "basal" metazoan taxa like sponges would allow an early formation of fast hyperpolarizing neurotransmission in metazoans. Gene duplications within the CCC subfamilies in vertebrates (in particular, KCCs, NKCCs, and NCCs) lend further evidence to the 2R hypothesis of two rounds of genome duplication at the base of the vertebrate lineage, especially in concert with our syntenic cluster analyses. This increased number of KCCs, NKCCs, and NCCs isoforms facilitates their further, important subfunctionalization in the vertebrate lineage.

Beschreibung: Standing genetic variation provides a rich reservoir of potentially useful mutations facilitating the adaptation to novel environments. Experimental evolution studies have demonstrated that rapid and strong phenotypic responses to selection can also be obtained in the laboratory. When combined with the next-generation sequencing technology, these experiments promise to identify the individual loci contributing to adaption. Nevertheless, until now, very little is known about the design of such evolve & resequencing (E&R) studies. Here, we use forward simulations of entire genomes to evaluate different experimental designs that aim to maximize the power to detect selected variants. We show that low linkage disequilibrium in the starting population, population size, duration of the experiment, and the number of replicates are the key factors in determining the power and accuracy of E&R studies. Furthermore, replication of E&R is more important for detecting the targets of selection than increasing the population size. Using an optimized design, beneficial loci with a selective advantage as low as s = 0.005 can be identified at the nucleotide level. Even when a large number of loci are selected simultaneously, up to 56% can be reliably detected without incurring large numbers of false positives. Our computer simulations suggest that, with an adequate experimental design, E&R studies are a powerful tool to identify adaptive mutations from standing genetic variation and thereby provide an excellent means to analyze the trajectories of selected alleles in evolving populations.

Beschreibung: Levels of selective constraint vary among proteins. Although strong constraint on a protein is often attributed to its functional importance, evolutionary rate may also be limited if a protein is fragile, such that a large proportion of amino acid replacements reduce its fitness. To determine the relative contributions of essentiality and fragility to selective constraint, we compared relationships of selection against nonsense mutations ( s non ) and selection against missense mutations ( s mis ) to protein sequence conservation ( K a ). As expected, s non is greater than s mis ; however, the correlation between s mis and K a is nearly three times stronger than the correlation between s non and K a . Moreover, examination of relationships to gene expression level, tissue specificity, and number of protein–protein interactions shows that s mis is more strongly correlated than s non to all three measures of biological function. Thus, our analysis reveals that slowly evolving proteins are under strong selective constraint primarily because they are fragile, and that this association likely exists because allowing a protein to function improperly, rather than removing it from a biological network, can negatively affect the functions of other molecules it interacts with and their downstream products.

Beschreibung: Although evolutionary studies of gene function often rely on RNA interference, the ideal approach would use reverse genetics to create null mutations for cross-species comparisons and forward genetics to identify novel genes in each species. We have used transcription activator-like effector nucleases (TALENs) to facilitate both approaches in Caenorhabditis nematodes. First, by combining golden gate cloning and TALEN technology, we can induce frameshifting mutations in any gene. Second, by combining this approach with bioinformatics we can predict and create the resources needed for forward genetic analysis in species like Caenorhabditis briggsae . Although developing genetic model organisms used to require years to isolate marker mutations, balancers, and tools, with TALENs, these reagents can now be produced in months. Furthermore, the analysis of nonsense mutants in related model organisms allows a directed approach for making these markers and tools. When used together, these methods could simplify the adaptation of other organisms for forward and reverse genetics.

Beschreibung: The gene regulatory network of floral organ cell fate specification of Arabidopsis thaliana is a robust developmental regulatory module. Although such finding was proposed to explain the overall conservation of floral organ types and organization among angiosperms, it has not been confirmed that the network components are conserved at the molecular level among flowering plants. Using the genomic data that have accumulated, we address the conservation of the genes involved in this network and the forces that have shaped its evolution during the divergence of angiosperms. We recovered the network gene homologs for 18 species of flowering plants spanning nine families. We found that all the genes are highly conserved with no evidence of positive selection. We studied the sequence conservation features of the genes in the context of their known biological function and the strength of the purifying selection acting upon them in relation to their placement within the network. Our results suggest an association between protein length and sequence conservation, evolutionary rates, and functional category. On the other hand, we found no significant correlation between the strength of purifying selection and gene placement. Our results confirm that the studied robust developmental regulatory module has been subjected to strong functional constraints. However, unlike previous studies, our results do not support the notion that network topology plays a major role in constraining evolutionary rates. We speculate that the dynamical functional role of genes within the network and not just its connectivity could play an important role in constraining evolution.

Beschreibung: The high frequency of alternative splicing among the serine/arginine-rich (SR) family of proteins in plants has been linked to important roles in gene regulation during development and in response to environmental stress. In this article, we have searched and manually annotated all the SR proteins in the genomes of maize and sorghum. The experimental validation of gene structure by reverse transcription-polymerase chain reaction (RT-PCR) analysis revealed, with few exceptions, that SR genes produced multiple isoforms of transcripts by alternative splicing. Despite sharing high structural similarity and conserved positions of the introns, the profile of alternative splicing diverged significantly between maize and sorghum for the vast majority of SR genes. These include many transcript isoforms discovered by RT-PCR and not represented in extant expressed sequence tag (EST) collection. However, we report the occurrence of various maize and sorghum SR mRNA isoforms that display evolutionary conservation of splicing events with their homologous SR genes in Arabidopsis and moss. Our data also indicate an important role of both 5' and 3' untranslated regions in the regulation of SR gene expression. These observations have potentially important implications for the processes of evolution and adaptation of plants to land.

Beschreibung: The nuclear genomes of euglenids contain three types of introns: conventional spliceosomal introns, nonconventional introns for which a splicing mechanism is unknown (variable noncanonical borders, RNA secondary structure bringing together intron ends), and so-called intermediate introns, which combine features of conventional and nonconventional introns. Analysis of two genes, tubA and tubB , from 20 species of euglenids reveals contrasting distribution patterns of conventional and nonconventional introns—positions of conventional introns are conserved, whereas those of the nonconventional ones are unique to individual species or small groups of closely related taxa. Moreover, in the group of phototrophic euglenids, 11 events of conventional intron loss versus 15 events of nonconventional intron gain were identified. A comparison of all nonconventional intron sequences highlighted the most conserved elements in their sequence and secondary structure. Our results led us to put forward two hypotheses. 1) The first one posits that mutational changes in intron sequence could lead to a change in their excision mechanism—intermediate introns would then be a transitional form between the conventional and nonconventional introns. 2) The second hypothesis concerns the origin of nonconventional introns—because of the presence of inverted repeats near their ends, insertion of MITE-like transposon elements is proposed as a possible source of new introns.

Beschreibung: Geraniaceae plastid genomes are highly rearranged, and each of the four genera already sequenced in the family has a distinct genome organization. This study reports plastid genome sequences of six additional species, Francoa sonchifolia , Melianthus villosus , and Viviania marifolia from Geraniales, and Pelargonium alternans , California macrophylla , and Hypseocharis bilobata from Geraniaceae. These genome sequences, combined with previously published species, provide sufficient taxon sampling to reconstruct the ancestral plastid genome organization of Geraniaceae and the rearrangements unique to each genus. The ancestral plastid genome of Geraniaceae has a 4 kb inversion and a reduced, Pelargonium -like small single copy region. Our ancestral genome reconstruction suggests that a few minor rearrangements occurred in the stem branch of Geraniaceae followed by independent rearrangements in each genus. The genomic comparison demonstrates that a series of inverted repeat boundary shifts and inversions played a major role in shaping genome organization in the family. The distribution of repeats is strongly associated with breakpoints in the rearranged genomes, and the proportion and the number of large repeats (〉20 bp and 〉60 bp) are significantly correlated with the degree of genome rearrangements. Increases in the degree of plastid genome rearrangements are correlated with the acceleration in nonsynonymous substitution rates (d N ) but not with synonymous substitution rates (d S ). Possible mechanisms that might contribute to this correlation, including DNA repair system and selection, are discussed.

Beschreibung: Many mitochondrial and plastid protein complexes contain subunits that are encoded in different genomes. In animals, nuclear-encoded mitochondrial proteins often exhibit rapid sequence evolution, which has been hypothesized to result from selection for mutations that compensate for changes in interacting subunits encoded in mutation-prone animal mitochondrial DNA. To test this hypothesis, we analyzed nuclear genes encoding cytosolic and organelle ribosomal proteins in flowering plants. The model angiosperm genus Arabidopsis exhibits low organelle mutation rates, typical of most plants. Nevertheless, we found that (nuclear-encoded) subunits of organelle ribosomes in Arabidopsis have higher amino acid sequence polymorphism and divergence than their counterparts in cytosolic ribosomes, suggesting that organelle ribosomes experience relaxed functional constraint. However, the observed difference between organelle and cytosolic ribosomes was smaller than in animals and could be partially attributed to rapid evolution in N-terminal organelle-targeting peptides that are not involved in ribosome function. To test the role of organelle mutation more directly, we used transcriptomic data from an angiosperm genus ( Silene ) with highly variable rates of organelle genome evolution. We found that Silene species with unusually fast-evolving mitochondrial and plastid DNA exhibited increased amino acid sequence divergence in ribosomal proteins targeted to the organelles but not in those that function in cytosolic ribosomes. Overall, these findings support the hypothesis that rapid organelle genome evolution has selected for compensatory mutations in nuclear-encoded proteins. We conclude that coevolution between interacting subunits encoded in different genomic compartments within the eukaryotic cell is an important determinant of variation in rates of protein sequence evolution.

Beschreibung: Studies of the Neanderthal and Denisovan genomes demonstrate archaic hominin introgression in Eurasians. Here, we present evidence of Neanderthal introgression within the chromosome 3p21.31 region, occurring with a high frequency in East Asians (ranging from 49.4% to 66.5%) and at a low frequency in Europeans. We also detected a signal of strong positive selection in this region only in East Asians. Our data indicate that likely candidate targets of selection include rs12488302-T and its associated alleles—among which four are nonsynonymous, including rs35455589-G in HYAL2 , a gene related to the cellular response to ultraviolet-B irradiation. Furthermore, suggestive evidence supports latitude-dependent selection, implicating a role of ultraviolet-B. Interestingly, the distribution of rs35455589-G suggests that this allele was lost during the exodus of ancestors of modern Eurasians from Africa and reintroduced to Eurasians from Neanderthals.

Beschreibung: Defining homologous genes is important in many evolutionary studies but raises obvious issues. Some of these issues are conceptual and stem from our assumptions of how a gene evolves, others are practical, and depend on the algorithmic decisions implemented in existing software. Therefore, to make progress in the study of homology, both ontological and epistemological questions must be considered. In particular, defining homologous genes cannot be solely addressed under the classic assumptions of strong tree thinking, according to which genes evolve in a strictly tree-like fashion of vertical descent and divergence and the problems of homology detection are primarily methodological. Gene homology could also be considered under a different perspective where genes evolve as "public goods," subjected to various introgressive processes. In this latter case, defining homologous genes becomes a matter of designing models suited to the actual complexity of the data and how such complexity arises, rather than trying to fit genetic data to some a priori tree-like evolutionary model, a practice that inevitably results in the loss of much information. Here we show how important aspects of the problems raised by homology detection methods can be overcome when even more fundamental roots of these problems are addressed by analyzing public goods thinking evolutionary processes through which genes have frequently originated. This kind of thinking acknowledges distinct types of homologs, characterized by distinct patterns, in phylogenetic and nonphylogenetic unrooted or multirooted networks. In addition, we define "family resemblances" to include genes that are related through intermediate relatives, thereby placing notions of homology in the broader context of evolutionary relationships. We conclude by presenting some payoffs of adopting such a pluralistic account of homology and family relationship, which expands the scope of evolutionary analyses beyond the traditional, yet relatively narrow focus allowed by a strong tree-thinking view on gene evolution.

Beschreibung: Carnivorous Lentibulariaceae exhibit the most sophisticated implementation of the carnivorous syndrome in plants. Their unusual lifestyle coincides with distinct genomic peculiarities such as the smallest angiosperm nuclear genomes and extremely high nucleotide substitution rates across all genomic compartments. Here, we report the complete plastid genomes from each of the three genera Pinguicula , Utricularia , and Genlisea , and investigate plastome-wide changes in their molecular evolution as the carnivorous syndrome unfolds. We observe a size reduction by up to 9% mostly due to the independent loss of genes for the plastid NAD(P)H dehydrogenase and altered proportions of plastid repeat DNA, as well as a significant plastome-wide increase of substitution rates and microstructural changes. Protein-coding genes across all gene classes show a disproportional elevation of nonsynonymous substitutions, particularly in Utricularia and Genlisea . Significant relaxation of purifying selection relative to noncarnivores occurs in the plastid-encoded fraction of the photosynthesis ATP synthase complex, the photosystem I, and in several other photosynthesis and metabolic genes. Shifts in selective regimes also affect housekeeping genes including the plastid-encoded polymerase, for which evidence for relaxed purifying selection was found once during the transition to carnivory, and a second time during the diversification of the family. Lentibulariaceae significantly exhibit enhanced rates of nucleotide substitution in most of the 130 noncoding regions. Various factors may underlie the observed patterns of relaxation of purifying selection and substitution rate increases, such as reduced net photosynthesis rates, alternative paths of nutrient uptake (including organic carbon), and impaired DNA repair mechanisms.

Beschreibung: Standard protein phylogenetic models use fixed rate matrices of amino acid interchange derived from analyses of large databases. Differences between the stationary amino acid frequencies of these rate matrices from those of a data set of interest are typically adjusted for by matrix multiplication that converts the empirical rate matrix to an exchangeability matrix which is then postmultiplied by the amino acid frequencies in the alignment. The result is a time-reversible rate matrix with stationary amino acid frequencies equal to the data set frequencies. On the basis of population genetics principles, we develop an amino acid substitution-selection model that parameterizes the fitness of an amino acid as the logarithm of the ratio of the frequency of the amino acid to the frequency of the same amino acid under no selection. The model gives rise to a different sequence of matrix multiplications to convert an empirical rate matrix to one that has stationary amino acid frequencies equal to the data set frequencies. We incorporated the substitution-selection model with an improved amino acid class frequency mixture (cF) model to partially take into account site-specific amino acid frequencies in the phylogenetic models. We show that 1) the selection models fit data significantly better than corresponding models without selection for most of the 21 test data sets; 2) both cF and cF selection models favored the phylogenetic trees that were inferred under current sophisticated models and methods for three difficult phylogenetic problems (the positions of microsporidia and breviates in eukaryote phylogeny and the position of the root of the angiosperm tree); and 3) for data simulated under site-specific residue frequencies, the cF selection models estimated trees closer to the generating trees than a standard model or cF without selection. We also explored several ways of estimating amino acid frequencies under neutral evolution that are required for these selection models. By better modeling the amino acid substitution process, the cF selection models will be valuable for phylogenetic inference and evolutionary studies.

Beschreibung: Gonadotropin-releasing hormone (GnRH) is a critical reproductive regulator in vertebrates. Homologous peptides are also found in invertebrates, with a variety of characterized functions. In the amphioxus, an invertebrate that provides the best model for the transition to vertebrates, four GnRH receptors (GnRHRs) were previously described, but their native ligands were not identified. Using a more sensitive search methodology with hidden Markov models, we identified the first GnRH-like peptide confirmed in the amphioxus Branchiostoma floridae. This peptide specifically activated one of the four GnRHRs. Although the primary structure of this peptide was divergent from any previously isolated GnRH peptide, the minimal conserved residues found in all other GnRH superfamily members were retained. The peptide was immunolocalized in proximity of the central canal of the anterior nerve cord, a region where other neuropeptides and receptors have been found. Additionally, the amphioxus GnRH-like gene was positioned in a locus surrounded by syntenic homologs of the human GnRH paralogon. The amphioxus GnRH-like peptide, with its distinct primary structure, activated a receptor with equal potency to multiple ligands that span the GnRH superfamily.

Beschreibung: Among land plants, angiosperms have the structurally most labile mitochondrial (mt) genomes. In contrast, the so-called early land plants (e.g., mosses) seem to have completely static mt chromosomes. We assembled the complete mt genomes from 12 mosses spanning the moss tree of life, to assess 1) the phylogenetic depth of the conserved mt gene content and order and 2) the correlation between scattered sequence repeats and gene order lability in land plants. The mt genome of most mosses is approximately 100 kb in size, and thereby the smallest among land plants. Based on divergence time estimates, moss mt genome structure has remained virtually frozen for 350 My, with only two independent gene losses and a single gene relocation detected across the macroevolutionary tree. This is the longest period of mt genome stasis demonstrated to date in a plant lineage. The complete lack of intergenic repeat sequences, considered to be essential for intragenomic recombinations, likely accounts for the evolutionary stability of moss mt genomes.

Beschreibung: The major families of chromatin remodelers have been conserved throughout eukaryotic evolution. Because they play broad, pleiotropic roles in gene regulation, it was not known if their functions could change rapidly. Here, we show that major alterations in the use of chromatin remodelers are possible, because the nucleosome remodeling factor (NURF) complex has acquired a unique role in the sperm/oocyte decision of the nematode Caenorhabditis briggsae . First, lowering the activity of C. briggsae NURF-1 or ISW-1, the core components of the NURF complex, causes germ cells to become oocytes rather than sperm. This observation is based on the analysis of weak alleles and null mutations that were induced with TALENs and on RNA interference. Second, qRT–polymerase chain reaction data show that the C. briggsae NURF complex promotes the expression of Cbr-fog-1 and Cbr-fog-3 , two genes that control the sperm/oocyte decision. This regulation occurs in the third larval stage and affects the expression of later spermatogenesis genes. Third, double mutants reveal that the NURF complex and the transcription factor TRA-1 act independently on Cbr-fog-1 and Cbr-fog-3 . TRA-1 binds both promoters, and computer analyses predict that these binding sites are buried in nucleosomes, so we suggest that the NURF complex alters chromatin structure to allow TRA-1 access to Cbr-fog-1 and Cbr-fog-3 . Finally, lowering NURF activity by mutation or RNA interference does not affect this trait in other nematodes, including the sister species C. nigoni , so it must have evolved recently. We conclude that altered chromatin remodeling could play an important role in evolutionary change.

Beschreibung: The captive genetic management of threatened species strives to preserve genetic diversity and avoid inbreeding to ensure populations remain available, healthy, and viable for future reintroduction. Determining and responding to the genetic status of captive populations is therefore paramount to these programs. Here, we genotyped 19 microsatellite loci for 240 captive giant pandas ( Ailuropoda melanoleuca ) (~64% of the captive population) from four breeding centers, Wolong (WL), Chengdu (CD), Louguantai (LGT), and Beijing (BJ), and analyzed 655 bp of mitochondrial DNA control region sequence for 220 of these animals. High levels of genetic diversity and low levels of inbreeding were estimated in the breeding centers, indicating that the captive population is genetically healthy and deliberate further genetic input from wild animals is unnecessary. However, the LGT population faces a higher risk of inbreeding, and significant genetic structure was detected among breeding centers, with LGT–CD and WL–BJ clustering separately. Based on these findings, we highlight that: 1) the LGT population should be managed as an independent captive population to resemble the genetic distinctness of their Qinling Mountain origins; 2) exchange between CD and WL should be encouraged because of similar wild founder sources; 3) the selection of captive individuals for reintroduction should consider their geographic origin, genetic background, and genetic contribution to wild populations; and 4) combining our molecular genetic data with existing pedigree data will better guide giant panda breeding and further reduce inbreeding into the future.

Beschreibung: Gene regulatory networks (GRNs) describe the progression of transcriptional states that take a single-celled zygote to a multicellular organism. It is well documented that GRNs can evolve extensively through mutations to cis -regulatory modules (CRMs). Transcription factor proteins that bind these CRMs may also evolve to produce novelty. Coding changes are considered to be rarer, however, because transcription factors are multifunctional and hence are more constrained to evolve in ways that will not produce widespread detrimental effects. Recent technological advances have unearthed a surprising variation in DNA-binding abilities, such that individual transcription factors may recognize both a preferred primary motif and an additional secondary motif. This provides a source of modularity in function. Here, we demonstrate that orthologous transcription factors can also evolve a changed preference for a secondary binding motif, thereby offering an unexplored mechanism for GRN evolution. Using protein-binding microarray, surface plasmon resonance, and in vivo reporter assays, we demonstrate an important difference in DNA-binding preference between Tbrain protein orthologs in two species of echinoderms, the sea star, Patiria miniata , and the sea urchin, Strongylocentrotus purpuratus . Although both orthologs recognize the same primary motif, only the sea star Tbr also has a secondary binding motif. Our in vivo assays demonstrate that this difference may allow for greater evolutionary change in timing of regulatory control. This uncovers a layer of transcription factor binding divergence that could exist for many pairs of orthologs. We hypothesize that this divergence provides modularity that allows orthologous transcription factors to evolve novel roles in GRNs through modification of binding to secondary sites.

Beschreibung: MicroRNAs (miRNAs) are a class of endogenous small noncoding RNAs that regulate gene expression either by degrading target mRNAs or by suppressing protein translation. miRNAs have been found to be involved in many biological processes, such as development, differentiation, and growth. However, the evolution of miRNA regulatory functions and networks has not been well studied. In this study, we conducted a cross-species analysis to study the evolution of cardiac miRNAs and their regulatory functions and networks. We found that conserved cardiac miRNA target genes have maintained highly conserved cardiac functions. Additionally, most of cardiac miRNA target genes in human with annotations of cardiac functions evolved from the corresponding homologous targets, which are also involved in heart development-related functions. On the basis of these results, we investigated the functional evolution of cardiac miRNAs and presented a functional evolutionary map. From this map, we identified the evolutionary time at which the cardiac miRNAs became involved in heart development or function and found that the biological processes of heart development evolved earlier than those of heart functions, for example, heart contraction/relaxation or cardiac hypertrophy. Our study of the evolution of the cardiac miRNA regulatory networks revealed the emergence of new regulatory functional branches during evolution. Furthermore, we discovered that early evolved cardiac miRNA target genes tend to participate in the early stages of heart development. This study sheds light on the evolution of developmental features of genes regulated by cardiac miRNAs.

Beschreibung: Calcium signaling is one of the most extensively employed signal transduction mechanisms in life. As life evolved into increasingly complex organisms, Ca 2+ acquired more extensive and varied functions. Here, we compare genes encoding proteins that govern Ca 2+ entry and exit across cells or organelles within organisms of early eukaryotic evolution into fungi, plants, and animals. Recent phylogenomics analyses reveal a complex Ca 2+ signaling machinery in the apusozoan protist Thecamonas trahens , a putative unicellular progenitor of Opisthokonta. We compare T. trahens Ca 2+ signaling to that in a marine bikont protist, Aurantiochytrium limacinum , and demonstrate the conservation of key Ca 2+ signaling molecules in the basally diverging alga Cyanophora paradoxa . Particularly, our findings reveal the conservation of the CatSper channel complex in Au. limacinum and C. paradoxa , suggesting that the CatSper complex likely originated from an ancestral Ca 2+ signaling machinery at the root of early eukaryotic evolution prior to the unikont/bikont split.

Beschreibung: Polymerases are essential for life, being responsible for replication, transcription, and the repair of nucleic acid molecules. Those that share a right-hand-shaped fold and catalytic site structurally similar to the DNA polymerase I of Escherichia coli may catalyze RNA- or DNA-dependent RNA polymerization, reverse transcription, or DNA replication in eukarya, archaea, bacteria, and their viruses. We have applied novel computational methods for structure-based clustering and phylogenetic analyses of this functionally diverse polymerase superfamily, which currently comprises six families. We identified a structural core common to all right-handed polymerases, composed of 57 amino acid residues, harboring two positionally and chemically conserved residues, the catalytic aspartates. The structural conservation within each of the six families is considerable, for example, the structural core shared by family Y DNA polymerases covers over 90% of the polymerase domain of the Sulfolobus solfataricus Dpo4. Our phylogenetic analyses propose an early separation of RNA-dependent polymerases that use primers from those that are primer-independent. Furthermore, the exchange of polymerase genes between viruses and their hosts is evident. Because of this horizontal gene transfer, the phylogeny of polymerases does not always reflect the evolutionary history of the corresponding organisms.

Beschreibung: Phylogenetic analyses of molecular data require a quantitative model for how sequences evolve. Traditionally, the details of the site-specific selection that governs sequence evolution are not known a priori, making it challenging to create evolutionary models that adequately capture the heterogeneity of selection at different sites. However, recent advances in high-throughput experiments have made it possible to quantify the effects of all single mutations on gene function. I have previously shown that such high-throughput experiments can be combined with knowledge of underlying mutation rates to create a parameter-free evolutionary model that describes the phylogeny of influenza nucleoprotein far better than commonly used existing models. Here, I extend this work by showing that published experimental data on TEM-1 beta-lactamase (Firnberg E, Labonte JW, Gray JJ, Ostermeier M. 2014. A comprehensive, high-resolution map of a gene’s fitness landscape. Mol Biol Evol. 31:1581–1592) can be combined with a few mutation rate parameters to create an evolutionary model that describes beta-lactamase phylogenies much better than most common existing models. This experimentally informed evolutionary model is superior even for homologs that are substantially diverged (about 35% divergence at the protein level) from the TEM-1 parent that was the subject of the experimental study. These results suggest that experimental measurements can inform phylogenetic evolutionary models that are applicable to homologs that span a substantial range of sequence divergence.

Beschreibung: In eukaryotic cells, identical proteins can be located in more than a single subcellular compartment, a phenomenon termed dual targeting. We hypothesized that dual-targeted proteins should be more evolutionary conserved than exclusive mitochondrial proteins, due to separate selective pressures administered by the different compartments to maintain the functions associated with the protein sequences. We employed codon usage bias, propensity for gene loss, phylogenetic relationships, conservation analysis at the DNA level, and gene expression, to test our hypothesis. Our findings indicate that, indeed, dual-targeted proteins are significantly more conserved than their exclusively targeted counterparts. We then used this trait of gene conservation, together with previously identified traits of dual-targeted proteins (such as protein net charge and mitochondrial targeting sequence strength) to 1) create, for the first time (due to addition of conservation parameters), a tool for the prediction of dual-targeted mitochondrial proteins based on protein and mRNA sequences, and 2) show that molecular mechanisms involving one versus two translation products are not correlated with specific dual-targeting parameters. Finally, we discuss what evolutionary pressure maintains protein dual targeting in eukaryotes and deduce, as we initially hypothesized, that it is the discrete functions of these proteins in the different subcellular compartments, regardless of their dual-targeting mechanism.

Beschreibung: Reliable estimates of the rate at which DNA accumulates mutations (the substitution rate) are crucial for our understanding of the evolution and past demography of virtually any species. In humans, there are considerable uncertainties around these rates, with substantial variation among recent published estimates. Substitution rates have traditionally been estimated by associating dated events to the root (e.g., the divergence between humans and chimpanzees) or to internal nodes in a phylogenetic tree (e.g., first entry into the Americas). The recent availability of ancient mitochondrial DNA sequences allows for a more direct calibration by assigning the age of the sequenced samples to the tips within the human phylogenetic tree. But studies also vary greatly in the methodology employed and in the sequence panels analyzed, making it difficult to tease apart the causes for the differences between previous estimates. To clarify this issue, we compiled a comprehensive data set of 350 ancient and modern human complete mitochondrial DNA genomes, among which 146 were generated for the purpose of this study and estimated substitution rates using calibrations based both on dated nodes and tips. Our results demonstrate that, for the same data set, estimates based on individual dated tips are far more consistent with each other than those based on nodes and should thus be considered as more reliable.

Beschreibung: Understanding the demographic history of populations and species is a central issue in evolutionary biology and molecular ecology. In this work, we develop a maximum-likelihood method for the inference of past changes in population size from microsatellite allelic data. Our method is based on importance sampling of gene genealogies, extended for new mutation models, notably the generalized stepwise mutation model (GSM). Using simulations, we test its performance to detect and characterize past reductions in population size. First, we test the estimation precision and confidence intervals coverage properties under ideal conditions, then we compare the accuracy of the estimation with another available method ( MSVAR ) and we finally test its robustness to misspecification of the mutational model and population structure. We show that our method is very competitive compared with alternative ones. Moreover, our implementation of a GSM allows more accurate analysis of microsatellite data, as we show that the violations of a single step mutation assumption induce very high bias toward false contraction detection rates. However, our simulation tests also showed some limits, which most importantly are large computation times for strong disequilibrium scenarios and a strong influence of some form of unaccounted population structure. This inference method is available in the latest implementation of the MIGRAINE software package.

Beschreibung: Haplotype-based scans to detect natural selection are useful to identify recent or ongoing positive selection in genomes. As both real and simulated genomic data sets grow larger, spanning thousands of samples and millions of markers, there is a need for a fast and efficient implementation of these scans for general use. Here, we present selscan , an efficient multithreaded application that implements Extended Haplotype Homozygosity (EHH), Integrated Haplotype Score (iHS), and Cross-population EHH (XPEHH). selscan accepts phased genotypes in multiple formats, including TPED, and performs extremely well on both simulated and real data and over an order of magnitude faster than existing available implementations. It calculates iHS on chromosome 22 (22,147 loci) across 204 CEU haplotypes in 353 s on one thread (33 s on 16 threads) and calculates XPEHH for the same data relative to 210 YRI haplotypes in 578 s on one thread (52 s on 16 threads). Source code and binaries (Windows, OSX, and Linux) are available at https://github.com/szpiech/selscan .

Beschreibung: Gene loss is one of the main drivers in the evolution of genomes and species. The demonstration that a gene has been lost by pseudogenization is truly complete when one finds the pseudogene in the orthologous genomic region with respect to active genes in other species. In some cases, the identification of such orthologous loci is not possible because of chromosomal rearrangements or if the gene of interest has not yet been sequenced. This question is particularly important in the case of birds because the genomes of avian species possess only about 15,000 predicted genes, in comparison with 20,000 in mammals. Yet, gene loss raises the question of which functions are affected by the changes in gene counts. We describe a systematic approach that makes it possible to demonstrate gene loss in the chicken genome even if a pseudogene has not been found. By using phylogenetic and synteny analysis in vertebrates, genome-wide comparisons between the chicken genome and expressed sequence tags, RNAseq data analysis, statistical analysis of the chicken genome, and radiation hybrid mapping, we show that resistin, TNF α, and PAI-1 ( SERPINE1 ), three genes encoding adipokines inhibiting insulin sensitivity, have been lost in chicken and zebra finch genomes. Moreover, omentin, a gene encoding an adipokine that enhances insulin sensitivity, has also been lost in the chicken genome. Overall, only one adipokine inhibiting insulin sensitivity and five adipokines enhancing insulin sensitivity are still present in the chicken genome. These genetic differences between mammals and chicken, given the functions of the genes in mammals, would have dramatic consequences on chicken endocrinology, leading to novel equilibriums especially in the regulation of energy metabolism, insulin sensitivity, as well as appetite and reproduction.

Beschreibung: Y chromosomes, with their reduced effective population size, lack of recombination, and male-limited transmission, present a unique collection of constraints for the operation of natural selection. Male-limited transmission may greatly increase the efficacy of selection for male-beneficial mutations, but the reduced effective size also inflates the role of random genetic drift. Together, these defining features of the Y chromosome are expected to influence rates and patterns of molecular evolution on the Y as compared with X-linked or autosomal loci. Here, we use sequence data from 11 genes in 9 Drosophila species to gain insight into the efficacy of natural selection on the Drosophila Y relative to the rest of the genome. Drosophila is an ideal system for assessing the consequences of Y-linkage for molecular evolution in part because the gene content of Drosophila Y chromosomes is highly dynamic, with orthologous genes being Y-linked in some species whereas autosomal in others. Our results confirm the expectation that the efficacy of natural selection at weakly selected sites is reduced on the Y chromosome. In contrast, purifying selection on the Y chromosome for strongly deleterious mutations does not appear to be compromised. Finally, we find evidence of recurrent positive selection for 4 of the 11 genes studied here. Our results thus highlight the variable nature of the mode and impact of natural selection on the Drosophila Y chromosome.

Beschreibung: MicroRNAs (miRNAs) are endogenous RNA molecules that regulate gene expression posttranscriptionally. To date, the emergence of miRNAs and their patterns of sequence evolution have been analyzed in great detail. However, the extent to which miRNA expression levels have evolved over time, the role different evolutionary forces play in shaping these changes, and whether this variation in miRNA expression can reveal the interplay between miRNAs and mRNAs remain poorly understood. This is especially true for miRNA expressed during key developmental transitions. Here, we assayed miRNA expression levels immediately before (≥18BPF [18 h before puparium formation]) and after (PF) the increase in the hormone ecdysone responsible for triggering metamorphosis. We did so in four strains of Drosophila melanogaster and two closely related species. In contrast to their sequence conservation, approximately 25% of miRNAs analyzed showed significant within-species variation in male expression levels at ≥18BPF and/or PF. Additionally, approximately 33% showed modifications in their pattern of expression bias between developmental timepoints . A separate analysis of the ≥18BPF and PF stages revealed that changes in miRNA abundance accumulate linearly over evolutionary time at PF but not at ≥18BPF. Importantly, ≥18BPF-enriched miRNAs showed the greatest variation in expression levels both within and between species, so are the less likely to evolve under stabilizing selection. Functional attributes, such as expression ubiquity, appeared more tightly associated with lower levels of miRNA expression polymorphism at PF than at ≥18BPF. Furthermore, ≥18BPF- and PF-enriched miRNAs showed opposite patterns of covariation in expression with mRNAs, which denoted the type of regulatory relationship between miRNAs and mRNAs. Collectively, our results show contrasting patterns of functional divergence associated with miRNA expression levels during Drosophila ontogeny.

Beschreibung: SVA (SINE-R-VNTR- Alu ) elements constitute the youngest family of composite non-LTR retrotransposons in hominoid primates. The sequence of their assembly, however, remains unclear. Recently, a second family of VNTR-containing composites, LAVA (L1- Alu -VNTR- Alu ), has been identified in gibbons. We now report the existence of two additional VNTR composite families, PVA ( PTGR2 -VNTR- Alu ) and FVA (FRAM-VNTR- Alu ), in the genome of Nomascus leucogenys . Like LAVA, they share the 5'- Alu -like region and VNTR with SVA, but differ at their 3'-ends. The 3'-end of PVA comprises part of the PTGR2 gene, whereas FVA is characterized by the presence of a partial FRAM element in its 3'-domain. Splicing could be identified as the mechanism of acquisition of the variant 3'-ends in all four families of VNTR composites. SVAs have been shown to be mobilized by the L1 protein machinery in trans . A critical role in this process has been ascribed to their 5'-hexameric repeat/ Alu -like region. The Alu -like region displays specific features in each of the VNTR composite families/subfamilies with characteristic deletions found in the evolutionary younger subfamilies. Using reciprocal exchanges between SVA_E and PVA/FVA elements, we demonstrate that the structure, not the presence of the (CCCTCT) n / Alu -like region determines mobilization capacity. Combination of LAVA and SVA_E domains does not yield any active elements—suggesting the use of different combinations of host factors for the two major groups of VNTR composites. Finally, we demonstrate that the LAVA 3'-L1ME5 fragment attenuates mobilization capacity.

Beschreibung: Small self-cleaving RNAs, such as the paradigmatic Hammerhead ribozyme (HHR), have been recently found widespread in DNA genomes across all kingdoms of life. In this work, we found that new HHR variants are preserved in the ancient family of Penelope -like elements (PLEs), a group of eukaryotic retrotransposons regarded as exceptional for encoding telomerase-like retrotranscriptases and spliceosomal introns. Our bioinformatic analysis revealed not only the presence of minimalist HHRs in the two flanking repeats of PLEs but also their massive and widespread occurrence in metazoan genomes. The architecture of these ribozymes indicates that they may work as dimers, although their low self-cleavage activity in vitro suggests the requirement of other factors in vivo. In plants, however, PLEs show canonical HHRs, whereas fungi and protist PLEs encode ribozyme variants with a stable active conformation as monomers. Overall, our data confirm the connection of self-cleaving RNAs with eukaryotic retroelements and unveil these motifs as a significant fraction of the encoded information in eukaryotic genomes.

Beschreibung: In air-breathing vertebrates, the physiologically optimal blood-O 2 affinity is jointly determined by the prevailing partial pressure of atmospheric O 2 , the efficacy of pulmonary O 2 transfer, and internal metabolic demands. Consequently, genetic variation in the oxygenation properties of hemoglobin (Hb) may be subject to spatially varying selection in species with broad elevational distributions. Here we report the results of a combined functional and evolutionary analysis of Hb polymorphism in the rufous-collared sparrow ( Zonotrichia capensis ), a species that is continuously distributed across a steep elevational gradient on the Pacific slope of the Peruvian Andes. We integrated a population genomic analysis that included all postnatally expressed Hb genes with functional studies of naturally occurring Hb variants, as well as recombinant Hb (rHb) mutants that were engineered through site-directed mutagenesis. We identified three clinally varying amino acid polymorphisms: Two in the α A -globin gene, which encodes the α-chain subunits of the major HbA isoform, and one in the α D -globin gene, which encodes the α-chain subunits of the minor HbD isoform. We then constructed and experimentally tested single- and double-mutant rHbs representing each of the alternative α A -globin genotypes that predominate at different elevations. Although the locus-specific patterns of altitudinal differentiation suggested a history of spatially varying selection acting on Hb polymorphism, the experimental tests demonstrated that the observed amino acid mutations have no discernible effect on respiratory properties of the HbA or HbD isoforms. These results highlight the importance of experimentally validating the hypothesized effects of genetic changes in protein function to avoid the pitfalls of adaptive storytelling.

Beschreibung: The 20 protein-coding amino acids are found in proteomes with different relative abundances. The most abundant amino acid, leucine, is nearly an order of magnitude more prevalent than the least abundant amino acid, cysteine. Amino acid metabolic costs differ similarly, constraining their incorporation into proteins. On the other hand, a diverse set of protein sequences is necessary to build functional proteomes. Here, we present a simple model for a cost-diversity trade-off postulating that natural proteomes minimize amino acid metabolic flux while maximizing sequence entropy. The model explains the relative abundances of amino acids across a diverse set of proteomes. We found that the data are remarkably well explained when the cost function accounts for amino acid chemical decay. More than 100 organisms reach comparable solutions to the trade-off by different combinations of proteome cost and sequence diversity. Quantifying the interplay between proteome size and entropy shows that proteomes can get optimally large and diverse.

Beschreibung: Orthology inference is central to phylogenomic analyses. Phylogenomic data sets commonly include transcriptomes and low-coverage genomes that are incomplete and contain errors and isoforms. These properties can severely violate the underlying assumptions of orthology inference with existing heuristics. We present a procedure that uses phylogenies for both homology and orthology assignment. The procedure first uses similarity scores to infer putative homologs that are then aligned, constructed into phylogenies, and pruned of spurious branches caused by deep paralogs, misassembly, frameshifts, or recombination. These final homologs are then used to identify orthologs. We explore four alternative tree-based orthology inference approaches, of which two are new. These accommodate gene and genome duplications as well as gene tree discordance. We demonstrate these methods in three published data sets including the grape family, Hymenoptera, and millipedes with divergence times ranging from approximately 100 to over 400 Ma. The procedure significantly increased the completeness and accuracy of the inferred homologs and orthologs. We also found that data sets that are more recently diverged and/or include more high-coverage genomes had more complete sets of orthologs. To explicitly evaluate sources of conflicting phylogenetic signals, we applied serial jackknife analyses of gene regions keeping each locus intact. The methods described here can scale to over 100 taxa. They have been implemented in python with independent scripts for each step, making it easy to modify or incorporate them into existing pipelines. All scripts are available from https://bitbucket.org/yangya/phylogenomic_dataset_construction .

Beschreibung: Recombination between double-stranded DNA molecules is a key genetic process which occurs in a wide variety of organisms. Usually, crossing-over (CO) occurs during meiosis between genotypes with 98.0–99.9% sequence identity, because within-population nucleotide diversity only rarely exceeds 2%. However, some species are hypervariable and it is unclear how CO can occur between genotypes with less than 90% sequence identity. Here, we study CO in Schizophyllum commune , a hypervariable cosmopolitan basidiomycete mushroom, a frequently encountered decayer of woody substrates. We crossed two haploid individuals, from the United States and from Russia, and obtained genome sequences for their 17 offspring. The average genetic distance between the parents was 14%, making it possible to study CO at very high resolution. We found reduced levels of linkage disequilibrium between loci flanking the CO sites indicating that they are mostly confined to hotspots of recombination. Furthermore, CO events preferentially occurred in regions under stronger negative selection, in particular within exons that showed reduced levels of nucleotide diversity. Apparently, in hypervariable species CO must avoid regions of higher divergence between the recombining genomes due to limitations imposed by the mismatch repair system, with regions under strong negative selection providing the opportunity for recombination. These patterns are opposite to those observed in a number of less variable species indicating that population genomics of hypervariable species may reveal novel biological phenomena.

Beschreibung: Blood coagulation occurs through a cascade of enzymes and cofactors that produces a fibrin clot, while otherwise maintaining hemostasis. The 11 human coagulation factors (FG, FII–FXIII) have been identified across all vertebrates, suggesting that they emerged with the first vertebrates around 500 Ma. Human FVIII, FIX, and FXI are associated with thousands of disease-causing mutations. Here, we evaluated the strength of selective pressures on the 14 genes coding for the 11 factors during vertebrate evolution, and compared these with human mutations in FVIII, FIX, and FXI. Positive selection was identified for fibrinogen (FG), FIII, FVIII, FIX, and FX in the mammalian Primates and Laurasiatheria and the Sauropsida (reptiles and birds). This showed that the coagulation system in vertebrates was under strong selective pressures, perhaps to adapt against blood-invading pathogens. The comparison of these results with disease-causing mutations reported in FVIII, FIX, and FXI showed that the number of disease-causing mutations, and the probability of positive selection were inversely related to each other. It was concluded that when a site was under positive selection, it was less likely to be associated with disease-causing mutations. In contrast, sites under negative selection were more likely to be associated with disease-causing mutations and be destabilizing. A residue-by-residue comparison of the FVIII, FIX, and FXI sequence alignments confirmed this. This improved understanding of evolutionary changes in FVIII, FIX, and FXI provided greater insight into disease-causing mutations, and better assessments of the codon sites that may be mutated in applications of gene therapy.

Beschreibung: In Drosophila, chromosomes have been extensively reorganized during evolution, with most rearrangements affecting the gene order in chromosomal elements but not their gene content. The level of reorganization and the evidence for breakpoint reuse vary both between and within elements. The subito gene stands out as a gene involved in multiple rearrangements both because of its active single-gene transposition and because it is the nearest gene to diverse rearrangements breakpoints. Indeed, subito has undergone three single-gene transpositions and it is the nearest gene to the breakpoints of other single-gene transpositions and of two chromosomal inversions. Given that subito is involved in meiosis and therefore active in the female germ line, the high number of nearby fixed breakages might be related among others to the presumed high accessibility of the subito region to the machinery associated with double-strand breaks repair. A second important contributor would be the reduced and simple regulatory region of subito , which would imply that a fraction of the rearrangements originating from subito nearby breakages would have not affected either its pattern or timing of expression and would have, thus, not resulted in reduced fitness.

Beschreibung: The role of gene duplication in generating new genes and novel functions is well recognized and is exemplified by the digestion-related protein lysozyme. In ruminants, duplicated chicken-type lysozymes facilitate the degradation of symbiotic bacteria in the foregut. Chicken-type lysozyme has also been reported to show chitinase-like activity, yet no study has examined the molecular evolution of lysozymes in species that specialize on eating insects. Insectivorous bats number over 900 species, and lysozyme expression in the mouths of some of these species is associated with the ingestion of insect cuticle, suggesting a chitinase role. Here, we show that chicken-type lysozyme has undergone multiple duplication events in a major family of insect-eating bats (Vespertilionidae) and that new duplicates have undergone molecular adaptation. Examination of duplicates from two insectivorous bats— Pipistrellus abramus and Scotophilus kuhlii —indicated that the new copy was highly expressed in the tongue, whereas the other one was less tissue-specific. Functional assays applied to pipistrelle lysozymes confirmed that, of the two copies, the tongue duplicate was more efficient at breaking down glycol chitin, a chitin derivative. These results suggest that the evolution of lysozymes in vespertilionid bats has likely been driven in part by natural selection for insectivory.

Beschreibung: The time required to transcribe genes with long primary transcripts may limit their ability to be expressed in cells with short mitotic cycles, a phenomenon termed intron delay. As such short cycles are a hallmark of the earliest stages of insect development, we tested the impact of intron delay on the Drosophila developmental transcriptome. We find that long zygotically expressed genes show substantial delay in expression relative to their shorter counterparts, which is not observed for maternally deposited transcripts. Patterns of RNA-seq coverage along transcripts show that this delay is consistent with their inability to completely transcribe long transcripts, but not with transcriptional initiation-based regulatory control. We further show that highly expressed zygotic genes maintain compact transcribed regions across the Drosophila phylogeny, allowing conservation of embryonic expression patterns. We propose that the physical constraints of intron delay affect patterns of expression and the evolution of gene structure of a substantial portion of the Drosophila transcriptome .

Beschreibung: Ralstonia solanacearum , the causal agent of a lethal bacterial wilt plant disease, infects an unusually wide range of hosts. These hosts can further be split into plants where R. solanacearum is known to cause disease (original hosts) and those where this bacterium can grow asymptomatically (distant hosts). Moreover, this pathogen is able to adapt to many plants as supported by field observations reporting emergence of strains with enlarged pathogenic properties. To investigate the genetic bases of host adaptation, we conducted evolution experiments by serial passages of a single clone of the pathogen on three original and two distant hosts over 300 bacterial generations and then analyzed the whole-genome of nine evolved clones. Phenotypic analysis of the evolved clones showed that the pathogen can increase its fitness on both original and distant hosts although the magnitude of fitness increase was greater on distant hosts. Only few genomic modifications were detected in evolved clones compared with the ancestor but parallel evolutionary changes in two genes were observed in independent evolved populations. Independent mutations in the regulatory gene efpR were selected for in three populations evolved on beans, a distant host. Reverse genetic approaches confirmed that these mutations were associated with fitness gain on bean plants. This work provides a first step toward understanding the within-host evolutionary dynamics of R. solanacearum during infection and identifying bacterial genes subjected to in planta selection. The discovery of EfpR as a determinant conditioning host adaptation of the pathogen illustrates how experimental evolution coupled with whole-genome sequencing is a potent tool to identify novel molecular players involved in central life-history traits.

Beschreibung: The origin of novel genes and their evolutionary fates are long-standing questions in evolutionary biology. These questions become more complicated for genes conserved across various lineages, such as TRIM5 , an antiretroviral restriction factor and a retrovirus capsid sensor in immune signaling. TRIM5 has been subjected to numerous pathogenic challenges and undergone dynamic evolution, making it an excellent example for studying gene diversification. Previous studies among several species showed that TRIM5 gained genetic and functional novelty in a lineage-specific manner, either through gene duplication or a cyclophilin A retrotransposing into the TRIM5 locus, creating the gene fusion known as TRIM5-Cyclophilin A ( TRIMCyp ). To date, the general pattern of TRIM5 across the mammalian lineage remains elusive. In this study, we surveyed 36 mammalian genomes to verify a potentially novel TRIM5 pattern that uniquely seems to have occurred in tree shrews ( Tupaia belangeri ), and found that both gene duplication and retrotransposition worked jointly to form a specific TRIM5/TRIMCyp cluster not found among other mammals. Evolutionary analyses showed that tree shrew TRIMCyp ( tsTRIMCyp ) originated independently in comparison with previously reported TRIMCyps and underwent strong positive selection, whereas no signal of positive selection was detected for other tree shrew TRIM5 ( tsTRIM5 ) genes. Functional assay results suggest a functional divergence between tsTRIMCyp and its closest paralog TRIM5-4 , likely reflecting different fates under diverse evolutionary forces. These findings present a rare example of novel gene origination resulting from a combination of gene duplication, retrotransposition, and exon shuffling processes, providing a new paradigm to study genetic innovations and evolutionary fates of duplicated genes.

Beschreibung: Chelicerata represents one of the oldest groups of arthropods, with a fossil record extending to the Cambrian, and is sister group to the remaining extant arthropods, the mandibulates. Attempts to resolve the internal phylogeny of chelicerates have achieved little consensus, due to marked discord in both morphological and molecular hypotheses of chelicerate phylogeny. The monophyly of Arachnida, the terrestrial chelicerates, is generally accepted, but has garnered little support from molecular data, which have been limited either in breadth of taxonomic sampling or in depth of sequencing. To address the internal phylogeny of this group, we employed a phylogenomic approach, generating transcriptomic data for 17 species in combination with existing data, including two complete genomes. We analyzed multiple data sets containing up to 1,235,912 sites across 3,644 loci, using alternative approaches to optimization of matrix composition. Here, we show that phylogenetic signal for the monophyly of Arachnida is restricted to the 500 slowest-evolving genes in the data set. Accelerated evolutionary rates in Acariformes, Pseudoscorpiones, and Parasitiformes potentially engender long-branch attraction artifacts, yielding nonmonophyly of Arachnida with increasing support upon incrementing the number of concatenated genes. Mutually exclusive hypotheses are supported by locus groups of variable evolutionary rate, revealing significant conflicts in phylogenetic signal. Analyses of gene-tree discordance indicate marked incongruence in relationships among chelicerate orders, whereas derived relationships are demonstrably robust. Consistently recovered and supported relationships include the monophyly of Chelicerata, Euchelicerata, Tetrapulmonata, and all orders represented by multiple terminals. Relationships supported by subsets of slow-evolving genes include Ricinulei + Solifugae; a clade comprised of Ricinulei, Opiliones, and Solifugae; and a clade comprised of Tetrapulmonata, Scorpiones, and Pseudoscorpiones. We demonstrate that outgroup selection without regard for branch length distribution exacerbates long-branch attraction artifacts and does not mitigate gene-tree discordance, regardless of high gene representation for outgroups that are model organisms. Arachnopulmonata (new name) is proposed for the clade comprising Scorpiones + Tetrapulmonata (previously named Pulmonata).

Beschreibung: Many organisms survive fluctuating and extreme environmental conditions by manifesting multiple distinct phenotypes during adulthood by means of developmental processes that enable phenotypic plasticity. We report on the discovery of putative plasticity-enabling genes that are involved in transforming the gill of the euryhaline teleost fish, Fundulus heteroclitus , from its freshwater to its seawater gill-type, a process that alters both morphology and function. Gene expression that normally enables osmotic plasticity is inhibited by arsenic. Gene sets defined by antagonistic interactions between arsenic and salinity show reduced transcriptional variation among individual fish, suggesting unusually accurate and precise regulatory control of these genes, consistent with the hypothesis that they participate in a canalized developmental response. We observe that natural selection acts to preserve canalized gene expression in populations of killifish that are most tolerant to abrupt salinity change and that these populations show the least variability in their transcription of genes enabling plasticity of the gill. We found that genes participating in this highly canalized and conserved plasticity-enabling response had significantly fewer and less complex associations with transcriptional regulators than genes that respond only to arsenic or salinity. Collectively these findings, which are drawn from the relationships between environmental challenge, plasticity, and canalization among populations, suggest that the selective processes that facilitate phenotypic plasticity do so by targeting the regulatory networks that gives rise to the response. These findings also provide a generalized, conceptual framework of how genes might interact with the environment and evolve toward the development of plastic traits.

Beschreibung: Detecting positive selection in species with heterogeneous habitats and complex demography is notoriously difficult and prone to statistical biases. The model plant Arabidopsis thaliana exemplifies this problem: In spite of the large amounts of data, little evidence for classic selective sweeps has been found. Moreover, many aspects of the demography are unclear, which makes it hard to judge whether the few signals are indeed signs of selection, or false positives caused by demographic events. Here, we focus on Swedish A. thaliana and we find that the demography can be approximated as a two-population model. Careful analysis of the data shows that such a two island model is characterized by a very old split time that significantly predates the last glacial maximum followed by secondary contact with strong migration. We evaluate selection based on this demography and find that this secondary contact model strongly affects the power to detect sweeps. Moreover, it affects the power differently for northern Sweden (more false positives) as compared with southern Sweden (more false negatives). However, even when the demographic history is accounted for, sweep signals in northern Sweden are stronger than in southern Sweden, with little or no positional overlap. Further simulations including the complex demography and selection confirm that this is not compatible with global selection acting on both populations, and thus can be taken as evidence for local selection within subpopulations of Swedish A. thaliana . This study demonstrates the necessity of combining demographic analyses and sweep scans for the detection of selection, particularly when selection acts predominantly local.

Beschreibung: Members of a gene family expressed in a single species often experience common selection pressures. Consequently, the molecular basis of complex adaptations may be expected to involve parallel evolutionary changes in multiple paralogs. Here, we use bacterial artificial chromosome library scans to investigate the evolution of the voltage-gated sodium channel (Na v ) family in the garter snake Thamnophis sirtalis , a predator of highly toxic Taricha newts. Newts possess tetrodotoxin (TTX), which blocks Na v ’s, arresting action potentials in nerves and muscle. Some Thamnophis populations have evolved resistance to extremely high levels of TTX. Previous work has identified amino acid sites in the skeletal muscle sodium channel Na v 1.4 that confer resistance to TTX and vary across populations. We identify parallel evolution of TTX resistance in two additional Na v paralogs, Na v 1.6 and 1.7, which are known to be expressed in the peripheral nervous system and should thus be exposed to ingested TTX. Each paralog contains at least one TTX-resistant substitution identical to a substitution previously identified in Na v 1.4. These sites are fixed across populations, suggesting that the resistant peripheral nerves antedate resistant muscle. In contrast, three sodium channels expressed solely in the central nervous system (Na v 1.1–1.3) showed no evidence of TTX resistance, consistent with protection from toxins by the blood–brain barrier. We also report the exon–intron structure of six Na v paralogs, the first such analysis for snake genes. Our results demonstrate that the molecular basis of adaptation may be both repeatable across members of a gene family and predictable based on functional considerations.

Beschreibung: Overcoming the stress of starvation is one of an organism’s most challenging phenotypic responses. Those organisms that frequently survive the challenge, by virtue of their fitness, will have evolved genomes that are shaped by their specific environments. Understanding this genotype–environment–phenotype relationship at a deep level will require quantitative predictive models of the complex molecular systems that link these aspects of an organism’s existence. Here, we treat one of the most fundamental molecular systems, protein synthesis, and the amino acid biosynthetic pathways involved in the stringent response to starvation. These systems face an inherent logical dilemma: Building an amino acid biosynthetic pathway to synthesize its product—the cognate amino acid of the pathway—may require that very amino acid when it is no longer available. To study this potential "catch-22," we have created a generic model of amino acid biosynthesis in response to sudden starvation. Our mathematical analysis and computational results indicate that there are two distinctly different outcomes: Partial recovery to a new steady state, or full system failure. Moreover, the cell’s fate is dictated by the cognate bias, the number of cognate amino acids in the corresponding biosynthetic pathway relative to the average number of that amino acid in the proteome. We test these implications by analyzing the proteomes of over 1,800 sequenced microbes, which reveals statistically significant evidence of low cognate bias, a genetic trait that would avoid the biosynthetic quandary. Furthermore, these results suggest that the pattern of cognate bias, which is readily derived by genome sequencing, may provide evolutionary clues to an organism’s natural environment.

Beschreibung: Gene loss, gain, and transfer play an important role in shaping the genomes of all organisms; however, the interplay of these processes in isolated populations, such as in obligate intracellular bacteria, is less understood. Despite a general trend towards genome reduction in these microbes, our phylogenomic analysis of the phylum Chlamydiae revealed that within the family Parachlamydiaceae, gene family expansions have had pronounced effects on gene content. We discovered that the largest gene families within the phylum are the result of rapid gene birth-and-death evolution. These large gene families are comprised of members harboring eukaryotic-like ubiquitination-related domains, such as F-box and BTB-box domains, marking the largest reservoir of these proteins found among bacteria. A heterologous type III secretion system assay suggests that these proteins function as effectors manipulating the host cell. The large disparity in copy number of members in these families between closely related organisms suggests that nonadaptive processes might contribute to the evolution of these gene families. Gene birth-and-death evolution in concert with genomic drift might represent a previously undescribed mechanism by which isolated bacterial populations diversify.

Beschreibung: The Ryukyu Islands are located to the southwest of the Japanese archipelago. Archaeological evidence has revealed the existence of prehistoric cultural differentiation between the northern Ryukyu islands of Amami and Okinawa, and the southern Ryukyu islands of Miyako and Yaeyama. To examine a genetic subdivision in the Ryukyu Islands, we conducted genome-wide single nucleotide polymorphism typing of inhabitants from the Okinawa Islands, the Miyako Islands, and the Yaeyama Islands. Principal component and cluster analyses revealed genetic differentiation among the island groups, especially between Okinawa and Miyako. No genetic affinity was observed between aboriginal Taiwanese and any of the Ryukyu populations. The genetic differentiation observed between the inhabitants of the Okinawa Islands and the Miyako Islands is likely to have arisen due to genetic drift rather than admixture with people from neighboring regions. Based on the observed genetic differences, the divergence time between the inhabitants of Okinawa and Miyako islands was dated to the Holocene. These findings suggest that the Pleistocene inhabitants, whose bones have been found on the southern Ryukyu Islands, did not make a major genetic contribution, if any, to the present-day inhabitants of the southern Ryukyu Islands.

Beschreibung: Estimating phylogenetic trees from sequence data is an extremely challenging and important statistical task. Within the maximum-likelihood paradigm, the best tree is a point estimate. To determine how strongly the data support such an evolutionary scenario, a hypothesis testing methodology is required. To this end, the Kishino–Hasegawa (KH) test was developed to determine whether one topology is significantly more supported by the sequence data than another one. This test and its derivatives are widely used in phylogenetics and phylogenomics. Here, we show that the KH test is biased in the presence of alignment error and can lead to erroneous conclusions. Using simulations we demonstrated that due to alignment errors the KH test often rejects one of the competing topologies, even though both topologies are equally supported by the data. Specifically, we show that the KH test favors the guide tree used to align the analyzed sequences. Further, branch length optimization renders the test too conservative. We propose two possible corrections for these biases. First, we evaluated the impact of removing unreliable alignment columns and found out that it decreases the bias at the cost of substantially reducing the test’s power. Second, we developed a parametric test that entirely abolishes the biases without data filtering. This test incorporates the alignment construction step into the test’s hypothesis, thus removing the above guide tree effect. We extend this methodology for the case of multiple-topology comparisons and demonstrate the applicability of the new methodology on an exemplary data set.

Beschreibung: Studies of natural selection, followed by functional validation, are shedding light on understanding of genetic mechanisms underlying human evolution and adaptation. Classic methods for detecting selection, such as the integrated haplotype score (iHS) and Fay and Wu’s H statistic, are useful for candidate gene searching underlying positive selection. These methods, however, have limited capability to localize causal variants in selection target regions. In this study, we developed a novel method based on conditional coalescent tree to detect recent positive selection by counting unbalanced mutations on coalescent gene genealogies. Extensive simulation studies revealed that our method is more robust than many other approaches against biases due to various demographic effects, including population bottleneck, expansion, or stratification, while not sacrificing its power. Furthermore, our method demonstrated its superiority in localizing causal variants from massive linked genetic variants. The rate of successful localization was about 20–40% higher than that of other state-of-the-art methods on simulated data sets. On empirical data, validated functional causal variants of four well-known positive selected genes were all successfully localized by our method, such as ADH1B , MCM6 , APOL1 , and HBB . Finally, the computational efficiency of this new method was much higher than that of iHS implementations, that is, 24–66 times faster than the REHH package, and more than 10,000 times faster than the original iHS implementation. These magnitudes make our method suitable for applying on large sequencing data sets. Software can be downloaded from https://github.com/wavefancy/scct .

Beschreibung: Modern sequencing technology now allows biologists to collect the entirety of molecular evidence for reconstructing evolutionary trees. We introduce a novel, user-friendly software package engineered for conducting state-of-the-art Bayesian tree inferences on data sets of arbitrary size. Our software introduces a nonblocking parallelization of Metropolis-coupled chains, modifications for efficient analyses of data sets comprising thousands of partitions and memory saving techniques. We report on first experiences with Bayesian inferences at the whole-genome level using the SuperMUC supercomputer and simulated data.