Genome-wide DNA polymorphism analyses using VariScan

BACKGROUND: DNA sequence polymorphisms analysis can provide valuable information on the evolutionary forces shaping nucleotide variation, and provides an insight into the functional significance of genomic regions. The recent ongoing genome projects will radically improve our capabilities to detect specific genomic regions shaped by natural selection. Current available methods and software, however, are unsatisfactory for such genome-wide analysis. RESULTS: We have developed methods for the analysis of DNA sequence polymorphisms at the genome-wide scale. These methods, which have been tested on a coalescent-simulated and actual data files from mouse and human, have been implemented in the VariScan software package version 2.0. Additionally, we have also incorporated a graphical-user interface. The main features of this software are: i) exhaustive population-genetic analyses including those based on the coalescent theory; ii) analysis adapted to the shallow data generated by the high-throughput genome projects; iii) use of genome annotations to conduct a comprehensive analyses separately for different functional regions; iv) identification of relevant genomic regions by the sliding-window and wavelet-multiresolution approaches; v) visualization of the results integrated with current genome annotations in commonly available genome browsers. CONCLUSION: VariScan is a powerful and flexible suite of software for the analysis of DNA polymorphisms. The current version implements new algorithms, methods, and capabilities, providing an important tool for an exhaustive exploratory analysis of genome-wide DNA polymorphism data.

Linkage disequilibrium between a melanin-based colour polymorphism and tail length in the barn owl

A linkage disequilibrium between sexually selected and life history traits can be explained by three mutually non-exclusive mechanisms. Genes coding for two traits may be located close on the same chromosome, genes responsible for variation in one of the trait may pleiotropically alter the other, and non-random pairing with respect to two traits may generate a non-physical linkage disequilibrium between their genes. Knowledge of which of these three mechanisms is responsible for a covariation between two traits is of interest to understand why differently ornamented individuals differ in several phenotypic aspects. In Switzerland, barn owls Tyto alba mate randomly with respect to a colour polymorphism generating a large range of variants between reddish-brown and white, males being lighter coloured than females. Several studies have shown that plumage coloration is not neutral with respect to some life history components. To test whether coloration is genetically associated with body size, partial cross-fostering experiments were performed by exchanging some hatchlings between nests. These experiments showed that darker biological fathers produce longer-tailed offspring. This sex-specific pattern is consistent with the hypothesis of non-physical linkage disequilibrium. In line with this hypothesis, darker coloured males were mated with longer-tailed females, whereas female coloration was not associated with tail length of their mate. The finding that dark nestlings had a longer tail than their pale siblings also supports the physical linkage and pleiotropy hypotheses. Therefore, non-random pairing can generate or strengthen a genetic covariation between a secondary sexual character and a morphological trait.

Inference of chromosomal inversion dynamics from Pool-Seq data in natural and laboratory populations of Drosophila melanogaster.

Sequencing of pools of individuals (Pool-Seq) represents a reliable and cost-effective approach for estimating genome-wide SNP and transposable element insertion frequencies. However, Pool-Seq does not provide direct information on haplotypes so that, for example, obtaining inversion frequencies has not been possible until now. Here, we have developed a new set of diagnostic marker SNPs for seven cosmopolitan inversions in Drosophila melanogaster that can be used to infer inversion frequencies from Pool-Seq data. We applied our novel marker set to Pool-Seq data from an experimental evolution study and from North American and Australian latitudinal clines. In the experimental evolution data, we find evidence that positive selection has driven the frequencies of In(3R)C and In(3R)Mo to increase over time. In the clinal data, we confirm the existence of frequency clines for In(2L)t, In(3L)P and In(3R)Payne in both North America and Australia and detect a previously unknown latitudinal cline for In(3R)Mo in North America. The inversion markers developed here provide a versatile and robust tool for characterizing inversion frequencies and their dynamics in Pool-Seq data from diverse D. melanogaster populations.

An 18-kDa translocator protein (TSPO) polymorphism explains differences in binding affinity of the PET radioligand PBR28

[(11)C]PBR28 binds the 18-kDa Translocator Protein (TSPO) and is used in positron emission tomography (PET) to detect microglial activation. However, quantitative interpretations of signal are confounded by large interindividual variability in binding affinity, which displays a trimodal distribution compatible with a codominant genetic trait. Here, we tested directly for an underlying genetic mechanism to explain this. Binding affinity of PBR28 was measured in platelets isolated from 41 human subjects and tested for association with polymorphisms in TSPO and genes encoding other proteins in the TSPO complex. Complete agreement was observed between the TSPO Ala147Thr genotype and PBR28 binding affinity phenotype (P value=3.1 x 10(-13)). The TSPO Ala147Thr polymorphism predicts PBR28 binding affinity in human platelets. As all second-generation TSPO PET radioligands tested hitherto display a trimodal distribution in binding affinity analogous to PBR28, testing for this polymorphism may allow quantitative interpretation of TSPO PET studies with these radioligands.

Association of genetic variation with systolic and diastolic blood pressure among African Americans: the Candidate Gene Association Resource study.

The prevalence of hypertension in African Americans (AAs) is higher than in other US groups; yet, few have performed genome-wide association studies (GWASs) in AA. Among people of European descent, GWASs have identified genetic variants at 13 loci that are associated with blood pressure. It is unknown if these variants confer susceptibility in people of African ancestry. Here, we examined genome-wide and candidate gene associations with systolic blood pressure (SBP) and diastolic blood pressure (DBP) using the Candidate Gene Association Resource (CARe) consortium consisting of 8591 AAs. Genotypes included genome-wide single-nucleotide polymorphism (SNP) data utilizing the Affymetrix 6.0 array with imputation to 2.5 million HapMap SNPs and candidate gene SNP data utilizing a 50K cardiovascular gene-centric array (ITMAT-Broad-CARe [IBC] array). For Affymetrix data, the strongest signal for DBP was rs10474346 (P= 3.6 × 10(-8)) located near GPR98 and ARRDC3. For SBP, the strongest signal was rs2258119 in C21orf91 (P= 4.7 × 10(-8)). The top IBC association for SBP was rs2012318 (P= 6.4 × 10(-6)) near SLC25A42 and for DBP was rs2523586 (P= 1.3 × 10(-6)) near HLA-B. None of the top variants replicated in additional AA (n = 11 882) or European-American (n = 69 899) cohorts. We replicated previously reported European-American blood pressure SNPs in our AA samples (SH2B3, P= 0.009; TBX3-TBX5, P= 0.03; and CSK-ULK3, P= 0.0004). These genetic loci represent the best evidence of genetic influences on SBP and DBP in AAs to date. More broadly, this work supports that notion that blood pressure among AAs is a trait with genetic underpinnings but also with significant complexity.

Two high throughput technologies to detect segmental aneuploidies identify new Williams-Beuren syndrome patients with atypical deletions.

OBJECTIVE: To develop and compare two new technologies for diagnosing a contiguous gene syndrome, the Williams-Beuren syndrome (WBS). METHODS: The first proposed method, named paralogous sequence quantification (PSQ), is based on the use of paralogous sequences located on different chromosomes and quantification of specific mismatches present at these loci using pyrosequencing technology. The second exploits quantitative real time polymerase chain reaction (QPCR) to assess the relative quantity of an analysed locus. RESULTS: A correct and unambiguous diagnosis was obtained for 100% of the analysed samples with either technique (n = 165 and n = 155, respectively). These methods allowed the identification of two patients with atypical deletions in a cohort of 182 WBS patients. Both patients presented with mild facial anomalies, mild mental retardation with impaired visuospatial cognition, supravalvar aortic stenosis, and normal growth indices. These observations are consistent with the involvement of GTF2IRD1 or GTF2I in some of the WBS facial features. CONCLUSIONS: Both PSQ and QPCR are robust, easy to interpret, and simple to set up. They represent a competitive alternative for the diagnosis of segmental aneuploidies in clinical laboratories. They have advantages over fluorescence in situ hybridisation or microsatellites/SNP genotyping for detecting short segmental aneuploidies as the former is costly and labour intensive while the latter depends on the informativeness of the polymorphisms.

Interleukin-1 receptor antagonist gene (IL-1RN) polymorphism is a predictive factor of clinical pregnancy after IVF.

BACKGROUND: Only 25% of IVF transfer cycles lead to a clinical pregnancy, calling for continued technical progress but also more in depth analysis of patients' individual characteristics. The interleukin-1 (IL-1) system and matrix metalloproteinases (MMPs) are strongly implicated in embryo implantation. The genes coding for IL-1Ra (gene symbol IL-1RN), IL-1beta, MMP2 and MMP9 bear functional polymorphisms. We analysed the maternal genetic profile at these polymorphic sites in IVF patients, to determine possible correlations with IVF outcome. METHODS: One hundred and sixty women undergoing an IVF cycle were enrolled and a buccal smear was obtained. The presence of IL-1RN variable number of tandem repeats and IL-1B + 3953, MMP2-1306 and MMP9-1562 single nucleotide substitutions were determined. Patients were divided into pregnancy failures (119), biochemical pregnancies (8) and clinical pregnancies (33). RESULTS: There was a 40% decrease in IL-1RN*2 allele frequency (P = 0.024) and a 45% decrease in IL-1RN*2 carrier status in the clinical pregnancy group as compared to the pregnancy failure group (P = 0.017). This decrease was still statistically significant after a multivariate logistic regression analysis. The likelihood of a clinical pregnancy was decreased accordingly in IL-1RN*2 carriers: odds ratio = 0.349, 95% confidence interval = 0.2-0.8, P = 0.017. The IL-1B, MMP2 and MMP9 polymorphisms showed no correlation with IVF outcome. CONCLUSIONS: IL-1RN*2 allele carriage is associated with a poor prognosis of achieving a pregnancy after IVF.

Identification and validation of copy number variants using SNP genotyping arrays from a large clinical cohort.

BACKGROUND: Genotypes obtained with commercial SNP arrays have been extensively used in many large case-control or population-based cohorts for SNP-based genome-wide association studies for a multitude of traits. Yet, these genotypes capture only a small fraction of the variance of the studied traits. Genomic structural variants (GSV) such as Copy Number Variation (CNV) may account for part of the missing heritability, but their comprehensive detection requires either next-generation arrays or sequencing. Sophisticated algorithms that infer CNVs by combining the intensities from SNP-probes for the two alleles can already be used to extract a partial view of such GSV from existing data sets. RESULTS: Here we present several advances to facilitate the latter approach. First, we introduce a novel CNV detection method based on a Gaussian Mixture Model. Second, we propose a new algorithm, PCA merge, for combining copy-number profiles from many individuals into consensus regions. We applied both our new methods as well as existing ones to data from 5612 individuals from the CoLaus study who were genotyped on Affymetrix 500K arrays. We developed a number of procedures in order to evaluate the performance of the different methods. This includes comparison with previously published CNVs as well as using a replication sample of 239 individuals, genotyped with Illumina 550K arrays. We also established a new evaluation procedure that employs the fact that related individuals are expected to share their CNVs more frequently than randomly selected individuals. The ability to detect both rare and common CNVs provides a valuable resource that will facilitate association studies exploring potential phenotypic associations with CNVs. CONCLUSION: Our new methodologies for CNV detection and their evaluation will help in extracting additional information from the large amount of SNP-genotyping data on various cohorts and use this to explore structural variants and their impact on complex traits.

LDL-cholesterol concentrations: a genome-wide association study.

BACKGROUND: LDL cholesterol has a causal role in the development of cardiovascular disease. Improved understanding of the biological mechanisms that underlie the metabolism and regulation of LDL cholesterol might help to identify novel therapeutic targets. We therefore did a genome-wide association study of LDL-cholesterol concentrations. METHODS: We used genome-wide association data from up to 11,685 participants with measures of circulating LDL-cholesterol concentrations across five studies, including data for 293 461 autosomal single nucleotide polymorphisms (SNPs) with a minor allele frequency of 5% or more that passed our quality control criteria. We also used data from a second genome-wide array in up to 4337 participants from three of these five studies, with data for 290,140 SNPs. We did replication studies in two independent populations consisting of up to 4979 participants. Statistical approaches, including meta-analysis and linkage disequilibrium plots, were used to refine association signals; we analysed pooled data from all seven populations to determine the effect of each SNP on variations in circulating LDL-cholesterol concentrations. FINDINGS: In our initial scan, we found two SNPs (rs599839 [p=1.7x10(-15)] and rs4970834 [p=3.0x10(-11)]) that showed genome-wide statistical association with LDL cholesterol at chromosomal locus 1p13.3. The second genome screen found a third statistically associated SNP at the same locus (rs646776 [p=4.3x10(-9)]). Meta-analysis of data from all studies showed an association of SNPs rs599839 (combined p=1.2x10(-33)) and rs646776 (p=4.8x10(-20)) with LDL-cholesterol concentrations. SNPs rs599839 and rs646776 both explained around 1% of the variation in circulating LDL-cholesterol concentrations and were associated with about 15% of an SD change in LDL cholesterol per allele, assuming an SD of 1 mmol/L. INTERPRETATION: We found evidence for a novel locus for LDL cholesterol on chromosome 1p13.3. These results potentially provide insight into the biological mechanisms that underlie the regulation of LDL cholesterol and might help in the discovery of novel therapeutic targets for cardiovascular disease.

Genome-wide association for abdominal subcutaneous and visceral adipose reveals a novel locus for visceral fat in women.

Body fat distribution, particularly centralized obesity, is associated with metabolic risk above and beyond total adiposity. We performed genome-wide association of abdominal adipose depots quantified using computed tomography (CT) to uncover novel loci for body fat distribution among participants of European ancestry. Subcutaneous and visceral fat were quantified in 5,560 women and 4,997 men from 4 population-based studies. Genome-wide genotyping was performed using standard arrays and imputed to ~2.5 million Hapmap SNPs. Each study performed a genome-wide association analysis of subcutaneous adipose tissue (SAT), visceral adipose tissue (VAT), VAT adjusted for body mass index, and VAT/SAT ratio (a metric of the propensity to store fat viscerally as compared to subcutaneously) in the overall sample and in women and men separately. A weighted z-score meta-analysis was conducted. For the VAT/SAT ratio, our most significant p-value was rs11118316 at LYPLAL1 gene (p = 3.1 × 10E-09), previously identified in association with waist-hip ratio. For SAT, the most significant SNP was in the FTO gene (p = 5.9 × 10E-08). Given the known gender differences in body fat distribution, we performed sex-specific analyses. Our most significant finding was for VAT in women, rs1659258 near THNSL2 (p = 1.6 × 10-08), but not men (p = 0.75). Validation of this SNP in the GIANT consortium data demonstrated a similar sex-specific pattern, with observed significance in women (p = 0.006) but not men (p = 0.24) for BMI and waist circumference (p = 0.04 [women], p = 0.49 [men]). Finally, we interrogated our data for the 14 recently published loci for body fat distribution (measured by waist-hip ratio adjusted for BMI); associations were observed at 7 of these loci. In contrast, we observed associations at only 7/32 loci previously identified in association with BMI; the majority of overlap was observed with SAT. Genome-wide association for visceral and subcutaneous fat revealed a SNP for VAT in women. More refined phenotypes for body composition and fat distribution can detect new loci not previously uncovered in large-scale GWAS of anthropometric traits.

Polymorphisms in toll-like receptor 4 and toll-like receptor 9 influence viral load in a seroincident cohort of HIV-1-infected individuals.

OBJECTIVES: Toll-like receptors (TLRs) are innate immune sensors that are integral to resisting chronic and opportunistic infections. Mounting evidence implicates TLR polymorphisms in susceptibilities to various infectious diseases, including HIV-1. We investigated the impact of TLR single nucleotide polymorphisms (SNPs) on clinical outcome in a seroincident cohort of HIV-1-infected volunteers. DESIGN: We analyzed TLR SNPs in 201 antiretroviral treatment-naive HIV-1-infected volunteers from a longitudinal seroincident cohort with regular follow-up intervals (median follow-up 4.2 years, interquartile range 4.4). Participants were stratified into two groups according to either disease progression, defined as peripheral blood CD4(+) T-cell decline over time, or peak and setpoint viral load. METHODS: Haplotype tagging SNPs from TLR2, TLR3, TLR4, and TLR9 were detected by mass array genotyping, and CD4(+) T-cell counts and viral load measurements were determined prior to antiretroviral therapy initiation. The association of TLR haplotypes with viral load and rapid progression was assessed by multivariate regression models using age and sex as covariates. RESULTS: Two TLR4 SNPs in strong linkage disequilibrium [1063 A/G (D299G) and 1363 C/T (T399I)] were more frequent among individuals with high peak viral load compared with low/moderate peak viral load (odds ratio 6.65, 95% confidence interval 2.19-20.46, P < 0.001; adjusted P = 0.002 for 1063 A/G). In addition, a TLR9 SNP previously associated with slow progression was found less frequently among individuals with high viral setpoint compared with low/moderate setpoint (odds ratio 0.29, 95% confidence interval 0.13-0.65, P = 0.003, adjusted P = 0.04). CONCLUSION: This study suggests a potentially new role for TLR4 polymorphisms in HIV-1 peak viral load and confirms a role for TLR9 polymorphisms in disease progression.

Comparative genomic and phylogeographic analysis of Mycobacterium leprae.

Reductive evolution and massive pseudogene formation have shaped the 3.31-Mb genome of Mycobacterium leprae, an unculturable obligate pathogen that causes leprosy in humans. The complete genome sequence of M. leprae strain Br4923 from Brazil was obtained by conventional methods (6x coverage), and Illumina resequencing technology was used to obtain the sequences of strains Thai53 (38x coverage) and NHDP63 (46x coverage) from Thailand and the United States, respectively. Whole-genome comparisons with the previously sequenced TN strain from India revealed that the four strains share 99.995% sequence identity and differ only in 215 polymorphic sites, mainly SNPs, and by 5 pseudogenes. Sixteen interrelated SNP subtypes were defined by genotyping both extant and extinct strains of M. leprae from around the world. The 16 SNP subtypes showed a strong geographical association that reflects the migration patterns of early humans and trade routes, with the Silk Road linking Europe to China having contributed to the spread of leprosy.

Alterations of the 5'untranslated region of SLC16A12 lead to age-related cataract.

PURPOSE. Knowledge of genetic factors predisposing to age-related cataract is very limited. The aim of this study was to identify DNA sequences that either lead to or predispose for this disease. METHODS. The candidate gene SLC16A12, which encodes a solute carrier of the monocarboxylate transporter family, was sequenced in 484 patients with cataract (134 with juvenile cataract, 350 with age-related cataract) and 190 control subjects. Expression studies included luciferase reporter assay and RT-PCR experiments. RESULTS. One patient with age-related cataract showed a novel heterozygous mutation (c.-17A>G) in the 5'untranslated region (5'UTR). This mutation is in cis with the minor G-allele of the single nucleotide polymorphism (SNP) rs3740030 (c.-42T/G), also within the 5'UTR. Using a luciferase reporter assay system, a construct with the patient's haplotype caused a significant upregulation of luciferase activity. In comparison, the SNP G-allele alone promoted less activity, but that amount was still significantly higher than the amount of the common T-allele. Analysis of SLC16A12 transcripts in surrogate tissue demonstrated striking allele-specific differences causing 5'UTR heterogeneity with respect to sequence and quantity. These differences in gene expression were mirrored in an allele-specific predisposition to age-related cataract, as determined in a Swiss population (odds ratio approximately 2.2; confidence intervals, 1.23-4.3). CONCLUSIONS. The monocarboxylate transporter SLC16A12 may contribute to age-related cataract. Sequences within the 5'UTR modulate translational efficiency with pathogenic consequences.

FTO genotype is associated with phenotypic variability of body mass index.

There is evidence across several species for genetic control of phenotypic variation of complex traits, such that the variance among phenotypes is genotype dependent. Understanding genetic control of variability is important in evolutionary biology, agricultural selection programmes and human medicine, yet for complex traits, no individual genetic variants associated with variance, as opposed to the mean, have been identified. Here we perform a meta-analysis of genome-wide association studies of phenotypic variation using ∼170,000 samples on height and body mass index (BMI) in human populations. We report evidence that the single nucleotide polymorphism (SNP) rs7202116 at the FTO gene locus, which is known to be associated with obesity (as measured by mean BMI for each rs7202116 genotype), is also associated with phenotypic variability. We show that the results are not due to scale effects or other artefacts, and find no other experiment-wise significant evidence for effects on variability, either at loci other than FTO for BMI or at any locus for height. The difference in variance for BMI among individuals with opposite homozygous genotypes at the FTO locus is approximately 7%, corresponding to a difference of ∼0.5 kilograms in the standard deviation of weight. Our results indicate that genetic variants can be discovered that are associated with variability, and that between-person variability in obesity can partly be explained by the genotype at the FTO locus. The results are consistent with reported FTO by environment interactions for BMI, possibly mediated by DNA methylation. Our BMI results for other SNPs and our height results for all SNPs suggest that most genetic variants, including those that influence mean height or mean BMI, are not associated with phenotypic variance, or that their effects on variability are too small to detect even with samples sizes greater than 100,000.

Enzyme polymorphism and clinical variability of diseases: a study of diabetes mellitus.

We investigated possible relations among four common neonatal manifestations of diabetic pregnancy (macrosomia, hypoglycemia, hypocalcemia, jaundice) and four enzyme polymorphisms (PGM1, ADA, AK1, ACP1 in a sample of infants born of diabetic mothers. The pattern of associations observed between the two sets of variables is consistent with known differences in enzymatic activity within phenotypes of each system, suggesting that low enzymatic activity may have unfavorable effects on fetal development and on adaptability of the neonate to the extrauterine environment, Some of the polymorphic enzymes studied influence fetal growth in normal pregnancy as well. Analysis of relations between genetic polymorphisms and the clinical pattern of common diseases may provide a better understanding of the genetic basis of the clinical variability of diseases within and between human populations.