Blood pressure loci identified with a gene-centric array.

Raised blood pressure (BP) is a major risk factor for cardiovascular disease. Previous studies have identified 47 distinct genetic variants robustly associated with BP, but collectively these explain only a few percent of the heritability for BP phenotypes. To find additional BP loci, we used a bespoke gene-centric array to genotype an independent discovery sample of 25,118 individuals that combined hypertensive case-control and general population samples. We followed up four SNPs associated with BP at our p < 8.56 × 10(-7) study-specific significance threshold and six suggestively associated SNPs in a further 59,349 individuals. We identified and replicated a SNP at LSP1/TNNT3, a SNP at MTHFR-NPPB independent (r(2) = 0.33) of previous reports, and replicated SNPs at AGT and ATP2B1 reported previously. An analysis of combined discovery and follow-up data identified SNPs significantly associated with BP at p < 8.56 × 10(-7) at four further loci (NPR3, HFE, NOS3, and SOX6). The high number of discoveries made with modest genotyping effort can be attributed to using a large-scale yet targeted genotyping array and to the development of a weighting scheme that maximized power when meta-analyzing results from samples ascertained with extreme phenotypes, in combination with results from nonascertained or population samples. Chromatin immunoprecipitation and transcript expression data highlight potential gene regulatory mechanisms at the MTHFR and NOS3 loci. These results provide candidates for further study to help dissect mechanisms affecting BP and highlight the utility of studying SNPs and samples that are independent of those studied previously even when the sample size is smaller than that in previous studies.

Intronic variants in the NFKB1 gene may influence hearing forecast in patients with unilateral sensorineural hearing loss in Meniere's disease.

Meniere's disease is an episodic vestibular syndrome associated with sensorineural hearing loss (SNHL) and tinnitus. Patients with MD have an elevated prevalence of several autoimmune diseases (rheumatoid arthritis, systemic lupus erythematosus, ankylosing spondylitis and psoriasis), which suggests a shared autoimmune background. Functional variants of several genes involved in the NF-κB pathway, such as REL, TNFAIP3, NFKB1 and TNIP1, have been associated with two or more immune-mediated diseases and allelic variations in the TLR10 gene may influence bilateral affectation and clinical course in MD. We have genotyped 716 cases of MD and 1628 controls by using the ImmunoChip, a high-density genotyping array containing 186 autoimmune loci, to explore the association of immune system related-loci with sporadic MD. Although no single nucleotide polymorphism (SNP) reached a genome-wide significant association (p<10(-8)), we selected allelic variants in the NF-kB pathway for further analyses to evaluate the impact of these SNPs in the clinical outcome of MD in our cohort. None of the selected SNPs increased susceptibility for MD in patients with uni or bilateral SNHL. However, two potential regulatory variants in the NFKB1 gene (rs3774937 and rs4648011) were associated with a faster hearing loss progression in patients with unilateral SNHL. So, individuals with unilateral MD carrying the C allele in rs3774937 or G allele in rs4648011 had a shorter mean time to reach hearing stage 3 (>40 dB HL) (log-rank test, corrected p values were p = 0.009 for rs3774937 and p = 0.003 for rs4648011, respectively). No variants influenced hearing in bilateral MD. Our data support that the allelic variants rs3774937 and rs4648011 can modify hearing outcome in patients with MD and unilateral SNHL.

First insights into the transcriptome and development of new genomic tools of a widespread circum-Mediterranean tree species, Pinus halepensis Mill.

Aleppo pine (Pinus halepensis Mill.) is a relevant conifer species for studying adaptive responses to drought and fire regimes in the Mediterranean region. In this study, we performed Illumina next-generation sequencing of two phenotypically divergent Aleppo pine accessions with the aims of (i) characterizing the transcriptome through Illumina RNA-Seq on trees phenotypically divergent for adaptive traits linked to fire adaptation and drought, (ii) performing a functional annotation of the assembled transcriptome, (iii) identifying genes with accelerated evolutionary rates, (iv) studying the expression levels of the annotated genes and (v) developing gene-based markers for population genomic and association genetic studies. The assembled transcriptome consisted of 48,629 contigs and covered about 54.6 Mbp. The comparison of Aleppo pine transcripts to Picea sitchensis protein-coding sequences resulted in the detection of 34,014 SNPs across species, with a Ka /Ks average value of 0.216, suggesting that the majority of the assembled genes are under negative selection. Several genes were differentially expressed across the two pine accessions with contrasted phenotypes, including a glutathione-s-transferase, a cellulose synthase and a cobra-like protein. A large number of new markers (3334 amplifiable SSRs and 28,236 SNPs) have been identified which should facilitate future population genomics and association genetics in this species. A 384-SNP Oligo Pool Assay for genotyping with the Illumina VeraCode technology has been designed which showed an high overall SNP conversion rate (76.6%). Our results showed that Illumina next-generation sequencing is a valuable technology to obtain an extensive overview on whole transcriptomes of nonmodel species with large genomes.

Gain-of-function haplotype in the epithelial calcium channel TRPV6 is a risk factor for renal calcium stone formation.

The rate-limiting step of dietary calcium absorption in the intestine requires the brush border calcium entry channel TRPV6. The TRPV6 gene was completely sequenced in 170 renal calcium stone patients. The frequency of an ancestral TRPV6 haplotype consisting of three non-synonymous polymorphisms (C157R, M378V, M681T) was significantly higher (P = 0.039) in calcium stone formers (8.4%; derived = 502, ancestral = 46) compared to non-stone-forming individuals (5.4%; derived = 645, ancestral = 37). Mineral metabolism was investigated on four different calcium regimens: (i) free-choice diet, (ii) low calcium diet, (iii) fasting and (iv) after a 1 g oral calcium load. When patients homozygous for the derived haplotype were compared with heterozygous patients, no differences were found with respect to the plasma concentrations of 1,25-vitamin D, PTH and calcium, and the urinary excretion of calcium. In one stone-forming patient, the ancestral haplotype was found to be homozygous. This patient had absorptive hypercalciuria. We therefore expressed the ancestral protein (157R+378V+681T) in Xenopus oocytes and found a significantly enhanced calcium permeability when tested by a (45)Ca(2+) uptake assay (7.11 +/- 1.93 versus 3.61 +/- 1.01 pmol/min/oocyte for ancestral versus derived haplotype, P < 0.01). These results suggest that the ancestral gain-of-function haplotype in TRPV6 plays a role in calcium stone formation in certain forms of absorptive hypercalciuria.

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.

IL28B alleles associated with poor hepatitis C virus (HCV) clearance protect against inflammation and fibrosis in patients infected with non-1 HCV genotypes.

Genetic polymorphisms near IL28B are associated with spontaneous and treatment-induced clearance of hepatitis C virus (HCV), two processes that require the appropriate activation of the host immune responses. Intrahepatic inflammation is believed to mirror such activation, but its relationship with IL28B polymorphisms has yet to be fully appreciated. We analyzed the association of IL28B polymorphisms with histological and follow-up features in 2335 chronically HCV-infected Caucasian patients. Assessable phenotypes before any antiviral treatment included necroinflammatory activity (n = 1,098), fibrosis (n = 1,527), fibrosis progression rate (n = 1,312), and hepatocellular carcinoma development (n = 1,915). Associations of alleles with the phenotypes were evaluated by univariate analysis and multivariate logistic regression, accounting for all relevant covariates. The rare G allele at IL28B marker rs8099917-previously shown to be at risk of treatment failure-was associated with lower activity (P = 0.04), lower fibrosis (P = 0.02) with a trend toward lower fibrosis progression rate (P = 0.06). When stratified according to HCV genotype, most significant associations were observed in patients infected with non-1 genotypes (P = 0.003 for activity, P = 0.001 for fibrosis, and P = 0.02 for fibrosis progression rate), where the odds ratio of having necroinflammation or rapid fibrosis progression for patients with IL28B genotypes TG or GG versus TT were 0.48 (95% confidence intervals 0.30-0.78) and 0.56 (0.35-0.92), respectively. IL28B polymorphisms were not predictive of the development of hepatocellular carcinoma. CONCLUSION: In chronic hepatitis C, IL28B variants associated with poor response to interferon therapy may predict slower fibrosis progression, especially in patients infected with non-1 HCV genotypes.

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.

Meta-analysis of gene-level associations for rare variants based on single-variant statistics.

Meta-analysis of genome-wide association studies (GWASs) has led to the discoveries of many common variants associated with complex human diseases. There is a growing recognition that identifying "causal" rare variants also requires large-scale meta-analysis. The fact that association tests with rare variants are performed at the gene level rather than at the variant level poses unprecedented challenges in the meta-analysis. First, different studies may adopt different gene-level tests, so the results are not compatible. Second, gene-level tests require multivariate statistics (i.e., components of the test statistic and their covariance matrix), which are difficult to obtain. To overcome these challenges, we propose to perform gene-level tests for rare variants by combining the results of single-variant analysis (i.e., p values of association tests and effect estimates) from participating studies. This simple strategy is possible because of an insight that multivariate statistics can be recovered from single-variant statistics, together with the correlation matrix of the single-variant test statistics, which can be estimated from one of the participating studies or from a publicly available database. We show both theoretically and numerically that the proposed meta-analysis approach provides accurate control of the type I error and is as powerful as joint analysis of individual participant data. This approach accommodates any disease phenotype and any study design and produces all commonly used gene-level tests. An application to the GWAS summary results of the Genetic Investigation of ANthropometric Traits (GIANT) consortium reveals rare and low-frequency variants associated with human height. The relevant software is freely available.

Association study of 182 candidate genes in anorexia nervosa.

We performed association studies with 5,151 SNPs that were judged as likely candidate genetic variations conferring susceptibility to anorexia nervosa (AN) based on location under reported linkage peaks, previous results in the literature (182 candidate genes), brain expression, biological plausibility, and estrogen responsivity. We employed a case-control design that tested each SNP individually as well as haplotypes derived from these SNPs in 1,085 case individuals with AN diagnoses and 677 control individuals. We also performed separate association analyses using three increasingly restrictive case definitions for AN: all individuals with any subtype of AN (All AN: n = 1,085); individuals with AN with no binge eating behavior (AN with No Binge Eating: n = 687); and individuals with the restricting subtype of AN (Restricting AN: n = 421). After accounting for multiple comparisons, there were no statistically significant associations for any individual SNP or haplotype block with any definition of illness. These results underscore the importance of large samples to yield appropriate power to detect genotypic differences in individuals with AN and also motivate complementary approaches involving Genome-Wide Association (GWA) studies, Copy Number Variation (CNV) analyses, sequencing-based rare variant discovery assays, and pathway-based analysis in order to make up for deficiencies in traditional candidate gene approaches to AN.

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.

Common genetic variation and the control of HIV-1 in humans.

To extend the understanding of host genetic determinants of HIV-1 control, we performed a genome-wide association study in a cohort of 2,554 infected Caucasian subjects. The study was powered to detect common genetic variants explaining down to 1.3% of the variability in viral load at set point. We provide overwhelming confirmation of three associations previously reported in a genome-wide study and show further independent effects of both common and rare variants in the Major Histocompatibility Complex region (MHC). We also examined the polymorphisms reported in previous candidate gene studies and fail to support a role for any variant outside of the MHC or the chemokine receptor cluster on chromosome 3. In addition, we evaluated functional variants, copy-number polymorphisms, epistatic interactions, and biological pathways. This study thus represents a comprehensive assessment of common human genetic variation in HIV-1 control in Caucasians.

Association of the hemochromatosis gene with pazopanib-induced transaminase elevation in renal cell carcinoma.

BACKGROUND & AIMS: Pazopanib has demonstrated clinical benefit in patients with advanced renal cell carcinoma (RCC) and is generally well tolerated. However, transaminase elevations have commonly been observed. This 2-stage study sought to identify genetic determinants of alanine transaminase (ALT) elevations in pazopanib-treated white patients with RCC.¦METHODS: Data from two separate clinical studies were used to examine the association of genetic polymorphisms with maximum on-treatment ALT levels.¦RESULTS: Of 6852 polymorphisms in 282 candidate genes examined in an exploratory dataset of 115 patients, 92 polymorphisms in 40 genes were significantly associated with ALT elevation (p<0.01). Two markers (rs2858996 and rs707889) in the HFE gene, which are not yet known to be associated with hemochromatosis, showed evidence for replication. Because of multiple comparisons, there was a 12% likelihood the replication occurred by chance. These two markers demonstrated strong linkage disequilibrium (r(2)=0.99). In the combined dataset, median (25-75th percentile) maximum ALT values were 1.2 (0.7-1.9), 1.1 (0.8-2.5), and 5.4 (1.9-7.6)×ULN for rs2858996 GG (n=148), GT (n=82), and TT (n=1 2) genotypes, respectively. All 12 TT patients had a maximum ALT>ULN, and 8 (67%) had ALT≥3×ULN. The odds ratio (95% CI) for ALT≥3×ULN for TT genotype was 39.7 (2.2-703.7) compared with other genotypes. As a predictor of ALT≥3×ULN, the TT genotype had a negative predictive value of 0.83 and positive predictive value of 0.67. No TT patients developed liver failure.¦CONCLUSIONS: The rs2858996/rs707889 polymorphisms in the HFE gene may be associated with reversible ALT elevation in pazo-panib-treated patients with RCC.