dbSNP: the NCBI database of genetic variation

In response to a need for a general catalog of genome variation to address the large-scale sampling designs required by association studies, gene mapping and evolutionary biology, the National Center for Biotechnology Information (NCBI) has established the dbSNP database [S.T.Sherry, M.Ward and K.Sirotkin (1999) Genome Res., 9, 677–679]. Submissions to dbSNP will be integrated with other sources of information at NCBI such as GenBank, PubMed, LocusLink and the Human Genome Project data. The complete contents of dbSNP are available to the public at website: http://www.ncbi.nlm.nih.gov/SNP. The complete contents of dbSNP can also be downloaded in multiple formats via anonymous FTP at ftp://ncbi.nlm.nih.gov/snp/.

Global survey of genetic variation in CCR5, RANTES, and MIP-1α: Impact on the epidemiology of the HIV-1 pandemic

Expression of CC chemokine receptor 5 (CCR5), the major coreceptor for HIV-1 cell entry, and its ligands (e.g., RANTES and MIP-1α) is widely regarded as central to the pathogenesis of HIV-1 infection. By surveying nearly 3,000 HIV+ and HIV− individuals from worldwide populations for polymorphisms in the genes encoding RANTES, MIP-1α, and CCR5, we show that the evolutionary histories of human populations have had a significant impact on the distribution of variation in these genes, and that this may be responsible, in part, for the heterogeneous nature of the epidemiology of the HIV-1 pandemic. The varied distribution of RANTES haplotypes (AC, GC, and AG) associated with population-specific HIV-1 transmission- and disease-modifying effects is a striking example. Homozygosity for the AC haplotype was associated with an increased risk of acquiring HIV-1 as well as accelerated disease progression in European Americans, but not in African Americans. Yet, the prevalence of the ancestral AC haplotype is high in individuals of African origin, but substantially lower in non-Africans. In a Japanese cohort, AG-containing RANTES haplotype pairs were associated with a delay in disease progression; however, we now show that their contribution to HIV-1 pathogenesis and epidemiology in other parts of the world is negligible because the AG haplotype is infrequent in non-Far East Asians. Thus, the varied distribution of RANTES, MIP-1α, and CCR5 haplotype pairs and their population-specific phenotypic effects on HIV-1 susceptibility and disease progression results in a complex pattern of biological determinants of HIV-1 epidemiology. These findings have important implications for the design, assessment, and implementation of effective HIV-1 intervention and prevention strategies.

Genetic variation in vulnerability to the behavioral effects of neonatal hippocampal damage in rats.

We explored how two independent variables, one genetic (i.e., specific rat strains) and another environmental (i.e., a developmental excitotoxic hippocampal lesion), contribute to phenotypic variation. Sprague-Dawley (SD), Fischer 344 (F344), and Lewis rats underwent two grades of neonatal excitotoxic damage: small and large ventral hippocampal (SVH and LVH) lesions. Locomotion was tested before puberty [postnatal day 35 (P35)] and after puberty (P56) following exposure to a novel environment or administration of amphetamine. The behavioral effects were strain- and lesion-specific. As shown previously, SD rats with LVH lesions displayed enhanced spontaneous and amphetamine-induced locomotion as compared with controls at P56, but not at P35. SVH lesions in SD rats had no effect at any age. In F344 rats with LVH lesions, enhanced spontaneous and amphetamine-induced locomotion appeared early (P35) and was exaggerated at P56. SVH lesions in F344 rats resulted in a pattern of effects analogous to LVH lesions in SD rats--i.e., postpubertal onset of hyperlocomotion (P56). In Lewis rats, LVH lesions had no significant effect on novelty- or amphetamine-induced locomotion at any age. These data show that the degree of genetic predisposition and the extent of early induced hippocampal defect contribute to the particular pattern of behavioral outcome. These results may have implications for modeling interactions of genetic and environmental factors involved in schizophrenia, a disorder characterized by phenotypic heterogeneity, genetic predisposition, a developmental hippocampal abnormality, and vulnerability to environmental stress.

Hospitalization risk due to respiratory illness associated with genetic variation at IFITM3 in patients with influenza A(H1N1)pdm09 infection: a case-control study

Recent studies suggest an association between the Interferon Inducible Transmembrane 3 (IFITM3) rs12252 variant and the course of influenza infection. However, it is not clear whether the reported association relates to influenza infection severity. The aim of this study was to estimate the hospitalization risk associated with this variant in Influenza Like Illness (ILI) patients during the H1N1 pandemic influenza. A case-control genetic association study was performed, using nasopharyngeal/oropharyngeal swabs collected during the H1N1 pandemic influenza. Laboratory diagnosis of influenza infection was performed by RT-PCR, the IFITM3 rs12252 was genotyped by RFLP and tested for association with hospitalization. Conditional logistic regression was performed to calculate the confounder-adjusted odds ratio of hospitalization associated with IFITM3 rs12252. We selected 312 ILI cases and 624 matched non-hospitalized controls. Within ILI Influenza A(H1N1)pdm09 positive patients, no statistical significant association was found between the variant and the hospitalization risk (Adjusted OR: 0.73 (95%CI: 0.33–1.50)). Regarding ILI Influenza A(H1N1)pdm09 negative patients, CT/CC genotype carriers had a higher risk of being hospitalized than patients with TT genotype (Adjusted OR: 2.54 (95%CI: 1.54–4.19)). The IFITM3 rs12252 variant was associated with respiratory infection hospitalization but not specifically in patients infected with Influenza A(H1N1)pdm09.

Genetic variation and inference of demographic histories in non-model species

Both long-term environmental changes such as those driven by the glacial cycles and more recent anthropogenic impacts have had major effects on the past demography in wild organisms. Within species, these changes are reflected in the amount and distribution of neutral genetic variation. In this thesis, mitochondrial and microsatellite DNA was analysed to investigate how environmental and anthropogenic factors have affected genetic diversity and structure in four ecologically different animal species. Paper I describes the post-glacial recolonisation history of the speckled-wood butterfly (Pararge aegeria) in Northern Europe. A decrease in genetic diversity with latitude and a marked population structure were uncovered, consistent with a hypothesis of repeated founder events during the postglacial recolonisation. Moreover, Approximate Bayesian Computation analyses indicate that the univoltine populations in Scandinavia and Finland originate from recolonisations along two routes, one on each side of the Baltic. Paper II aimed to investigate how past sea-level rises affected the population history of the convict surgeonfish (Acanthurus triostegus) in the Indo-Pacific. Assessment of the species’ demographic history suggested a population expansion that occurred approximately at the end of the last glaciation. Moreover, the results demonstrated an overall lack of phylogeographic structure, probably due to the high dispersal rates associated with the species’ pelagic larval stage. Populations at the species’ eastern range margin were significantly differentiated from other populations, which likely is a consequence of their geographic isolation. In Paper III, we assessed the effect of human impact on the genetic variation of European moose (Alces alces) in Sweden. Genetic analyses revealed a spatial structure with two genetic clusters, one in northern and one in southern Sweden, which were separated by a narrow transition zone. Moreover, demographic inference suggested a recent population bottleneck. The inferred timing of this bottleneck coincided with a known reduction in population size in the 19th and early 20th century due to high hunting pressure. In Paper IV, we examined the effect of an indirect but well-described human impact, via environmental toxic chemicals (PCBs), on the genetic variation of Eurasian otters (Lutra lutra) in Sweden. Genetic clustering assignment revealed differentiation between otters in northern and southern Sweden, but also in the Stockholm region. ABC analyses indicated a decrease in effective population size in both northern and southern Sweden. Moreover, comparative analyses of historical and contemporary samples demonstrated a more severe decline in genetic diversity in southern Sweden compared to northern Sweden, in agreement with the levels of PCBs found.

Genetic variation of the scleractinian coral Stylophora pistillata, from western Pacific reefs

Most scleractinian coral species are widely distributed across the tropical and subtropical Indo-Pacific. However, the genetic connectivity between populations of corals separated by large distances (thousands of kilometers) is not well known. We analyzed variability in the nucleotide sequence of the internal transcribed spacer-1 (ITS-1) of the nuclear ribosomal gene unit in the ubiquitous coral Stylophora pistillata, across the western Pacific Ocean. Eight populations from Japan, Malaysia, and the northern and southern Great Barrier Reef (GBR) were studied. Phylogenetic analyses and analysis of molecular variance (AMOVA) clearly revealed that there is panmixia among these coral populations. AMOVA showed that ITS-1 sequence variability was greater within populations (78.37%) than among populations (12.06%). These patterns strongly suggest high levels of connectivity across the species' latitudinal distribution range in the western Pacific, as is seen in many marine invertebrates.

Genetic variation for carbon isotope discrimination in sunflower: Association with transpiration efficiency and evidence for cytoplasmic inheritance

Plants accumulate isotopes of carbon at different rates because of discrimination against C-13 relative to C-12. In plants that fix carbon by the C-3 pathway, the amount of discrimination correlates negatively with transpiration efficiency (TE) where TE is the amount of dry matter accumulated per unit water transpired. Therefore, carbon isotope discrimination (Delta) has become a useful tool for selecting genotypes with improved TE and performance in dry environments. Surveys of 161 sunflower (Helianthus spp.) genotypes of diverse origin revealed a large and unprecedented range of genetic variation for Delta (19.5-23.8parts per thousand). A strong negative genetic correlation (r(g)) between TE and Delta (r(g) = -0.87, P < 0.001) was observed in glasshouse studies. Gas exchange measurements of field grown plants indicated that Delta was strongly correlated with stomatal conductance to water vapor (g), (r(g) 0.64, P < 0.01), and the ratio of net assimilation rate (A) to g, (r(g) = 0.86, P < 0.001), an instantaneous measure of TE. Genotype CMSHA89MAX1 had the lowest TE (and highest Delta) of all genotypes tested in these studies and low yields in hybrid combination. Backcrossing studies showed that the TE of this genotype was due to an adverse effect of the MAX1 cytoplasm, which was inherited from the diploid perennial H. maximiliani Schrader. Overall, these studies suggested that there is an excellent opportunity for breeders to develop sunflower germplasm with improved TE. This can be achieved, in part, by avoiding cytoplasms such as the MAX1 cytoplasm.

The effects of habitat fragmentation due to forestry plantation establishment on the demography and genetic variation of a marsupial carnivore, Antechinus agilis

We conducted a demographic and genetic study to investigate the effects of fragmentation due to the establishment of an exotic softwood plantation on populations of a small marsupial carnivore, the agile antechinus (Antechinus agilis), and the factors influencing the persistence of those populations in the fragmented habitat. The first aspect of the study was a descriptive analysis of patch occupancy and population size, in which we found a patch occupancy rate of 70% among 23 sites in the fragmented habitat compared to 100% among 48 sites with the same habitat characteristics in unfragmented habitat. Mark-recapture analyses yielded most-likely population size estimates of between 3 and 85 among the 16 occupied patches in the fragmented habitat. Hierarchical partitioning and model selection were used to identify geographic and habitat-related characteristics that influence patch occupancy and population size. Patch occupancy was primarily influenced by geographic isolation and habitat quality (vegetation basal area). The variance in population size among occupied sites was influenced primarily by forest type (dominant Eucalyptus species) and, to a lesser extent, by patch area and topographic context (gully sites had larger populations). A comparison of the sex ratios between the samples from the two habitat contexts revealed a significant deficiency of males in the fragmented habitat. We hypothesise that this is due to male-biased dispersal in an environment with increased dispersal-associated mortality. The population size and sex ratio data were incorporated into a simulation study to estimate the proportion of genetic diversity that would have been lost over the known timescale since fragmentation if the patch populations had been totally isolated. The observed difference in genetic diversity (gene diversity and allelic richness at microsatellite and mitochondrial markers) between 16 fragmented and 12 unfragmented sites was extremely low and inconsistent with the isolation of the patch populations. Our results show that although the remnant habitat patches comprise approximately 2% of the study area, they can support non-isolated populations. However, the distribution of agile antechinus populations in the fragmented system is dependent on habitat quality and patch connectivity. (C) 2004 Elsevier Ltd. All rights reserved.

A framework for the study of genetic variation in migratory behaviour

Evolutionary change results from selection acting on genetic variation. For migration to be successful, many different aspects of an animal's physiology and behaviour need to function in a co-coordinated way. Changes in one migratory trait are therefore likely to be accompanied by changes in other migratory and life-history traits. At present, we have some knowledge of the pressures that operate at the various stages of migration, but we know very little about the extent of genetic variation in various aspects of the migratory syndrome. As a consequence, our ability to predict which species is capable of what kind of evolutionary change, and at which rate, is limited. Here, we review how our evolutionary understanding of migration may benefit from taking a quantitative-genetic approach and present a framework for studying the causes of phenotypic variation. We review past research, that has mainly studied single migratory traits in captive birds, and discuss how this work could be extended to study genetic variation in the wild and to account for genetic correlations and correlated selection. In the future, reaction-norm approaches may become very important, as they allow the study of genetic and environmental effects on phenotypic expression within a single framework, as well as of their interactions. We advocate making more use of repeated measurements on single individuals to study the causes of among-individual variation in the wild, as they are easier to obtain than data on relatives and can provide valuable information for identifying and selecting traits. This approach will be particularly informative if it involves systematic testing of individuals under different environmental conditions. We propose extending this research agenda by using optimality models to predict levels of variation and covariation among traits and constraints. This may help us to select traits in which we might expect genetic variation, and to identify the most informative environmental axes. We also recommend an expansion of the passerine model, as this model does not apply to birds, like geese, where cultural transmission of spatio-temporal information is an important determinant of migration patterns and their variation.

Interpreting Human Genetic Variation through Genomics and In Vivo Models

Improvements in genomic technology, both in the increased speed and reduced cost of sequencing, have expanded the appreciation of the abundance of human genetic variation. However the sheer amount of variation, as well as the varying type and genomic content of variation, poses a challenge in understanding the clinical consequence of a single mutation. This work uses several methodologies to interpret the observed variation in the human genome, and presents novel strategies for the prediction of allele pathogenicity.

Using the zebrafish model system as an in vivo assay of allele function, we identified a novel driver of Bardet-Biedl Syndrome (BBS) in CEP76. A combination of targeted sequencing of 785 cilia-associated genes in a cohort of BBS patients and subsequent in vivo functional assays recapitulating the human phenotype gave strong evidence for the role of CEP76 mutations in the pathology of an affected family. This portion of the work demonstrated the necessity of functional testing in validating disease-associated mutations, and added to the catalogue of known BBS disease genes.

Further study into the role of copy-number variations (CNVs) in a cohort of BBS patients showed the significant contribution of CNVs to disease pathology. Using high-density array comparative genomic hybridization (aCGH) we were able to identify pathogenic CNVs as small as several hundred bp. Dissection of constituent gene and in vivo experiments investigating epistatic interactions between affected genes allowed for an appreciation of several paradigms by which CNVs can contribute to disease. This study revealed that the contribution of CNVs to disease in BBS patients is much higher than previously expected, and demonstrated the necessity of consideration of CNV contribution in future (and retrospective) investigations of human genetic disease.

Finally, we used a combination of comparative genomics and in vivo complementation assays to identify second-site compensatory modification of pathogenic alleles. These pathogenic alleles, which are found compensated in other species (termed compensated pathogenic deviations [CPDs]), represent a significant fraction (from 3 – 10%) of human disease-associated alleles. In silico pathogenicity prediction algorithms, a valuable method of allele prioritization, often misrepresent these alleles as benign, leading to omission of possibly informative variants in studies of human genetic disease. We created a mathematical model that was able to predict CPDs and putative compensatory sites, and functionally showed in vivo that second-site mutation can mitigate the pathogenicity of disease alleles. Additionally, we made publically available an in silico module for the prediction of CPDs and modifier sites.

These studies have advanced the ability to interpret the pathogenicity of multiple types of human variation, as well as made available tools for others to do so as well.

Dissecting the Functional Impacts of Non-Coding Genetic Variation

A large proportion of the variation in traits between individuals can be attributed to variation in the nucleotide sequence of the genome. The most commonly studied traits in human genetics are related to disease and disease susceptibility. Although scientists have identified genetic causes for over 4,000 monogenic diseases, the underlying mechanisms of many highly prevalent multifactorial inheritance disorders such as diabetes, obesity, and cardiovascular disease remain largely unknown. Identifying genetic mechanisms for complex traits has been challenging because most of the variants are located outside of protein-coding regions, and determining the effects of such non-coding variants remains difficult. In this dissertation, I evaluate the hypothesis that such non-coding variants contribute to human traits and diseases by altering the regulation of genes rather than the sequence of those genes. I will specifically focus on studies to determine the functional impacts of genetic variation associated with two related complex traits: gestational hyperglycemia and fetal adiposity. At the genomic locus associated with maternal hyperglycemia, we found that genetic variation in regulatory elements altered the expression of the HKDC1 gene. Furthermore, we demonstrated that HKDC1 phosphorylates glucose in vitro and in vivo, thus demonstrating that HKDC1 is a fifth human hexokinase gene. At the fetal-adiposity associated locus, we identified variants that likely alter VEPH1 expression in preadipocytes during differentiation. To make such studies of regulatory variation high-throughput and routine, we developed POP-STARR, a novel high throughput reporter assay that can empirically measure the effects of regulatory variants directly from patient DNA. By combining targeted genome capture technologies with STARR-seq, we assayed thousands of haplotypes from 760 individuals in a single experiment. We subsequently used POP-STARR to identify three key features of regulatory variants: that regulatory variants typically have weak effects on gene expression; that the effects of regulatory variants are often coordinated with respect to disease-risk, suggesting a general mechanism by which the weak effects can together have phenotypic impact; and that nucleotide transversions have larger impacts on enhancer activity than transitions. Together, the findings presented here demonstrate successful strategies for determining the regulatory mechanisms underlying genetic associations with human traits and diseases, and value of doing so for driving novel biological discovery.