A mega-analysis of genome-wide association studies for major depressive disorder.

Prior genome-wide association studies (GWAS) of major depressive disorder (MDD) have met with limited success. We sought to increase statistical power to detect disease loci by conducting a GWAS mega-analysis for MDD. In the MDD discovery phase, we analyzed more than 1.2 million autosomal and X chromosome single-nucleotide polymorphisms (SNPs) in 18 759 independent and unrelated subjects of recent European ancestry (9240 MDD cases and 9519 controls). In the MDD replication phase, we evaluated 554 SNPs in independent samples (6783 MDD cases and 50 695 controls). We also conducted a cross-disorder meta-analysis using 819 autosomal SNPs with P<0.0001 for either MDD or the Psychiatric GWAS Consortium bipolar disorder (BIP) mega-analysis (9238 MDD cases/8039 controls and 6998 BIP cases/7775 controls). No SNPs achieved genome-wide significance in the MDD discovery phase, the MDD replication phase or in pre-planned secondary analyses (by sex, recurrent MDD, recurrent early-onset MDD, age of onset, pre-pubertal onset MDD or typical-like MDD from a latent class analyses of the MDD criteria). In the MDD-bipolar cross-disorder analysis, 15 SNPs exceeded genome-wide significance (P<5 × 10(-8)), and all were in a 248 kb interval of high LD on 3p21.1 (chr3:52 425 083-53 822 102, minimum P=5.9 × 10(-9) at rs2535629). Although this is the largest genome-wide analysis of MDD yet conducted, its high prevalence means that the sample is still underpowered to detect genetic effects typical for complex traits. Therefore, we were unable to identify robust and replicable findings. We discuss what this means for genetic research for MDD. The 3p21.1 MDD-BIP finding should be interpreted with caution as the most significant SNP did not replicate in MDD samples, and genotyping in independent samples will be needed to resolve its status.

Association of adiposity genetic variants with menarche timing in 92,105 women of European descent.

Obesity is of global health concern. There are well-described inverse relationships between female pubertal timing and obesity. Recent genome-wide association studies of age at menarche identified several obesity-related variants. Using data from the ReproGen Consortium, we employed meta-analytical techniques to estimate the associations of 95 a priori and recently identified obesity-related (body mass index (weight (kg)/height (m)(2)), waist circumference, and waist:hip ratio) single-nucleotide polymorphisms (SNPs) with age at menarche in 92,116 women of European descent from 38 studies (1970-2010), in order to estimate associations between genetic variants associated with central or overall adiposity and pubertal timing in girls. Investigators in each study performed a separate analysis of associations between the selected SNPs and age at menarche (ages 9-17 years) using linear regression models and adjusting for birth year, site (as appropriate), and population stratification. Heterogeneity of effect-measure estimates was investigated using meta-regression. Six novel associations of body mass index loci with age at menarche were identified, and 11 adiposity loci previously reported to be associated with age at menarche were confirmed, but none of the central adiposity variants individually showed significant associations. These findings suggest complex genetic relationships between menarche and overall obesity, and to a lesser extent central obesity, in normal processes of growth and development.

A genome-wide association study of resistance to HIV infection in highly exposed uninfected individuals with hemophilia A.

Human genetic variation contributes to differences in susceptibility to HIV-1 infection. To search for novel host resistance factors, we performed a genome-wide association study (GWAS) in hemophilia patients highly exposed to potentially contaminated factor VIII infusions. Individuals with hemophilia A and a documented history of factor VIII infusions before the introduction of viral inactivation procedures (1979-1984) were recruited from 36 hemophilia treatment centers (HTCs), and their genome-wide genetic variants were compared with those from matched HIV-infected individuals. Homozygous carriers of known CCR5 resistance mutations were excluded. Single nucleotide polymorphisms (SNPs) and inferred copy number variants (CNVs) were tested using logistic regression. In addition, we performed a pathway enrichment analysis, a heritability analysis, and a search for epistatic interactions with CCR5 Δ32 heterozygosity. A total of 560 HIV-uninfected cases were recruited: 36 (6.4%) were homozygous for CCR5 Δ32 or m303. After quality control and SNP imputation, we tested 1 081 435 SNPs and 3686 CNVs for association with HIV-1 serostatus in 431 cases and 765 HIV-infected controls. No SNP or CNV reached genome-wide significance. The additional analyses did not reveal any strong genetic effect. Highly exposed, yet uninfected hemophiliacs form an ideal study group to investigate host resistance factors. Using a genome-wide approach, we did not detect any significant associations between SNPs and HIV-1 susceptibility, indicating that common genetic variants of major effect are unlikely to explain the observed resistance phenotype in this population.

Genetic Analyses Reveal a Role for Vitamin D Insufficiency in HCV-Associated Hepatocellular Carcinoma Development.

BACKGROUND: Vitamin D insufficiency has been associated with the occurrence of various types of cancer, but causal relationships remain elusive. We therefore aimed to determine the relationship between genetic determinants of vitamin D serum levels and the risk of developing hepatitis C virus (HCV)-related hepatocellular carcinoma (HCC). METHODOLOGYPRINCIPAL FINDINGS: Associations between CYP2R1, GC, and DHCR7 genotypes that are determinants of reduced 25-hydroxyvitamin D (25[OH]D3) serum levels and the risk of HCV-related HCC development were investigated for 1279 chronic hepatitis C patients with HCC and 4325 without HCC, respectively. The well-known associations between CYP2R1 (rs1993116, rs10741657), GC (rs2282679), and DHCR7 (rs7944926, rs12785878) genotypes and 25(OH)D3 serum levels were also apparent in patients with chronic hepatitis C. The same genotypes of these single nucleotide polymorphisms (SNPs) that are associated with reduced 25(OH)D3 serum levels were found to be associated with HCV-related HCC (P = 0.07 [OR = 1.13, 95% CI = 0.99-1.28] for CYP2R1, P = 0.007 [OR = 1.56, 95% CI = 1.12-2.15] for GC, P = 0.003 [OR = 1.42, 95% CI = 1.13-1.78] for DHCR7; ORs for risk genotypes). In contrast, no association between these genetic variations and liver fibrosis progression rate (P>0.2 for each SNP) or outcome of standard therapy with pegylated interferon-α and ribavirin (P>0.2 for each SNP) was observed, suggesting a specific influence of the genetic determinants of 25(OH)D3 serum levels on hepatocarcinogenesis. CONCLUSIONSSIGNIFICANCE: Our data suggest a relatively weak but functionally relevant role for vitamin D in the prevention of HCV-related hepatocarcinogenesis.

Genome-wide meta-analysis of common variant differences between men and women.

The male-to-female sex ratio at birth is constant across world populations with an average of 1.06 (106 male to 100 female live births) for populations of European descent. The sex ratio is considered to be affected by numerous biological and environmental factors and to have a heritable component. The aim of this study was to investigate the presence of common allele modest effects at autosomal and chromosome X variants that could explain the observed sex ratio at birth. We conducted a large-scale genome-wide association scan (GWAS) meta-analysis across 51 studies, comprising overall 114 863 individuals (61 094 women and 53 769 men) of European ancestry and 2 623 828 common (minor allele frequency >0.05) single-nucleotide polymorphisms (SNPs). Allele frequencies were compared between men and women for directly-typed and imputed variants within each study. Forward-time simulations for unlinked, neutral, autosomal, common loci were performed under the demographic model for European populations with a fixed sex ratio and a random mating scheme to assess the probability of detecting significant allele frequency differences. We do not detect any genome-wide significant (P < 5 × 10(-8)) common SNP differences between men and women in this well-powered meta-analysis. The simulated data provided results entirely consistent with these findings. This large-scale investigation across ~115 000 individuals shows no detectable contribution from common genetic variants to the observed skew in the sex ratio. The absence of sex-specific differences is useful in guiding genetic association study design, for example when using mixed controls for sex-biased traits.

Variant within the promoter region of the CHRNA3 gene associated with FTN dependence is not related to self-reported willingness to quit smoking.

INTRODUCTION: Common variation in the CHRNA5-CHRNA3-CHRNB4 gene region is robustly associated with smoking quantity. Conversely, the association between one of the most significant single nucleotide polymorphisms (SNPs; rs1051730 within the CHRNA3 gene) with perceived difficulty or willingness to quit smoking among current smokers is unknown. METHODS: Cross-sectional study including current smokers, 502 women, and 552 men. Heaviness of smoking index (HSI), difficulty, attempting, and intention to quit smoking were assessed by questionnaire. RESULTS: The rs1051730 SNP was associated with increased HSI (age, gender, and education-adjusted mean ± SE: 2.6 ± 0.1, 2.2 ± 0.1, and 2.0 ± 0.1 for AA, AG, and GG genotypes, respectively, p < .01). Multivariate logistic regression adjusting for gender, age, education, leisure-time physical activity, and personal history of cardiovascular or lung disease showed rs1051730 to be associated with higher smoking dependence (odds ratio [OR] and 95% CI for each additional A-allele: 1.38 [1.11-1.72] for smoking more than 20 cigarette equivalents/day; 1.31 [1.00-1.71] for an HSI ≥5 and 1.32 [1.05-1.65] for smoking 5 min after waking up) and borderline associated with difficulty to quit (OR = 1.29 [0.98-1.70]), but this relationship was no longer significant after adjusting for nicotine dependence. Also, no relationship was found with willingness (OR = 1.03 [0.85-1.26]), attempt (OR = 1.00 [0.83-1.20]), or preparation (OR = 0.95 [0.38-2.38]) to quit. Similar findings were obtained for other SNPs, but their effect on nicotine dependence was no longer significant after adjusting for rs1051730. Conclusions: These data confirm the effect of rs1051730 on nicotine dependence but failed to find any relationship with difficulty, willingness, and motivation to quit.

Multi-ethnic fine-mapping of 14 central adiposity loci.

The Genetic Investigation of Anthropometric Traits (GIANT) consortium identified 14 loci in European Ancestry (EA) individuals associated with waist-to-hip ratio (WHR) adjusted for body mass index. These loci are wide and narrowing the signals remains necessary. Twelve of 14 loci identified in GIANT EA samples retained strong associations with WHR in our joint EA/individuals of African Ancestry (AA) analysis (log-Bayes factor >6.1). Trans-ethnic analyses at five loci (TBX15-WARS2, LYPLAL1, ADAMTS9, LY86 and ITPR2-SSPN) substantially narrowed the signals to smaller sets of variants, some of which are in regions that have evidence of regulatory activity. By leveraging varying linkage disequilibrium structures across different populations, single-nucleotide polymorphisms (SNPs) with strong signals and narrower credible sets from trans-ethnic meta-analysis of central obesity provide more precise localizations of potential functional variants and suggest a possible regulatory role. Meta-analysis results for WHR were obtained from 77 167 EA participants from GIANT and 23 564 AA participants from the African Ancestry Anthropometry Genetics Consortium. For fine mapping we interrogated SNPs within ± 250 kb flanking regions of 14 previously reported index SNPs from loci discovered in EA populations by performing trans-ethnic meta-analysis of results from the EA and AA meta-analyses. We applied a Bayesian approach that leverages allelic heterogeneity across populations to combine meta-analysis results and aids in fine-mapping shared variants at these locations. We annotated variants using information from the ENCODE Consortium and Roadmap Epigenomics Project to prioritize variants for possible functionality.

Genome-wide re-sequencing of multidrug-resistant Mycobacterium leprae Airaku-3.

Genotyping and molecular characterization of drug resistance mechanisms in Mycobacterium leprae enables disease transmission and drug resistance trends to be monitored. In the present study, we performed genome-wide analysis of Airaku-3, a multidrug-resistant strain with an unknown mechanism of resistance to rifampicin. We identified 12 unique non-synonymous single-nucleotide polymorphisms (SNPs) including two in the transporter-encoding ctpC and ctpI genes. In addition, two SNPs were found that improve the resolution of SNP-based genotyping, particularly for Venezuelan and South East Asian strains of M. leprae.

TRAIL/TRAIL Receptor System and Susceptibility to Multiple Sclerosis

The TNF-related apoptosis inducing ligand (TRAIL)/TRAIL receptor system participates in crucial steps in immune cell activation or differentiation. It is able to inhibit proliferation and activation of T cells and to induce apoptosis of neurons and oligodendrocytes, and seems to be implicated in autoimmune diseases. Thus, TRAIL and TRAIL receptor genes are potential candidates for involvement in susceptibility to multiple sclerosis (MS). To test whether single-nucleotide polymorphisms (SNPs) in the human genes encoding TRAIL, TRAILR-1, TRAILR-2, TRAILR-3 and TRAILR-4 are associated with MS susceptibility, we performed a candidate gene case-control study in the Spanish population. 59 SNPs in the TRAIL and TRAIL receptor genes were analysed in 628 MS patients and 660 controls, and validated in an additional cohort of 295 MS patients and 233 controls. Despite none of the SNPs withstood the highly conservative Bonferroni correction, three SNPs showing uncorrected p values<0.05 were successfully replicated: rs4894559 in TRAIL gene, p = 9.8×10(-4), OR = 1.34; rs4872077, in TRAILR-1 gene, p = 0.005, OR = 1.72; and rs1001793 in TRAILR-2 gene, p = 0.012, OR = 0.84. The combination of the alleles G/T/A in these SNPs appears to be associated with a reduced risk of developing MS (p = 2.12×10(-5), OR = 0.59). These results suggest that genes of the TRAIL/TRAIL receptor system exerts a genetic influence on MS.

A genome-wide association study reveals variants in ARL15 that influence adiponectin levels.

The adipocyte-derived protein adiponectin is highly heritable and inversely associated with risk of type 2 diabetes mellitus (T2D) and coronary heart disease (CHD). We meta-analyzed 3 genome-wide association studies for circulating adiponectin levels (n = 8,531) and sought validation of the lead single nucleotide polymorphisms (SNPs) in 5 additional cohorts (n = 6,202). Five SNPs were genome-wide significant in their relationship with adiponectin (P< or =5x10(-8)). We then tested whether these 5 SNPs were associated with risk of T2D and CHD using a Bonferroni-corrected threshold of P< or =0.011 to declare statistical significance for these disease associations. SNPs at the adiponectin-encoding ADIPOQ locus demonstrated the strongest associations with adiponectin levels (P-combined = 9.2x10(-19) for lead SNP, rs266717, n = 14,733). A novel variant in the ARL15 (ADP-ribosylation factor-like 15) gene was associated with lower circulating levels of adiponectin (rs4311394-G, P-combined = 2.9x10(-8), n = 14,733). This same risk allele at ARL15 was also associated with a higher risk of CHD (odds ratio [OR] = 1.12, P = 8.5x10(-6), n = 22,421) more nominally, an increased risk of T2D (OR = 1.11, P = 3.2x10(-3), n = 10,128), and several metabolic traits. Expression studies in humans indicated that ARL15 is well-expressed in skeletal muscle. These findings identify a novel protein, ARL15, which influences circulating adiponectin levels and may impact upon CHD risk.

Analysis of TNFAIP3, a feedback inhibitor of nuclear factor-κB and the neighbor intergenic 6q23 region in rheumatoid arthritis susceptibility

INTRODUCTION Genome-wide association studies of rheumatoid arthritis (RA) have identified an association of the disease with a 6q23 region devoid of genes. TNFAIP3, an RA candidate gene, flanks this region, and polymorphisms in both the TNFAIP3 gene and the intergenic region are associated with systemic lupus erythematosus. We hypothesized that there is a similar association with RA, including polymorphisms in TNFAIP3 and the intergenic region. METHODS To test this hypothesis, we selected tag-single nucleotide polymorphisms (SNPs) in both loci. They were analyzed in 1,651 patients with RA and 1,619 control individuals of Spanish ancestry. RESULTS Weak evidence of association was found both in the 6q23 intergenic region and in the TNFAIP3 locus. The rs582757 SNP and a common haplotype in the TNFAIP3 locus exhibited association with RA. In the intergenic region, two SNPs were associated, namely rs609438 and rs13207033. The latter was only associated in patients with anti-citrullinated peptide antibodies. Overall, statistical association was best explained by the interdependent contribution of SNPs from the two loci TNFAIP3 and the 6q23 intergenic region. CONCLUSIONS Our data are consistent with the hypothesis that several RA genetic factors exist in the 6q23 region, including polymorphisms in the TNFAIP3 gene, like that previously described for systemic lupus erythematosus.

Genome-wide linkage and association analyses to identify genes influencing adiponectin levels: the GEMS Study.

Adiponectin has a variety of metabolic effects on obesity, insulin sensitivity, and atherosclerosis. To identify genes influencing variation in plasma adiponectin levels, we performed genome-wide linkage and association scans of adiponectin in two cohorts of subjects recruited in the Genetic Epidemiology of Metabolic Syndrome Study. The genome-wide linkage scan was conducted in families of Turkish and southern European (TSE, n = 789) and Northern and Western European (NWE, N = 2,280) origin. A whole genome association (WGA) analysis (500K Affymetrix platform) was carried out in a set of unrelated NWE subjects consisting of approximately 1,000 subjects with dyslipidemia and 1,000 overweight subjects with normal lipids. Peak evidence for linkage occurred at chromosome 8p23 in NWE subjects (lod = 3.10) and at chromosome 3q28 near ADIPOQ, the adiponectin structural gene, in TSE subjects (lod = 1.70). In the WGA analysis, the single-nucleotide polymorphisms (SNPs) most strongly associated with adiponectin were rs3774261 and rs6773957 (P < 10(-7)). These two SNPs were in high linkage disequilibrium (r(2) = 0.98) and located within ADIPOQ. Interestingly, our fourth strongest region of association (P < 2 x 10(-5)) was to an SNP within CDH13, whose protein product is a newly identified receptor for high-molecular-weight species of adiponectin. Through WGA analysis, we confirmed previous studies showing SNPs within ADIPOQ to be strongly associated with variation in adiponectin levels and further observed these to have the strongest effects on adiponectin levels throughout the genome. We additionally identified a second gene (CDH13) possibly influencing variation in adiponectin levels. The impact of these SNPs on health and disease has yet to be determined.

Is the RET proto-oncogene involved in the pathogenesis of intestinal neuronal dysplasia type B?

Hirschsprung disease (HSCR) is defined by the absence of intramural ganglia of Meissner and Auerbach along variable lengths of the gastrointestinal tract. Intestinal neuronal dysplasia (IND) type B is characterized by the malformation of the parasympathetic submucous plexus of the gut. A connection appears to exist between these two enteric nervous system abnormalities. Due to the major role played by the RET proto-oncogene in HSCR, we sought to determine whether this gene was also related to INDB. dHPLC techniques were employed to screen the RET coding region in 23 patients presenting with INDB and 30 patients with a combined HSCR+INDB phenotype. In addition, eight RET single nucleotide polymorphisms (SNPs) were strategically selected and genotyped by TaqMan technology. The distribution of SNPs and haplotypes was compared among the different groups of patients (INDB, HSCR+INDB, HSCR) and the controls. We found several RET mutations in our patients and some differences in the distribution of the RET SNPs among the groups of study. Our results suggest an involvement of RET in the pathogenesis of intestinal INDB, although by different molecular mechanisms than those leading to HSCR. Further investigation is warranted to elucidate these precise mechanisms and to clarify the genetic nature of INDB.

Sequence analysis of coding DNA fragments of pfcrt and pfmdr-1 genes in Plasmodium falciparum isolates from Odisha, India

The global emergence and spread of malaria parasites resistant to antimalarial drugs is the major problem in malaria control. The genetic basis of the parasite's resistance to the antimalarial drug chloroquine (CQ) is well-documented, allowing for the analysis of field isolates of malaria parasites to address evolutionary questions concerning the origin and spread of CQ-resistance. Here, we present DNA sequence analyses of both the second exon of the Plasmodium falciparum CQ-resistance transporter (pfcrt) gene and the 5' end of the P. falciparum multidrug-resistance 1 (pfmdr-1) gene in 40 P. falciparum field isolates collected from eight different localities of Odisha, India. First, we genotyped the samples for the pfcrt K76T and pfmdr-1 N86Y mutations in these two genes, which are the mutations primarily implicated in CQ-resistance. We further analyzed amino acid changes in codons 72-76 of the pfcrt haplotypes. Interestingly, both the K76T and N86Y mutations were found to co-exist in 32 out of the total 40 isolates, which were of either the CVIET or SVMNT haplotype, while the remaining eight isolates were of the CVMNK haplotype. In total, eight nonsynonymous single nucleotide polymorphisms (SNPs) were observed, six in the pfcrt gene and two in the pfmdr-1 gene. One poorly studied SNP in the pfcrt gene (A97T) was found at a high frequency in many P. falciparum samples. Using population genetics to analyze these two gene fragments, we revealed comparatively higher nucleotide diversity in the pfcrt gene than in the pfmdr-1 gene. Furthermore, linkage disequilibrium was found to be tight between closely spaced SNPs of the pfcrt gene. Finally, both the pfcrt and the pfmdr-1 genes were found to evolve under the standard neutral model of molecular evolution.

Polymorphism analysis of the CTLA-4 gene in paracoccidioidomycosis patients

The CTLA-4 protein is expressed in activated T cells and plays an essential role in the immune response through its regulatory effect on T cell activation. Polymorphisms of the CTLA-4 gene have been correlated with autoimmune, neoplastic and infectious illnesses. This work aimed to verify possible associations between single nucleotide polymorphisms (SNPs) in CTLA-4, -318C/T in the promoter and +49A/G in exon 1 and paracoccidioidomycosis (PCM) caused by Paracoccidioides brasiliensis. For this purpose, 66 chronic form PCM patients and 76 healthy controls had their allele, genotype and haplotype frequencies determined. The genetic admixture structure of the patients and controls was evaluated to eliminate ancestral bias. The comparison of frequencies indicated no significant differences between patients and controls that could link the SNPs to PCM. Groups were admixture matched with no difference observed in population ancestry inference, indicating that the absence of association between CTLA-4 polymorphisms and PCM could not be attributed to ancestral bias. This study showed that there was no association between the CTLA-4 SNPs -318 and +49 and the resistance or susceptibility to PCM.