Tandem repeat sequence variation as causative cis-eQTLs for protein-coding gene expression variation: the case of CSTB.

Association studies have revealed expression quantitative trait loci (eQTLs) for a large number of genes. However, the causative variants that regulate gene expression levels are generally unknown. We hypothesized that copy-number variation of sequence repeats contribute to the expression variation of some genes. Our laboratory has previously identified that the rare expansion of a repeat c.-174CGGGGCGGGGCG in the promoter region of the CSTB gene causes a silencing of the gene, resulting in progressive myoclonus epilepsy. Here, we genotyped the repeat length and quantified CSTB expression by quantitative real-time polymerase chain reaction in 173 lymphoblastoid cell lines (LCLs) and fibroblast samples from the GenCord collection. The majority of alleles contain either two or three copies of this repeat. Independent analysis revealed that the c.-174CGGGGCGGGGCG repeat length is strongly associated with CSTB expression (P = 3.14 × 10(-11)) in LCLs only. Examination of both genotyped and imputed single-nucleotide polymorphisms (SNPs) within 2 Mb of CSTB revealed that the dodecamer repeat represents the strongest cis-eQTL for CSTB in LCLs. We conclude that the common two or three copy variation is likely the causative cis-eQTL for CSTB expression variation. More broadly, we propose that polymorphic tandem repeats may represent the causative variation of a fraction of cis-eQTLs in the genome.

Common variants in the ATP2B1 gene are associated with susceptibility to hypertension: the Japanese Millennium Genome Project.

Hypertension is one of the most common complex genetic disorders. We have described previously 38 single nucleotide polymorphisms (SNPs) with suggestive association with hypertension in Japanese individuals. In this study we extend our previous findings by analyzing a large sample of Japanese individuals (n=14 105) for the most associated SNPs. We also conducted replication analyses in Japanese of susceptibility loci for hypertension identified recently from genome-wide association studies of European ancestries. Association analysis revealed significant association of the ATP2B1 rs2070759 polymorphism with hypertension (P=5.3×10(-5); allelic odds ratio: 1.17 [95% CI: 1.09 to 1.26]). Additional SNPs in ATP2B1 were subsequently genotyped, and the most significant association was with rs11105378 (odds ratio: 1.31 [95% CI: 1.21 to 1.42]; P=4.1×10(-11)). Association of rs11105378 with hypertension was cross-validated by replication analysis with the Global Blood Pressure Genetics consortium data set (odds ratio: 1.13 [95% CI: 1.05 to 1.21]; P=5.9×10(-4)). Mean adjusted systolic blood pressure was highly significantly associated with the same SNP in a meta-analysis with individuals of European descent (P=1.4×10(-18)). ATP2B1 mRNA expression levels in umbilical artery smooth muscle cells were found to be significantly different among rs11105378 genotypes. Seven SNPs discovered in published genome-wide association studies were also genotyped in the Japanese population. In the combined analysis with replicated 3 genes, FGF5 rs1458038, CYP17A1, rs1004467, and CSK rs1378942, odds ratio of the highest risk group was 2.27 (95% CI: 1.65 to 3.12; P=4.6×10(-7)) compared with the lower risk group. In summary, this study confirmed common genetic variation in ATP2B1, as well as FGF5, CYP17A1, and CSK, to be associated with blood pressure levels and risk of hypertension.

Depressive disorder moderates the effect of the FTO gene on body mass index.

There is evidence that obesity-related disorders are increased among people with depression. Variation in the FTO (fat mass and obesity associated) gene has been shown to contribute to common forms of human obesity. This study aimed to investigate the genetic influence of polymorphisms in FTO in relation to body mass index (BMI) in two independent samples of major depressive disorder (MDD) cases and controls. We analysed 88 polymorphisms in the FTO gene in a clinically ascertained sample of 2442 MDD cases and 809 controls (Radiant Study). In all, 8 of the top 10 single-nucleotide polymorphisms (SNPs) showing the strongest associations with BMI were followed-up in a population-based cohort (PsyCoLaus Study) consisting of 1292 depression cases and 1690 controls. Linear regression analyses of the FTO variants and BMI yielded 10 SNPs significantly associated with increased BMI in the depressive group but not the control group in the Radiant sample. The same pattern was found in the PsyCoLaus sample. We found a significant interaction between genotype and affected status in relation to BMI for seven SNPs in Radiant (P<0.0057), with PsyCoLaus giving supportive evidence for five SNPs (P-values between 0.03 and 0.06), which increased in significance when the data were combined in a meta-analysis. This is the first study investigating FTO and BMI within the context of MDD, and the results indicate that having a history of depression moderates the effect of FTO on BMI. This finding suggests that FTO is involved in the mechanism underlying the association between mood disorders and obesity.

Genetic vulnerability factor for schizophrenia: association with glutamate cysteine ligase modifier gene and abnormal enzyme activity

Background: We previously reported in schizophrenia patients a decreased level of glutathione ([GSH]), the principal non-protein antioxidant and redox regulator, both in cerebrospinal-fluid and prefrontal cortex. To identify possible genetic causation, we studied genes involved in GSH metabolism. Methods: Genotyping: mass spectrometry analysis of polymerase chain reaction (PCR) amplified DNA fragments purified from peripheral blood. Gene expression: real-time PCR of total RNA isolated from fibroblast cultures derived from skin of patients (DSM-IV) and healthy controls (DIGS). Results: Case-control association study of single nucleotide polymorphisms (SNP) from the GSH key synthesizing enzyme glutamate-cysteine-ligase (GCL) modifier subunit (GCLM) was performed in two populations: Swiss (patients/controls: 40/31) and Danish (349/348). We found a strong association of SNP rs2301022 in GCLM gene (Danish: c2=3.2; P=0.001 after correction for multiple testing). Evidence for GCLM as a risk factor was confirmed in linkage study of NIMH families. Moreover, we observed a decrease in GCLM mRNA levels in patient fibroblasts, consistently with the association study. Interestingly, Dalton and collaborators reported in GCLM knock-out mice an increased feedback inhibition of GCL activity, resulting in 60% decrease of brain [GSH], a situation analogous to patients. These mice also exhibited an increased sensitivity to oxidative stress. Similarly, under oxidative stress conditions, GCL enzymatic activity was also decreased in patient fibroblasts. Conclusions: These results at the genetic and functional levels, combined with observations that GSH deficient models reveal morphological, electrophysiological, and behavioral anomalies analogous to those observed in patients, suggest that GCLM allelic variant is a vulnerability factor for schizophrenia.

IL28B single nucleotide polymorphisms in the treatment of hepatitis C.

Recent genome-wide association studies (GWAS) have identified genetic variations near the IL28B gene which are strongly associated with spontaneous and treatment-induced clearance of hepatitis C virus (HCV) infection. Protective IL28B variations are strongly associated with on-treatment viral kinetics and approximately 2-fold increased sustained virologic response (SVR) rates in HCV genotype 1 and 4 patients. In HCV genotype 1 patients, IL28B variations were shown to be the strongest pre-treatment predictor of virologic response. In the treatment of HCV genotype 2 and 3 infected patients, IL28B variations play only a minor role. Preliminary data indicate that IL28B variations are also associated with treatment outcome of regimens, including directly acting antiviral (DAA) agents, though their impact seems to be attenuated compared to standard treatment. Here, we review these important findings and discuss possible implications for clinical decision making in the treatment of HCV infection.

Contribution of CYP2B6 alleles in explaining extreme (S)-methadone plasma levels: a CYP2B6 gene resequencing study.

BACKGROUND: (S)-Methadone, metabolized mainly by CYP2B6, shows a wide interindividual variability in its pharmacokinetics and pharmacodynamics. METHODS: Resequencing of the CYP2B6 gene was performed in 12 and 35 selected individuals with high (S)-methadone plasma exposure and low (S)-methadone plasma exposure, respectively, from a previously described cohort of 276 patients undergoing methadone maintenance treatment. Selected genetic polymorphisms were then analyzed in the complete cohort. RESULTS: The rs35303484 (*11; c136A>G; M46V) polymorphism was overrepresented in the high (S)-methadone level group, whereas the rs3745274 (*9; c516G>T; Q172H), rs2279344 (c822+183G>A), and rs8192719 (c1294+53C>T) polymorphisms were underrepresented in the low (S)-methadone level group, suggesting an association with decreased CYP2B6 activity. Conversely, the rs3211371 (*5; c1459C>T; R487C) polymorphism was overrepresented in the low-level group, indicating an increased CYP2B6 activity. A higher allele frequency was found in the high-level group compared with the low-level group for rs3745274 (*9; c516G>T; Q172H), rs2279343 (*4; c785A>G; K262R) (together representing CYP2B6*6), rs8192719 (c1294+53C>T), and rs2279344 (c822+183G>A), suggesting their involvement in decreased CYP2B6 activity. These results should be replicated in larger independent cohorts. CONCLUSION: Known genetic polymorphisms in CYP2B6 contribute toward explaining extreme (S)-methadone plasma levels observed in a cohort of patients following methadone maintenance treatment.

Interleukin-1 receptor antagonist gene (IL-1RN) polymorphism is a predictive factor of clinical pregnancy after IVF.

BACKGROUND: Only 25% of IVF transfer cycles lead to a clinical pregnancy, calling for continued technical progress but also more in depth analysis of patients' individual characteristics. The interleukin-1 (IL-1) system and matrix metalloproteinases (MMPs) are strongly implicated in embryo implantation. The genes coding for IL-1Ra (gene symbol IL-1RN), IL-1beta, MMP2 and MMP9 bear functional polymorphisms. We analysed the maternal genetic profile at these polymorphic sites in IVF patients, to determine possible correlations with IVF outcome. METHODS: One hundred and sixty women undergoing an IVF cycle were enrolled and a buccal smear was obtained. The presence of IL-1RN variable number of tandem repeats and IL-1B + 3953, MMP2-1306 and MMP9-1562 single nucleotide substitutions were determined. Patients were divided into pregnancy failures (119), biochemical pregnancies (8) and clinical pregnancies (33). RESULTS: There was a 40% decrease in IL-1RN*2 allele frequency (P = 0.024) and a 45% decrease in IL-1RN*2 carrier status in the clinical pregnancy group as compared to the pregnancy failure group (P = 0.017). This decrease was still statistically significant after a multivariate logistic regression analysis. The likelihood of a clinical pregnancy was decreased accordingly in IL-1RN*2 carriers: odds ratio = 0.349, 95% confidence interval = 0.2-0.8, P = 0.017. The IL-1B, MMP2 and MMP9 polymorphisms showed no correlation with IVF outcome. CONCLUSIONS: IL-1RN*2 allele carriage is associated with a poor prognosis of achieving a pregnancy after IVF.

Meta-analysis of gene expression profiles in breast cancer: toward a unified understanding of breast cancer subtyping and prognosis signatures.

INTRODUCTION: Breast cancer subtyping and prognosis have been studied extensively by gene expression profiling, resulting in disparate signatures with little overlap in their constituent genes. Although a previous study demonstrated a prognostic concordance among gene expression signatures, it was limited to only one dataset and did not fully elucidate how the different genes were related to one another nor did it examine the contribution of well-known biological processes of breast cancer tumorigenesis to their prognostic performance. METHOD: To address the above issues and to further validate these initial findings, we performed the largest meta-analysis of publicly available breast cancer gene expression and clinical data, which are comprised of 2,833 breast tumors. Gene coexpression modules of three key biological processes in breast cancer (namely, proliferation, estrogen receptor [ER], and HER2 signaling) were used to dissect the role of constituent genes of nine prognostic signatures. RESULTS: Using a meta-analytical approach, we consolidated the signatures associated with ER signaling, ERBB2 amplification, and proliferation. Previously published expression-based nomenclature of breast cancer 'intrinsic' subtypes can be mapped to the three modules, namely, the ER-/HER2- (basal-like), the HER2+ (HER2-like), and the low- and high-proliferation ER+/HER2- subtypes (luminal A and B). We showed that all nine prognostic signatures exhibited a similar prognostic performance in the entire dataset. Their prognostic abilities are due mostly to the detection of proliferation activity. Although ER- status (basal-like) and ERBB2+ expression status correspond to bad outcome, they seem to act through elevated expression of proliferation genes and thus contain only indirect information about prognosis. Clinical variables measuring the extent of tumor progression, such as tumor size and nodal status, still add independent prognostic information to proliferation genes. CONCLUSION: This meta-analysis unifies various results of previous gene expression studies in breast cancer. It reveals connections between traditional prognostic factors, expression-based subtyping, and prognostic signatures, highlighting the important role of proliferation in breast cancer prognosis.

Restricted gene flow at specific parts of the shrew genome in chromosomal hybrid zones.

The species and races of the shrews of the Sorex araneus group exhibit a broad range of chromosomal polymorphisms. European taxa of this group are parapatric and form contact or hybrid zones that span an extraordinary variety of situations, ranging from absolute genetic isolation to almost free gene flow. This variety seems to depend for a large part on the chromosome composition of populations, which are primarily differentiated by various Robertsonian fusions of a subset of acrocentric chromosomes. Previous studies suggested that chromosomal rearrangements play a causative role in the speciation process. In such models, gene flow should be more restricted for markers on chromosomes involved in rearrangements than on chromosomes common in both parent species. In the present study, we address the possibility of such differential gene flow in the context of two genetically very similar but karyotypically different hybrid zones between species of the S. araneus group using microsatellite loci mapped to the chromosome arm level. Interspecific genetic structure across rearranged chromosomes was in general larger than across common chromosomes. However, the difference between the two classes of chromosomes was only significant in the hybrid zone where the complexity of hybrids is expected to be larger. These differences did not distinguish populations within species. Therefore, the rearranged chromosomes appear to affect the reproductive barrier between karyotypic species, although the strength of this effect depends on the complexity of the hybrids produced.

Contrasting diffusion of Quaternary gene pools across Europe: the case of the arctic-alpine Gentiana nivalis L. (Gentianaceae).

The fate of European arctic-alpine species during Pleistocene climatic oscillations still remains debated. Did these cold-adapted species invade much of the continental steppe or did they remain restricted to warmer slopes of inner mountain massifs? To examine this question, we investigated the phylogeography of Gentiana nivalis, a typical European arctic-alpine plant species. Genome fingerprinting analyses revealed that four genetic pools are actually unevenly distributed across the continent. One cluster covers almost all mountain massifs as well as northern areas, and thus coincides with a scenario of past distribution covering a large part of the European glacial steppe. In contrast, the three other lineages are strongly restricted spatially to western, central, and eastern Alps, respectively, thus arguing towards a scenario of in situ glacial survival. The coexistence of lineages with such contrasting demographic histories in Europe challenges our classical view of refugia and corroborates several hypotheses of biogeographers from the twentieth century.

The vitamin D receptor gene bAt (CCA) haplotype impairs the response to pegylated-interferon/ribavirin-based therapy in chronic hepatitis C patients.

BACKGROUND: Chronic hepatitis C infection is a major cause of end-stage liver disease. Therapy outcome is influenced by 25-OH vitamin D deficiency. To further address this observation, our study investigates the impact of the vitamin D receptor (NR1I1) haplotype and combined effects of plasma vitamin D levels in a well-described cohort of hepatitis C patients. METHODS: A total of 155 chronic hepatitis C patients were recruited from the Swiss Hepatitis C Cohort Study for NR1I1 genotyping and plasma 25-OH vitamin D level measurement. NR1I1 genotype data and combined effects of plasma 25-OH vitamin D level were analysed regarding therapy response (sustained virological response). RESULTS: A strong association was observed between therapy non-response and the NR1I1 CCA (bAt) haplotype consisting of rs1544410 (BsmI) C, rs7975232 (ApaI) C and rs731236 (TaqI) A alleles. Of the HCV patients carrying the CCA haplotype, 50.3% were non-responders (odds ratio [OR] 1.69, 95% CI 1.07, 2.67; P=0.028). A similar association was observed for the combinational CCCCAA genotype (OR 2.94, 95% CI 1.36, 6.37; P=0.007). The combinational CCCCAA genotype was confirmed as an independent risk factor for non-response in multivariate analysis (OR 2.50, 95% CI 1.07, 5.87; P=0.034). Analysing combined effects, a significant impact of low 25-OH vitamin D levels on sustained virological response were only seen in patients with the unfavourable NR1I1 CCA (bAt) haplotype (OR for non-SVR 3.55; 95% CI 1.005, 12.57; P=0.049). CONCLUSIONS: NR1I1 vitamin D receptor polymorphisms influence response to pegylated-interferon/ribavirin-based therapy in chronic hepatitis C and exert an additive genetic predisposition to previously described low 25-OH vitamin D serum levels.

Diet and gene interactions influence the skeletal response to polyunsaturated fatty acids.

Diets rich in omega-3s have been thought to prevent both obesity and osteoporosis. However, conflicting findings are reported, probably as a result of gene by nutritional interactions. Peroxisome proliferator-activated receptor-gamma (PPARγ) is a nuclear receptor that improves insulin sensitivity but causes weight gain and bone loss. Fish oil is a natural agonist for PPARγ and thus may exert its actions through the PPARγ pathway. We examined the role of PPARγ in body composition changes induced by a fish or safflower oil diet using two strains of C57BL/6J (B6); i.e. B6.C3H-6T (6T) congenic mice created by backcrossing a small locus on Chr 6 from C3H carrying 'gain of function' polymorphisms in the Pparγ gene onto a B6 background, and C57BL/6J mice. After 9months of feeding both diets to female mice, body weight, percent fat and leptin levels were less in mice fed the fish oil vs those fed safflower oil, independent of genotype. At the skeletal level, fish oil preserved vertebral bone mineral density (BMD) and microstructure in B6 but not in 6T mice. Moreover, fish oil consumption was associated with an increase in bone marrow adiposity and a decrease in BMD, cortical thickness, ultimate force and plastic energy in femur of the 6T but not the B6 mice. These effects paralleled an increase in adipogenic inflammatory and resorption markers in 6T but not B6. Thus, compared to safflower oil, fish oil (high ratio omega-3/-6) prevents weight gain, bone loss, and changes in trabecular microarchitecture in the spine with age. These beneficial effects are absent in mice with polymorphisms in the Pparγ gene (6T), supporting the tenet that the actions of n-3 fatty acids on bone microstructure are likely to be genotype dependent. Thus caution must be used in interpreting dietary intervention trials with skeletal endpoints in mice and in humans.

Interleukin-1 receptor antagonist gene polymorphism and circulating levels of human immunodeficiency virus type 1 RNA in Brazilian women.

Interleukin-1 receptor antagonist (IL-1ra) gene polymorphisms in 83 human immunodeficiency virus (HIV)-seropositive women were evaluated. Fourteen of the subjects (16.9%) were homozygous for IL-1ra allele 2 (IL-1RN*2). These women had a lower median level of HIV RNA than did women homozygous for allele 1 (IL-1RN*1) (P = 0.01) or heterozygous for both alleles (P = 0.04). Among 46 subjects not receiving antiretroviral treatment, HIV levels were also reduced in IL-1RN*2 homozygous individuals (P < 0.05). There was no relation between IL-1ra alleles and CD4 levels.

Absence of association between specific common variants of the obesity-related FTO gene and psychological and behavioral eating disorder phenotypes.

Extensive population-based genome-wide association studies have identified an association between the FTO gene and BMI; however, the mechanism of action is still unknown. To determine whether FTO may influence weight regulation through psychological and behavioral factors, seven single-nucleotide polymorphisms (SNPs) of the FTO gene were genotyped in 1,085 individuals with anorexia nervosa (AN) and 677 healthy weight controls from the international Price Foundation Genetic Studies of Eating Disorders. Each SNP was tested in association with eating disorder phenotypes and measures that have previously been associated with eating behavior pathology: trait anxiety, harm-avoidance, novelty seeking, impulsivity, obsessionality, compulsivity, and concern over mistakes. After appropriate correction for multiple comparisons, no significant associations between individual FTO gene SNPs and eating disorder phenotypes or related eating behavior pathology were identified in cases or controls. Thus, this study found no evidence that FTO gene variants associated with weight regulation in the general population are associated with eating disorder phenotypes in AN participants or matched controls. © 2011 Wiley-Liss, Inc.

Polymorphisms in Toll-like receptor 9 influence the clinical course of HIV-1 infection.

BACKGROUND: The clinical course of HIV-1 infection is highly variable among individuals, at least in part as a result of genetic polymorphisms in the host. Toll-like receptors (TLRs) have a key role in innate immunity and mutations in the genes encoding these receptors have been associated with increased or decreased susceptibility to infections. OBJECTIVES: To determine whether single-nucleotide polymorphisms (SNPs) in TLR2-4 and TLR7-9 influenced the natural course of HIV-1 infection. METHODS: Twenty-eight SNPs in TLRs were analysed in HAART-naive HIV-positive patients from the Swiss HIV Cohort Study. The SNPs were detected using Sequenom technology. Haplotypes were inferred using an expectation-maximization algorithm. The CD4 T cell decline was calculated using a least-squares regression. Patients with a rapid CD4 cell decline, less than the 15th percentile, were defined as rapid progressors. The risk of rapid progression associated with SNPs was estimated using a logistic regression model. Other candidate risk factors included age, sex and risk groups (heterosexual, homosexual and intravenous drug use). RESULTS: Two SNPs in TLR9 (1635A/G and +1174G/A) in linkage disequilibrium were associated with the rapid progressor phenotype: for 1635A/G, odds ratio (OR), 3.9 [95% confidence interval (CI),1.7-9.2] for GA versus AA and OR, 4.7 (95% CI,1.9-12.0) for GG versus AA (P = 0.0008). CONCLUSION: Rapid progression of HIV-1 infection was associated with TLR9 polymorphisms. Because of its potential implications for intervention strategies and vaccine developments, additional epidemiological and experimental studies are needed to confirm this association.