Genome-wide association study identifies two loci strongly affecting transferrin glycosylation.

Polysaccharide sidechains attached to proteins play important roles in cell-cell and receptor-ligand interactions. Variation in the carbohydrate component has been extensively studied for the iron transport protein transferrin, because serum levels of the transferrin isoforms asialotransferrin + disialotransferrin (carbohydrate-deficient transferrin, CDT) are used as biomarkers of excessive alcohol intake. We conducted a genome-wide association study to assess whether genetic factors affect CDT concentration in serum. CDT was measured in three population-based studies: one in Switzerland (CoLaus study, n = 5181) and two in Australia (n = 1509, n = 775). The first cohort was used as the discovery panel and the latter ones served as replication. Genome-wide single-nucleotide polymorphism (SNP) typing data were used to identify loci with significant associations with CDT as a percentage of total transferrin (CDT%). The top three SNPs in the discovery panel (rs2749097 near PGM1 on chromosome 1, and missense polymorphisms rs1049296, rs1799899 in TF on chromosome 3) were successfully replicated , yielding genome-wide significant combined association with CDT% (P = 1.9 × 10(-9), 4 × 10(-39), 5.5 × 10(-43), respectively) and explain 5.8% of the variation in CDT%. These allelic effects are postulated to be caused by variation in availability of glucose-1-phosphate as a precursor of the glycan (PGM1), and variation in transferrin (TF) structure.

A genomewide association study of smoking relapse in four European population-based samples.

OBJECTIVES: Genomewide association studies (GWAS) have identified clear evidence of genetic markers for nicotine dependence. Other smoking phenotypes have been tested, but the results are less consistent. The tendency to relapse versus the ability to maintain long-term abstinence has received little attention in genetic studies; thus, our aim was to provide a better biological understanding of this phenotype through the identification of genetic loci associated with smoking relapse. METHODS: We carried out a GWAS on data from two European population-based collections, including a total of 835 cases (relapsers) and 990 controls (abstainers). Top-ranked findings from the discovery phase were tested for replication in two additional independent European population-based cohorts. RESULTS: Of the seven top markers from the discovery phase, none were consistently associated with smoking relapse across all samples and none reached genomewide significance. A single-nucleotide polymorphism rs1008509, within the Xylosyltransferase II (XYLT2) gene, was suggestively associated with smoking relapse in the discovery phase (β=-0.504; P=5.6E-06) and in the first replication sample (ALSPAC) (β=-0.27; P=0.004; n=1932), but not in the second sample (KORA) (β=0.19; P=0.138; n=912). We failed to identify an association between loci implicated previously in other smoking phenotypes and smoking relapse. CONCLUSION: Although no genomewide significant findings emerged from this study, we found that loci implicated in other smoking phenotypes were not associated with smoking relapse, which suggests that the neurobiology of smoking relapse and long-term abstinence may be distinct from biological mechanisms implicated in the development of nicotine dependence.

Quantitative sequence analysis of FBN1 premature termination codons provides evidence for incomplete NMD in leukocytes.

We improved, evaluated, and used Sanger sequencing for quantification of single nucleotide polymorphism (SNP) variants in transcripts and gDNA samples. This improved assay resulted in highly reproducible relative allele frequencies (e.g., for a heterozygous gDNA 50.0+/-1.4%, and for a missense mutation-bearing transcript 46.9+/-3.7%) with a lower detection limit of 3-9%. It provided excellent accuracy and linear correlation between expected and observed relative allele frequencies. This sequencing assay, which can also be used for the quantification of copy number variations (CNVs), methylations, mosaicisms, and DNA pools, enabled us to analyze transcripts of the FBN1 gene in fibroblasts and blood samples of patients with suspected Marfan syndrome not only qualitatively but also quantitatively. We report a total of 18 novel and 19 known FBN1 sequence variants leading to a premature termination codon (PTC), 26 of which we analyzed by quantitative sequencing both at gDNA and cDNA levels. The relative amounts of PTC-containing FBN1 transcripts in fresh and PAXgene-stabilized blood samples were significantly higher (33.0+/-3.9% to 80.0+/-7.2%) than those detected in affected fibroblasts with inhibition of nonsense-mediated mRNA decay (NMD) (11.0+/-2.1% to 25.0+/-1.8%), whereas in fibroblasts without NMD inhibition no mutant alleles could be detected. These results provide evidence for incomplete NMD in leukocytes and have particular importance for RNA-based analyses not only in FBN1 but also in other genes.

Association Study of Common Genetic Variants and HIV-1 Acquisition in 6,300 Infected Cases and 7,200 Controls.

Multiple genome-wide association studies (GWAS) have been performed in HIV-1 infected individuals, identifying common genetic influences on viral control and disease course. Similarly, common genetic correlates of acquisition of HIV-1 after exposure have been interrogated using GWAS, although in generally small samples. Under the auspices of the International Collaboration for the Genomics of HIV, we have combined the genome-wide single nucleotide polymorphism (SNP) data collected by 25 cohorts, studies, or institutions on HIV-1 infected individuals and compared them to carefully matched population-level data sets (a list of all collaborators appears in Note S1 in Text S1). After imputation using the 1,000 Genomes Project reference panel, we tested approximately 8 million common DNA variants (SNPs and indels) for association with HIV-1 acquisition in 6,334 infected patients and 7,247 population samples of European ancestry. Initial association testing identified the SNP rs4418214, the C allele of which is known to tag the HLA-B*57:01 and B*27:05 alleles, as genome-wide significant (p = 3.6×10(-11)). However, restricting analysis to individuals with a known date of seroconversion suggested that this association was due to the frailty bias in studies of lethal diseases. Further analyses including testing recessive genetic models, testing for bulk effects of non-genome-wide significant variants, stratifying by sexual or parenteral transmission risk and testing previously reported associations showed no evidence for genetic influence on HIV-1 acquisition (with the exception of CCR5Δ32 homozygosity). Thus, these data suggest that genetic influences on HIV acquisition are either rare or have smaller effects than can be detected by this sample size.

Altered brain-derived neurotrophic factor blood levels and gene variability are associated with anorexia and bulimia

Murine models and association studies in eating disorder (ED) patients have shown a role for the brain-derived neurotrophic factor (BDNF) in eating behavior. Some studies have shown association of BDNF -270C/T single-nucleotide polymorphism (SNP) with bulimia nervosa (BN), while BDNF Val66Met variant has been shown to be associated with both BN and anorexia nervosa (AN). To further test the role of this neurotrophin in humans, we screened 36 SNPs in the BDNF gene and tested for their association with ED and plasma BDNF levels as a quantitative trait. We performed a family-based association study in 106 ED nuclear families and analyzed BDNF blood levels in 110 ED patients and in 50 sib pairs discordant for ED. The rs7124442T/rs11030102C/rs11030119G haplotype was found associated with high BDNF levels (mean BDNF TCG haplotype carriers = 43.6 ng/ml vs. mean others 23.0 ng/ml, P = 0.016) and BN (Z = 2.64; P recessive = 0.008), and the rs7934165A/270T haplotype was associated with AN (Z =-2.64; P additive = 0.008). The comparison of BDNF levels in 50 ED discordant sib pairs showed elevated plasma BDNF levels for the ED group (mean controls = 41.0 vs. mean ED = 52.7; P = 0.004). Our data strongly suggest that altered BDNF levels modulated by BDNF gene variability are associated with the susceptibility to ED, providing physiological evidence that BDNF plays a role in the development of AN and BN, and strongly arguing for its involvement in eating behavior and body weight regulation.

Decay of linkage disequilibrium within genes across HGDP-CEPH human samples: most population isolates do not show increased LD

Background: It is well known that the pattern of linkage disequilibrium varies between human populations, with remarkable geographical stratification. Indirect association studies routinely exploit linkage disequilibrium around genes, particularly in isolated populations where it is assumed to be higher. Here, we explore both the amount and the decay of linkage disequilibrium with physical distance along 211 gene regions, most of them related to complex diseases, across 39 HGDP-CEPH population samples, focusing particularly on the populations defined as isolates. Within each gene region and population we use r2 between all possible single nucleotide polymorphism (SNP) pairs as a measure of linkage disequilibrium and focus on the proportion of SNP pairs with r2 greater than 0.8.Results: Although the average r2 was found to be significantly different both between and within continental regions, a much higher proportion of r2 variance could be attributed to differences between continental regions (2.8% vs. 0.5%, respectively). Similarly, while the proportion of SNP pairs with r2 > 0.8 was significantly different across continents for all distance classes, it was generally much more homogenous within continents, except in the case of Africa and the Americas. The only isolated populations with consistently higher LD in all distance classes with respect to their continent are the Kalash (Central South Asia) and the Surui (America). Moreover, isolated populations showed only slightly higher proportions of SNP pairs with r2 > 0.8 per gene region than non-isolated populations in the same continent. Thus, the number of SNPs in isolated populations that need to be genotyped may be only slightly less than in non-isolates. Conclusion: The "isolated population" label by itself does not guarantee a greater genotyping efficiency in association studies, and properties other than increased linkage disequilibrium may make these populations interesting in genetic epidemiology.

Impact of genetic SLC28 transporter and ITPA variants on ribavirin serum level, hemoglobin drop and therapeutic response in patients with HCV infection.

BACKGROUND & AIMS: In the last decade, pegylated interferon-α (PegIFN-α) plus ribavirin (RBV) was the standard treatment of chronic hepatitis C for genotype 1, and it remains the standard for genotypes 2 and 3. Recent studies reported associations between RBV-induced anemia and genetic polymorphisms of concentrative nucleoside transporters such as CNT3 (encoded by SLC28A3) and inosine triphosphatase (encoded by ITPA). We aimed at studying genetic determinants of RBV kinetics, efficacy and treatment-associated anemia. METHODS: We included 216 patients from two Swiss study cohorts (61% HCV genotype 1, 39% genotypes 2 or 3). Patients were analyzed for SLC28A2 single nucleotide polymorphism (SNP) rs11854484, SLC28A3 rs56350726, and SLC28A3 rs10868138 as well as ITPA SNPs rs1127354 and rs7270101, and followed for treatment-associated hemoglobin changes and sustained virological response (SVR). In 67 patients, RBV serum levels were additionally measured during treatment. RESULTS: Patients with SLC28A2 rs11854484 genotype TT had higher dosage- and body weight-adjusted RBV levels than those with genotypes TC or CC (p=0.02 and p=0.06 at weeks 4 and 8, respectively). ITPA SNP rs1127354 was associated with hemoglobin drop ≥3 g/dl during treatment, in genotype (relative risk (RR)=2.1, 95% CI 1.3-3.5) as well as allelic analyses (RR=2.0, 95%CI 1.2-3.4). SLC28A3 rs56350726 was associated with SVR in genotype (RR=2.2; 95% CI 1.1-4.3) as well as allelic analyses (RR=2.0, 95% CI 1.1-3.4). CONCLUSIONS: The newly identified association between RBV serum levels and SLC28A2 rs11854484 genotype, as well as the replicated association of ITPA and SLC28A3 genetic polymorphisms with RBV-induced anemia and treatment response, may support individualized treatment of chronic hepatitis C and warrant further investigation in larger studies.

Stratification by smoking status reveals an association of CHRNA5-A3-B4 genotype with body mass index in never smokers.

We previously used a single nucleotide polymorphism (SNP) in the CHRNA5-A3-B4 gene cluster associated with heaviness of smoking within smokers to confirm the causal effect of smoking in reducing body mass index (BMI) in a Mendelian randomisation analysis. While seeking to extend these findings in a larger sample we found that this SNP is associated with 0.74% lower body mass index (BMI) per minor allele in current smokers (95% CI -0.97 to -0.51, P = 2.00 × 10(-10)), but also unexpectedly found that it was associated with 0.35% higher BMI in never smokers (95% CI +0.18 to +0.52, P = 6.38 × 10(-5)). An interaction test confirmed that these estimates differed from each other (P = 4.95 × 10(-13)). This difference in effects suggests the variant influences BMI both via pathways unrelated to smoking, and via the weight-reducing effects of smoking. It would therefore be essentially undetectable in an unstratified genome-wide association study of BMI, given the opposite association with BMI in never and current smokers. This demonstrates that novel associations may be obscured by hidden population sub-structure. Stratification on well-characterized environmental factors known to impact on health outcomes may therefore reveal novel genetic associations.

Genetic Loci Associated with C-reactive Protein Levels and Risk of Coronary Heart Disease.

CONTEXT: Plasma levels of C-reactive protein (CRP) are independently associated with risk of coronary heart disease, but whether CRP is causally associated with coronary heart disease or merely a marker of underlying atherosclerosis is uncertain. OBJECTIVE: To investigate association of genetic loci with CRP levels and risk of coronary heart disease. DESIGN, SETTING, AND PARTICIPANTS: We first carried out a genome-wide association (n = 17,967) and replication study (n = 13,615) to identify genetic loci associated with plasma CRP concentrations. Data collection took place between 1989 and 2008 and genotyping between 2003 and 2008. We carried out a mendelian randomization study of the most closely associated single-nucleotide polymorphism (SNP) in the CRP locus and published data on other CRP variants involving a total of 28,112 cases and 100,823 controls, to investigate the association of CRP variants with coronary heart disease. We compared our finding with that predicted from meta-analysis of observational studies of CRP levels and risk of coronary heart disease. For the other loci associated with CRP levels, we selected the most closely associated SNP for testing against coronary heart disease among 14,365 cases and 32,069 controls. MAIN OUTCOME MEASURE: Risk of coronary heart disease. RESULTS: Polymorphisms in 5 genetic loci were strongly associated with CRP levels (% difference per minor allele): SNP rs6700896 in LEPR (-14.8%; 95% confidence interval [CI], -17.6% to -12.0%; P = 6.2 x 10(-22)), rs4537545 in IL6R (-11.5%; 95% CI, -14.4% to -8.5%; P = 1.3 x 10(-12)), rs7553007 in the CRP locus (-20.7%; 95% CI, -23.4% to -17.9%; P = 1.3 x 10(-38)), rs1183910 in HNF1A (-13.8%; 95% CI, -16.6% to -10.9%; P = 1.9 x 10(-18)), and rs4420638 in APOE-CI-CII (-21.8%; 95% CI, -25.3% to -18.1%; P = 8.1 x 10(-26)). Association of SNP rs7553007 in the CRP locus with coronary heart disease gave an odds ratio (OR) of 0.98 (95% CI, 0.94 to 1.01) per 20% lower CRP level. Our mendelian randomization study of variants in the CRP locus showed no association with coronary heart disease: OR, 1.00; 95% CI, 0.97 to 1.02; per 20% lower CRP level, compared with OR, 0.94; 95% CI, 0.94 to 0.95; predicted from meta-analysis of the observational studies of CRP levels and coronary heart disease (z score, -3.45; P < .001). SNPs rs6700896 in LEPR (OR, 1.06; 95% CI, 1.02 to 1.09; per minor allele), rs4537545 in IL6R (OR, 0.94; 95% CI, 0.91 to 0.97), and rs4420638 in the APOE-CI-CII cluster (OR, 1.16; 95% CI, 1.12 to 1.21) were all associated with risk of coronary heart disease. CONCLUSION: The lack of concordance between the effect on coronary heart disease risk of CRP genotypes and CRP levels argues against a causal association of CRP with coronary heart disease.

ZNRD1 (zinc ribbon domain-containing 1) is a host cellular factor that influences HIV-1 replication and disease progression.

BACKGROUND: Human immunodeficiency virus (HIV) takes advantage of multiple host proteins to support its own replication. The gene ZNRD1 (zinc ribbon domain-containing 1) has been identified as encoding a potential host factor that influenced disease progression in HIV-positive individuals in a genomewide association study and also significantly affected HIV replication in a large-scale in vitro short interfering RNA (siRNA) screen. Genes and polymorphisms identified by large-scale analysis need to be followed up by means of functional assays and resequencing efforts to more precisely map causal genes. METHODS: Genotyping and ZNRD1 gene resequencing for 208 HIV-positive subjects (119 who experienced long-term nonprogression [LTNP] and 89 who experienced normal disease progression) was done by either TaqMan genotyping assays or direct sequencing. Genetic association analysis was performed with the SNPassoc package and Haploview software. siRNA and short hairpin RNA (shRNA) specifically targeting ZNRD1 were used to transiently or stably down-regulate ZNRD1 expression in both lymphoid and nonlymphoid cells. Cells were infected with X4 and R5 HIV strains, and efficiency of infection was assessed by reporter gene assay or p24 assay. RESULTS: Genetic association analysis found a strong statistically significant correlation with the LTNP phenotype (single-nucleotide polymorphism rs1048412; [Formula: see text]), independently of HLA-A10 influence. siRNA-based functional analysis showed that ZNRD1 down-regulation by siRNA or shRNA impaired HIV-1 replication at the transcription level in both lymphoid and nonlymphoid cells. CONCLUSION: Genetic association analysis unequivocally identified ZNRD1 as an independent marker of LTNP to AIDS. Moreover, in vitro experiments pointed to viral transcription as the inhibited step. Thus, our data strongly suggest that ZNRD1 is a host cellular factor that influences HIV-1 replication and disease progression in HIV-positive individuals.

Estudio Farmacogenético en Pacientes con Cáncer de Pulmón No Microcítico Avanzado tratados con Quimioterapia basada en Platino

a partir de ADN genómico obtenido de las células nucleadas de sangre periférica de 103 pacientes con Cáncer de Pulmón No Microcítico (CPNM) avanzado tratados con quimioterapia basada en platino, hemos analizado la asociación entre supervivencia y cinco SNPs (Single Nucleotide Polymorphism) pertenecientes a dos grupos de genes: i) de la via metabólica del ácido fólico (Timidilato Sintetasa (TS), Metil-tetrahidrofolato Reductasa (MTHFR) y, ii) de la vía de reparación del ADN (Excision repair cross-complemeting group 1 (ERCC1) y Xeroderma pigmentosum group D (XPD).

Global Genomic Diversity of Human Papillomavirus 6 Based on 724 Isolates and 190 Complete Genome Sequences.

Human papillomavirus type 6 (HPV6) is the major etiological agent of anogenital warts and laryngeal papillomas and has been included in both the quadrivalent and nonavalent prophylactic HPV vaccines. This study investigated the global genomic diversity of HPV6, using 724 isolates and 190 complete genomes from six continents, and the association of HPV6 genomic variants with geographical location, anatomical site of infection/disease, and gender. Initially, a 2,800-bp E5a-E5b-L1-LCR fragment was sequenced from 492/530 (92.8%) HPV6-positive samples collected for this study. Among them, 130 exhibited at least one single nucleotide polymorphism (SNP), indel, or amino acid change in the E5a-E5b-L1-LCR fragment and were sequenced in full. A global alignment and maximum likelihood tree of 190 complete HPV6 genomes (130 fully sequenced in this study and 60 obtained from sequence repositories) revealed two variant lineages, A and B, and five B sublineages: B1, B2, B3, B4, and B5. HPV6 (sub)lineage-specific SNPs and a 960-bp representative region for whole-genome-based phylogenetic clustering within the L2 open reading frame were identified. Multivariate logistic regression analysis revealed that lineage B predominated globally. Sublineage B3 was more common in Africa and North and South America, and lineage A was more common in Asia. Sublineages B1 and B3 were associated with anogenital infections, indicating a potential lesion-specific predilection of some HPV6 sublineages. Females had higher odds for infection with sublineage B3 than males. In conclusion, a global HPV6 phylogenetic analysis revealed the existence of two variant lineages and five sublineages, showing some degree of ethnogeographic, gender, and/or disease predilection in their distribution. IMPORTANCE: This study established the largest database of globally circulating HPV6 genomic variants and contributed a total of 130 new, complete HPV6 genome sequences to available sequence repositories. Two HPV6 variant lineages and five sublineages were identified and showed some degree of association with geographical location, anatomical site of infection/disease, and/or gender. We additionally identified several HPV6 lineage- and sublineage-specific SNPs to facilitate the identification of HPV6 variants and determined a representative region within the L2 gene that is suitable for HPV6 whole-genome-based phylogenetic analysis. This study complements and significantly expands the current knowledge of HPV6 genetic diversity and forms a comprehensive basis for future epidemiological, evolutionary, functional, pathogenicity, vaccination, and molecular assay development studies.

Integration of genome-wide association studies with biological knowledge identifies six novel genes related to kidney function.

In conducting genome-wide association studies (GWAS), analytical approaches leveraging biological information may further understanding of the pathophysiology of clinical traits. To discover novel associations with estimated glomerular filtration rate (eGFR), a measure of kidney function, we developed a strategy for integrating prior biological knowledge into the existing GWAS data for eGFR from the CKDGen Consortium. Our strategy focuses on single nucleotide polymorphism (SNPs) in genes that are connected by functional evidence, determined by literature mining and gene ontology (GO) hierarchies, to genes near previously validated eGFR associations. It then requires association thresholds consistent with multiple testing, and finally evaluates novel candidates by independent replication. Among the samples of European ancestry, we identified a genome-wide significant SNP in FBXL20 (P = 5.6 × 10(-9)) in meta-analysis of all available data, and additional SNPs at the INHBC, LRP2, PLEKHA1, SLC3A2 and SLC7A6 genes meeting multiple-testing corrected significance for replication and overall P-values of 4.5 × 10(-4)-2.2 × 10(-7). Neither the novel PLEKHA1 nor FBXL20 associations, both further supported by association with eGFR among African Americans and with transcript abundance, would have been implicated by eGFR candidate gene approaches. LRP2, encoding the megalin receptor, was identified through connection with the previously known eGFR gene DAB2 and extends understanding of the megalin system in kidney function. These findings highlight integration of existing genome-wide association data with independent biological knowledge to uncover novel candidate eGFR associations, including candidates lacking known connections to kidney-specific pathways. The strategy may also be applicable to other clinical phenotypes, although more testing will be needed to assess its potential for discovery in general.

Genome-wide arrays in routine diagnostics of hematological malignancies.

Over the last three decades, cytogenetic analysis of malignancies has become an integral part of disease evaluation and prediction of prognosis or responsiveness to therapy. In most diagnostic laboratories, conventional karyotyping, in conjunction with targeted fluorescence in situ hybridization analysis, is routinely performed to detect recurrent aberrations with prognostic implications. However, the genetic complexity of cancer cells requires a sensitive genome-wide analysis, enabling the detection of small genomic changes in a mixed cell population, as well as of regions of homozygosity. The advent of comprehensive high-resolution genomic tools, such as molecular karyotyping using comparative genomic hybridization or single-nucleotide polymorphism microarrays, has overcome many of the limitations of traditional cytogenetic techniques and has been used to study complex genomic lesions in, for example, leukemia. The clinical impact of the genomic copy-number and copy-neutral alterations identified by microarray technologies is growing rapidly and genome-wide array analysis is evolving into a diagnostic tool, to better identify high-risk patients and predict patients' outcomes from their genomic profiles. Here, we review the added clinical value of an array-based genome-wide screen in leukemia, and discuss the technical challenges and an interpretation workflow in applying arrays in the acquired cytogenetic diagnostic setting.

In vitro whole-genome analysis identifies a susceptibility locus for HIV-1.

Advances in large-scale analysis of human genomic variability provide unprecedented opportunities to study the genetic basis of susceptibility to infectious agents. We report here the use of an in vitro system for the identification of a locus on HSA8q24.3 associated with cellular susceptibility to HIV-1. This locus was mapped through quantitative linkage analysis using cell lines from multigeneration families, validated in vitro, and followed up by two independent association studies in HIV-positive individuals. Single nucleotide polymorphism rs2572886, which is associated with cellular susceptibility to HIV-1 in lymphoblastoid B cells and in primary T cells, was also associated with accelerated disease progression in one of two cohorts of HIV-1-infected patients. Biological analysis suggests a role of the rs2572886 region in the regulation of the LY6 family of glycosyl-phosphatidyl-inositol (GPI)-anchored proteins. Genetic analysis of in vitro cellular phenotypes provides an attractive approach for the discovery of susceptibility loci to infectious agents.