Population Analysis of Large Copy Number Variants and Hotspots of Human Genetic Disease

Copy number variants (CNVs) contribute to human genetic and phenotypic diversity. However, the distribution of larger CNVs in the general population remains largely unexplored. We identify large variants in ~2500 individuals by using Illumina SNP data, with an emphasis on “hotspots” prone to recurrent mutations. We find variants larger than 500 kb in 5%–10% of individuals and variants greater than 1 Mb in 1%–2%. In contrast to previous studies, we find limited evidence for stratification of CNVs in geographically distinct human populations. Importantly, our sample size permits a robust distinction between truly rare and polymorphic but low-frequency copy number variation. We find that a significant fraction of individual CNVs larger than 100 kb are rare and that both gene density and size are strongly anticorrelated with allele frequency. Thus, although large CNVs commonly exist in normal individuals, which suggests that size alone can not be used as a predictor of pathogenicity, such variation is generally deleterious. Considering these observations, we combine our data with published CNVs from more than 12,000 individuals contrasting control and neurological disease collections. This analysis identifies known disease loci and highlights additional CNVs (e.g., 3q29, 16p12, and 15q25.2) for further investigation. This study provides one of the first analyses of large, rare (0.1%–1%) CNVs in the general population, with insights relevant to future analyses of genetic disease.     

梨小食心虫微卫星标记的扩增稳定性及多态性分析

【目的】筛选适合我国梨小食心虫Grapholita molesta种群遗传学研究的微卫星位点,并依据所筛选的微卫星位点进行梨小食心虫地理种群的遗传多样性分析。【方法】利用欧洲梨小食心虫和苹果蠹蛾Cydia pomonella种群中已报道的11个微卫星位点, 分析各位点在我国12个种群257头梨小食心虫样本中的扩增稳定性,再进行其多态性分析,筛选适合的位点,然后进行种群遗传多态性分析。【结果】在分析的11个微卫星位点中, 位点Gm01, Gm03, Gm04和Cyd15无法稳定扩增; 位点Gm05扩增成功率较低, 位点Gm07遗传多态性较低; 而位点Gm02, Gm06, Gm08, Gm09和Gm10等扩增效果稳定且遗传多态性丰富。这5个稳定扩增的微卫星位点平均等位基因数量（N_A）为7.417~12.500, 平均观察杂合度（H_o）为0.366~0.655, 平均期望杂合度（H_e）为0.642~0.846, 多态信息含量（PIC）为0.800~0.935。【结论】本研究成功筛选出位点Gm02, Gm06, Gm08, Gm09和Gm10等5个微卫星位点。基于这5个微卫星位点标记的结果显示, 山东和陕西不同梨小食心虫地理种群均具有丰富的遗传多样性。 这5个位点可以适用于我国梨小食心虫种群的进一步遗传分析研究。     

Predicting signatures of "synthetic associations" and "natural associations" from empirical patterns of human genetic variation

Genome-wide association studies (GWAS) have in recent years discovered thousands of associated markers for hundreds of phenotypes. However, associated loci often only explain a relatively small fraction of heritability and the link between association and causality has yet to be uncovered for most loci. Rare causal variants have been suggested as one scenario that may partially explain these shortcomings. Specifically, Dickson et al. recently reported simulations of rare causal variants that lead to association signals of common, tag single nucleotide polymorphisms, dubbed "synthetic associations". However, an open question is what practical implications synthetic associations have for GWAS. Here, we explore the signatures exhibited by such "synthetic associations" and their implications based on patterns of genetic variation observed in human populations, thus accounting for human evolutionary history -a force disregarded in previous simulation studies. This is made possible by human population genetic data from HapMap 3 consisting of both resequencing and array-based genotyping data for the same set of individuals from multiple populations. We report that synthetic associations tend to be further away from the underlying risk alleles compared to "natural associations" (i.e. associations due to underlying common causal variants), but to a much lesser extent than previously predicted, with both the age and the effect size of the risk allele playing a part in this phenomenon. We find that while a synthetic association has a lower probability of capturing causal variants within its linkage disequilibrium block, sequencing around the associated variant need not extend substantially to have a high probability of capturing at least one causal variant. We also show that the minor allele frequency of synthetic associations is lower than of natural associations for most, but not all, loci that we explored. Finally, we find the variance in associated allele frequency to be a potential indicator of synthetic associations.     

Insights into recently fragmented vole populations from combined genetic and demographic data

We combined demographic and genetic data to evaluate the effects of habitat fragmentation on the population structure of the California red-backed vole (Clethrionomys californicus). We analysed variation in the mitochondrial DNA (mtDNA) control region and five nuclear microsatellite loci in small samples collected from two forest fragments and an unfragmented control site in 1990-91. We intensively sampled the same forest fragments and two different control sites in 1998 and 1999. Vole abundances fluctuated greatly at sizes below 50 individuals per fragment. Fragment populations had significantly lower mtDNA allelic diversity than controls, but not nuclear heterozygosity or numbers of alleles. The use of only trapping and/or mtDNA marker data would imply that fragment populations are at least partially isolated and vulnerable to inbreeding depression. In contrast, the abundance estimates combined with microsatellite data show that small fragment populations must be linked to nearby forests by high rates of migration. These results provide evidence for the usefulness of combining genetic and demographic data to understand nonequilibrium population structure in recently fragmented habitats.     

The Power of Gene-Based Rare Variant Methods to Detect Disease-Associated Variation and Test Hypotheses About Complex Disease

Genome and exome sequencing in large cohorts enables characterization of the role of rare variation in complex diseases. Success in this endeavor, however, requires investigators to test a diverse array of genetic hypotheses which differ in the number, frequency and effect sizes of underlying causal variants. In this study, we evaluated the power of gene-based association methods to interrogate such hypotheses, and examined the implications for study design. We developed a flexible simulation approach, using 1000 Genomes data, to (a) generate sequence variation at human genes in up to 10K case-control samples, and (b) quantify the statistical power of a panel of widely used gene-based association tests under a variety of allelic architectures, locus effect sizes, and significance thresholds. For loci explaining ~1% of phenotypic variance underlying a common dichotomous trait, we find that all methods have low absolute power to achieve exome-wide significance (~5-20% power at α=2.5×10^-6) in 3K individuals; even in 10K samples, power is modest (~60%). The combined application of multiple methods increases sensitivity, but does so at the expense of a higher false positive rate. MiST, SKAT-O, and KBAC have the highest individual mean power across simulated datasets, but we observe wide architecture-dependent variability in the individual loci detected by each test, suggesting that inferences about disease architecture from analysis of sequencing studies can differ depending on which methods are used. Our results imply that tens of thousands of individuals, extensive functional annotation, or highly targeted hypothesis testing will be required to confidently detect or exclude rare variant signals at complex disease loci.     

Microsatellite variation among divergent populations of stalk-eyed flies, genus Cyrtodiopsis

Microsatellite primers are often developed in one species and used to assess neutral variability in related species. Such analyses may be confounded by ascertainment bias (i.e. a decline in amplification success and allelic variability with increasing genetic distance from the source of the microsatellites). In addition, other factors, such as the size of the microsatellite, whether it consists of perfect or interrupted tandem repeats, and whether it is autosomal or X-linked, can affect variation. To test the relative importance of these factors on microsatellite variation, we examine patterns of amplification and allelic diversity in 52 microsatellite loci amplified from five individuals in each of six populations of Cyrtodiopsis stalk-eyed flies that range from 2.2 % to 11.2% mitochondrial DNA sequence divergence from the population used for microsatellite development. We find that amplification success and most measures of allelic diversity declined with genetic distance from the source population, in some cases an order of magnitude faster than in birds or mammals. The median and range of the repeat array length did not decline with genetic distance. In addition, for loci on the X chromosome, we find evidence of lower observed heterozygosity compared with loci on autosomes. The differences in variability between X-linked and autosomal loci are not adequately explained by differences in effective population sizes of the chromosomes. We suggest, instead, that periodic selection events associated with X-chromosome meiotic drive, which is present in many of these populations, reduces X-linked variation.     

A Genome-Wide Survey of Genetic Variation in Gorillas Using Reduced Representation Sequencing

All non-human great apes are endangered in the wild, and it is therefore important to gain an understanding of their demography and genetic diversity. Whole genome assembly projects have provided an invaluable foundation for understanding genetics in all four genera, but to date genetic studies of multiple individuals within great ape species have largely been confined to mitochondrial DNA and a small number of other loci. Here, we present a genome-wide survey of genetic variation in gorillas using a reduced representation sequencing approach, focusing on the two lowland subspecies. We identify 3,006,670 polymorphic sites in 14 individuals: 12 western lowland gorillas (Gorilla gorilla gorilla) and 2 eastern lowland gorillas (Gorilla beringei graueri). We find that the two species are genetically distinct, based on levels of heterozygosity and patterns of allele sharing. Focusing on the western lowland population, we observe evidence for population substructure, and a deficit of rare genetic variants suggesting a recent episode of population contraction. In western lowland gorillas, there is an elevation of variation towards telomeres and centromeres on the chromosomal scale. On a finer scale, we find substantial variation in genetic diversity, including a marked reduction close to the major histocompatibility locus, perhaps indicative of recent strong selection there. These findings suggest that despite their maintaining an overall level of genetic diversity equal to or greater than that of humans, population decline, perhaps associated with disease, has been a significant factor in recent and long-term pressures on wild gorilla populations.     

Genome-wide Characterization of Shared and Distinct Genetic Components that Influence Blood Lipid Levels in Ethnically Diverse Human Populations

Blood lipid concentrations are heritable risk factors associated with atherosclerosis and cardiovascular diseases. Lipid traits exhibit considerable variation among populations of distinct ancestral origin as well as between individuals within a population. We performed association analyses to identify genetic loci influencing lipid concentrations in African American and Hispanic American women in the Women’s Health Initiative SNP Health Association Resource. We validated one African-specific high-density lipoprotein cholesterol locus at CD36 as well as 14 known lipid loci that have been previously implicated in studies of European populations. Moreover, we demonstrate striking similarities in genetic architecture (loci influencing the trait, direction and magnitude of genetic effects, and proportions of phenotypic variation explained) of lipid traits across populations. In particular, we found that a disproportionate fraction of lipid variation in African Americans and Hispanic Americans can be attributed to genomic loci exhibiting statistical evidence of association in Europeans, even though the precise genes and variants remain unknown. At the same time, we found substantial allelic heterogeneity within shared loci, characterized both by population-specific rare variants and variants shared among multiple populations that occur at disparate frequencies. The allelic heterogeneity emphasizes the importance of including diverse populations in future genetic association studies of complex traits such as lipids; furthermore, the overlap in lipid loci across populations of diverse ancestral origin argues that additional knowledge can be gleaned from multiple populations.     

Drift,selection and adaptive variation in small populations of a threatened rattlesnake

An important goal of conservation genetics is to determine if the viability of small populations is reduced by a loss of adaptive variation due to genetic drift. Here, we assessed the impact of drift and selection on direct measures of adaptive variation (toxin loci encoding venom proteins) in the eastern massasauga rattlesnake (Sistrurus catenatus), a threatened reptile that exists in small isolated populations. We estimated levels of individual polymorphism in 46 toxin loci and 1,467 control loci across 12 populations of this species, and compared the results with patterns of selection on the same loci following speciation of S. catenatus and its closest relative, the western massasauga (S. tergeminus). Multiple lines of evidence suggest that both drift and selection have had observable impacts on standing adaptive variation. In support of drift effects, we found little evidence for selection on toxin variation within populations and a significant positive relationship between current levels of adaptive variation and long‐ and short‐term estimates of effective population size. However, we also observed levels of directional selection on toxin loci among populations that are broadly similar to patterns predicted from interspecific selection analyses that pre‐date the effects of recent drift, and that functional variation in these loci persists despite small short‐term effective sizes. This suggests that much of the adaptive variation present in populations may represent an example of “drift debt,” a nonequilibrium state where present‐day levels of variation overestimate the amount of functional genetic diversity present in future populations.     

Aberrant Gene Expression in Humans

Gene expression as an intermediate molecular phenotype has been a focus of research interest. In particular, studies of expression quantitative trait loci (eQTL) have offered promise for understanding gene regulation through the discovery of genetic variants that explain variation in gene expression levels. Existing eQTL methods are designed for assessing the effects of common variants, but not rare variants. Here, we address the problem by establishing a novel analytical framework for evaluating the effects of rare or private variants on gene expression. Our method starts from the identification of outlier individuals that show markedly different gene expression from the majority of a population, and then reveals the contributions of private SNPs to the aberrant gene expression in these outliers. Using population-scale mRNA sequencing data, we identify outlier individuals using a multivariate approach. We find that outlier individuals are more readily detected with respect to gene sets that include genes involved in cellular regulation and signal transduction, and less likely to be detected with respect to the gene sets with genes involved in metabolic pathways and other fundamental molecular functions. Analysis of polymorphic data suggests that private SNPs of outlier individuals are enriched in the enhancer and promoter regions of corresponding aberrantly-expressed genes, suggesting a specific regulatory role of private SNPs, while the commonly-occurring regulatory genetic variants (i.e., eQTL SNPs) show little evidence of involvement. Additional data suggest that non-genetic factors may also underlie aberrant gene expression. Taken together, our findings advance a novel viewpoint relevant to situations wherein common eQTLs fail to predict gene expression when heritable, rare inter-individual variation exists. The analytical framework we describe, taking into consideration the reality of differential phenotypic robustness, may be valuable for investigating complex traits and conditions.     

Population structure and conservation genetic assessment of the endangered Pugnose Shiner,Notropis anogenus

The Pugnose Shiner is a small minnow with a fragmented distribution across the Great Lakes and Upper Mississippi River in North America. The species is listed federally as endangered in Canada, and in the United States its status varies by state, from Special Concern to Endangered (as well as Extirpated). We conducted a thorough genetic assessment of the Pugnose Shiner using both microsatellite loci and mitochondrial DNA collected for samples across the species range. Our results indicate high levels of population differentiation suggesting restricted dispersal, in some cases at very small geographical scales. We also found strong evidence of small effective population sizes and one case of a genetic bottleneck. Although significant range-wide genetic variation was observed in both microsatellite loci and mitochondrial DNA, the species is best characterized as a single evolutionarily significant unit for conservation purposes.     

Molecular evolution of the phytochrome gene family in sorghum: changing rates of synonymous and replacement evolution

The photoreceptor phytochromes, encoded by a small gene family, are responsible for controlling the expression of a number of light-responsive genes and photomorphogenic events, including agronomically important phenotypes such as flowering time and shade-avoidance behavior. The understanding and control of flowering time are particularly important goals in sorghum cultivar development for diverse environments, and naturally occurring variation in the phytochrome genes might prove useful in breeding programs. Also of interest is whether variation observed at the phytochrome loci in domesticated sorghum, or in particular races, is a result of human selection. Population genetic studies can reveal evidence of such selection in patterns of polymorphism and divergence. In this study we report a population genetic analysis of the PHY gene family in Sorghum bicolor (L.) Moench in a diverse panel including both cultivated and wild accessions. We show that the level of nucleotide variation in all gene family members is about half the average for this species, consistent with purifying selection acting on these loci. However, the rate of amino acid substitution is accelerated at PHYC compared to the other two loci. In comparisons to a closely related sorghum species, PHYC shows a pattern of intermediate frequency amino acid changes that differ from the patterns observed in comparisons across longer evolutionary distances. There is also a departure from expected patterns of polymorphism and divergence at synonymous sites in PHYC, although the data do not fit a simple model of directional or diversifying selection. Cultivated sorghum has a level of variation similar to that of wild relatives (ssp. verticilliflorum), but many polymorphisms are subspecies-specific, including several amino acid variants.     

The effects of sample size on population genetic diversity estimates in song sparrows Melospiza melodia

To empirically determine the effects of sample size on commonly used measures of average genetic diversity, we genotyped 200 song sparrows Melospiza melodia from two populations, one genetically depauperate (n=100) and the other genetically diverse (n=100), using eight microsatellite loci. These genotypes were used to randomly create 10,000 datasets of differing sizes (5 to 50) for each population to determine what the effects of sample size might be on several estimates of genetic diversity (number of alleles per locus, average observed heterozygosity, and unbiased average expected heterozygosity) in natural populations of conservation concern. We found that at small sample sizes of 5 to 10 individuals, estimates of unbiased heterozygosity outperformed those based on observed heterozygosity or allelic diversity for both low- and high-diversity populations. We also found that when comparing across populations in which different numbers of individuals were sampled, rarefaction provided a useful way to compare estimates of allelic diversity. We recommend that standard errors should be reported for all diversity estimators, especially when sample sizes are small. We also recommend that at least 20 to 30 individuals be sampled in microsatellite studies that assess genetic diversity when working in a population that has an unknown level of diversity. However, research on critically endangered populations (where large sample sizes are impossible or extremely difficult to obtain) should include measures of genetic diversity even if sample sizes are less than ideal. These estimates can be useful in assessing the genetic diversity of the population.     

Allozyme variation and differentiation in African and Indian rhinoceroses

We studied variation at 25 to 31 allozymic loci in African and Asian rhinoceroses. Four taxa in three genera were examined: African Ceratotherium simum simum (northern white rhinoceros), C. s. cottoni (southern white rhinoceros), Diceros bicornis (black rhinoceros), and Rhinoceros unicornis (Indian rhinoceros). Extremely small amounts of intraspecific variation were observed in sample sizes of 2 to 10 presumably unrelated individuals per taxon: P = .00-.10, H = 0.00-0.02. We examined demographic bottlenecks and sampling errors as possible reasons for the low levels of detectable variation. The very small intraspecific genetic distance (D = 0.005) between the two living white rhinoceros subspecies is far less than the distance that has been reported for other mammal subspecies. The mean D value of 0.32 +/- 0.11 between the two African genera was also less than expected given the divergence time of greater than 7 million years suggested by the fossil record. Rhinoceroses may be evolving more slowly at the structural gene loci than are some other mammal groups. The estimate of D = 1.05 +/- 0.24 for the African-Indian split supports this idea, as the lineage diverged at least 26 million years ago. Our results contribute to the currently available scientific information on which management decisions aimed toward saving endangered rhinoceroses should be based.     

Demographic Genetics of Brown Trout (Salmo Trutta) and Estimation of Effective Population Size from Temporal Change of Allele Frequencies

We studied temporal allele frequency shifts over 15 years and estimated the genetically effective size of four natural populations of brown trout (Salmo trutta L.) on the basis of the variation at 14 polymorphic allozyme loci. The allele frequency differences between consecutive cohorts were significant in all four populations. There were no indications of natural selection, and we conclude that random genetic drift is the most likely cause of temporal allele frequency shifts at the loci examined. Effective population sizes were estimated from observed allele frequency shifts among cohorts, taking into consideration the demographic characteristics of each population. The estimated effective sizes of the four populations range from 52 to 480 individuals, and we conclude that the effective size of natural brown trout populations may differ considerably among lakes that are similar in size and other apparent characteristics. In spite of their different effective sizes all four populations have similar levels of genetic variation (average heterozygosity) indicating that excessive loss of genetic variability has been retarded, most likely because of gene flow among neighboring populations.     

Partitioning of genetic variation across the genome using multimarker methods in a wild bird population

The underlying basis of genetic variation in quantitative traits, in terms of the number of causal variants and the size of their effects, is largely unknown in natural populations. The expectation is that complex quantitative trait variation is attributable to many, possibly interacting, causal variants, whose effects may depend upon the sex, age and the environment in which they are expressed. A recently developed methodology in animal breeding derives a value of relatedness among individuals from high‐density genomic marker data, to estimate additive genetic variance within livestock populations. Here, we adapt and test the effectiveness of these methods to partition genetic variation for complex traits across genomic regions within ecological study populations where individuals have varying degrees of relatedness. We then apply this approach for the first time to a natural population and demonstrate that genetic variation in wing length in the great tit (Parus major) reflects contributions from multiple genomic regions. We show that a polygenic additive mode of gene action best describes the patterns observed, and we find no evidence of dosage compensation for the sex chromosome. Our results suggest that most of the genomic regions that influence wing length have the same effects in both sexes. We found a limited amount of genetic variance in males that is attributed to regions that have no effects in females, which could facilitate the sexual dimorphism observed for this trait. Although this exploratory work focuses on one complex trait, the methodology is generally applicable to any trait for any laboratory or wild population, paving the way for investigating sex‐, age‐ and environment‐specific genetic effects and thus the underlying genetic architecture of phenotype in biological study systems.     

Inferring population history and demography using microsatellites, mitochondrial DNA, and major histocompatibility complex (MHC) genes

Microsatellites and mitochondrial DNA (mtDNA) have traditionally been used in population genetics because of their variability and presumed neutrality, whereas genes of the major histocompatibility complex (MHC) are increasingly of interest because strong selective pressures shape their standing variation. Despite the potential for MHC genes, microsatellites, and mtDNA sequences to complement one another in deciphering population history and demography, the three are rarely used in tandem. Here we report on MHC, microsatellite, and mtDNA variability in a single large population of the eastern tiger salamander (Ambystoma tigrinum tigrinum). We use the mtDNA mismatch distribution and, on microsatellite data, the imbalance index and bottleneck tests to infer aspects of population history and demography. Haplotype and allelic variation was high at all loci surveyed, and heterozygosity was high at the nuclear loci. We find concordance among neutral molecular markers that suggests our study population originated from post-Pleistocene expansions of multiple, fragmented sources that shared few migrants. Differences in N(e) estimates derived from haploid and diploid genetic markers are potentially attributable to secondary contact among source populations that experienced rapid mtDNA divergence and comparatively low levels of nuclear DNA divergence. We find strong evidence of natural selection acting on MHC genes and estimate long-term effective population sizes (N(e)) that are very large, making small selection intensities significant evolutionary forces in this population.     

Genetic Analysis of Variation in Human Meiotic Recombination

The number of recombination events per meiosis varies extensively among individuals. This recombination phenotype differs between female and male, and also among individuals of each gender. In this study, we used high-density SNP genotypes of over 2,300 individuals and their offspring in two datasets to characterize recombination landscape and to map the genetic variants that contribute to variation in recombination phenotypes. We found six genetic loci that are associated with recombination phenotypes. Two of these (RNF212 and an inversion on chromosome 17q21.31) were previously reported in the Icelandic population, and this is the first replication in any other population. Of the four newly identified loci (KIAA1462, PDZK1, UGCG, NUB1), results from expression studies provide support for their roles in meiosis. Each of the variants that we identified explains only a small fraction of the individual variation in recombination. Notably, we found different sequence variants associated with female and male recombination phenotypes, suggesting that they are regulated by different genes. Characterization of genetic variants that influence natural variation in meiotic recombination will lead to a better understanding of normal meiotic events as well as of non-disjunction, the primary cause of pregnancy loss.     

Analysis of quantitative trait loci that influence animal behavior

Behavioral differences between inbred strains of mice and rats have a genetic basis that can now be dissected using quantitative trait locus (QTL) analysis. Over the last 10 years, a large number of genetic loci that influence behavior have been mapped. In this article I review what that information has revealed about the genetic architecture of behavior. I show that most behaviors are influenced by QTL of small effect, each contributing to less than 10% of the variance of a behavioral trait. The small effect of each QTL on behavioral variation suggests that the mutational spectrum is different from that which results in Mendelian disorders. Regions of DNA should be appropriately prioritized to find the molecular variants, for instance by looking at sequences that control the level of gene expression rather than variants in coding regions. While the number of allelic loci that can contribute to a trait is large, this is not necessarily the case: the analysis of selected strains shows that a remarkably small number of QTL can explain the bulk of the genetic variation in behavior. I conclude by arguing that genetic mapping has more to offer than a starting point for positional cloning projects. With advances in multivariate analyses, mapping can also test hypotheses about the psychological processes that give rise to behavioral variation.