Trans genomic capture and sequencing of primate exomes reveals new targets of positive selection

Patterns of positive selection in seven ant genomes

The evolution of ants is marked by remarkable adaptations that allowed the development of very complex social systems. To identify how ant-specific adaptations are associated with patterns of molecular evolution, we searched for signs of positive selection on amino-acid changes in proteins. We identified 24 functional categories of genes which were enriched for positively selected genes in the ant lineage. We also reanalyzed genome-wide data sets in bees and flies with the same methodology to check whether positive selection was specific to ants or also present in other insects. Notably, genes implicated in immunity were enriched for positively selected genes in the three lineages, ruling out the hypothesis that the evolution of hygienic behaviors in social insects caused a major relaxation of selective pressure on immune genes. Our scan also indicated that genes implicated in neurogenesis and olfaction started to undergo increased positive selection before the evolution of sociality in Hymenoptera. Finally, the comparison between these three lineages allowed us to pinpoint molecular evolution patterns that were specific to the ant lineage. In particular, there was ant-specific recurrent positive selection on genes with mitochondrial functions, suggesting that mitochondrial activity was improved during the evolution of this lineage. This might have been an important step toward the evolution of extreme lifespan that is a hallmark of ants.

Microcephaly genes evolved adaptively throughout the evolution of eutherian mammals

BACKGROUND

Genes associated with the neurodevelopmental disorder microcephaly display a strong signature of adaptive evolution in primates. Comparative data suggest a link between selection on some of these loci and the evolution of primate brain size. Whether or not either positive selection or this phenotypic association are unique to primates is unclear, but recent studies in cetaceans suggest at least two microcephaly genes evolved adaptively in other large brained mammalian clades.

RESULTS

Here we analyse the evolution of seven microcephaly loci, including three recently identified loci, across 33 eutherian mammals. We find extensive evidence for positive selection having acted on the majority of these loci not just in primates but also across non-primate mammals. Furthermore, the patterns of selection in major mammalian clades are not significantly different. Using phylogenetically corrected comparative analyses, we find that the evolution of two microcephaly loci, ASPM and CDK5RAP2, are correlated with neonatal brain size in Glires and Euungulata, the two most densely sampled non-primate clades.

CONCLUSIONS

Together with previous results, this suggests that ASPM and CDK5RAP2 may have had a consistent role in the evolution of brain size in mammals. Nevertheless, several limitations of currently available data and gene-phenotype tests are discussed, including sparse sampling across large evolutionary distances, averaging gene-wide rates of evolution, potential phenotypic variation and evolutionary reversals. We discuss the implications of our results for studies of the genetic basis of brain evolution, and explicit tests of gene-phenotype hypotheses.

Comparative primate genomics: emerging patterns of genome content and dynamics

Advances in genome sequencing technologies have created new opportunities for comparative primate genomics. Genome assemblies have been published for various primate species, and analyses of several others are underway. Whole-genome assemblies for the great apes provide remarkable new information about the evolutionary origins of the human genome and the processes involved. Genomic data for macaques and other non-human primates offer valuable insights into genetic similarities and differences among species that are used as models for disease-related research. This Review summarizes current knowledge regarding primate genome content and dynamics, and proposes a series of goals for the near future.

Selection on plant male function genes identifies candidates for reproductive isolation of yellow monkeyflowers

Understanding the genetic basis of reproductive isolation promises insight into speciation and the origins of biological diversity. While progress has been made in identifying genes underlying barriers to reproduction that function after fertilization (post-zygotic isolation), we know much less about earlier acting pre-zygotic barriers. Of particular interest are barriers involved in mating and fertilization that can evolve extremely rapidly under sexual selection, suggesting they may play a prominent role in the initial stages of reproductive isolation. A significant challenge to the field of speciation genetics is developing new approaches for identification of candidate genes underlying these barriers, particularly among non-traditional model systems. We employ powerful proteomic and genomic strategies to study the genetic basis of conspecific pollen precedence, an important component of pre-zygotic reproductive isolation among yellow monkeyflowers (Mimulus spp.) resulting from male pollen competition. We use isotopic labeling in combination with shotgun proteomics to identify more than 2,000 male function (pollen tube) proteins within maternal reproductive structures (styles) of M. guttatus flowers where pollen competition occurs. We then sequence array-captured pollen tube exomes from a large outcrossing population of M. guttatus, and identify those genes with evidence of selective sweeps or balancing selection consistent with their role in pollen competition. We also test for evidence of positive selection on these genes more broadly across yellow monkeyflowers, because a signal of adaptive divergence is a common feature of genes causing reproductive isolation. Together the molecular evolution studies identify 159 pollen tube proteins that are candidate genes for conspecific pollen precedence. Our work demonstrates how powerful proteomic and genomic tools can be readily adapted to non-traditional model systems, allowing for genome-wide screens towards the goal of identifying the molecular basis of genetically complex traits.

Human testis-specific genes are under relaxed negative selection

Recent studies have suggested that selective forces and constraints acting on genes varied during human evolution depending on the organ in which they are expressed. To gain insight into the evolution of organ determined negative selection forces, we compared the non-synonymous SNP diversity of genes expressed in different organs. Based on a HAPMAP dataset, we determined for each SNP its frequency in 11 human populations and, in each case, predicted whether or not the change it produces is deleterious. We have shown that, for all organs under study, SNPs predicted to be deleterious are present at a significantly lower frequency than SNPs predicted to be tolerated. However, testis-specific genes contain a higher proportion of deleterious SNPs than other organs. This study shows that negative selection is acting on the whole human genome, but that the action of negative selection is relaxed on testis-specific genes. This result adds to and expands the hypothesis of a recent evolutionary change in the human male reproductive system and its behavior.

A model-based analysis of GC-biased gene conversion in the human and chimpanzee genomes

GC-biased gene conversion (gBGC) is a recombination-associated process that favors the fixation of G/C alleles over A/T alleles. In mammals, gBGC is hypothesized to contribute to variation in GC content, rapidly evolving sequences, and the fixation of deleterious mutations, but its prevalence and general functional consequences remain poorly understood. gBGC is difficult to incorporate into models of molecular evolution and so far has primarily been studied using summary statistics from genomic comparisons. Here, we introduce a new probabilistic model that captures the joint effects of natural selection and gBGC on nucleotide substitution patterns, while allowing for correlations along the genome in these effects. We implemented our model in a computer program, called phastBias, that can accurately detect gBGC tracts about 1 kilobase or longer in simulated sequence alignments. When applied to real primate genome sequences, phastBias predicts gBGC tracts that cover roughly 0.3% of the human and chimpanzee genomes and account for 1.2% of human-chimpanzee nucleotide differences. These tracts fall in clusters, particularly in subtelomeric regions; they are enriched for recombination hotspots and fast-evolving sequences; and they display an ongoing fixation preference for G and C alleles. They are also significantly enriched for disease-associated polymorphisms, suggesting that they contribute to the fixation of deleterious alleles. The gBGC tracts provide a unique window into historical recombination processes along the human and chimpanzee lineages. They supply additional evidence of long-term conservation of megabase-scale recombination rates accompanied by rapid turnover of hotspots. Together, these findings shed new light on the evolutionary, functional, and disease implications of gBGC. The phastBias program and our predicted tracts are freely available.

Interpreting patterns of DNA sequence variation in the genomes of closely related species is critically important for understanding the causes and functional effects of nucleotide substitutions. Classical models describe patterns of substitution in terms of the fundamental forces of mutation, recombination, neutral drift, and natural selection. However, an entirely separate force, called GC-biased gene conversion (gBGC), also appears to have an important influence on substitution patterns in many species. gBGC is a recombination-associated evolutionary process that favors the fixation of strong (G/C) over weak (A/T) alleles. In mammals, gBGC is thought to promote variation in GC content, rapidly evolving sequences, and the fixation of deleterious mutations. However, its genome-wide influence remains poorly understood, in part because, it is difficult to incorporate gBGC into statistical models of evolution. In this paper, we describe a new evolutionary model that jointly describes the effects of selection and gBGC and apply it to the human and chimpanzee genomes. Our genome-wide predictions of gBGC tracts indicate that gBGC has been an important force in recent human evolution. Our publicly available computer program, called phastBias, and our genome-wide predictions will enable other researchers to consider gBGC in their analyses.

Targeted enrichment: maximizing orthologous gene comparisons across deep evolutionary time

Estimated phylogenies of evolutionarily diverse taxa will be well supported and more likely to be historically accurate when the analysis contains large amounts of data–many genes sequenced across many taxa. Inferring such phylogenies for non-model organisms is challenging given limited resources for whole-genome sequencing. We take advantage of genomic data from a single species to test the limits of hybridization-based enrichment of hundreds of exons across frog species that diverged up to 250 million years ago. Enrichment success for a given species depends greatly on the divergence time between it and the reference species, and the resulting alignment contains a significant proportion of missing data. However, our alignment generates a well-supported phylogeny of frogs, suggesting that this technique is a practical solution towards resolving relationships across deep evolutionary time.

Evidence for polygenic adaptation to pathogens in the human genome

Most approaches aiming at finding genes involved in adaptive events have focused on the detection of outlier loci, which resulted in the discovery of individually "significant" genes with strong effects. However, a collection of small effect mutations could have a large effect on a given biological pathway that includes many genes, and such a polygenic mode of adaptation has not been systematically investigated in humans. We propose here to evidence polygenic selection by detecting signals of adaptation at the pathway or gene set level instead of analyzing single independent genes. Using a gene-set enrichment test to identify genome-wide signals of adaptation among human populations, we find that most pathways globally enriched for signals of positive selection are either directly or indirectly involved in immune response. We also find evidence for long-distance genotypic linkage disequilibrium, suggesting functional epistatic interactions between members of the same pathway. Our results show that past interactions with pathogens have elicited widespread and coordinated genomic responses, and suggest that adaptation to pathogens can be considered as a primary example of polygenic selection.

Social environment influences the relationship between genotype and gene expression in wild baboons

Variation in the social environment can have profound effects on survival and reproduction in wild social mammals. However, we know little about the degree to which these effects are influenced by genetic differences among individuals, and conversely, the degree to which social environmental variation mediates genetic reaction norms. To better understand these relationships, we investigated the potential for dominance rank, social connectedness and group size to modify the effects of genetic variation on gene expression in the wild baboons of the Amboseli basin. We found evidence for a number of gene-environment interactions (GEIs) associated with variation in the social environment, encompassing social environments experienced in adulthood as well as persistent effects of early life social environment. Social connectedness, maternal dominance rank and group size all interacted with genotype to influence gene expression in at least one sex, and either in early life or in adulthood. These results suggest that social and behavioural variation, akin to other factors such as age and sex, can impact the genotype-phenotype relationship. We conclude that GEIs mediated by the social environment are important in the evolution and maintenance of individual differences in wild social mammals, including individual differences in responses to social stressors.

Dissecting structural and nucleotide genome-wide variation in inbred Iberian pigs

BACKGROUND

In contrast to international pig breeds, the Iberian breed has not been admixed with Asian germplasm. This makes it an important model to study both domestication and relevance of Asian genes in the pig. Besides, Iberian pigs exhibit high meat quality as well as appetite and propensity to obesity. Here we provide a genome wide analysis of nucleotide and structural diversity in a reduced representation library from a pool (n=9 sows) and shotgun genomic sequence from a single sow of the highly inbred Guadyerbas strain. In the pool, we applied newly developed tools to account for the peculiarities of these data.

RESULTS

A total of 254,106 SNPs in the pool (79.6 Mb covered) and 643,783 in the Guadyerbas sow (1.47 Gb covered) were called. The nucleotide diversity (1.31x10-3 per bp in autosomes) is very similar to that reported in wild boar. A much lower than expected diversity in the X chromosome was confirmed (1.79x10-4 per bp in the individual and 5.83x10-4 per bp in the pool). A strong (0.70) correlation between recombination and variability was observed, but not with gene density or GC content. Multicopy regions affected about 4% of annotated pig genes in their entirety, and 2% of the genes partially. Genes within the lowest variability windows comprised interferon genes and, in chromosome X, genes involved in behavior like HTR2C or MCEP2. A modified Hudson-Kreitman-Aguadé test for pools also indicated an accelerated evolution in genes involved in behavior, as well as in spermatogenesis and in lipid metabolism.

CONCLUSIONS

This work illustrates the strength of current sequencing technologies to picture a comprehensive landscape of variability in livestock species, and to pinpoint regions containing genes potentially under selection. Among those genes, we report genes involved in behavior, including feeding behavior, and lipid metabolism. The pig X chromosome is an outlier in terms of nucleotide diversity, which suggests selective constraints. Our data further confirm the importance of structural variation in the species, including Iberian pigs, and allowed us to identify new paralogs for known gene families.

Duplicate abalone egg coat proteins bind sperm lysin similarly, but evolve oppositely, consistent with molecular mimicry at fertilization

Sperm and egg proteins constitute a remarkable paradigm in evolutionary biology: despite their fundamental role in mediating fertilization (suggesting stasis), some of these molecules are among the most rapidly evolving ones known, and their divergence can lead to reproductive isolation. Because of strong selection to maintain function among interbreeding individuals, interacting fertilization proteins should also exhibit a strong signal of correlated divergence among closely related species. We use evidence of such molecular co-evolution to target biochemical studies of fertilization in North Pacific abalone (Haliotis spp.), a model system of reproductive protein evolution. We test the evolutionary rates (d_N/d_S) of abalone sperm lysin and two duplicated egg coat proteins (VERL and VEZP14), and find a signal of co-evolution specific to ZP-N, a putative sperm binding motif previously identified by homology modeling. Positively selected residues in VERL and VEZP14 occur on the same face of the structural model, suggesting a common mode of interaction with sperm lysin. We test this computational prediction biochemically, confirming that the ZP-N motif is sufficient to bind lysin and that the affinities of VERL and VEZP14 are comparable. However, we also find that on phylogenetic lineages where lysin and VERL evolve rapidly, VEZP14 evolves slowly, and vice versa. We describe a model of sexual conflict that can recreate this pattern of anti-correlated evolution by assuming that VEZP14 acts as a VERL mimic, reducing the intensity of sexual conflict and slowing the co-evolution of lysin and VERL.

Interacting sperm and egg proteins must co-evolve to maintain compatibility at fertilization, so their divergence among species should be correlated—lineages with rapidly evolving sperm proteins should have rapidly evolving egg proteins. We use this expectation to target biochemical studies of fertilization in a model system (abalone). We study a discrete functional domain (ZP-N) found in a pair of duplicated egg coat proteins, and we find the ZP-N motif from both proteins bind sperm lysin (a protein important for sperm passage of the egg coat) in a similar fashion. ZP-N is a feature of vertebrate and invertebrate egg coat proteins, as well as yeast mating recognition proteins, demonstrating its broad significance in sexual reproduction. Unexpectedly, we find that the ZP-N motifs of VEZP14 and VERL exhibit inverse patterns of co-evolution with lysin, suggesting that these duplicates may have opposite functions in fertilization. Using computer simulations, we model a novel explanation for this pattern whereby VEZP14 acts as a decoy of VERL in order to decrease the effective amount of sperm lysin and slow the rate of fertilization. Such molecular mimicry could complement other well-established fertilization blocks that females use to control rates of fertilization and limit polyspermy.

Comparative population genomics of the ejaculate in humans and the great apes

The rapid molecular evolution of reproductive genes is nearly ubiquitous across animals, yet the selective forces and functional targets underlying this divergence remain poorly understood. Humans and closely related species of great apes show strongly divergent mating systems, providing a powerful system to investigate the influence of sperm competition on the evolution of reproductive genes. This is complemented by detailed information on male reproductive biology and unparalleled genomic resources in humans. Here, we have used custom microarrays to capture and sequence 285 genes encoding proteins present in the ejaculate as well as 101 randomly selected control genes in 21 gorillas, 20 chimpanzees, 20 bonobos, and 20 humans. In total, we have generated >25× average genomic coverage per individual for over 1 million target base pairs. Our analyses indicate high levels of evolutionary constraint across much of the ejaculate combined with more rapid evolution of genes involved in immune defense and proteolysis. We do not find evidence for appreciably more positive selection along the lineage leading to bonobos and chimpanzees, although this would be predicted given more intense sperm competition in these species. Rather, the extent of positive and negative selection depended more on the effective population sizes of the species. Thus, general patterns of male reproductive protein evolution among apes and humans depend strongly on gene function but not on inferred differences in the intensity of sperm competition among extant species.

Protein evolution in two co-occurring types of Symbiodinium: an exploration into the genetic basis of thermal tolerance in Symbiodinium clade D

Background

The symbiosis between reef-building corals and photosynthetic dinoflagellates (Symbiodinium) is an integral part of the coral reef ecosystem, as corals are dependent on Symbiodinium for the majority of their energy needs. However, this partnership is increasingly at risk due to changing climatic conditions. It is thought that functional diversity within Symbiodinium may allow some corals to rapidly adapt to different environments by changing the type of Symbiodinium with which they partner; however, very little is known about the molecular basis of the functional differences among symbiont groups. One group of Symbiodinium that is hypothesized to be important for the future of reefs is clade D, which, in general, seems to provide the coral holobiont (i.e., coral host and associated symbiont community) with elevated thermal tolerance. Using high-throughput sequencing data from field-collected corals we assembled, de novo, draft transcriptomes for Symbiodinium clades C and D. We then explore the functional basis of thermal tolerance in clade D by comparing rates of coding sequence evolution among the four clades of Symbiodinium most commonly found in reef-building corals (A-D).

Results

We are able to highlight a number of genes and functional categories as candidates for involvement in the increased thermal tolerance of clade D. These include a fatty acid desaturase, molecular chaperones and proteins involved in photosynthesis and the thylakoid membrane. We also demonstrate that clades C and D co-occur within most of the sampled colonies of Acropora hyacinthus, suggesting widespread potential for this coral species to acclimatize to changing thermal conditions via ‘shuffling’ the proportions of these two clades from within their current symbiont communities.

Conclusions

Transcriptome-wide analysis confirms that the four main Symbiodinium clades found within corals exhibit extensive evolutionary divergence (18.5-27.3% avg. pairwise nucleotide difference). Despite these evolutionary distinctions, many corals appear to host multiple clades simultaneously, which may allow for rapid acclimatization to changing environmental conditions. This study provides a first step toward understanding the molecular basis of functional differences between Symbiodinium clades by highlighting a number of genes with signatures consistent with positive selection along the thermally tolerant clade D lineage.

Positive selection in the major histocompatibility complex class III region of cynomolgus macaques (Macaca fascicularis) of the Philippines origin

The cynomolgus macaque (Macaca fascicularis) is a model of choice among primates for the study of local adaptation processes because of its mixed and wide insular and continental distribution. In a previous study, by using 12 markers [5 microsatellites located in the major histocompatibility complex (MHC) region and 7 outside MHC], we have detected a signal of positive selection on the microsatellite DRACA located inside the gene DRA. In order to refine the location of this signal of positive selection in the MHC region, we studied the genetic diversity of 36 markers (18 microsatellites spread across the MHC region and 18 autosomal microsatellites outside MHC) in a sample of 254 individuals from four populations (Vietnam, Java, the Philippines, and Mauritius). We estimated for each locus the deviation of F(st) from a neutral model by using two methods based on contrasted demographic models. The two approaches showed a signal of positive selection in the MHC class III region that is much more significant than the one previously reported for the marker DRACA which could have been influenced by a hitchhiking effect due to its proximity with the class III region.

Exome RNA sequencing reveals rare and novel alternative transcripts

RNA sequencing has become an important method to perform hypothesis-free characterization of global gene expression. One of the limitations of RNA sequencing is that most sequence reads represent highly expressed transcripts, whereas low level transcripts are challenging to detect. To combine the benefits of traditional expression arrays with the advantages of RNA sequencing, we have used whole exome enrichment prior to sequencing of total RNA. We show that whole exome capture can be successfully applied to cDNA to study the transcriptional landscape in human tissues. By introducing the exome enrichment step, we are able to identify transcripts present at very low levels, which are below the level of detection in conventional RNA sequencing. Although the enrichment increases the ability to detect presence of transcripts, it also lowers the accuracy of quantification of expression levels. Our results yield a large number of novel exons and splice isoforms, suggesting that conventional RNA sequencing methods only detect a small fraction of the full transcript diversity. We propose that whole exome enrichment of RNA is a suitable strategy for genome-wide discovery of novel transcripts, alternative splice variants and fusion genes.

Transcriptome-based exon capture enables highly cost-effective comparative genomic data collection at moderate evolutionary scales

BACKGROUND

To date, exon capture has largely been restricted to species with fully sequenced genomes, which has precluded its application to lineages that lack high quality genomic resources. We developed a novel strategy for designing array-based exon capture in chipmunks (Tamias) based on de novo transcriptome assemblies. We evaluated the performance of our approach across specimens from four chipmunk species.

RESULTS

We selectively targeted 11,975 exons (~4 Mb) on custom capture arrays, and enriched over 99% of the targets in all libraries. The percentage of aligned reads was highly consistent (24.4-29.1%) across all specimens, including in multiplexing up to 20 barcoded individuals on a single array. Base coverage among specimens and within targets in each species library was uniform, and the performance of targets among independent exon captures was highly reproducible. There was no decrease in coverage among chipmunk species, which showed up to 1.5% sequence divergence in coding regions. We did observe a decline in capture performance of a subset of targets designed from a much more divergent ground squirrel genome (30 My), however, over 90% of the targets were also recovered. Final assemblies yielded over ten thousand orthologous loci (~3.6 Mb) with thousands of fixed and polymorphic SNPs among species identified.

CONCLUSIONS

Our study demonstrates the potential of a transcriptome-enabled, multiplexed, exon capture method to create thousands of informative markers for population genomic and phylogenetic studies in non-model species across the tree of life.

Targeted capture of homoeologous coding and noncoding sequence in polyploid cotton

Targeted sequence capture is a promising technology in many areas in biology. These methods enable efficient and relatively inexpensive sequencing of hundreds to thousands of genes or genomic regions from many more individuals than is practical using whole-genome sequencing approaches. Here, we demonstrate the feasibility of target enrichment using sequence capture in polyploid cotton. To capture and sequence both members of each gene pair (homeologs) of wild and domesticated Gossypium hirsutum, we created custom hybridization probes to target 1000 genes (500 pairs of homeologs) using information from the cotton transcriptome. Two widely divergent samples of G. hirsutum were hybridized to four custom NimbleGen capture arrays containing probes for targeted genes. We show that the two coresident homeologs in the allopolyploid nucleus were efficiently captured with high coverage. The capture efficiency was similar between the two accessions and independent of whether the samples were multiplexed. A significant amount of flanking, nontargeted sequence (untranslated regions and introns) was also captured and sequenced along with the targeted exons. Intraindividual heterozygosity is low in both wild and cultivated Upland cotton, as expected from the high level of inbreeding in natural G. hirsutum and bottlenecks accompanying domestication. In addition, levels of heterozygosity appeared asymmetrical with respect to genome (A_T or D_T) in cultivated cotton. The approach used here is general, scalable, and may be adapted for many different research inquiries involving polyploid plant genomes.

Targeted sequence capture as a powerful tool for evolutionary analysis

Next-generation sequencing technologies (NGS) have revolutionized biological research by significantly increasing data generation while simultaneously decreasing the time to data output. For many ecologists and evolutionary biologists, the research opportunities afforded by NGS are substantial; even for taxa lacking genomic resources, large-scale genome-level questions can now be addressed, opening up many new avenues of research. While rapid and massive sequencing afforded by NGS increases the scope and scale of many research objectives, whole genome sequencing is often unwarranted and unnecessarily complex for specific research questions. Recently developed targeted sequence enrichment, coupled with NGS, represents a beneficial strategy for enhancing data generation to answer questions in ecology and evolutionary biology. This marriage of technologies offers researchers a simple method to isolate and analyze a few to hundreds, or even thousands, of genes or genomic regions from few to many samples in a relatively efficient and effective manner. These strategies can be applied to questions at both the infra- and interspecific levels, including those involving parentage, gene flow, divergence, phylogenetics, reticulate evolution, and many more. Here we provide a brief overview of targeted sequence enrichment, and emphasize the power of this technology to increase our ability to address a wide range of questions of interest to ecologists and evolutionary biologists, particularly for those working with taxa for which few genomic resources are available.