Proteomic and Metabolomic Analyses of Mitochondrial Complex I-deficient Mouse Model Generated by Spontaneous B2 Short Interspersed Nuclear Element (SINE) Insertion into NADH Dehydrogenase (Ubiquinone) Fe-S Protein 4 (Ndufs4) Gene.

Eukaryotic cells generate energy in the form of ATP, through a network of mitochondrial complexes and electron carriers known as the oxidative phosphorylation system. In mammals, mitochondrial complex I (CI) is the largest component of this system, comprising 45 different subunits encoded by mitochondrial and nuclear DNA. Humans diagnosed with mutations in the gene NDUFS4, encoding a nuclear DNA-encoded subunit of CI (NADH dehydrogenase ubiquinone Fe-S protein 4), typically suffer from Leigh syndrome, a neurodegenerative disease with onset in infancy or early childhood. Mitochondria from NDUFS4 patients usually lack detectable NDUFS4 protein and show a CI stability/assembly defect. Here, we describe a recessive mouse phenotype caused by the insertion of a transposable element into Ndufs4, identified by a novel combined linkage and expression analysis. Designated Ndufs4(fky), the mutation leads to aberrant transcript splicing and absence of NDUFS4 protein in all tissues tested of homozygous mice. Physical and behavioral symptoms displayed by Ndufs4(fky/fky) mice include temporary fur loss, growth retardation, unsteady gait, and abnormal body posture when suspended by the tail. Analysis of CI in Ndufs4(fky/fky) mice using blue native PAGE revealed the presence of a faster migrating crippled complex. This crippled CI was shown to lack subunits of the "N assembly module", which contains the NADH binding site, but contained two assembly factors not present in intact CI. Metabolomic analysis of the blood by tandem mass spectrometry showed increased hydroxyacylcarnitine species, implying that the CI defect leads to an imbalanced NADH/NAD(+) ratio that inhibits mitochondrial fatty acid β-oxidation.

IFN stimulated gene expression in the liver is a better predictor of treatment response in chronic hepatitis c than the IL28b (IFN lambda 3) genotype

Background: Therapy of chronic hepatitis C (CHC) with pegIFNa/ribavirin achieves sustained virologic response (SVR) in ~55%. Pre-activation of the endogenous interferon system in the liver is associated non-response (NR). Recently, genome-wide association studies described associations of allelic variants near the IL28B (IFNλ3) gene with treatment response and with spontaneous clearance of the virus. We investigated if the IL28B genotype determines the constitutive expression of IFN stimulated genes (ISGs) in the liver of patients with CHC. Methods: We genotyped 93 patients with CHC for 3 IL28B single nucleotide polymorphisms (SNPs, rs12979860, rs8099917, rs12980275), extracted RNA from their liver biopsies and quantified the expression of IL28B and of 8 previously identified classifier genes which discriminate between SVR and NR (IFI44L, RSAD2, ISG15, IFI22, LAMP3, OAS3, LGALS3BP and HTATIP2). Decision tree ensembles in the form of a random forest classifier were used to calculate the relative predictive power of these different variables in a multivariate analysis. Results: The minor IL28B allele (bad risk for treatment response) was significantly associated with increased expression of ISGs, and, unexpectedly, with decreased expression of IL28B. Stratification of the patients into SVR and NR revealed that ISG expression was conditionally independent from the IL28B genotype, i.e. there was an increased expression of ISGs in NR compared to SVR irrespective of the IL28B genotype. The random forest feature score (RFFS) identified IFI27 (RFFS = 2.93), RSAD2 (1.88) and HTATIP2 (1.50) expression and the HCV genotype (1.62) as the strongest predictors of treatment response. ROC curves of the IL28B SNPs showed an AUC of 0.66 with an error rate (ERR) of 0.38. A classifier with the 3 best classifying genes showed an excellent test performance with an AUC of 0.94 and ERR of 0.15. The addition of IL28B genotype information did not improve the predictive power of the 3-gene classifier. Conclusions: IL28B genotype and hepatic ISG expression are conditionally independent predictors of treatment response in CHC. There is no direct link between altered IFNλ3 expression and pre-activation of the endogenous system in the liver. Hepatic ISG expression is by far the better predictor for treatment response than IL28B genotype.

Genome-wide association study of metabolic traits reveals novel gene-metabolite-disease links.

Metabolic traits are molecular phenotypes that can drive clinical phenotypes and may predict disease progression. Here, we report results from a metabolome- and genome-wide association study on (1)H-NMR urine metabolic profiles. The study was conducted within an untargeted approach, employing a novel method for compound identification. From our discovery cohort of 835 Caucasian individuals who participated in the CoLaus study, we identified 139 suggestively significant (P<5×10(-8)) and independent associations between single nucleotide polymorphisms (SNP) and metabolome features. Fifty-six of these associations replicated in the TasteSensomics cohort, comprising 601 individuals from São Paulo of vastly diverse ethnic background. They correspond to eleven gene-metabolite associations, six of which had been previously identified in the urine metabolome and three in the serum metabolome. Our key novel findings are the associations of two SNPs with NMR spectral signatures pointing to fucose (rs492602, P = 6.9×10(-44)) and lysine (rs8101881, P = 1.2×10(-33)), respectively. Fine-mapping of the first locus pinpointed the FUT2 gene, which encodes a fucosyltransferase enzyme and has previously been associated with Crohn's disease. This implicates fucose as a potential prognostic disease marker, for which there is already published evidence from a mouse model. The second SNP lies within the SLC7A9 gene, rare mutations of which have been linked to severe kidney damage. The replication of previous associations and our new discoveries demonstrate the potential of untargeted metabolomics GWAS to robustly identify molecular disease markers.

Single-cell gene-expression profiling reveals qualitatively distinct CD8 T cells elicited by different gene-based vaccines.

CD8 T cells play a key role in mediating protective immunity against selected pathogens after vaccination. Understanding the mechanism of this protection is dependent upon definition of the heterogeneity and complexity of cellular immune responses generated by different vaccines. Here, we identify previously unrecognized subsets of CD8 T cells based upon analysis of gene-expression patterns within single cells and show that they are differentially induced by different vaccines. Three prime-boost vector combinations encoding HIV Env stimulated antigen-specific CD8 T-cell populations of similar magnitude, phenotype, and functionality. Remarkably, however, analysis of single-cell gene-expression profiles enabled discrimination of a majority of central memory (CM) and effector memory (EM) CD8 T cells elicited by the three vaccines. Subsets of T cells could be defined based on their expression of Eomes, Cxcr3, and Ccr7, or Klrk1, Klrg1, and Ccr5 in CM and EM cells, respectively. Of CM cells elicited by DNA prime-recombinant adenoviral (rAd) boost vectors, 67% were Eomes(-) Ccr7(+) Cxcr3(-), in contrast to only 7% and 2% stimulated by rAd5-rAd5 or rAd-LCMV, respectively. Of EM cells elicited by DNA-rAd, 74% were Klrk1(-) Klrg1(-)Ccr5(-) compared with only 26% and 20% for rAd5-rAd5 or rAd5-LCMV. Definition by single-cell gene profiling of specific CM and EM CD8 T-cell subsets that are differentially induced by different gene-based vaccines will facilitate the design and evaluation of vaccines, as well as enable our understanding of mechanisms of protective immunity.

Promoter rearrangements cause species-specific hepatic regulation of the glyoxylate reductase/hydroxypyruvate reductase gene by the peroxisome proliferator-activated receptor alpha.

In liver, the glyoxylate cycle contributes to two metabolic functions, urea and glucose synthesis. One of the key enzymes in this pathway is glyoxylate reductase/hydroxypyruvate reductase (GRHPR) whose dysfunction in human causes primary hyperoxaluria type 2, a disease resulting in oxalate accumulation and formation of kidney stones. In this study, we provide evidence for a transcriptional regulation by the peroxisome proliferator-activated receptor alpha (PPARalpha) of the mouse GRHPR gene in liver. Mice fed with a PPARalpha ligand or in which PPARalpha activity is enhanced by fasting increase their GRHPR gene expression via a peroxisome proliferator response element located in the promoter region of the gene. Consistent with these observations, mice deficient in PPARalpha present higher plasma levels of oxalate in comparison with their wild type counterparts. As expected, the administration of a PPARalpha ligand (Wy-14,643) reduces the plasma oxalate levels. Surprisingly, this effect is also observed in null mice, suggesting a PPARalpha-independent action of the compound. Despite a high degree of similarity between the transcribed region of the human and mouse GRHPR gene, the human promoter has been dramatically reorganized, which has resulted in a loss of PPARalpha regulation. Overall, these data indicate a species-specific regulation by PPARalpha of GRHPR, a key gene of the glyoxylate cycle.

Clinical and genetic investigation of a large Tunisian family with complete achromatopsia: identification of a new nonsense mutation in GNAT2 gene.

Complete achromatopsia is a rare autosomal recessive disease associated with CNGA3, CNGB3, GNAT2 and PDE6C mutations. This retinal disorder is characterized by complete loss of color discrimination due to the absence or alteration of the cones function. The purpose of the present study was the clinical and the genetic characterization of achromatopsia in a large consanguineous Tunisian family. Ophthalmic evaluation included a full clinical examination, color vision testing and electroretinography. Linkage analysis using microsatellite markers flanking CNGA3, CNGB3, GNAT2 and PDE6C genes was performed. Mutations were screened by direct sequencing. A total of 12 individuals were diagnosed with congenital complete achromatopsia. They are members of six nuclear consanguineous families belonging to the same large consanguineous family. Linkage analysis revealed linkage to GNAT2. Mutational screening of GNAT2 revealed three intronic variations c.119-69G>C, c.161+66A>T and c.875-31G>C that co-segregated with a novel mutation p.R313X. An identical GNAT2 haplotype segregating with this mutation was identified, indicating a founder mutation. All patients were homozygous for the p.R313X mutation. This is the first report of the clinical and genetic investigation of complete achromatopsia in North Africa and the largest family with recessive achromatopsia involving GNAT2; thus, providing a unique opportunity for genotype-phenotype correlation for this extremely rare condition.

Combined clonotyping and gene signatures of individual tumor-reactive CD8 T-cells define their functional relevance following peptide vaccination in melanoma patients

Novel cancer vaccines are capableto efficiently induce and boost humantumor antigen specific T-cells. However,the properties of these CD8T-cells are only partially characterized.For in depth investigation ofT-cells following Melan-A/MART-1peptide vaccination in melanoma patients,we conducted a detailed prospectivestudy at the single cell level.We first sorted individual human naiveand effector CD8 T-cells from peripheralblood by flow cytometry, andtested a modified RT-PCR protocolincluding a global amplification ofexpressed mRNAs to obtain sufficientcDNAfromsingle cells.We successfullydetected the expression ofseveral specific genes of interest evendown to 106-fold dilution (equivalentto 10-5 cell). We then analyzed tumor-specific effector memory (EM)CD8T-cell subpopulations ex vivo, assingle cells from vaccinated melanomapatients. To elucidate the hallmarksof effective immunity the genesignatures were defined by a panel ofgenes related to effector functions(e.g. IFN-, granzyme B, perforin),and individual clonotypes were identifiedaccording to the expression ofdistinct T-cell receptors (TCR). Usingthis novel single cell analysis approach,we observed that T-cell differentiationis clonotype dependent,with a progressive restriction in TCRBV clonotype diversity from EMCD28pos to EMCD28neg subsets. However,the effector function gene imprintingis clonotype-independent,but dependent on differentiation,since it correlates with the subset oforigin (EMCD28pos or EMCD28neg). We also conducted a detailedcomparative analysis after vaccinationwith natural vs. analog Melan-Apeptide. We found that the peptideused for vaccination determines thefunctional outcome of individualT-cell clonotypes, with native peptideinducing more potent effector functions.Yet, selective clonotypic expansionwith differentiation was preservedregardless of the peptide usedfor vaccination. In summary, the exvivo single cell RT-PCR approach ishighly sensitive and efficient, andrepresents a reliable and powerfultool to refine our current view of molecularprocesses taking place duringT-cell differentiation.

Variation in novel exons (RACEfrags) of the MECP2 gene in Rett syndrome patients and controls.

The study of transcription using genomic tiling arrays has lead to the identification of numerous additional exons. One example is the MECP2 gene on the X chromosome; using 5'RACE and RT-PCR in human tissues and cell lines, we have found more than 70 novel exons (RACEfrags) connecting to at least one annotated exon.. We sequenced all MECP2-connected exons and flanking sequences in 3 groups: 46 patients with the Rett syndrome and without mutations in the currently annotated exons of the MECP2 and CDKL5 genes; 32 patients with the Rett syndrome and identified mutations in the MECP2 gene; 100 control individuals from the same geoethnic group. Approximately 13 kb were sequenced per sample, (2.4 Mb of DNA resequencing). A total of 75 individuals had novel rare variants (mostly private variants) but no statistically significant difference was found among the 3 groups. These results suggest that variants in the newly discovered exons may not contribute to Rett syndrome. Interestingly however, there are about twice more variants in the novel exons than in the flanking sequences (44 vs. 21 for approximately 1.3 Mb sequenced for each class of sequences, p=0.0025). Thus the evolutionary forces that shape these novel exons may be different than those of neighboring sequences.

EGFR signalling as a negative regulator of Notch1 gene transcription and function in proliferating keratinocytes and cancer.

The Notch1 gene has an important role in mammalian cell-fate decision and tumorigenesis. Upstream control mechanisms for transcription of this gene are still poorly understood. In a chemical genetics screen for small molecule activators of Notch signalling, we identified epidermal growth factor receptor (EGFR) as a key negative regulator of Notch1 gene expression in primary human keratinocytes, intact epidermis and skin squamous cell carcinomas (SCCs). The underlying mechanism for negative control of the Notch1 gene in human cells, as well as in a mouse model of EGFR-dependent skin carcinogenesis, involves transcriptional suppression of p53 by the EGFR effector c-Jun. Suppression of Notch signalling in cancer cells counteracts the differentiation-inducing effects of EGFR inhibitors while, at the same time, synergizing with these compounds in induction of apoptosis. Thus, our data reveal a key role of EGFR signalling in the negative regulation of Notch1 gene transcription, of potential relevance for combinatory approaches for cancer therapy.

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.

New gene functions in megakaryopoiesis and platelet formation.

Platelets are the second most abundant cell type in blood and are essential for maintaining haemostasis. Their count and volume are tightly controlled within narrow physiological ranges, but there is only limited understanding of the molecular processes controlling both traits. Here we carried out a high-powered meta-analysis of genome-wide association studies (GWAS) in up to 66,867 individuals of European ancestry, followed by extensive biological and functional assessment. We identified 68 genomic loci reliably associated with platelet count and volume mapping to established and putative novel regulators of megakaryopoiesis and platelet formation. These genes show megakaryocyte-specific gene expression patterns and extensive network connectivity. Using gene silencing in Danio rerio and Drosophila melanogaster, we identified 11 of the genes as novel regulators of blood cell formation. Taken together, our findings advance understanding of novel gene functions controlling fate-determining events during megakaryopoiesis and platelet formation, providing a new example of successful translation of GWAS to function.

CUBN is a gene locus for albuminuria.

Identification of genetic risk factors for albuminuria may alter strategies for early prevention of CKD progression, particularly among patients with diabetes. Little is known about the influence of common genetic variants on albuminuria in both general and diabetic populations. We performed a meta-analysis of data from 63,153 individuals of European ancestry with genotype information from genome-wide association studies (CKDGen Consortium) and from a large candidate gene study (CARe Consortium) to identify susceptibility loci for the quantitative trait urinary albumin-to-creatinine ratio (UACR) and the clinical diagnosis microalbuminuria. We identified an association between a missense variant (I2984V) in the CUBN gene, which encodes cubilin, and both UACR (P = 1.1 × 10(-11)) and microalbuminuria (P = 0.001). We observed similar associations among 6981 African Americans in the CARe Consortium. The associations between this variant and both UACR and microalbuminuria were significant in individuals of European ancestry regardless of diabetes status. Finally, this variant associated with a 41% increased risk for the development of persistent microalbuminuria during 20 years of follow-up among 1304 participants with type 1 diabetes in the prospective DCCT/EDIC Study. In summary, we identified a missense CUBN variant that associates with levels of albuminuria in both the general population and in individuals with diabetes.

Glycogen synthase 2 is a novel target gene of peroxisome proliferator-activated receptors.

Glycogen synthase 2 (Gys-2) is the ratelimiting enzyme in the storage of glycogen in liver and adipose tissue, yet little is known about regulation of Gys-2 transcription. The peroxisome proliferator-activated receptors (PPARs) are transcription factors involved in the regulation of lipid and glucose metabolism and might be hypothesized to govern glycogen synthesis as well. Here, we show that Gys-2 is a direct target gene of PPARalpha, PPARbeta/delta and PPARgamma. Expression of Gys-2 is significantly reduced in adipose tissue of PPARalpha-/-, PPARbeta/delta-/- and PPARgamma+/- mice. Furthermore, synthetic PPARbeta/delta, and gamma agonists markedly up-regulate Gys-2 mRNA and protein expression in mouse 3T3-L1 adipocytes. In liver, PPARalpha deletion leads to decreased glycogen levels in the refed state, which is paralleled by decreased expression of Gys-2 in fasted and refed state. Two putative PPAR response elements (PPREs) were identified in the mouse Gys-2 gene: one in the upstream promoter (DR-1prom) and one in intron 1 (DR-1int). It is shown that DR-1int is the response element for PPARs, while DR-1prom is the response element for Hepatic Nuclear Factor 4 alpha (HNF4alpha). In adipose tissue, which does not express HNF4alpha, DR-1prom is occupied by PPARbeta/delta and PPARgamma, yet binding does not translate into transcriptional activation of Gys-2. Overall, we conclude that mouse Gys-2 is a novel PPAR target gene and that transactivation by PPARs and HNF4alpha is mediated by two distinct response elements.

Effects of (p)ppGpp on the progression of the cell cycle of Caulobacter crescentus.

Bacteria must control the progression of their cell cycle in response to nutrient availability. This regulation can be mediated by guanosine tetra- or pentaphosphate [(p)ppGpp], which are synthesized by enzymes of the RelA/SpoT homologue (Rsh) family, particularly under starvation conditions. Here, we study the effects of (p)ppGpp on the cell cycle of Caulobacter crescentus, an oligotrophic bacterium with a dimorphic life cycle. C. crescentus divides asymmetrically, producing a motile swarmer cell that cannot replicate its chromosome and a sessile stalked cell that is replication competent. The swarmer cell rapidly differentiates into a stalked cell in appropriate conditions. An artificial increase in the levels of (p)ppGpp in nonstarved C. crescentus cells was achieved by expressing a truncated relA gene from Escherichia coli, encoding a constitutively active (p)ppGpp synthetase. By combining single-cell microscopy, flow cytometry approaches, and swarming assays, we show that an increase in the intracellular concentration of (p)ppGpp is sufficient to slow down the swarmer-to-stalked cell differentiation process and to delay the initiation of chromosome replication. We also present evidence that the intracellular levels of two master regulators of the cell cycle of C. crescentus, DnaA and CtrA, are modulated in response to (p)ppGpp accumulation, even in the absence of actual starvation. CtrA proteolysis and DnaA synthesis seem indirectly inhibited by (p)ppGpp accumulation. By extending the life span of the motile nonreproductive swarmer cell and thus promoting dispersal and foraging functions over multiplication under starvation conditions, (p)ppGpp may play a central role in the ecological adaptation of C. crescentus to nutritional stresses.

Mirror extreme BMI phenotypes associated with gene dosage at the chromosome 16p11.2 locus.

Both obesity and being underweight have been associated with increased mortality. Underweight, defined as a body mass index (BMI) ≤ 18.5 kg per m(2) in adults and ≤ -2 standard deviations from the mean in children, is the main sign of a series of heterogeneous clinical conditions including failure to thrive, feeding and eating disorder and/or anorexia nervosa. In contrast to obesity, few genetic variants underlying these clinical conditions have been reported. We previously showed that hemizygosity of a ∼600-kilobase (kb) region on the short arm of chromosome 16 causes a highly penetrant form of obesity that is often associated with hyperphagia and intellectual disabilities. Here we show that the corresponding reciprocal duplication is associated with being underweight. We identified 138 duplication carriers (including 132 novel cases and 108 unrelated carriers) from individuals clinically referred for developmental or intellectual disabilities (DD/ID) or psychiatric disorders, or recruited from population-based cohorts. These carriers show significantly reduced postnatal weight and BMI. Half of the boys younger than five years are underweight with a probable diagnosis of failure to thrive, whereas adult duplication carriers have an 8.3-fold increased risk of being clinically underweight. We observe a trend towards increased severity in males, as well as a depletion of male carriers among non-medically ascertained cases. These features are associated with an unusually high frequency of selective and restrictive eating behaviours and a significant reduction in head circumference. Each of the observed phenotypes is the converse of one reported in carriers of deletions at this locus. The phenotypes correlate with changes in transcript levels for genes mapping within the duplication but not in flanking regions. The reciprocal impact of these 16p11.2 copy-number variants indicates that severe obesity and being underweight could have mirror aetiologies, possibly through contrasting effects on energy balance.