The yeast two-hybrid (Y2H) system is the most widely applied methodology for systematic protein-protein interaction (PPI) screening and the generation of comprehensive interaction networks. We developed a novel Y2H interaction screening procedure using DNA microarrays for high-throughput quantitative PPI detection. Applying a global pooling and selection scheme to a large collection of human open reading frames, proof-of-principle Y2H interaction screens were performed for the human neurodegenerative disease proteins huntingtin and ataxin-1. Using systematic controls for unspecific Y2H results and quantitative benchmarking, we identified and scored a large number of known and novel partner proteins for both huntingtin and ataxin-1. Moreover, we show that this parallelized screening procedure and the global inspection of Y2H interaction data are uniquely suited to define specific PPI patterns and their alteration by disease-causing mutations in huntingtin and ataxin-1. This approach takes advantage of the specificity and flexibility of DNA microarrays and of the existence of solid-related statistical methods for the analysis of DNA microarray data, and allows a quantitative approach toward interaction screens in human and in model organisms.

Alternative splicing, polyadenylation of pre-messenger RNA molecules and differential promoter usage can produce a variety of transcript isoforms whose respective expression levels are regulated in time and space, thus contributing specific biological functions. However, the repertoire of mammalian alternative transcripts and their regulation are still poorly understood. Second-generation sequencing is now opening unprecedented routes to address the analysis of entire transcriptomes. Here, we developed methods that allow the prediction and quantification of alternative isoforms derived solely from exon expression levels in RNA-Seq data. These are based on an explicit statistical model and enable the prediction of alternative isoforms within or between conditions using any known gene annotation, as well as the relative quantification of known transcript structures. Applying these methods to a human RNA-Seq dataset, we validated a significant fraction of the predictions by RT-PCR. Data further showed that these predictions correlated well with information originating from junction reads. A direct comparison with exon arrays indicated improved performances of RNA-Seq over microarrays in the prediction of skipped exons. Altogether, the set of methods presented here comprehensively addresses multiple aspects of alternative isoform analysis. The software is available as an open-source R-package called Solas at http://cmb.molgen.mpg.de/2ndGenerationSequencing/Solas/.

Down syndrome (DS; trisomy 21) is the most common genetic cause of mental retardation in the human population and key molecular networks dysregulated in DS are still unknown. Many different experimental techniques have been applied to analyse the effects of dosage imbalance at the molecular and phenotypical level, however, currently no integrative approach exists that attempts to extract the common information.

The computational prediction of alternative splicing from high-throughput sequencing data is inherently difficult and necessitates robust statistical measures because the differential splicing signal is overlaid by influencing factors such as gene expression differences and simultaneous expression of multiple isoforms amongst others. In this work we describe ARH-seq, a discovery tool for differential splicing in case-control studies that is based on the information-theoretic concept of entropy. ARH-seq works on high-throughput sequencing data and is an extension of the ARH method that was originally developed for exon microarrays. We show that the method has inherent features, such as independence of transcript exon number and independence of differential expression, what makes it particularly suited for detecting alternative splicing events from sequencing data. In order to test and validate our workflow we challenged it with publicly available sequencing data derived from human tissues and conducted a comparison with eight alternative computational methods. In order to judge the performance of the different methods we constructed a benchmark data set of true positive splicing events across different tissues agglomerated from public databases and show that ARH-seq is an accurate, computationally fast and high-performing method for detecting differential splicing events.

Type 2 diabetes (T2D) has a strong genetic component. Most of the gene variants driving the pathogenesis of T2D seem to target pancreatic β-cell function. To identify novel gene variants acting at early stage of the disease, we analyzed whole transcriptome data to identify differential expression (DE) and alternative exon splicing (AS) transcripts in pancreatic islets collected from two metabolically diverse mouse strains at 6 weeks of age after three weeks of high-fat-diet intervention. Our analysis revealed 1218 DE and 436 AS genes in islets from NZO/Hl vs C3HeB/FeJ. Whereas some of the revealed genes present well-established markers for β-cell failure, such as Cd36 or Aldh1a3, we identified numerous DE/AS genes that have not been described in context with β-cell function before. The gene Lgals2, previously associated with human T2D development, was DE as well as AS and localizes in a quantitative trait locus (QTL) for blood glucose on Chr.15 that we reported recently in our N2(NZOxC3H) population. In addition, pathway enrichment analysis of DE and AS genes showed an overlap of only half of the revealed pathways, indicating that DE and AS in large parts influence different pathways in T2D development. PPARG and adipogenesis pathways, two well-established metabolic pathways, were overrepresented for both DE and AS genes, probably as an adaptive mechanism to cope for increased cellular stress. Our results provide guidance for the identification of novel T2D candidate genes and demonstrate the presence of numerous AS transcripts possibly involved in islet function and maintenance of glucose homeostasis.

We analyzed the transcriptional signatures of mouse bone marrow-derived macrophages at different times after infection with promastigotes of the protozoan parasite Leishmania major. Ingenuity Pathway Analysis revealed that the macrophage metabolic pathways including carbohydrate and lipid metabolisms were among the most altered pathways at later time points of infection. Indeed, L. major promastiogtes induced increased mRNA levels of the glucose transporter and almost all of the genes associated with glycolysis and lactate dehydrogenase, suggesting a shift to anaerobic glycolysis. On the other hand, L. major promastigotes enhanced the expression of scavenger receptors involved in the uptake of Low-Density Lipoprotein (LDL), inhibited the expression of genes coding for proteins regulating cholesterol efflux, and induced the synthesis of triacylglycerides. These data suggested that Leishmania infection disturbs cholesterol and triglycerides homeostasis and may lead to cholesterol accumulation and foam cell formation. Using Filipin and Bodipy staining, we showed cholesterol and triglycerides accumulation in infected macrophages. Moreover, Bodipy-positive lipid droplets accumulated in close proximity to parasitophorous vacuoles, suggesting that intracellular L. major may take advantage of these organelles as high-energy substrate sources. While the effect of infection on cholesterol accumulation and lipid droplet formation was independent on parasite development, our data indicate that anaerobic glycolysis is actively induced by L. major during the establishment of infection.

Multiple functional genomics data for complex human diseases have been published and made available by researchers worldwide. The main goal of these studies is the detailed analysis of a particular aspect of the disease. Complementary, meta-analysis approaches try to extract supersets of disease genes and interaction networks by integrating and combining these individual studies using statistical approaches.

To nominate novel disease genes for obesity and type 2 diabetes (T2D), we recently generated two mouse backcross populations of the T2D-susceptible New Zealand Obese (NZO/HI) mouse strain and two genetically different, lean and T2D-resistant strains, 129P2/OlaHsd and C3HeB/FeJ. Comparative linkage analysis of our two female backcross populations identified seven novel body fat-associated quantitative trait loci (QTL). Only the locus Nbw14 (NZO body weight on chromosome 14) showed linkage to obesity-related traits in both backcross populations, indicating that the causal gene variant is likely specific for the NZO strain as NZO allele carriers in both crosses displayed elevated body weight and fat mass. To identify candidate genes for Nbw14, we used a combined approach of gene expression and haplotype analysis to filter for NZO-specific gene variants in gonadal white adipose tissue, defined as the main QTL-target tissue. Only two genes, Arl11 and Sgcg, fulfilled our candidate criteria. In addition, expression QTL analysis revealed cis-signals for both genes within the Nbw14 locus. Moreover, retroviral overexpression of Sgcg in 3T3-L1 adipocytes resulted in increased insulin-stimulated glucose uptake. In humans, mRNA levels of SGCG correlated with body mass index and body fat mass exclusively in diabetic subjects, suggesting that SGCG may present a novel marker for metabolically unhealthy obesity. In conclusion, our comparative-cross analysis could substantially improve the mapping resolution of the obesity locus Nbw14. Future studies will throw light on the mechanism by which Sgcg may protect from the development of obesity.

Type 2 diabetes (T2D) represents a multifactorial metabolic disease with a strong genetic predisposition. Despite elaborate efforts in identifying the genetic variants determining individual susceptibility towards T2D, the majority of genetic factors driving disease development remain poorly understood. With the aim to identify novel T2D risk genes we previously generated an N2 outcross population using the two inbred mouse strains New Zealand obese (NZO) and C3HeB/FeJ (C3H). A linkage study performed in this population led to the identification of the novel T2D-associated quantitative trait locus (QTL) Nbg15 (NZO blood glucose on chromosome 15, Logarithm of odds (LOD) 6.6). In this study we used a combined approach of positional cloning, gene expression analyses and in silico predictions of DNA polymorphism on gene/protein function to dissect the genetic variants linking Nbg15 to the development of T2D. Moreover, we have generated congenic strains that associated the distal sublocus of Nbg15 to mechanisms altering pancreatic beta cell function. In this sublocus, Cbx6, Fam135b and Kdelr3 were nominated as potential causative genes associated with the Nbg15 driven effects. Moreover, a putative mutation in the Kdelr3 gene from NZO was identified, negatively influencing adaptive responses associated with pancreatic beta cell death and induction of endoplasmic reticulum stress. Importantly, knockdown of Kdelr3 in cultured Min6 beta cells altered insulin granules maturation and pro-insulin levels, pointing towards a crucial role of this gene in islets function and T2D susceptibility.