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

Background: T reatment o f chronic hepatitis C i s evolving, a nd direct acting antivirals ( DAAs) are now a dded to p egylated interferon-α ( Peg- INF-α) and ribavirin (RBV) for the treatment o f hepatitis C v irus ( HCV) genotype 1 infection. DAAs c ause d ifferent side effects and can even worsen RBV induced hemolytic anemia. T herefore, identifying host genetic d eterminants of R BV bioavailability and therapeutic e fficacy will remain crucial for individualized treatment. Recent d ata showed associations between R BV induced h emolytic anemia and genetic polymorphisms o f concentrative nucleoside transporters s uch as C NT3 (SLC28A3) and i nosine t riphosphatase (ITPA). T o analyze t he association of genetic variants of SLC28 transporters and ITPA with RBV induced hemolytic anemia and treatment o utcome. Methods: I n our study, 173 patients f rom t he S wiss Hepatitis C C ohort Study and 2 2 patients from Swiss Association for the Study of the Liver study 24 (61% HCV g enotype 1, 3 9% genotypes 2 o r 3) were analyzed for SLC28A2 single nucleotide p olymorphism (SNP) rs11854484, SLC28A3 rs56350726 and SLC28A3 rs10868138 as well as ITPA SNPs rs1127354 and rs7270101. RBV serum levels during treatment were measured in 49 patients. Results: SLC28A2 r s11854484 genotype TT was associated with significantly higher dosage- and body weight-adjusted RBV levels as compared to genotypes TC and CC (p=0.04 and p=0.02 at weeks 4 and 8, respectively). ITPA SNPs rs1127354 and rs7270101 were associated with h emolytic a nemia both in genotype as w ell as i n allelic a nalyses. SLC28A3 rs56350726 genotype TT (vs. AT/AA, RR=2.1; 95% CI 1.1-4.1) as well as the T allele (vs. A; RR=1.8, 95% CI 1.1-3.2) were associated with increased SVR rates. The combined analysis of overall ITPA activity and SLC28 v ariants together revealed n o significant a dditive effects on either treatment-related anemia or SVR. Conclusions: T he newly identified association between RBV serum levels a nd SLC28A2 rs11854484 genotype as well as the replicated association of ITPA and SLC28A3 g enetic p olymorphisms w ith RBV induced hemolytic anemia and treatment r esponse underpin the need for further studies on host genetic d eterminants of R BV bioavailability and therapeutic e fficacy f or individualized treatment of chronic hepatitis C.

Social chromosome variants differentially affect queen determination and the survival of workers in the fire ant Solenopsis invicta.

Intraspecific variation in social organization is common, yet the underlying causes are rarely known. An exception is the fire ant Solenopsis invicta in which the existence of two distinct forms of social colony organization is under the control of the two variants of a pair of social chromosomes, SB and Sb. Colonies containing exclusively SB/SB workers accept only one single queen and she must be SB/SB. By contrast, when colonies contain more than 10% of SB/Sb workers, they accept several queens but only SB/Sb queens. The variants of the social chromosome are associated with several additional important phenotypic differences, including the size, fecundity and dispersal strategies of queens, aggressiveness of workers, and sperm count in males. However, little is known about whether social chromosome variants affect fitness in other life stages. Here, we perform experiments to determine whether differential selection occurs during development and in adult workers. We find evidence that the Sb variant of the social chromosome increases the likelihood of female brood to develop into queens and that adult SB/Sb workers, the workers that cull SB/SB queens, are overrepresented in comparison to SB/SB workers. This demonstrates that supergenes such as the social chromosome can have complex effects on phenotypes at various stages of development.

Rare genomic structural variants in complex disease: lessons from the replication of associations with obesity.

The limited ability of common variants to account for the genetic contribution to complex disease has prompted searches for rare variants of large effect, to partly explain the 'missing heritability'. Analyses of genome-wide genotyping data have identified genomic structural variants (GSVs) as a source of such rare causal variants. Recent studies have reported multiple GSV loci associated with risk of obesity. We attempted to replicate these associations by similar analysis of two familial-obesity case-control cohorts and a population cohort, and detected GSVs at 11 out of 18 loci, at frequencies similar to those previously reported. Based on their reported frequencies and effect sizes (OR≥25), we had sufficient statistical power to detect the large majority (80%) of genuine associations at these loci. However, only one obesity association was replicated. Deletion of a 220 kb region on chromosome 16p11.2 has a carrier population frequency of 2×10(-4) (95% confidence interval [9.6×10(-5)-3.1×10(-4)]); accounts overall for 0.5% [0.19%-0.82%] of severe childhood obesity cases (P = 3.8×10(-10); odds ratio = 25.0 [9.9-60.6]); and results in a mean body mass index (BMI) increase of 5.8 kg.m(-2) [1.8-10.3] in adults from the general population. We also attempted replication using BMI as a quantitative trait in our population cohort; associations with BMI at or near nominal significance were detected at two further loci near KIF2B and within FOXP2, but these did not survive correction for multiple testing. These findings emphasise several issues of importance when conducting rare GSV association, including the need for careful cohort selection and replication strategy, accurate GSV identification, and appropriate correction for multiple testing and/or control of false discovery rate. Moreover, they highlight the potential difficulty in replicating rare CNV associations across different populations. Nevertheless, we show that such studies are potentially valuable for the identification of variants making an appreciable contribution to complex disease.

Physical activity attenuates the influence of FTO variants on obesity risk: a meta-analysis of 218,166 adults and 19,268 children.

BACKGROUND: The FTO gene harbors the strongest known susceptibility locus for obesity. While many individual studies have suggested that physical activity (PA) may attenuate the effect of FTO on obesity risk, other studies have not been able to confirm this interaction. To confirm or refute unambiguously whether PA attenuates the association of FTO with obesity risk, we meta-analyzed data from 45 studies of adults (n = 218,166) and nine studies of children and adolescents (n = 19,268). METHODS AND FINDINGS: All studies identified to have data on the FTO rs9939609 variant (or any proxy [r(2)>0.8]) and PA were invited to participate, regardless of ethnicity or age of the participants. PA was standardized by categorizing it into a dichotomous variable (physically inactive versus active) in each study. Overall, 25% of adults and 13% of children were categorized as inactive. Interaction analyses were performed within each study by including the FTO×PA interaction term in an additive model, adjusting for age and sex. Subsequently, random effects meta-analysis was used to pool the interaction terms. In adults, the minor (A-) allele of rs9939609 increased the odds of obesity by 1.23-fold/allele (95% CI 1.20-1.26), but PA attenuated this effect (p(interaction)  = 0.001). More specifically, the minor allele of rs9939609 increased the odds of obesity less in the physically active group (odds ratio  = 1.22/allele, 95% CI 1.19-1.25) than in the inactive group (odds ratio  = 1.30/allele, 95% CI 1.24-1.36). No such interaction was found in children and adolescents. CONCLUSIONS: The association of the FTO risk allele with the odds of obesity is attenuated by 27% in physically active adults, highlighting the importance of PA in particular in those genetically predisposed to obesity.

Immunochip SNP array identifies novel genetic variants conferring susceptibility to candidaemia.

Candidaemia is the fourth most common cause of bloodstream infection, with a high mortality rate of up to 40%. Identification of host genetic factors that confer susceptibility to candidaemia may aid in designing adjunctive immunotherapeutic strategies. Here we hypothesize that variation in immune genes may predispose to candidaemia. We analyse 118,989 single-nucleotide polymorphisms (SNPs) across 186 loci known to be associated with immune-mediated diseases in the largest candidaemia cohort to date of 217 patients of European ancestry and a group of 11,920 controls. We validate the significant associations by comparison with a disease-matched control group. We observe significant association between candidaemia and SNPs in the CD58 (P = 1.97 × 10(-11); odds ratio (OR) = 4.68), LCE4A-C1orf68 (P = 1.98 × 10(-10); OR = 4.25) and TAGAP (P = 1.84 × 10(-8); OR = 2.96) loci. Individuals carrying two or more risk alleles have an increased risk for candidaemia of 19.4-fold compared with individuals carrying no risk allele. We identify three novel genetic risk factors for candidaemia, which we subsequently validate for their role in antifungal host defence.

Syntaxin1a variants lacking an N-peptide or bearing the LE mutation bind to Munc18a in a closed conformation.

In neurons, soluble N-ethylmaleimide-sensitive factor attachment receptor (SNARE) proteins drive the fusion of synaptic vesicles to the plasma membrane through the formation of a four-helix SNARE complex. Members of the Sec1/Munc18 protein family regulate membrane fusion through interactions with the syntaxin family of SNARE proteins. The neuronal protein Munc18a interacts with a closed conformation of the SNARE protein syntaxin1a (Syx1a) and with an assembled SNARE complex containing Syx1a in an open conformation. The N-peptide of Syx1a (amino acids 1-24) has been implicated in the transition of Munc18a-bound Syx1a to Munc18a-bound SNARE complex, but the underlying mechanism is not understood. Here we report the X-ray crystal structures of Munc18a bound to Syx1a with and without its native N-peptide (Syx1aΔN), along with small-angle X-ray scattering (SAXS) data for Munc18a bound to Syx1a, Syx1aΔN, and Syx1a L165A/E166A (LE), a mutation thought to render Syx1a in a constitutively open conformation. We show that all three complexes adopt the same global structure, in which Munc18a binds a closed conformation of Syx1a. We also identify a possible structural connection between the Syx1a N-peptide and SNARE domain that might be important for the transition of closed-to-open Syx1a in SNARE complex assembly. Although the role of the N-peptide in Munc18a-mediated SNARE complex assembly remains unclear, our results demonstrate that the N-peptide and LE mutation have no effect on the global conformation of the Munc18a-Syx1a complex.

Response to antiretroviral treatment in HIV-1-infected individuals with allelic variants of the multidrug resistance transporter 1: a pharmacogenetics study.

BACKGROUND:HIV-1-infected patients vary considerably by their response to antiretroviral treatment, drug concentrations in plasma, toxic events, and rate of immune recovery. This variability could have a genetic basis. We did a pharmacogenetics study to analyse the association between response to antiretroviral treatment and allelic variants of several genes. METHODS:In 123 patients, we did PCR analyses of the gene for the multidrug-resistance transporter (MDR1), which codes for P-glycoprotein, of genes coding for isoenzymes of cytochrome P450, CYP3A4, CYP3A5, CYP2D6, and CYP2C19, and of the gene for the chemokine receptor CCR5. We measured concentrations in plasma of the antiretroviral agents efavirenz and nelfinavir by high-performance liquid-chromatography, and measured levels of P-glycoprotein expression, CD4-cell count, and HIV-1 viraemia. FINDINGS: Median drug concentrations in patients with the MDR1 3435 TT, CT, and CC genotypes were at the 30th, 50th, and 75th percentiles, respectively (p=0.0001). In patients with CYP2D6 extensive-metaboliser or poor-metaboliser alleles, median drug concentrations were at percentiles 45 and 62.5, respectively (p=0.04). Patients with the MDR1 TT genotype 6 months after starting treatment had a greater rise in CD4-cell count (257 cells/microL) than patients with the CT (165 cells/microL) and CC (121 cells/microL) genotype (p=0.0048), and the best recovery of naïve CD4-cells. INTERPRETATION:The polymorphism MDR1 3435 C/T predicts immune recovery after initiation of antiretroviral treatment. This finding suggests that P-glycoprotein has an important role in admittance of antiretroviral drugs to restricted compartments in vivo.

A genome-wide association study of early menopause and the combined impact of identified variants.

Early menopause (EM) affects up to 10% of the female population, reducing reproductive lifespan considerably. Currently, it constitutes the leading cause of infertility in the western world, affecting mainly those women who postpone their first pregnancy beyond the age of 30 years. The genetic aetiology of EM is largely unknown in the majority of cases. We have undertaken a meta-analysis of genome-wide association studies (GWASs) in 3493 EM cases and 13 598 controls from 10 independent studies. No novel genetic variants were discovered, but the 17 variants previously associated with normal age at natural menopause as a quantitative trait (QT) were also associated with EM and primary ovarian insufficiency (POI). Thus, EM has a genetic aetiology which overlaps variation in normal age at menopause and is at least partly explained by the additive effects of the same polygenic variants. The combined effect of the common variants captured by the single nucleotide polymorphism arrays was estimated to account for ∼30% of the variance in EM. The association between the combined 17 variants and the risk of EM was greater than the best validated non-genetic risk factor, smoking.

Caffeine intake and CYP1A2 variants associated with high caffeine intake protect non-smokers from hypertension.

The 15q24.1 locus, including CYP1A2, is associated with blood pressure (BP). The CYP1A2 rs762551 C allele is associated with lower CYP1A2 enzyme activity. CYP1A2 metabolizes caffeine and is induced by smoking. The association of caffeine consumption with hypertension remains controversial. We explored the effects of CYP1A2 variants and CYP1A2 enzyme activity on BP, focusing on caffeine as the potential mediator of CYP1A2 effects. Four observational (n = 16 719) and one quasi-experimental studies (n = 106) including European adults were conducted. Outcome measures were BP, caffeine intake, CYP1A2 activity and polymorphisms rs762551, rs1133323 and rs1378942. CYP1A2 variants were associated with hypertension in non-smokers, but not in smokers (CYP1A2-smoking interaction P = 0.01). Odds ratios (95% CIs) for hypertension for rs762551 CC, CA and AA genotypes were 1 (reference), 0.78 (0.59-1.02) and 0.66 (0.50-0.86), respectively, P = 0.004. Results were similar for the other variants. Higher CYP1A2 activity was linearly associated with lower BP after quitting smoking (P = 0.049 and P = 0.02 for systolic and diastolic BP, respectively), but not while smoking. In non-smokers, the CYP1A2 variants were associated with higher reported caffeine intake, which in turn was associated with lower odds of hypertension and lower BP (P = 0.01). In Mendelian randomization analyses using rs1133323 as instrument, each cup of caffeinated beverage was negatively associated with systolic BP [-9.57 (-16.22, -2.91) mmHg]. The associations of CYP1A2 variants with BP were modified by reported caffeine intake. These observational and quasi-experimental results strongly support a causal role of CYP1A2 in BP control via caffeine intake.

Orosomucoid system: 17 additional orosomucoid variants and proposal for a new nomenclature.

There are two forms of orosomucoid (ORM) in the sera of most individuals. They are encoded by two separate but closely linked loci, ORM1 and ORM2. A number of variants have been identified in various populations. Duplication and nonexpression are also observed in some populations. Thus, the ORM system is very complicated and its nomenclature is very confusing. In order to propose a new nomenclature, ORM variants detected by several laboratories have been compared and characterized by isoelectric focusing (IEF) followed by immunoprinting. A total of 57 different alleles including 17 new ones were identified. The 27 alleles were assigned to the ORM1 locus, and the others to the ORM2 locus. The designations ORM*F1, ORM1*F2, ORM1*S and ORM2*M were adopted for the four common alleles instead of ORM1*1, ORM1*3, ORM1*2 and ORM2*1 (ORM2*A), respectively. The variants were designated alpha numerically according to their relative mobilities after IEF in a pH gradient of 4.5-5.4 with Triton X-100 and glycerol. For the duplicated genes a prefix is added to a combined name of two alleles, e.g. ORM1*dB9S. Silent alleles were named ORM1*Q0 and ORM2*Q0 conventionally. In addition, the effects of diseases to ORM band patterns after IEF are also discussed.

Genetic variants in novel pathways influence blood pressure and cardiovascular disease risk.

Blood pressure is a heritable trait influenced by several biological pathways and responsive to environmental stimuli. Over one billion people worldwide have hypertension (≥140 mm Hg systolic blood pressure or  ≥90 mm Hg diastolic blood pressure). Even small increments in blood pressure are associated with an increased risk of cardiovascular events. This genome-wide association study of systolic and diastolic blood pressure, which used a multi-stage design in 200,000 individuals of European descent, identified sixteen novel loci: six of these loci contain genes previously known or suspected to regulate blood pressure (GUCY1A3-GUCY1B3, NPR3-C5orf23, ADM, FURIN-FES, GOSR2, GNAS-EDN3); the other ten provide new clues to blood pressure physiology. A genetic risk score based on 29 genome-wide significant variants was associated with hypertension, left ventricular wall thickness, stroke and coronary artery disease, but not kidney disease or kidney function. We also observed associations with blood pressure in East Asian, South Asian and African ancestry individuals. Our findings provide new insights into the genetics and biology of blood pressure, and suggest potential novel therapeutic pathways for cardiovascular disease prevention.

Long-term efficacy of ciliary muscle gene transfer of three sFlt-1 variants in a rat model of laser-induced choroidal neovascularization.

Inhibition of vascular endothelial growth factor (VEGF) has become the standard of care for patients presenting with wet age-related macular degeneration. However, monthly intravitreal injections are required for optimal efficacy. We have previously shown that electroporation enabled ciliary muscle gene transfer results in sustained protein secretion into the vitreous for up to 9 months. Here, we evaluated the long-term efficacy of ciliary muscle gene transfer of three soluble VEGF receptor-1 (sFlt-1) variants in a rat model of laser-induced choroidal neovascularization (CNV). All three sFlt-1 variants significantly diminished vascular leakage and neovascularization as measured by fluorescein angiography (FA) and flatmount choroid at 3 weeks. FA and infracyanine angiography demonstrated that inhibition of CNV was maintained for up to 6 months after gene transfer of the two shortest sFlt-1 variants. Throughout, clinical efficacy was correlated with sustained VEGF neutralization in the ocular media. Interestingly, treatment with sFlt-1 induced a 50% downregulation of VEGF messenger RNA levels in the retinal pigment epithelium and the choroid. We demonstrate for the first time that non-viral gene transfer can achieve a long-term reduction of VEGF levels and efficacy in the treatment of CNV.Gene Therapy advance online publication, 27 June 2013; doi:10.1038/gt.2013.36.

Relevance of biallelic versus monoallelic TNFRSF13B mutations in distinguishing disease-causing from risk-increasing TNFRSF13B variants in antibody deficiency syndromes.

TNFRSF13B encodes transmembrane activator and calcium modulator and cyclophilin ligand interactor (TACI), a B cell- specific tumor necrosis factor (TNF) receptor superfamily member. Both biallelic and monoallelic TNFRSF13B mutations were identified in patients with common variable immunodeficiency disorders. The genetic complexity and variable clinical presentation of TACI deficiency prompted us to evaluate the genetic, immunologic, and clinical condition in 50 individuals with TNFRSF13B alterations, following screening of 564 unrelated patients with hypogammaglobulinemia. We identified 13 new sequence variants. The most frequent TNFRSF13B variants (C104R and A181E; n=39; 6.9%) were also present in a heterozygous state in 2% of 675 controls. All patients with biallelic mutations had hypogammaglobulinemia and nearly all showed impaired binding to a proliferation-inducing ligand (APRIL). However, the majority (n=41; 82%) of the pa-tients carried monoallelic changes in TNFRSF13B. Presence of a heterozygous mutation was associated with antibody deficiency (P< .001, relative risk 3.6). Heterozygosity for the most common mutation, C104R, was associated with disease (P< .001, relative risk 4.2). Furthermore, heterozygosity for C104R was associated with low numbers of IgD(-)CD27(+) B cells (P= .019), benign lymphoproliferation (P< .001), and autoimmune complications (P= .001). These associations indicate that C104R heterozygosity increases the risk for common variable immunodeficiency disorders and influences clinical presentation.

Analysis of stop-gain and frameshift variants in human innate immunity genes.

Loss-of-function variants in innate immunity genes are associated with Mendelian disorders in the form of primary immunodeficiencies. Recent resequencing projects report that stop-gains and frameshifts are collectively prevalent in humans and could be responsible for some of the inter-individual variability in innate immune response. Current computational approaches evaluating loss-of-function in genes carrying these variants rely on gene-level characteristics such as evolutionary conservation and functional redundancy across the genome. However, innate immunity genes represent a particular case because they are more likely to be under positive selection and duplicated. To create a ranking of severity that would be applicable to innate immunity genes we evaluated 17,764 stop-gain and 13,915 frameshift variants from the NHLBI Exome Sequencing Project and 1,000 Genomes Project. Sequence-based features such as loss of functional domains, isoform-specific truncation and nonsense-mediated decay were found to correlate with variant allele frequency and validated with gene expression data. We integrated these features in a Bayesian classification scheme and benchmarked its use in predicting pathogenic variants against Online Mendelian Inheritance in Man (OMIM) disease stop-gains and frameshifts. The classification scheme was applied in the assessment of 335 stop-gains and 236 frameshifts affecting 227 interferon-stimulated genes. The sequence-based score ranks variants in innate immunity genes according to their potential to cause disease, and complements existing gene-based pathogenicity scores. Specifically, the sequence-based score improves measurement of functional gene impairment, discriminates across different variants in a given gene and appears particularly useful for analysis of less conserved genes.