Six Novel Loci Associated with Circulating VEGF Levels Identified by a Meta-analysis of Genome-Wide Association Studies

Vascular endothelial growth factor (VEGF) is an angiogenic and neurotrophic factor, secreted by endothelial cells, known to impact various physiological and disease processes from cancer to cardiovascular disease and to be pharmacologically modifiable. We sought to identify novel loci associated with circulating VEGF levels through a genome-wide association meta-analysis combining data from European-ancestry individuals and using a dense variant map from 1000 genomes imputation panel. Six discovery cohorts including 13,312 samples were analyzed, followed by in-silico and de-novo replication studies including an additional 2,800 individuals. A total of 10 genome-wide significant variants were identified at 7 loci. Four were novel loci (5q14.3, 10q21.3, 16q24.2 and 18q22.3) and the leading variants at these loci were rs114694170 (MEF2C, P = 6.79x10^-13), rs74506613 (JMJD1C, P = 1.17x10^-19), rs4782371 (ZFPM1, P = 1.59x10^-9) and rs2639990 (ZADH2, P = 1.72x10^-8), respectively. We also identified two new independent variants (rs34528081, VEGFA, P = 1.52x10^-18; rs7043199, VLDLR-AS1, P = 5.12x10^-14) at the 3 previously identified loci and strengthened the evidence for the four previously identified SNPs (rs6921438, LOC100132354, P = 7.39x10^-1467; rs1740073, C6orf223, P = 2.34x10^-17; rs6993770, ZFPM2, P = 2.44x10^-60; rs2375981, KCNV2, P = 1.48x10^-100). These variants collectively explained up to 52% of the VEGF phenotypic variance. We explored biological links between genes in the associated loci using Ingenuity Pathway Analysis that emphasized their roles in embryonic development and function. Gene set enrichment analysis identified the ERK5 pathway as enriched in genes containing VEGF associated variants. eQTL analysis showed, in three of the identified regions, variants acting as both cis and trans eQTLs for multiple genes. Most of these genes, as well as some of those in the associated loci, were involved in platelet biogenesis and functionality, suggesting the importance of this process in regulation of VEGF levels. This work also provided new insights into the involvement of genes implicated in various angiogenesis related pathologies in determining circulating VEGF levels. The understanding of the molecular mechanisms by which the identified genes affect circulating VEGF levels could be important in the development of novel VEGF-related therapies for such diseases.     

A Comprehensive Analysis of Common and Rare Variants to Identify Adiposity Loci in Hispanic Americans: The IRAS Family Study (IRASFS)

Obesity is growing epidemic affecting 35% of adults in the United States. Previous genome-wide association studies (GWAS) have identified numerous loci associated with obesity. However, the majority of studies have been completed in Caucasians focusing on total body measures of adiposity. Here we report the results from genome-wide and exome chip association studies focusing on total body measures of adiposity including body mass index (BMI), percent body fat (PBF) and measures of fat deposition including waist circumference (WAIST), waist-hip ratio (WHR), subcutaneous adipose tissue (SAT), and visceral adipose tissue (VAT) in Hispanic Americans (n_max = 1263) from the Insulin Resistance Atherosclerosis Family Study (IRASFS). Five SNPs from two novel loci attained genome-wide significance (P<5.00x10^-8) in IRASFS. A missense SNP in the isocitrate dehydrogenase 1 gene (IDH1) was associated with WAIST (rs34218846, MAF = 6.8%, P_DOM = 1.62x10^-8). This protein is postulated to play an important role in fat and cholesterol biosynthesis as demonstrated in cell and knock-out animal models. Four correlated intronic SNPs in the Zinc finger, GRF-type containing 1 gene (ZGRF1; SNP rs1471880, MAF = 48.1%, P_DOM = 1.00x10^-8) were strongly associated with WHR. The exact biological function of ZGRF1 and the connection with adiposity remains unclear. SNPs with p-values less than 5.00x10^-6 from IRASFS were selected for replication. Meta-analysis was computed across seven independent Hispanic-American cohorts (n_max = 4156) and the strongest signal was rs1471880 (P_DOM = 8.38x10^-6) in ZGRF1 with WAIST. In conclusion, a genome-wide and exome chip association study was conducted that identified two novel loci (IDH1 and ZGRF1) associated with adiposity. While replication efforts were inconclusive, when taken together with the known biology, IDH1 and ZGRF1 warrant further evaluation.     

A genome-wide significant association on chromosome 15 for congenital entropion in Swiss White Alpine sheep

Entropion is a known congenital disorder in sheep presumed to be heritable but no causative genetic variant has been reported. Affected lambs show a variable inward rolling of the lower eyelids leading to blindness in severe cases. In Switzerland, the Swiss White Alpine (SWA) breed showed a significantly higher prevalence for entropion than other breeds. A GWAS using 150 SWA sheep (90 affected lambs and 60 controls), based on 600k SNP data, revealed a genome-wide significant signal on chromosome 15. The 0.2 Mb associated region contains functional candidate genes, SMTNL1 and CTNND1. Pathogenic variants in human CTNND1 cause blepharocheilodontic syndrome 2, a rare disorder including eyelid anomalies, and SMTNL1 regulates contraction and relaxation of skeletal and smooth muscle. WGS of a single entropion-affected lamb revealed two private missense variants in SMTNL1 and CTNND1. Subsequent genotyping of both variants in 231 phenotyped SWA sheep was performed. The SMTNL1 variant p.(Asp452Asn) affects an evolutionary conserved residue within an important domain and represents a rare allele, which occurred also in controls. The p.(Glu943Lys) variant in CTNND1 represents a common variant unlikely to cause entropion as the mutant allele occurred more frequently in non-affected sheep. Therefore, we propose that these protein-changing variants are unlikely to explain the phenotype. Additionally, WGS of three further disconcordant pairs of full siblings was carried out but revealed no obvious causative variant. Finally, we conclude that entropion represents a more complex disease caused by different non-coding regulatory variants.     

Genome-wide association study of body weight in Australian Merino sheep reveals an orthologous region on OAR6 to human and bovine genomic regions affecting height and weight

Background

Body weight (BW) is an important trait for meat production in sheep. Although over the past few years, numerous quantitative trait loci (QTL) have been detected for production traits in cattle, few QTL studies have been reported for sheep, with even fewer on meat production traits. Our objective was to perform a genome-wide association study (GWAS) with the medium-density Illumina Ovine SNP50 BeadChip to identify genomic regions and corresponding haplotypes associated with BW in Australian Merino sheep.

Methods

A total of 1781 Australian Merino sheep were genotyped using the medium-density Illumina Ovine SNP50 BeadChip. Among the 53 862 single nucleotide polymorphisms (SNPs) on this array, 48 640 were used to perform a GWAS using a linear mixed model approach. Genotypes were phased with hsphase; to estimate SNP haplotype effects, linkage disequilibrium blocks were identified in the detected QTL region.

Results

Thirty-nine SNPs were associated with BW at a Bonferroni-corrected genome-wide significance threshold of 1 %. One region on sheep (Ovis aries) chromosome 6 (OAR6) between 36.15 and 38.56 Mb, included 13 significant SNPs that were associated with BW; the most significant SNP was OAR6_41936490.1 (P = 2.37 × 10⁻¹⁶) at 37.69 Mb with an allele substitution effect of 2.12 kg, which corresponds to 0.248 phenotypic standard deviations for BW. The region that surrounds this association signal on OAR6 contains three genes: leucine aminopeptidase 3 (LAP3), which is involved in the processing of the oxytocin precursor; NCAPG non-SMC condensin I complex, subunit G (NCAPG), which is associated with foetal growth and carcass size in cattle; and ligand dependent nuclear receptor corepressor-like (LCORL), which is associated with height in humans and cattle.

Conclusions

The GWAS analysis detected 39 SNPs associated with BW in sheep and a major QTL region was identified on OAR6. In several other mammalian species, regions that are syntenic with this region have been found to be associated with body size traits, which may reflect that the underlying biological mechanisms share a common ancestry. These findings should facilitate the discovery of causative variants for BW and contribute to marker-assisted selection.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0142-4) contains supplementary material, which is available to authorized users.     

Novel Genetic Locus Implicated for HIV-1 Acquisition with Putative Regulatory Links to HIV Replication and Infectivity: A Genome-Wide Association Study

Fifty percent of variability in HIV-1 susceptibility is attributable to host genetics. Thus identifying genetic associations is essential to understanding pathogenesis of HIV-1 and important for targeting drug development. To date, however, CCR5 remains the only gene conclusively associated with HIV acquisition. To identify novel host genetic determinants of HIV-1 acquisition, we conducted a genome-wide association study among a high-risk sample of 3,136 injection drug users (IDUs) from the Urban Health Study (UHS). In addition to being IDUs, HIV- controls were frequency-matched to cases on environmental exposures to enhance detection of genetic effects. We tested independent replication in the Women’s Interagency HIV Study (N=2,533). We also examined publicly available gene expression data to link SNPs associated with HIV acquisition to known mechanisms affecting HIV replication/infectivity. Analysis of the UHS nominated eight genetic regions for replication testing. SNP rs4878712 in FRMPD1 met multiple testing correction for independent replication (P=1.38x10^-4), although the UHS-WIHS meta-analysis p-value did not reach genome-wide significance (P=4.47x10^-7 vs. P<5.0x10^-8) Gene expression analyses provided promising biological support for the protective G allele at rs4878712 lowering risk of HIV: (1) the G allele was associated with reduced expression of FBXO10 (r=-0.49, P=6.9x10^-5); (2) FBXO10 is a component of the Skp1-Cul1-F-box protein E3 ubiquitin ligase complex that targets Bcl-2 protein for degradation; (3) lower FBXO10 expression was associated with higher BCL2 expression (r=-0.49, P=8x10^-5); (4) higher basal levels of Bcl-2 are known to reduce HIV replication and infectivity in human and animal in vitro studies. These results suggest new potential biological pathways by which host genetics affect susceptibility to HIV upon exposure for follow-up in subsequent studies.     

A Genome-Wide Association Study Uncovers a Genetic Locus Associated with Thoracic-to-Hip Ratio in Koreans

The thoracic-to-hip circumference ratio (THR) is an anthropometric marker recently described as a predictor of type 2 diabetes. In this study, we performed a genome-wide association study (GWAS) followed by confirmatory analyses to identify genetic markers associated with THR. A total of 7,240 Korean subjects (4,988 for the discovery stage and 2,252 for the confirmatory analyses) were recruited for this study, and genome-wide single nucleotide polymorphism (SNP) genotyping of the initial 4,988 individuals was performed using Affymetrix Human SNP array 5.0. Linear regression analysis was then performed to adjust for the effects of age, sex, and current diabetes medication status on the THR of the study subjects. In the initial discovery stage, there was a statistically nominal association between minor alleles of SNP markers on chromosomes 4, 8, 10, and 12, and THR changes (p < 5.0 × 10⁻⁶). The subsequent confirmatory analyses of these markers, however, only detected a significant association between two SNPs in the HECTD4 gene and decreased THRs. Notably, this association was detected in male (rs11066280: p = 1.14 × 10⁻²; rs2074356: p = 1.10 × 10⁻²), but not in female subjects. Meanwhile, the combined results from the two analyses (initial and confirmatory) indicated that minor alleles of these two intronic variants exhibited a significant genome-wide association with decreased THR in the male subjects (n = 3,155; rs11066280: effect size = −0.008624, p = 6.19 × 10⁻⁹; rs2074356: effect size = −0.008762, p = 1.89 × 10⁻⁸). Furthermore, minor alleles of these two SNPs exhibited protective effects on patients’ risks for developing type 2 diabetes. In conclusion, we have identified two genetic variations in HECTD4 that are associated with THR, particularly in men.     

Genome-Wide Association Study of Parity in Bangladeshi Women

Human fertility is a complex trait determined by gene-environment interactions in which genetic factors represent a significant component. To better understand inter-individual variability in fertility, we performed one of the first genome-wide association studies (GWAS) of common fertility phenotypes, lifetime number of pregnancies and number of children in a developing country population. The fertility phenotype data and DNA samples were obtained at baseline recruitment from individuals participating in a large prospective cohort study in Bangladesh. GWAS analyses of fertility phenotypes were conducted among 1,686 married women. One SNP on chromosome 4 was non-significantly associated with number of children at P <10^-7 and number of pregnancies at P <10^-6. This SNP is located in a region without a gene within 1 Mb. One SNP on chromosome 6 was non-significantly associated with extreme number of children at P <10^-6. The closest gene to this SNP is HDGFL1, a hepatoma-derived growth factor. When we excluded hormonal contraceptive users, a SNP on chromosome 5 was non-significantly associated at P <10^-5 for number of children and number of pregnancies. This SNP is located near C5orf64, an open reading frame, and ZSWIM6, a zinc ion binding gene. We also estimated the heritability of these phenotypes from our genotype data using GCTA (Genome-wide Complex Trait Analysis) for number of children (h_g ² = 0.149, SE = 0.24, p-value = 0.265) and number of pregnancies (h_g ² = 0.007, SE = 0.22, p-value = 0.487). Our genome-wide association study and heritability estimates of number of pregnancies and number of children in Bangladesh did not confer strong evidence of common variants for parity variation. However, our results suggest that future studies may want to consider the role of 3 notable SNPs in their analysis.     

A comparison in association and linkage genome-wide scans for alcoholism susceptibility genes using single-nucleotide polymorphisms

We conducted genome-wide linkage scans using both microsatellite and single-nucleotide polymorphism (SNP) markers. Regions showing the strongest evidence of linkage to alcoholism susceptibility genes were identified. Haplotype analyses using a sliding-window approach for SNPs in these regions were performed. In addition, we performed a genome-wide association scan using SNP data. SNPs in these regions with evidence of association (P 相似文献   

Identification of Non-HLA Genes Associated with Celiac Disease and Country-Specific Differences in a Large,International Pediatric Cohort

Objectives

There are significant geographical differences in the prevalence and incidence of celiac disease that cannot be explained by HLA alone. More than 40 loci outside of the HLA region have been associated with celiac disease. We investigated the roles of these non-HLA genes in the development of tissue transglutaminase autoantibodies (tTGA) and celiac disease in a large international prospective cohort study.

Methods

A total of 424,788 newborns from the US and European general populations and first-degree relatives with type 1 diabetes were screened for specific HLA genotypes. Of these, 21,589 carried 1 of the 9 HLA genotypes associated with increased risk for type 1 diabetes and celiac disease; we followed 8676 of the children in a 15 y prospective follow-up study. Genotype analyses were performed on 6010 children using the Illumina ImmunoChip. Levels of tTGA were measured in serum samples using radio-ligand binding assays; diagnoses of celiac disease were made based on persistent detection of tTGA and biopsy analysis. Data were analyzed using Cox proportional hazards analyses.

Results

We found 54 single-nucleotide polymorphisms (SNPs) in 5 genes associated with celiac disease (TAGAP, IL18R1, RGS21, PLEK, and CCR9) in time to celiac disease analyses (10⁻⁴>P>5.8x10⁻⁶). The hazard ratios (HR) for the SNPs with the smallest P values in each region were 1.59, 1.45, 2.23, 2.64, and 1.40, respectively. Outside of regions previously associated with celiac disease, we identified 10 SNPs in 8 regions that could also be associated with the disease (P<10⁻⁴). A SNP near PKIA (rs117128341, P = 6.5x10⁻⁸, HR = 2.8) and a SNP near PFKFB3 (rs117139146, P<2.8x10⁻⁷, HR = 4.9) reached the genome-wide association threshold in subjects from Sweden. Analyses of time to detection of tTGA identified 29 SNPs in 2 regions previously associated with celiac disease (CTLA4, P = 1.3x10⁻⁶, HR = 0.76 and LPP, P = 2.8x10⁻⁵, HR = .80) and 6 SNPs in 5 regions not previously associated with celiac disease (P<10⁻⁴); non-HLA genes are therefore involved in development of tTGA.

Conclusions

In conclusion, using a genetic analysis of a large international cohort of children, we associated celiac disease development with 5 non-HLA regions previously associated with the disease and 8 regions not previously associated with celiac disease. We identified 5 regions associated with development of tTGA. Two loci associated with celiac disease progression reached a genome-wide association threshold in subjects from Sweden.     

Genome-Wide Association Studies in Dogs and Humans Identify ADAMTS20 as a Risk Variant for Cleft Lip and Palate

Cleft lip with or without cleft palate (CL/P) is the most commonly occurring craniofacial birth defect. We provide insight into the genetic etiology of this birth defect by performing genome-wide association studies in two species: dogs and humans. In the dog, a genome-wide association study of 7 CL/P cases and 112 controls from the Nova Scotia Duck Tolling Retriever (NSDTR) breed identified a significantly associated region on canine chromosome 27 (unadjusted p=1.1 x 10^-13; adjusted p= 2.2 x 10^-3). Further analysis in NSDTR families and additional full sibling cases identified a 1.44 Mb homozygous haplotype (chromosome 27: 9.29 – 10.73 Mb) segregating with a more complex phenotype of cleft lip, cleft palate, and syndactyly (CLPS) in 13 cases. Whole-genome sequencing of 3 CLPS cases and 4 controls at 15X coverage led to the discovery of a frameshift mutation within ADAMTS20 (c.1360_1361delAA (p.Lys453Ilefs*3)), which segregated concordant with the phenotype. In a parallel study in humans, a family-based association analysis (DFAM) of 125 CL/P cases, 420 unaffected relatives, and 392 controls from a Guatemalan cohort, identified a suggestive association (rs10785430; p =2.67 x 10^-6) with the same gene, ADAMTS20. Sequencing of cases from the Guatemalan cohort was unable to identify a causative mutation within the coding region of ADAMTS20, but four coding variants were found in additional cases of CL/P. In summary, this study provides genetic evidence for a role of ADAMTS20 in CL/P development in dogs and as a candidate gene for CL/P development in humans.     

Genome-Wide Association and Trans-ethnic Meta-Analysis for Advanced Diabetic Kidney Disease: Family Investigation of Nephropathy and Diabetes (FIND)

Diabetic kidney disease (DKD) is the most common etiology of chronic kidney disease (CKD) in the industrialized world and accounts for much of the excess mortality in patients with diabetes mellitus. Approximately 45% of U.S. patients with incident end-stage kidney disease (ESKD) have DKD. Independent of glycemic control, DKD aggregates in families and has higher incidence rates in African, Mexican, and American Indian ancestral groups relative to European populations. The Family Investigation of Nephropathy and Diabetes (FIND) performed a genome-wide association study (GWAS) contrasting 6,197 unrelated individuals with advanced DKD with healthy and diabetic individuals lacking nephropathy of European American, African American, Mexican American, or American Indian ancestry. A large-scale replication and trans-ethnic meta-analysis included 7,539 additional European American, African American and American Indian DKD cases and non-nephropathy controls. Within ethnic group meta-analysis of discovery GWAS and replication set results identified genome-wide significant evidence for association between DKD and rs12523822 on chromosome 6q25.2 in American Indians (P = 5.74x10^-9). The strongest signal of association in the trans-ethnic meta-analysis was with a SNP in strong linkage disequilibrium with rs12523822 (rs955333; P = 1.31x10^-8), with directionally consistent results across ethnic groups. These 6q25.2 SNPs are located between the SCAF8 and CNKSR3 genes, a region with DKD relevant changes in gene expression and an eQTL with IPCEF1, a gene co-translated with CNKSR3. Several other SNPs demonstrated suggestive evidence of association with DKD, within and across populations. These data identify a novel DKD susceptibility locus with consistent directions of effect across diverse ancestral groups and provide insight into the genetic architecture of DKD.     

A GWAS Study on Liver Function Test Using eMERGE Network Participants

Introduction

Liver enzyme levels and total serum bilirubin are under genetic control and in recent years genome-wide population-based association studies have identified different susceptibility loci for these traits. We conducted a genome-wide association study in European ancestry participants from the Electronic Medical Records and Genomics (eMERGE) Network dataset of patient medical records with available genotyping data in order to identify genetic contributors to variability in serum bilirubin levels and other liver function tests and to compare the effects between adult and pediatric populations.

Methods

The process of whole genome imputation of eMERGE samples with standard quality control measures have been described previously. After removing missing data and outliers based on principal components (PC) analyses, 3294 samples from European ancestry were used for the GWAS study. The association between each single nucleotide polymorphism (SNP) and total serum bilirubin and other liver function tests was tested using linear regression, adjusting for age, gender, site, platform and ancestry principal components (PC).

Results

Consistent with previous results, a strong association signal has been detected for UGT1A gene cluster (best SNP rs887829, beta = 0.15, p = 1.30x10^-118) for total serum bilirubin level. Indeed, in this region more than 176 SNPs (or indels) had p<10⁻⁸ spanning 150Kb on the long arm of chromosome 2q37.1. In addition, we found a similar level of magnitude in a pediatric group (p = 8.26x10^-47, beta = 0.17). Further imputation using sequencing data as a reference panel revealed association of other markers including known TA7 repeat indels (rs8175347) (p = 9.78x10^-117) and rs111741722 (p = 5.41x10^-119) which were in proxy (r2 = 0.99) with rs887829. Among rare variants, two Asian subjects homozygous for coding SNP rs4148323 (G71R) were identified. Additional known effects for total serum bilirubin were also confirmed including organic anion transporters SLCO1B1-SLCO1B3, TDRP and ZMYND8 at FDR<0.05 with no gene-gene interaction effects. Phenome-wide association studies (PheWAS) suggest a protective effect of TA7 repeat against cerebrovascular disease in an adult cohort (OR = 0.75, p = 0.0008). Among other liver function tests, we also confirmed the previous effect of the ABO blood group locus for variation in serum alkaline phosphatase (rs579459, p = 9.44x10^-15).

Conclusions

Taken together, our data present interesting findings with strong confirmation of previous effects by simply using the eMERGE electronic health record phenotyping. In addition, our findings indicate that similar to the adult population, the UGT1A1 is the main locus responsible for normal variation of serum bilirubin in pediatric populations.     

Common Variants in CLDN2 and MORC4 Genes Confer Disease Susceptibility in Patients with Chronic Pancreatitis

A recent genome-wide association study (GWAS) identified association with variants in X-linked CLDN2 and MORC4, and PRSS1-PRSS2 loci with chronic pancreatitis (CP) in North American patients of European ancestry. We selected 9 variants from the reported GWAS and replicated the association with CP in Indian patients by genotyping 1807 unrelated Indians of Indo-European ethnicity, including 519 patients with CP and 1288 controls. The etiology of CP was idiopathic in 83.62% and alcoholic in 16.38% of 519 patients. Our study confirmed a significant association of 2 variants in CLDN2 gene (rs4409525—OR 1.71, P = 1.38 x 10^-09; rs12008279—OR 1.56, P = 1.53 x 10^-04) and 2 variants in MORC4 gene (rs12688220—OR 1.72, P = 9.20 x 10^-09; rs6622126—OR 1.75, P = 4.04x10^-05) in Indian patients with CP. We also found significant association at PRSS1-PRSS2 locus (OR 0.60; P = 9.92 x 10^-06) and SAMD12-TNFRSF11B (OR 0.49, 95% CI [0.31–0.78], P = 0.0027). A variant in the gene MORC4 (rs12688220) showed significant interaction with alcohol (OR for homozygous and heterozygous risk allele -14.62 and 1.51 respectively, P = 0.0068) suggesting gene-environment interaction. A combined analysis of the genes CLDN2 and MORC4 based on an effective risk allele score revealed a higher percentage of individuals homozygous for the risk allele in CP cases with 5.09 fold enhanced risk in individuals with 7 or more effective risk alleles compared with individuals with 3 or less risk alleles (P = 1.88 x 10^-14). Genetic variants in CLDN2 and MORC4 genes were associated with CP in Indian patients.     

A genome-wide association study of BMI in American Indians

Numerous studies have been done to understand genetic contributors to BMI, but only a limited number of studies have been done in nonwhite groups such as American Indians. A genome-wide association study (GWAS) for BMI was therefore performed in Pima Indians. BMI measurements from a longitudinal study of 1,120 Pima Indians and 454,194 single-nucleotide polymorphisms (SNPs) from the 1 million Affymetrix SNP panel were used (35% of SNPs were excluded due to minor allele frequency <0.05). Data included BMI measured at multiple examinations collected from 1965 to 2004, as well as the maximum BMI at one of these visits. General and within-family tests were performed using a maximum-likelihood based mixed model procedure. No SNP reached a genome-wide significance level (estimated at P < 4.94 × 10(-7)). For repeated measures analyses, the strongest associations for general and within-family tests mapped to two different regions on chromosome 6 (rs9342220 (P = 1.39 × 10(-6)) and rs7758764 (P = 2.51 × 10(-6)), respectively). For maximum BMI, the strongest association for the general tests mapped to chromosome 4 (rs17612333; P = 1.98 × 10(-6)) and to chromosome 3 (rs11127958; P = 1.53 × 10(-6)) for the within-family tests. Further analysis is important because only a few of these regions have been previously implicated in a GWAS and genetic susceptibility may differ by ethnicity.     

Linkage Analysis in Autoimmune Addison’s Disease: NFATC1 as a Potential Novel Susceptibility Locus

BackgroundAutoimmune Addison’s disease (AAD) is a rare, highly heritable autoimmune endocrinopathy. It is possible that there may be some highly penetrant variants which confer disease susceptibility that have yet to be discovered.MethodsDNA samples from 23 multiplex AAD pedigrees from the UK and Norway (50 cases, 67 controls) were genotyped on the Affymetrix SNP 6.0 array. Linkage analysis was performed using Merlin. EMMAX was used to carry out a genome-wide association analysis comparing the familial AAD cases to 2706 UK WTCCC controls. To explore some of the linkage findings further, a replication study was performed by genotyping 64 SNPs in two of the four linked regions (chromosomes 7 and 18), on the Sequenom iPlex platform in three European AAD case-control cohorts (1097 cases, 1117 controls). The data were analysed using a meta-analysis approach.ResultsIn a parametric analysis, applying a rare dominant model, loci on chromosomes 7, 9 and 18 had LOD scores >2.8. In a non-parametric analysis, a locus corresponding to the HLA region on chromosome 6, known to be associated with AAD, had a LOD score >3.0. In the genome-wide association analysis, a SNP cluster on chromosome 2 and a pair of SNPs on chromosome 6 were associated with AAD (P <5x10^-7). A meta-analysis of the replication study data demonstrated that three chromosome 18 SNPs were associated with AAD, including a non-synonymous variant in the NFATC1 gene.ConclusionThis linkage study has implicated a number of novel chromosomal regions in the pathogenesis of AAD in multiplex AAD families and adds further support to the role of HLA in AAD. The genome-wide association analysis has also identified a region of interest on chromosome 2. A replication study has demonstrated that the NFATC1 gene is worthy of future investigation, however each of the regions identified require further, systematic analysis.     

Common Variation at 1q24.1 (ALDH9A1) Is a Potential Risk Factor for Renal Cancer

So far six susceptibility loci for renal cell carcinoma (RCC) have been discovered by genome-wide association studies (GWAS). To identify additional RCC common risk loci, we performed a meta-analysis of published GWAS (totalling 2,215 cases and 8,566 controls of Western-European background) with imputation using 1000 Genomes Project and UK10K Project data as reference panels and followed up the most significant association signals [22 single nucleotide polymorphisms (SNPs) and 3 indels in eight genomic regions] in 383 cases and 2,189 controls from The Cancer Genome Atlas (TCGA). A combined analysis identified a promising susceptibility locus mapping to 1q24.1 marked by the imputed SNP rs3845536 (P _combined =2.30x10^-8). Specifically, the signal maps to intron 4 of the ALDH9A1 gene (aldehyde dehydrogenase 9 family, member A1). We further evaluated this potential signal in 2,461 cases and 5,081 controls from the International Agency for Research on Cancer (IARC) GWAS of RCC cases and controls from multiple European regions. In contrast to earlier findings no association was shown in the IARC series (P=0.94; P _combined =2.73x10^-5). While variation at 1q24.1 represents a potential risk locus for RCC, future replication analyses are required to substantiate our observation.     

Genetic Determinants of Pelvic Organ Prolapse among African American and Hispanic Women in the Women’s Health Initiative

Current evidence suggests a multifactorial etiology to pelvic organ prolapse (POP), including genetic predisposition. We conducted a genome-wide association study of POP in African American (AA) and Hispanic (HP) women from the Women’s Health Initiative Hormone Therapy study. Cases were defined as any POP (grades 1–3) or moderate/severe POP (grades 2–3), while controls had grade 0 POP. We performed race-specific multiple logistic regression analyses between SNPs imputed to 1000 genomes in relation to POP (grade 0 vs 1–3; grade 0 vs 2–3) adjusting for age at diagnosis, body mass index, parity, and genetic ancestry. There were 1274 controls and 1427 cases of any POP and 317 cases of moderate/severe POP. Although none of the analyses reached genome-wide significance (p<5x10^-8), we noted variants in several loci that met p<10⁻⁶. In race-specific analysis of grade 0 vs 2–3, intronic SNPs in the CPE gene (rs28573326, OR:2.14; 95% CI 1.62–2.83; p = 1.0x10^-7) were associated with POP in AAs, and SNPs in the gene AL132709.5 (rs1950626, OR:2.96; 95% CI 1.96–4.48, p = 2.6x10^-7) were associated with POP in HPs. Inverse variance fixed-effect meta-analysis of the race-specific results showed suggestive signals for SNPs in the DPP6 gene (rs11243354, OR:1.36; p = 4.2x10^-7) in the grade 0 vs 1–3 analyses and for SNPs around PGBD5 (rs740494, OR:2.17; p = 8.6x10^-7) and SHC3 (rs2209875, OR:0.60; p = 9.3x10^-7) in the grade 0 vs 2–3 analyses. While we did not identify genome-wide significant findings, we document several SNPs reaching suggestive statistical significance. Further interrogation of POP in larger minority samples is warranted.     

Paternal genetic variants and risk of obstructive heart defects: A parent-of-origin approach

Previous research on risk factors for obstructive heart defects (OHDs) focused on maternal and infant genetic variants, prenatal environmental exposures, and their potential interaction effects. Less is known about the role of paternal genetic variants or environmental exposures and risk of OHDs. We examined parent-of-origin effects in transmission of alleles in the folate, homocysteine, or transsulfuration pathway genes on OHD occurrence in offspring. We used data on 569 families of liveborn infants with OHDs born between October 1997 and August 2008 from the National Birth Defects Prevention Study to conduct a family-based case-only study. Maternal, paternal, and infant DNA were genotyped using an Illumina Golden Gate custom single nucleotide polymorphism (SNP) panel. Relative risks (RR), 95% confidence interval (CI), and likelihood ratio tests from log-linear models were used to estimate the parent-of-origin effect of 877 SNPs in 60 candidate genes in the folate, homocysteine, and transsulfuration pathways on the risk of OHDs. Bonferroni correction was applied for multiple testing. We identified 3 SNPs in the transsulfuration pathway and 1 SNP in the folate pathway that were statistically significant after Bonferroni correction. Among infants who inherited paternally-derived copies of the G allele for rs6812588 in the RFC1 gene, the G allele for rs1762430 in the MGMT gene, and the A allele for rs9296695 and rs4712023 in the GSTA3 gene, RRs for OHD were 0.11 (95% CI: 0.04, 0.29, P = 9.16x10^-7), 0.30 (95% CI: 0.17, 0.53, P = 9.80x10^-6), 0.34 (95% CI: 0.20, 0.57, P = 2.28x10^-5), and 0.34 (95% CI: 0.20, 0.58, P = 3.77x10^-5), respectively, compared to infants who inherited maternally-derived copies of the same alleles. We observed statistically significant decreased risk of OHDs among infants who inherited paternal gene variants involved in folate and transsulfuration pathways.