Genomic Targets of Positive Selection in Giant Mice from Gough Island.

A key challenge in understanding how organisms adapt to their environments is to identify the mutations and genes that make it possible. By comparing patterns of sequence variation to neutral predictions across genomes, the targets of positive selection can be located. We applied this logic to house mice that invaded Gough Island (GI), an unusual population that shows phenotypic and ecological hallmarks of selection. We used massively parallel short-read sequencing to survey the genomes of 14 GI mice. We computed a set of summary statistics to capture diverse aspects of variation across these genome sequences, used approximate Bayesian computation to reconstruct a null demographic model, and then applied machine learning to estimate the posterior probability of positive selection in each region of the genome. Using a conservative threshold, 1,463 5-kb windows show strong evidence for positive selection in GI mice but not in a mainland reference population of German mice. Disproportionate shares of these selection windows contain genes that harbor derived nonsynonymous mutations with large frequency differences. Over-represented gene ontologies in selection windows emphasize neurological themes. Inspection of genomic regions harboring many selection windows with high posterior probabilities pointed to genes with known effects on exploratory behavior and body size as potential targets. Some genes in these regions contain candidate adaptive variants, including missense mutations and/or putative regulatory mutations. Our results provide a genomic portrait of adaptation to island conditions and position GI mice as a powerful system for understanding the genetic component of natural selection.

A first genetic portrait of synaptonemal complex variation.

The synaptonemal complex (SC) is a proteinaceous scaffold required for synapsis and recombination between homologous chromosomes during meiosis. Although the SC has been linked to differences in genome-wide crossover rates, the genetic basis of standing variation in SC structure remains unknown. To investigate the possibility that recombination evolves through changes to the SC, we characterized the genetic architecture of SC divergence on two evolutionary timescales. Applying a novel digital image analysis technique to spermatocyte spreads, we measured total SC length in 9,532 spermatocytes from recombinant offspring of wild-derived mouse strains with differences in this fundamental meiotic trait. Using this large dataset, we identified the first known genomic regions involved in the evolution of SC length. Distinct loci affect total SC length divergence between and within subspecies, with the X chromosome contributing to both. Joint genetic analysis of MLH1 foci-immunofluorescent markers of crossovers-from the same spermatocytes revealed that two of the identified loci also confer differences in the genome-wide recombination rate. Causal mediation analysis suggested that one pleiotropic locus acts early in meiosis to designate crossovers prior to SC assembly, whereas a second locus primarily shapes crossover number through its effect on SC length. One genomic interval shapes the relationship between SC length and recombination rate, likely modulating the strength of crossover interference. Our findings pinpoint SC formation as a key step in the evolution of recombination and demonstrate the power of genetic mapping on standing variation in the context of the recombination pathway.

A genomic portrait of human microsatellite variation.

Rapid advances in DNA sequencing and genotyping technologies are beginning to reveal the scope and pattern of human genomic variation. Although single nucleotide polymorphisms (SNPs) have been intensively studied, the extent and form of variation at other types of molecular variants remain poorly understood. Polymorphism at the most variable loci in the human genome, microsatellites, has rarely been examined on a genomic scale without the ascertainment biases that attend typical genotyping studies. We conducted a genomic survey of variation at microsatellites with at least three perfect repeats by comparing two complete genome sequences, the Human Genome Reference sequence and the sequence of J. Craig Venter. The genomic proportion of polymorphic loci was 2.7%, much higher than the rate of SNP variation, with marked heterogeneity among classes of loci. The proportion of variable loci increased substantially with repeat number. Repeat lengths differed in levels of variation, with longer repeat lengths generally showing higher polymorphism at the same repeat number. Microsatellite variation was weakly correlated with regional SNP number, indicating modest effects of shared genealogical history. Reductions in variation were detected at microsatellites located in introns, in untranslated regions, in coding exons, and just upstream of transcription start sites, suggesting the presence of selective constraints. Our results provide new insights into microsatellite mutational processes and yield a preview of patterns of variation that will be obtained in genomic surveys of larger numbers of individuals.

Acute phase cytotoxic T lymphocyte escape is a hallmark of simian immunodeficiency virus infection.

Cytotoxic T-lymphocyte (CTL) responses peak coincident with the decline in acute HIV viremia. Despite two reports of CTL-resistant HIV variants emerging during acute infection, the contribution of acute CTL escape to HIV pathogenesis remains unclear. Difficulties inherent in studying acute HIV infection can be overcome by modeling virus-host interactions in SIV-infected rhesus macaques. We sequenced 21 complete simian immunodeficiency virus (SIV)mac239 genomes at four weeks post-infection to determine the extent of acute CTL escape. Here we show that viruses from 19 of 21 macaques escaped from CTLs during acute infection and that these escape-selecting CTLs were responsive to lower concentrations of peptide than other SIV-specific CTLs. Interestingly, CTLs that require low peptide concentrations for stimulation (high 'functional avidity') are particularly effective at controlling other viral infections. Our results suggest that acute viral escape from CTLs is a hallmark of SIV infection and that CTLs with high functional avidity can rapidly select for escape variants.

A genomewide comparison of population structure at STRPs and nearby SNPs in humans.

Patterns of population structure provide insights into evolutionary processes and help identify groups of individuals for genotype-phenotype association studies. With increasing availability of polymorphic molecular markers across genomes, the examination of population structure using large numbers of unlinked loci has become a common practice in evolutionary biology and human genetics. The two classes of molecular variation most widely used for this purpose, short tandem repeat polymorphisms (STRPs) and single-nucleotide polymorphisms (SNPs), differ in mutational properties expected to affect population structure. To measure the relative ability of these loci to describe population structure, we compared diversity at neighboring STRPs and SNPs from 720 genomic regions in the four populations that comprise the Human HapMap. Comparing loci from the same genomic regions allowed us to focus on the contribution of mutational differences (rather than variation in genealogical history) to disparities in population structure between STRPs and SNPs. Relative to average values for SNPs from the same regions, STRPs had lower F(st), but higher G(st)' and I(n) values. STRP-SNP correlations in population structure across genomic regions were statistically significant but weak in magnitude. Separate analyses by repeat type showed that these correlations were driven primarily by tetranucleotide and trinucleotide STRPs; measures of population structure at dinucleotides and SNPs were not significantly correlated. Pairwise comparisons among populations revealed effects of divergence time on differences in population structure between STRPs and SNPs. Collectively, these results confirm that individual STRPs can provide more information about population structure than individual SNPs, but suggest that the difference in structure at STRPs and SNPs depends on local genealogical history. Our study motivates theoretical comparisons of population structure at loci with different mutational properties.

Weak Correlation between Nucleotide Variation and Recombination Rate across the House Mouse Genome.

Positive selection and purifying selection reduce levels of variation at linked neutral loci. One consequence of these processes is that the amount of neutral diversity and the meiotic recombination rate are predicted to be positively correlated across the genome-a prediction met in some species but not others. To better document the prevalence of selection at linked sites, we used new and published whole-genome sequences to survey nucleotide variation in population samples of the western European house mouse (Mus musculus domesticus) from Germany, France, and Gough Island, a remote volcanic island in the south Atlantic. Correlations between sequence variation and recombination rates estimated independently from dense linkage maps were consistently very weak (ρ ≤ 0.06), though they exceeded conventional significance thresholds. This pattern persisted in comparisons between genomic regions with the highest and lowest recombination rates, as well as in models incorporating the density of transcribed sites, the density of CpG dinucleotides, and divergence between mouse and rat as covariates. We conclude that natural selection affects linked neutral variation in a restricted manner in the western European house mouse.

Giant Island Mice Exhibit Widespread Gene Expression Changes in Key Metabolic Organs.

Island populations repeatedly evolve extreme body sizes, but the genomic basis of this pattern remains largely unknown. To understand how organisms on islands evolve gigantism, we compared genome-wide patterns of gene expression in Gough Island mice, the largest wild house mice in the world, and mainland mice from the WSB/EiJ wild-derived inbred strain. We used RNA-seq to quantify differential gene expression in three key metabolic organs: gonadal adipose depot, hypothalamus, and liver. Between 4,000 and 8,800 genes were significantly differentially expressed across the evaluated organs, representing between 20% and 50% of detected transcripts, with 20% or more of differentially expressed transcripts in each organ exhibiting expression fold changes of at least 2×. A minimum of 73 candidate genes for extreme size evolution, including Irs1 and Lrp1, were identified by considering differential expression jointly with other data sets: 1) genomic positions of published quantitative trait loci for body weight and growth rate, 2) whole-genome sequencing of 16 wild-caught Gough Island mice that revealed fixed single-nucleotide differences between the strains, and 3) publicly available tissue-specific regulatory elements. Additionally, patterns of differential expression across three time points in the liver revealed that Arid5b potentially regulates hundreds of genes. Functional enrichment analyses pointed to cell cycling, mitochondrial function, signaling pathways, inflammatory response, and nutrient metabolism as potential causes of weight accumulation in Gough Island mice. Collectively, our results indicate that extensive gene regulatory evolution in metabolic organs accompanied the rapid evolution of gigantism during the short time house mice have inhabited Gough Island.

Endoderm and pancreatic islet lineage differentiation from human embryonic stem cells.

Human embryonic stem cells (HESCs) are a potential source of insulin-producing tissue for transplantation. Recent studies have begun to define factors that promote definitive endoderm formation from HESCs, but conditions permitting complete islet specification in vitro have not been described. Here, we study spontaneous differentiation of HESCs to definitive endoderm and pancreatic progenitor cells, and begin to determine which aspects of the protocol are required for this cell fate commitment. HESCs were differentiated in culture for up to 10 weeks, including an embryoid body (EB) formation step. Modifications to the protocol included elimination of the EB phase, varying initial cell cluster size when forming EBs, and addition of mesoderm-derived cells to EBs. Differentiated cells were analyzed by reverse transcription-polymerase chain reaction (RT-PCR) and immunohistochemistry. HESCs are capable of spontaneous differentiation to cells expressing the definitive endoderm and pancreatic progenitor markers Foxa2, Sox17, and Pdx1, and ultimately, some cells express islet endocrine hormones. This differentiation occurs to a much greater extent when an EB formation step is included. Increased expression of endoderm markers during and after EB formation also correlated strongly with the size of cell clusters used to start EBs, as well as the addition of mesoderm- derived embryonic cells. This study demonstrates that a subset of differentiated HESC progeny adopt an endoderm fate and exhibit the capacity for further pancreatic lineage specification in vitro. Basal conditions were established for examining factors that can commit HESC-derived endoderm cells to specific pancreatic lineages.

Tat-vaccinated macaques do not control simian immunodeficiency virus SIVmac239 replication.

The regulatory proteins of human immunodeficiency virus may represent important vaccine targets. Here we assessed the role of Tat-specific cytotoxic T lymphocytes (CTL) in controlling pathogenic simian immunodeficiency virus SIVmac239 replication after using a DNA-prime, vaccinia virus Ankara-boost vaccine regimen. Despite the induction of Tat-specific CTL, there was no significant reduction in either peak or viral set point compared to that of controls.

Effects of cytotoxic T lymphocytes (CTL) directed against a single simian immunodeficiency virus (SIV) Gag CTL epitope on the course of SIVmac239 infection.

Vaccine-induced cytotoxic T lymphocytes (CTL) have been implicated in the control of virus replication in simian immunodeficiency virus (SIV)-challenged and simian-human immunodeficiency virus-challenged macaques. Therefore, we wanted to test the impact that vaccine-induced CTL responses against an immunodominant Gag epitope might have in the absence of other immune responses. By themselves, these strong CTL responses failed to control SIVmac239 replication.

Major histocompatibility complex class I alleles associated with slow simian immunodeficiency virus disease progression bind epitopes recognized by dominant acute-phase cytotoxic-T-lymphocyte responses.

Certain major histocompatibility complex class I (MHC-I) alleles are associated with delayed disease progression in individuals infected with human immunodeficiency virus (HIV) and in macaques infected with simian immunodeficiency virus (SIV). However, little is known about the influence of these MHC alleles on acute-phase cellular immune responses. Here we follow 51 animals infected with SIV(mac)239 and demonstrate a dramatic association between Mamu-A*01 and -B*17 expression and slowed disease progression. We show that the dominant acute-phase cytotoxic T lymphocyte (CTL) responses in animals expressing these alleles are largely directed against two epitopes restricted by Mamu-A*01 and one epitope restricted by Mamu-B*17. One Mamu-A*01-restricted response (Tat(28-35)SL8) and the Mamu-B*17-restricted response (Nef(165-173)IW9) typically select for viral escape variants in early SIV(mac)239 infection. Interestingly, animals expressing Mamu-A*1 and -B*17 have less variation in the Tat(28-35)SL8 epitope during chronic infection than animals that express only Mamu-A*01. Our results show that MHC-I alleles that are associated with slow progression to AIDS bind epitopes recognized by dominant CTL responses during acute infection and underscore the importance of understanding CTL responses during primary HIV infection.