Exome-wide analysis of rare coding variation identifies novel associations with COPD and airflow limitation in MOCS3, IFIT3 and SERPINA12

Funding: British Women’s Heart and Health Study is funded by the Department of Health grant no. 90049 and the British Heart Foundation grant no. PG/09/022. British Regional Heart Study is supported by the British Heart Foundation (grant RG/ 13/16/30528). CB (COPDBEAT) received funding from the Medical Research Council UK (grant no. G0601369), CB (COPDBEAT) and AJW (UKCOPD) were supported by the National Institute for Health Research (NIHR Leicester Biomedical Research Unit). MB (COPDBEAT) received funding from the NIHR (grant no. PDF-2013-06-052). Hertfordshire Cohort Study received support from the Medical Research Council, Arthritis Research UK, the International Osteoporosis Foundation and the British Heart Foundation; NIHR Biomedical Research Centre in Nutrition, University of Southampton; NIHR Musculoskeletal Biomedical Research Unit, University of Oxford. Generation Scotland: Scottish Family Health Study is funded by the Chief Scientist Office, Scottish Government Health Directorates, grant number CZD/16/6 and the Scottish Funding Council grant HR03006. EU COPD Gene Scan is funded by the European Union, grant no. QLG1-CT-2001-01012. English Longitudinal Study of Aging is funded by the Institute of Aging, NIH grant No. AG1764406S1. GoDARTs is funded by the Wellcome Trust grants 072960, 084726 and 104970. MDT has been supported by MRC fellowship G0902313. UK Biobank Lung Exome Variant Evaluation study was funded by a Medical Research Council strategic award to MDT, IPH, DPS and LVW (MC_PC_12010)

Rare coding variants in the phospholipase D3 gene confer risk for Alzheimer's disease.

Genome-wide association studies (GWAS) have identified several risk variants for late-onset Alzheimer's disease (LOAD)1, 2. These common variants have replicable but small effects on LOAD risk and generally do not have obvious functional effects. Low-frequency coding variants, not detected by GWAS, are predicted to include functional variants with larger effects on risk. To identify low-frequency coding variants with large effects on LOAD risk, we carried out whole-exome sequencing (WES) in 14 large LOAD families and follow-up analyses of the candidate variants in several large LOAD case–control data sets. A rare variant in PLD3 (phospholipase D3; Val232Met) segregated with disease status in two independent families and doubled risk for Alzheimer’s disease in seven independent case–control series with a total of more than 11,000 cases and controls of European descent. Gene-based burden analyses in 4,387 cases and controls of European descent and 302 African American cases and controls, with complete sequence data for PLD3, reveal that several variants in this gene increase risk for Alzheimer’s disease in both populations. PLD3 is highly expressed in brain regions that are vulnerable to Alzheimer’s disease pathology, including hippocampus and cortex, and is expressed at significantly lower levels in neurons from Alzheimer’s disease brains compared to control brains. Overexpression of PLD3 leads to a significant decrease in intracellular amyloid-β precursor protein (APP) and extracellular Aβ42 and Aβ40 (the 42- and 40-residue isoforms of the amyloid-β peptide), and knockdown of PLD3 leads to a significant increase in extracellular Aβ42 and Aβ40. Together, our genetic and functional data indicate that carriers of PLD3 coding variants have a twofold increased risk for LOAD and that PLD3 influences APP processing. This study provides an example of how densely affected families may help to identify rare variants with large effects on risk for disease or other complex traits.

Comprehensive analysis of RET common and rare variant patientsin a series of Spanish Hirschsprung patients confirms a synergistic effect of both kinds of events

Background. RET is the major gene associated to Hirschsprung disease (HSCR) with differential contributions of its rare and common, coding and noncoding mutations to the multifactorial nature of this pathology. In the present study, we have performed a comprehensive study of our HSCR series evaluating the involvement of both RET rare variants (RVs) and common variants (CVs) in the context of the disease. Methods. RET mutational screening was performed by dHPLC and direct sequencing for the identification of RVs. In addition Taqman technology was applied for the genotyping of 3 RET CVs previously associated to HSCR, including a variant lying in an enhancer domain within RET intron 1 (rs2435357). Statistical analyses were performed using the SPSS v.17.0 to analyze the distribution of the variants. Results. Our results confirm the strongest association to HSCR for the "enhancer" variant, and demonstrate a significantly higher impact of it in male versus female patients. Integration of the RET RVs and CVs analysis showed that in 91.66% of cases with both kinds of mutational events, the enhancer allele is in trans with the allele bearing the RET RV. Conclusions. A gender effect exists on both the transmission and distribution of rare coding and common HSCR causing mutations. In addition, these RET CVs and RVs seem to act in a synergistic way leading to HSCR phenotype.

Allele and genotype frequencies of polymorphic cytochromes P4502D6, 2C19 and 2E1 in Aborigines from Western Australia

The polymorphisms of the important xenobiotic metabolizing enzymes CYP2D6, CYP2C19 and CYP2E1 have been studied extensively in a large number of populations and show significant heterogeneity in the frequency of different alleles/genotypes and in the prevalence of the extensive and poor metabolizer phenotypes, Understanding of inter-ethnic differences in genotypes is important in prediction of either beneficial or adverse effects from therapeutic agents and other xenobiotics. Since no data were available for Australian Aborigines, we investigated the frequencies of alleles and genotypes for CYP2D6, CYP2C19 and CYP2E1 in a population living in the far north of Western Australia. Because of its geographical isolation, this population can serve as a model to study the impact of evolutionary forces on the distribution of different alleles for xenobiotic metabolizing enzymes. Twelve CYP2D6 alleles were analysed, The wild-type allele *1 was the most frequent (85.8%) and the non-functional alleles (*4, *5, *16) had an overall frequency of less than 10%. Only one subject (0.4%) was a poor metabolizer for CYP2D6 because of the genotype *5/*5, For CYP2C19, the frequencies of the *1 (wild-type) and the non-functional (*2 and *3) alleles were 50.2%, 35.5% and 14.3%, respectively. The combined CYP2C19 genotypes (*2/*2, *2/*3 or *3/*3) correspond to a predicted frequency of 25.6% for the CYP2C19 poor metabolizer phenotype, For CYP2E1, only one subject had the rare c2 allele giving an overall allele frequency of 0.2%. For CYP2D6 and CYP2C19, allele frequencies and predicted phenotypes differed significantly from those for Caucasians but were similar to those for Orientals indicating a close relationship to East Asian populations. Differences between Aborigines and Orientals in allele frequencies for CYP2D6*10 and CYP2E1 c2 may have arisen through natural selection, or genetic drift, respectively, Pharmacogenetics 11:69-76 (C) 2001 Lippincott Williams & Wilkins.

Evolutionary and convergence stability for continuous phenotypes in finite populations derived from two-allele models.

The evolution of a quantitative phenotype is often envisioned as a trait substitution sequence where mutant alleles repeatedly replace resident ones. In infinite populations, the invasion fitness of a mutant in this two-allele representation of the evolutionary process is used to characterize features about long-term phenotypic evolution, such as singular points, convergence stability (established from first-order effects of selection), branching points, and evolutionary stability (established from second-order effects of selection). Here, we try to characterize long-term phenotypic evolution in finite populations from this two-allele representation of the evolutionary process. We construct a stochastic model describing evolutionary dynamics at non-rare mutant allele frequency. We then derive stability conditions based on stationary average mutant frequencies in the presence of vanishing mutation rates. We find that the second-order stability condition obtained from second-order effects of selection is identical to convergence stability. Thus, in two-allele systems in finite populations, convergence stability is enough to characterize long-term evolution under the trait substitution sequence assumption. We perform individual-based simulations to confirm our analytic results.

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.

IL28B polymorphisms leading to poor response to treatment are associated with low necroinflammatory activity and slow fibrosis progression in HCV genotype non-1-infected patients

Background and Aims: Genetic polymorphisms near IL28Bhave been associated with spontaneous and treatment-inducedclearance of hepatitis C virus (HCV). This is believed to proceed viathe appropriate activation of innate and adaptive immune responsestargeting infected hepatocytes. Intrahepatic inflammation is thereflection of the host cell immune response, but its relationshipwith IL28B polymorphisms has yet to be fully appreciated.Methods: We analyzed the association of IL28B polymorphismswith Metavir activity (≥1) and fibrosis scores (≥2) in 1114 HCVinfectedCaucasian patients enrolled in the Swiss Hepatitis C CohortStudy (629, 127, 268 and 110 infected with HCV genotype 1, 2, 3and 4, respectively). In a subgroup of 915 patients with an estimateddate of infection, the association between IL28B polymorphismsand fibrosis progression rate (FPR > median) was assessed. Singlenucleotide polymorphisms (SNPs) of interest were extracted froma dataset generated in a genome-wide association study and/orgenotyped by TaqMan assay. Associations of alleles with differentdegrees of activity and fibrosis were evaluated using an additivemodel of inheritance by multivariate logistic regression, accountingfor all relevant covariates.Results: The rare G allele at marker rs8099917 was associated withlower activity (P = 0.008) and fibrosis (P = 0.01), as well as slower FPR(P = 0.02). Most striking associations were observed among patientsinfected with non-1 genotypes (P = 0.002 for activity, P = 0.002 forfibrosis and P = 0.005 for FPR). In genotype 1-infected patients, theassociation with activity was observed only in the recessive model(P = 0.04), whereas other associations were not significant (P = 0.7for fibrosis and P = 0.4 for FPR).Conclusions: In chronic hepatitis C, IL28B polymorphisms linkedwith a poor virological response to therapy are also associated withreduced intrahepatic necroinflammation and slower liver diseaseprogression. These observations underscore the role played by thehost immune response in clearing HCV, especially in patients withHCV genotypes non-1.

Apolipoprotein CIII polymorphism and triglyceride levels of a Japanese population living in Southern Brazil

Apolipoprotein CIII (apo-CIII) participates in the regulation of triglyceride-rich lipoprotein metabolism. Several polymorphic sites have been detected within and around the apo-CIII gene. Here, we examined the relationship between apo-CIII SstI polymorphism (CC, CG, GG genotypes) and plasma triglyceride (TG) levels in a group of 159 Japanese individuals living in Southern Brazil. The sample was divided into a group of Japanese descendants (N = 51) with high TG (HTG; >200 mg/dL) and a group of Japanese descendants (N = 108) with normal TG (NTG; <200 mg/dL). TG and total cholesterol levels were analyzed by an enzymatic method using the Labtest-Diagnostic kit and high- and low-density lipoproteins by a direct method using the Labtest-Diagnostic kit and DiaSys Diagnostic System International kit, respectively. A 428-bp sequence of apo-CIII gene was amplified using oligonucleotide primers 5' GGT GAC CGA TGG CTT CAG TTC CCT GA 3' and 5' CAG AAG GTG GAT AGA GCG CTG GCC T 3'. The PCR products were digested with a restriction endonuclease SstI. Rare G allele was highly prevalent in our study population (0.416) compared to Caucasians (0.00-0.11). G allele was almost two times more prevalent in the HTG group compared to the NTG group (P < 0.001). The genotype distribution was consistent with the Hardy-Weinberg equilibrium. There was a significant association between rare G allele and HTG in Japanese individuals living in Southern Brazil as indicated by one-way ANOVA, P < 0.05.

Pleiotropic effects of TCF7L2 gene variants and its modulation in the metabolic syndrome: from the LIPGENE study

Aims/hypothesis: Variants of the TCF7L2 gene predict the development of type 2 diabetes mellitus (T2DM). We investigated the associations between gene variants of TCF7L2 and clinical features of the metabolic syndrome (MetS) (an entity often preceeding T2DM), and their interaction with non-genetic factors, including plasma saturated fatty acids (SFA) concentration and insulin resistance (IR). Methods: Fasting lipid profiles, insulin sensitivity, insulin secretion, anthropometrics, blood pressure and 10 gene variations of the TCF7L2 gene were determined in 450 subjects with MetS. Results: Several single nucleotide polymorphisms (SNP) showed phenotypic associations independent of SFA or IR. Carriers of the rare T allele of rs7903146, and of three other SNPs in linkage disequilibrium with rs7903146, had lower blood pressure and insulin secretion. High IR and the presence of the T-allele of rs7903146 acted synergistically to define those with reduced insulin secretion. Carriers of the minor allele of rs290481 exhibited an altered lipid profile, with increased plasma levels of apolipoprotein B, non-esterified fatty acids, cholesterol and apolipoprotein B in triglyceride rich lipoproteins, and LDL cholesterol. Carriers of the minor allele of rs11196224 that had higher plasma SFA levels showed elevated procoagulant/proinflammatory biomarkers, impaired insulin secretion and increased IR, whereas carriers of the minor allele of rs17685538 with high plasma SFA levels exhibited higher blood pressure. Conclusions/interpretation: SNP in the TCF7L2 gene are associated with differences in insulin secretion, blood pressure, blood lipids and coagulation in MetS patients, and may be modulated by SFA in plasma or IR.

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.

New transposons to generate GFP protein fusions in Candida albicans

Transposon elements are important tools for gene function analysis, for example they can be used to easily create genome-wide collections of insertion mutants. Transposons may also carry sequences coding for an epitope or fluorescent marker useful for protein expression and localization analysis. We have developed three new Tn5-based transposons that incorporate a GFP (green fluorescent protein) coding sequence to generate fusion proteins in the important fungal pathogen Candida albicans. Each transposon also contains the URA3 and Kan(R) genes for yeast and bacterial selection, respectively. After in vitro transposition, the insertional allele is transferred to the chromosomal locus by homologous recombination. Transposons Tn5-CaGFP and Tn5-CaGFP-URA3:FLIP can generate C-terminal truncated GFP fusions. A URA3 flipper recycling cassette was incorporated into the transposon Th5-CaGFP-UFRA3:FLIP. After the induction of Flip recombinase to excise the marker, the heterozygous strain is transformed again in order to obtain a GFP-tagged homozygous strains. In the Tn5-CaGFP-FL transposon the markers are flanked by a rare-cutting enzyme. After in vitro transposition into a plasmid-borne target gene, the markers are eliminated by restriction digestion and religation, resulting in a construct coding for full-length GFP-fusion proteins. This transposon can generate plasmid libraries of GFP insertions in proteins where N- or C-terminal tagging may alter localization. We tested our transposon system by mutagenizing the essential septin CDC3 gene. The results indicate that the Cdc3 C-terminal extension is important for correct septin filament assembly. The transposons described here provide a new system to obtain global gene expression and protein localization data in C. albicans. (c) 2008 Elsevier B.V. All rights reserved.

eRF3a/GSPT1 12-GGC allele increases the susceptibility for breast cancer development

It is now widely recognized that translation factors are involved in cancer development and that components of the translation machinery that are deregulated in cancer cells may become targets for cancer therapy. The eukaryotic Release Factor 3 (eRF3) is a GTPase that associates with eRF1 in a complex that mediates translation termination. eRF3a/GSPT1 first exon contains a (GGC)n expansion coding for proteins with different N-terminal extremities. Herein we show that the longer allele (12-GGC) is present in 5.1% (7/137) of the breast cancer patients analysed and is absent in the control population (0/135), corresponding to an increased risk for cancer development, as revealed by Odds Ratio analysis. mRNA quantification suggests that patients with the 12-GGC allele overexpress eRF3a/GSPT1 in tumor tissues relative to the normal adjacent tissues. However, using an in vivo assay for translation termination in HEK293 cells, we do not detect any difference in the activity of the eRF3a proteins encoded by the various eRF3a/GSPT1 alleles. Although the connection between the presence of eRF3a/GSPT1 12-GGC allele and tumorigenesis is still unknown, our data suggest that the presence of the 12-GGC allele provides a potential novel risk marker for various types of cancer.

Trypanosoma cruzi ribosomal protein S4: characterization of its coding locus, analysis of transcripts, and antigenicity of the protein

Two allelic genomic fragments containing ribosomal protein S4 encoding genes (rpS4) from Trypanosoma cruzi (CL-Brener strain) were isolated and characterized. One allele comprises two complete tandem repeats of a sequence encoding an rpS4 gene. In the other, only one rpS4 gene is found. Sequence comparison to the accessed data in the genome project database reveals that our two-copy allele corresponds to a variant haplotype. However, the deduced aminoacid sequence of all the gene copies is identical. The rpS4 transcripts processing sites were determined by comparison of genomic sequences with published cDNA data. The obtained sequence data demonstrates that rpS4 genes are expressed in epimastigotes, amastigotes, and trypomastigotes. A recombinant version of rpS4 was found to be an antigenic: it was recognized by 62.5% of the individuals with positive serology for T. cruzi and by 93.3% of patients with proven chronic chagasic disease.

Comprehensive analysis of NRG1 common and rare variants in Hirschsprung patients.

Hirschsprung disease (HSCR, OMIM 142623) is a developmental disorder characterized by the absence of ganglion cells along variable lengths of the distal gastrointestinal tract, which results in tonic contraction of the aganglionic gut segment and functional intestinal obstruction. The RET proto-oncogene is the major gene for HSCR with differential contributions of its rare and common, coding and noncoding mutations to the multifactorial nature of this pathology. Many other genes have been described to be associated with the pathology, as NRG1 gene (8p12), encoding neuregulin 1, which is implicated in the development of the enteric nervous system (ENS), and seems to contribute by both common and rare variants. Here we present the results of a comprehensive analysis of the NRG1 gene in the context of the disease in a series of 207 Spanish HSCR patients, by both mutational screening of its coding sequence and evaluation of 3 common tag SNPs as low penetrance susceptibility factors, finding some potentially damaging variants which we have functionally characterized. All of them were found to be associated with a significant reduction of the normal NRG1 protein levels. The fact that those mutations analyzed alter NRG1 protein would suggest that they would be related with HSCR disease not only in Chinese but also in a Caucasian population, which reinforces the implication of NRG1 gene in this pathology.

Truncating mutations in the last exon of NOTCH2 cause a rare skeletal disorder with osteoporosis.

Hajdu-Cheney syndrome is a rare autosomal dominant skeletal disorder with facial anomalies, osteoporosis and acro-osteolysis. We sequenced the exomes of six unrelated individuals with this syndrome and identified heterozygous nonsense and frameshift mutations in NOTCH2 in five of them. All mutations cluster to the last coding exon of the gene, suggesting that the mutant mRNA products escape nonsense-mediated decay and that the resulting truncated NOTCH2 proteins act in a gain-of-function manner.