A MULTILEVEL MODEL TO ADDRESS BATCH EFFECTS IN COPY NUMBER ESTIMATION USING SNP ARRAYS

Submicroscopic changes in chromosomal DNA copy number dosage are common and have been implicated in many heritable diseases and cancers. Recent high-throughput technologies have a resolution that permits the detection of segmental changes in DNA copy number that span thousands of basepairs across the genome. Genome-wide association studies (GWAS) may simultaneously screen for copy number-phenotype and SNP-phenotype associations as part of the analytic strategy. However, genome-wide array analyses are particularly susceptible to batch effects as the logistics of preparing DNA and processing thousands of arrays often involves multiple laboratories and technicians, or changes over calendar time to the reagents and laboratory equipment. Failure to adjust for batch effects can lead to incorrect inference and requires inefficient post-hoc quality control procedures that exclude regions that are associated with batch. Our work extends previous model-based approaches for copy number estimation by explicitly modeling batch effects and using shrinkage to improve locus-specific estimates of copy number uncertainty. Key features of this approach include the use of diallelic genotype calls from experimental data to estimate batch- and locus-specific parameters of background and signal without the requirement of training data. We illustrate these ideas using a study of bipolar disease and a study of chromosome 21 trisomy. The former has batch effects that dominate much of the observed variation in quantile-normalized intensities, while the latter illustrates the robustness of our approach to datasets where as many as 25% of the samples have altered copy number. Locus-specific estimates of copy number can be plotted on the copy-number scale to investigate mosaicism and guide the choice of appropriate downstream approaches for smoothing the copy number as a function of physical position. The software is open source and implemented in the R package CRLMM available at Bioconductor (http:www.bioconductor.org).

EFFICIENT EVALUATION OF RANKING PROCEDURES WHEN THE NUMBER OF UNITS IS LARGE WITH APPLICATION TO SNP IDENTIFICATION

Simulation-based assessment is a popular and frequently necessary approach to evaluation of statistical procedures. Sometimes overlooked is the ability to take advantage of underlying mathematical relations and we focus on this aspect. We show how to take advantage of large-sample theory when conducting a simulation using the analysis of genomic data as a motivating example. The approach uses convergence results to provide an approximation to smaller-sample results, results that are available only by simulation. We consider evaluating and comparing a variety of ranking-based methods for identifying the most highly associated SNPs in a genome-wide association study, derive integral equation representations of the pre-posterior distribution of percentiles produced by three ranking methods, and provide examples comparing performance. These results are of interest in their own right and set the framework for a more extensive set of comparisons.

Novel genetic risk markers for ulcerative colitis in the IL2/IL21 region are in epistasis with IL23R and suggest a common genetic background for ulcerative colitis and celiac disease

OBJECTIVES: Recently, a genome-wide association study showed that single-nucleotide polymorphisms (SNPs) in the chromosome 4q27 region containing IL2 and IL21 are associated with celiac disease. Given the increased prevalence of inflammatory bowel disease (IBD) among celiac disease patients, we investigated the possible involvement of these SNPs in IBD. METHODS: Five SNPs strongly associated with celiac disease within the KIAA1109/TENR/IL2/IL21 linkage disequilibrium block on chromosome 4q27 and one coding SNP within the IL21 gene were analyzed in a large German IBD cohort. The study population comprised a total of 2,948 Caucasian individuals, including 1,461 IBD patients (ulcerative colitis (UC): n=514, Crohn's disease (CD): n=947) and 1,487 healthy unrelated controls. RESULTS: Three of the five celiac disease risk markers had a protective effect on UC susceptibility, and this effect remained significant after correcting for multiple testing: rs6840978: P=0.0082, P(corr)=0.049, odds ratio (OR) 0.77, 95% confidence interval (CI) 0.63-0.93; rs6822844: P=0.0028, P(corr)=0.017, OR 0.73, 95% CI 0.59-0.90; rs13119723: P=0.0058, P(corr)=0.035, OR 0.75, 95% CI 0.61-0.92. A haplotype consisting of the six SNPs tested was markedly associated with UC susceptibility (P=0.0025, P(corr)=0.015, OR 0.72, 95% CI 0.58-0.89). Moreover, in UC, epistasis was observed between the IL23R SNP rs1004819 and three SNPs in the KIAA1109/TENR/IL2/IL21 block (rs13151961, rs13119723, and rs6822844). CONCLUSIONS: Similar to other autoimmune diseases such as celiac disease, rheumatoid arthritis, type 1 diabetes, Graves' disease, and psoriatic arthritis, genetic variation in the chromosome 4q27 region predisposes to UC, suggesting a common genetic background for these diseases.

Host genetic determinants of spontaneous hepatitis C clearance

Acute infection with the hepatitis C virus (HCV) induces a wide range of innate and adaptive immune responses. A total of 20-50% of acutely HCV-infected individuals permanently control the virus, referred to as 'spontaneous hepatitis C clearance', while the infection progresses to chronic hepatitis C in the majority of cases. Numerous studies have examined host genetic determinants of hepatitis C infection outcome and revealed the influence of genetic polymorphisms of human leukocyte antigens, killer immunoglobulin-like receptors, chemokines, interleukins and interferon-stimulated genes on spontaneous hepatitis C clearance. However, most genetic associations were not confirmed in independent cohorts, revealed opposing results in diverse populations or were limited by varying definitions of hepatitis C outcomes or small sample size. Coordinated efforts are needed in the search for key genetic determinants of spontaneous hepatitis C clearance that include well-conducted candidate genetic and genome-wide association studies, direct sequencing and follow-up functional studies.

Exploring the genetic link between RLS and ADHD

Attention deficit/hyperactivity disorder (ADHD) is a highly heritable neurodevelopmental disorder of childhood onset. Clinical and biological evidence points to shared common central nervous system (CNS) pathology of ADHD and restless legs syndrome (RLS). It was hypothesized that variants previously found to be associated with RLS in two large genome-wide association studies (GWA), will also be associated with ADHD. SNPs located in MEIS1 (rs2300478), BTBD9 (rs9296249, rs3923809, rs6923737), and MAP2K5 (rs12593813, rs4489954) as well as three SNPs tagging the identified haplotype in MEIS1 (rs6710341, rs12469063, rs4544423) were genotyped in a well characterized German sample of 224 families comprising one or more affected sibs (386 children) and both parents. We found no evidence for preferential transmission of the hypothesized variants to ADHD. Subsequent analyses elicited nominal significant association with haplotypes consisting of the three SNPs in BTBD9 (chi2 = 14.8, df = 7, nominal p = 0.039). According to exploratory post hoc analyses, the major contribution to this finding came from the A-A-A-haplotype with a haplotype-wise nominal p-value of 0.009. However, this result did not withstand correction for multiple testing. In view of our results, RLS risk alleles may have a lower effect on ADHD than on RLS or may not be involved in ADHD. The negative findings may additionally result from genetic heterogeneity of ADHD, i.e. risk alleles for RLS may only be relevant for certain subtypes of ADHD. Genes relevant to RLS remain interesting candidates for ADHD; particularly BTBD9 needs further study, as it has been related to iron storage, a potential pathophysiological link between RLS and certain subtypes of ADHD.

Accumulating mutations in series of haplotypes at the KIT and MITF loci are major determinants of white markings in Franches-Montagnes horses.

Coat color and pattern variations in domestic animals are frequently inherited as simple monogenic traits, but a number are known to have a complex genetic basis. While the analysis of complex trait data remains a challenge in all species, we can use the reduced haplotypic diversity in domestic animal populations to gain insight into the genomic interactions underlying complex phenotypes. White face and leg markings are examples of complex traits in horses where little is known of the underlying genetics. In this study, Franches-Montagnes (FM) horses were scored for the occurrence of white facial and leg markings using a standardized scoring system. A genome-wide association study (GWAS) was performed for several white patterning traits in 1,077 FM horses. Seven quantitative trait loci (QTL) affecting the white marking score with p-values p≤10(-4) were identified. Three loci, MC1R and the known white spotting genes, KIT and MITF, were identified as the major loci underlying the extent of white patterning in this breed. Together, the seven loci explain 54% of the genetic variance in total white marking score, while MITF and KIT alone account for 26%. Although MITF and KIT are the major loci controlling white patterning, their influence varies according to the basic coat color of the horse and the specific body location of the white patterning. Fine mapping across the MITF and KIT loci was used to characterize haplotypes present. Phylogenetic relationships among haplotypes were calculated to assess their selective and evolutionary influences on the extent of white patterning. This novel approach shows that KIT and MITF act in an additive manner and that accumulating mutations at these loci progressively increase the extent of white markings.

A mutation in the SUV39H2 gene in Labrador Retrievers with hereditary nasal parakeratosis (HNPK) provides insights into the epigenetics of keratinocyte differentiation.

Hereditary nasal parakeratosis (HNPK), an inherited monogenic autosomal recessive skin disorder, leads to crusts and fissures on the nasal planum of Labrador Retrievers. We performed a genome-wide association study (GWAS) using 13 HNPK cases and 23 controls. We obtained a single strong association signal on chromosome 2 (p(raw) = 4.4×10⁻¹⁴). The analysis of shared haplotypes among the 13 cases defined a critical interval of 1.6 Mb with 25 predicted genes. We re-sequenced the genome of one case at 38× coverage and detected 3 non-synonymous variants in the critical interval with respect to the reference genome assembly. We genotyped these variants in larger cohorts of dogs and only one was perfectly associated with the HNPK phenotype in a cohort of more than 500 dogs. This candidate causative variant is a missense variant in the SUV39H2 gene encoding a histone 3 lysine 9 (H3K9) methyltransferase, which mediates chromatin silencing. The variant c.972T>G is predicted to change an evolutionary conserved asparagine into a lysine in the catalytically active domain of the enzyme (p.N324K). We further studied the histopathological alterations in the epidermis in vivo. Our data suggest that the HNPK phenotype is not caused by hyperproliferation, but rather delayed terminal differentiation of keratinocytes. Thus, our data provide evidence that SUV39H2 is involved in the epigenetic regulation of keratinocyte differentiation ensuring proper stratification and tight sealing of the mammalian epidermis.

A frameshift mutation in the cubilin gene (CUBN) in Border Collies with Imerslund-Gräsbeck syndrome (selective cobalamin malabsorption).

Imerslund-Gräsbeck syndrome (IGS) or selective cobalamin malabsorption has been described in humans and dogs. IGS occurs in Border Collies and is inherited as a monogenic autosomal recessive trait in this breed. Using 7 IGS cases and 7 non-affected controls we mapped the causative mutation by genome-wide association and homozygosity mapping to a 3.53 Mb interval on chromosome 2. We re-sequenced the genome of one affected dog at ∼10× coverage and detected 17 non-synonymous variants in the critical interval. Two of these non-synonymous variants were in the cubilin gene (CUBN), which is known to play an essential role in cobalamin uptake from the ileum. We tested these two CUBN variants for association with IGS in larger cohorts of dogs and found that only one of them was perfectly associated with the phenotype. This variant, a single base pair deletion (c.8392delC), is predicted to cause a frameshift and premature stop codon in the CUBN gene. The resulting mutant open reading frame is 821 codons shorter than the wildtype open reading frame (p.Q2798Rfs*3). Interestingly, we observed an additional nonsense mutation in the MRC1 gene encoding the mannose receptor, C type 1, which was in perfect linkage disequilibrium with the CUBN frameshift mutation. Based on our genetic data and the known role of CUBN for cobalamin uptake we conclude that the identified CUBN frameshift mutation is most likely causative for IGS in Border Collies.

A COL11A2 mutation in Labrador retrievers with mild disproportionate dwarfism.

We describe a mild form of disproportionate dwarfism in Labrador Retrievers, which is not associated with any obvious health problems such as secondary arthrosis. We designate this phenotype as skeletal dysplasia 2 (SD2). It is inherited as a monogenic autosomal recessive trait with incomplete penetrance primarily in working lines of the Labrador Retriever breed. Using 23 cases and 37 controls we mapped the causative mutation by genome-wide association and homozygosity mapping to a 4.44 Mb interval on chromosome 12. We re-sequenced the genome of one affected dog at 30x coverage and detected 92 non-synonymous variants in the critical interval. Only two of these variants, located in the lymphotoxin A (LTA) and collagen alpha-2(XI) chain gene (COL11A2), respectively, were perfectly associated with the trait. Previously described COL11A2 variants in humans or mice lead to skeletal dysplasias and/or deafness. The dog variant associated with disproportionate dwarfism, COL11A2:c.143G>C or p.R48P, probably has only a minor effect on collagen XI function, which might explain the comparatively mild phenotype seen in our study. The identification of this candidate causative mutation thus widens the known phenotypic spectrum of COL11A2 mutations. We speculate that non-pathogenic COL11A2 variants might even contribute to the heritable variation in height.

A non-coding genomic duplication at the HMX1 locus is associated with crop ears in highland cattle.

Highland cattle with congenital crop ears have notches of variable size on the tips of both ears. In some cases, cartilage deformation can be seen and occasionally the external ears are shortened. We collected 40 cases and 80 controls across Switzerland. Pedigree data analysis confirmed a monogenic autosomal dominant mode of inheritance with variable expressivity. All affected animals could be traced back to a single common ancestor. A genome-wide association study was performed and the causative mutation was mapped to a 4 Mb interval on bovine chromosome 6. The H6 family homeobox 1 (HMX1) gene was selected as a positional and functional candidate gene. By whole genome re-sequencing of an affected Highland cattle, we detected 6 non-synonymous coding sequence variants and two variants in an ultra-conserved element at the HMX1 locus with respect to the reference genome. Of these 8 variants, only a non-coding 76 bp genomic duplication (g.106720058_106720133dup) located in the conserved region was perfectly associated with crop ears. The identified copy number variation probably results in HMX1 misregulation and possible gain-of-function. Our findings confirm the role of HMX1 during the development of the external ear. As it is sometimes difficult to phenotypically diagnose Highland cattle with slight ear notches, genetic testing can now be used to improve selection against this undesired trait.

Contribution of Genetic Background, Traditional Risk Factors, and HIV-Related Factors to Coronary Artery Disease Events in HIV-Positive Persons

Background Persons infected with human immunodeficiency virus (HIV) have increased rates of coronary artery disease (CAD). The relative contribution of genetic background, HIV-related factors, antiretroviral medications, and traditional risk factors to CAD has not been fully evaluated in the setting of HIV infection. Methods In the general population, 23 common single-nucleotide polymorphisms (SNPs) were shown to be associated with CAD through genome-wide association analysis. Using the Metabochip, we genotyped 1875 HIV-positive, white individuals enrolled in 24 HIV observational studies, including 571 participants with a first CAD event during the 9-year study period and 1304 controls matched on sex and cohort. Results A genetic risk score built from 23 CAD-associated SNPs contributed significantly to CAD (P = 2.9×10−4). In the final multivariable model, participants with an unfavorable genetic background (top genetic score quartile) had a CAD odds ratio (OR) of 1.47 (95% confidence interval [CI], 1.05–2.04). This effect was similar to hypertension (OR = 1.36; 95% CI, 1.06–1.73), hypercholesterolemia (OR = 1.51; 95% CI, 1.16–1.96), diabetes (OR = 1.66; 95% CI, 1.10–2.49), ≥1 year lopinavir exposure (OR = 1.36; 95% CI, 1.06–1.73), and current abacavir treatment (OR = 1.56; 95% CI, 1.17–2.07). The effect of the genetic risk score was additive to the effect of nongenetic CAD risk factors, and did not change after adjustment for family history of CAD. Conclusions In the setting of HIV infection, the effect of an unfavorable genetic background was similar to traditional CAD risk factors and certain adverse antiretroviral exposures. Genetic testing may provide prognostic information complementary to family history of CAD.

Multiple independent genetic factors at NOS1AP modulate the QT interval in a multi-ethnic population.

Extremes of electrocardiographic QT interval are associated with increased risk for sudden cardiac death (SCD); thus, identification and characterization of genetic variants that modulate QT interval may elucidate the underlying etiology of SCD. Previous studies have revealed an association between a common genetic variant in NOS1AP and QT interval in populations of European ancestry, but this finding has not been extended to other ethnic populations. We sought to characterize the effects of NOS1AP genetic variants on QT interval in the multi-ethnic population-based Dallas Heart Study (DHS, n = 3,072). The SNP most strongly associated with QT interval in previous samples of European ancestry, rs16847548, was the most strongly associated in White (P = 0.005) and Black (P = 3.6 x 10(-5)) participants, with the same direction of effect in Hispanics (P = 0.17), and further showed a significant SNP x sex-interaction (P = 0.03). A second SNP, rs16856785, uncorrelated with rs16847548, was also associated with QT interval in Blacks (P = 0.01), with qualitatively similar results in Whites and Hispanics. In a previously genotyped cohort of 14,107 White individuals drawn from the combined Atherosclerotic Risk in Communities (ARIC) and Cardiovascular Health Study (CHS) cohorts, we validated both the second locus at rs16856785 (P = 7.63 x 10(-8)), as well as the sex-interaction with rs16847548 (P = 8.68 x 10(-6)). These data extend the association of genetic variants in NOS1AP with QT interval to a Black population, with similar trends, though not statistically significant at P<0.05, in Hispanics. In addition, we identify a strong sex-interaction and the presence of a second independent site within NOS1AP associated with the QT interval. These results highlight the consistent and complex role of NOS1AP genetic variants in modulating QT interval.

Common noncoding UMOD gene variants induce salt-sensitive hypertension and kidney damage by increasing uromodulin expression

Hypertension and chronic kidney disease (CKD) are complex traits representing major global health problems1,2. Multiple genome-wide association studies have identified common variants in the promoter of the UMOD gene3–9, which encodes uromodulin, the major protein secreted in normal urine, that cause independent susceptibility to CKD and hypertension. Despite compelling genetic evidence for the association between UMOD risk variants and disease susceptibility in the general population, the underlying biological mechanism is not understood. Here, we demonstrate that UMOD risk variants increased UMOD expression in vitro and in vivo. Uromodulin overexpression in transgenic mice led to salt-sensitive hypertension and to the presence of age-dependent renal lesions similar to those observed in elderly individuals homozygous for UMOD promoter risk variants. The link between uromodulin and hypertension is due to activation of the renal sodium cotransporter NKCC2. We demonstrated the relevance of this mechanism in humans by showing that pharmacological inhibition of NKCC2 was more effective in lowering blood pressure in hypertensive patients who are homozygous for UMOD promoter risk variants than in other hypertensive patients. Our findings link genetic susceptibility to hypertension and CKD to the level of uromodulin expression and uromodulin’s effect on salt reabsorption in the kidney. These findings point to uromodulin as a therapeutic target for lowering blood pressure and preserving renal function.

Gene by BMI Interactions Influencing C-Reactive Protein Levels in European-Americans

C-Reactive Protein (CRP) is a biomarker indicating tissue damage, inflammation, and infection. High-sensitivity CRP (hsCRP) is an emerging biomarker often used to estimate an individual’s risk for future coronary heart disease (CHD). hsCRP levels falling below 1.00 mg/l indicate a low risk for developing CHD, levels ranging between 1.00 mg/l and 3.00 mg/l indicate an elevated risk, and levels exceeding 3.00 mg/l indicate high risk. Multiple Genome-Wide Association Studies (GWAS) have identified a number of genetic polymorphisms which influence CRP levels. SNPs implicated in such studies have been found in or near genes of interest including: CRP, APOE, APOC, IL-6, HNF1A, LEPR, and GCKR. A strong positive correlation has also been found to exist between CRP levels and BMI, a known risk factor for CHD and a state of chronic inflammation. We conducted a series of analyses designed to identify loci which interact with BMI to influence CRP levels in a subsample of European-Americans in the ARIC cohort. In a stratified GWA analysis, 15 genetic regions were identified as having significantly (p-value < 2.00*10-3) distinct effects on hsCRP levels between the two obesity strata: lean (18.50 kg/m2 < BMI < 24.99 kg/m2) and obese (BMI ≥ 30.00 kg/m2). A GWA analysis performed on all individuals combined (i.e. not a priori stratified for obesity status) with the inclusion of an additional parameter for BMI by gene interaction, identified 11 regions which interact with BMI to influence hsCRP levels. Two regions containing the genes GJA5 and GJA8 (on chromosome 1) and FBXO11 (on chromosome 2) were identified in both methods of analysis suggesting that these genes possibly interact with BMI to influence hsCRP levels. We speculate that atrial fibrillation (AF), age-related cataracts and the TGF-β pathway may be the biological processes influenced by the interaction of GJA5, GJA8 and FBXO11, respectively, with BMI to cause changes in hsCRP levels. Future studies should focus on the influence of gene x bmi interaction on AF, age-related cataracts and TGF-β.

THE EXPRESSION PROFILE OF A GWAS DETECTED GENE (NINJ2) IN THE BRAIN OF A MOUSE MODEL OF ISCHEMIC STROKE

Stroke is the third leading cause of death and a major debilitating disease in the United States. Multiple factors, including genetic factors, contribute to the development of the disease. Genome-wide association studies (GWAS) have contributed to the identification of genetic loci influencing risk for complex diseases, such as stroke. In 2010, a GWAS of incident stroke was performed in four large prospective cohorts from the USA and Europe and identified an association of two Single Nucleotide Polymorphisms (SNPs) on chromosome 12p13 with a greater risk of ischemic stroke in individuals of European and African-American ancestry. These SNPs are located 11 Kb upstream of the nerve injury-induced gene 2, Ninjurin2 (NINJ2), suggesting that this gene may be involved in stroke pathogenesis. NINJ2 is a cell adhesion molecule induced in the distal glial cells from a sciatic-nerve injury at 7-days after injury. In an effort to ascribe a possible role of NINJ2 in stroke, we have assessed changes in the level of gene and protein expression of NINJ2 following a time-course from a transiently induced middle cerebral artery ischemic stroke in mice brains. We report an increase in the gene expression of NINJ2 in the ischemic and peri-infarct (ipsilateral) cortical tissues at 7 and 14-days after stroke. We also report an increase in the protein expression of NINJ2 in the cortex of both the ipsilateral and contralateral cortical tissues at the same time-points. We conclude that the expression of NINJ2 is regulated by an ischemic stroke in the cortex and is consistent with NINJ2 being involved in the recovery time-points of stroke.