João Fadista

A systems genetics approach identifies genes and pathways for type 2 diabetes in human islets.

Close to 50 genetic loci have been associated with type 2 diabetes (T2D), but they explain only 15% of the heritability. In an attempt to identify additional T2D genes, we analyzed global gene expression in human islets from 63 donors. Using 48 genes located near T2D risk variants, we identified gene coexpression and protein-protein interaction networks that were strongly associated with islet insulin secretion and HbA(1c). We integrated our data to form a rank list of putative T2D genes, of which CHL1, LRFN2, RASGRP1, and PPM1K were validated in INS-1 cells to influence insulin secretion, whereas GPR120 affected apoptosis in islets. Expression variation of the top 20 genes explained 24% of the variance in HbA(1c) with no claim of the direction. The data present a global map of genes associated with islet dysfunction and demonstrate the value of systems genetics for the identification of genes potentially involved in T2D.

Altered DNA methylation and differential expression of genes influencing metabolism and inflammation in adipose tissue from subjects with type 2 diabetes

Genetics, epigenetics, and environment may together affect the susceptibility for type 2 diabetes (T2D). Our aim was to dissect molecular mechanisms underlying T2D using genome-wide expression and DNA methylation data in adipose tissue from monozygotic twin pairs discordant for T2D and independent case-control cohorts. In adipose tissue from diabetic twins, we found decreased expression of genes involved in oxidative phosphorylation; carbohydrate, amino acid, and lipid metabolism; and increased expression of genes involved in inflammation and glycan degradation. The most differentially expressed genes included ELOVL6, GYS2, FADS1, SPP1 (OPN), CCL18, and IL1RN. We replicated these results in adipose tissue from an independent case-control cohort. Several candidate genes for obesity and T2D (e.g., IRS1 and VEGFA) were differentially expressed in discordant twins. We found a heritable contribution to the genome-wide DNA methylation variability in twins. Differences in methylation between monozygotic twin pairs discordant for T2D were subsequently modest. However, 15,627 sites, representing 7,046 genes including PPARG, KCNQ1, TCF7L2, and IRS1, showed differential DNA methylation in adipose tissue from unrelated subjects with T2D compared with control subjects. A total of 1,410 of these sites also showed differential DNA methylation in the twins discordant for T2D. For the differentially methylated sites, the heritability estimate was 0.28. We also identified copy number variants (CNVs) in monozygotic twin pairs discordant for T2D. Taken together, subjects with T2D exhibit multiple transcriptional and epigenetic changes in adipose tissue relevant to the development of the disease.

Global genomic and transcriptomic analysis of human pancreatic islets reveals novel genes influencing glucose metabolism

Significance We provide a comprehensive catalog of novel genetic variants influencing gene expression and metabolic phenotypes in human pancreatic islets. The data also show that the path from genetic variation (SNP) to gene expression is more complex than hitherto often assumed, and that we need to consider that genetic variation can also influence function of a gene by influencing exon usage or splice isoforms (sQTL), allelic imbalance, RNA editing, and expression of noncoding RNAs, which in turn can influence expression of target genes. Genetic variation can modulate gene expression, and thereby phenotypic variation and susceptibility to complex diseases such as type 2 diabetes (T2D). Here we harnessed the potential of DNA and RNA sequencing in human pancreatic islets from 89 deceased donors to identify genes of potential importance in the pathogenesis of T2D. We present a catalog of genetic variants regulating gene expression (eQTL) and exon use (sQTL), including many long noncoding RNAs, which are enriched in known T2D-associated loci. Of 35 eQTL genes, whose expression differed between normoglycemic and hyperglycemic individuals, siRNA of tetraspanin 33 (TSPAN33), 5′-nucleotidase, ecto (NT5E), transmembrane emp24 protein transport domain containing 6 (TMED6), and p21 protein activated kinase 7 (PAK7) in INS1 cells resulted in reduced glucose-stimulated insulin secretion. In addition, we provide a genome-wide catalog of allelic expression imbalance, which is also enriched in known T2D-associated loci. Notably, allelic imbalance in paternally expressed gene 3 (PEG3) was associated with its promoter methylation and T2D status. Finally, RNA editing events were less common in islets than previously suggested in other tissues. Taken together, this study provides new insights into the complexity of gene regulation in human pancreatic islets and better understanding of how genetic variation can influence glucose metabolism.

Copy number variation in the bovine genome

BackgroundCopy number variations (CNVs), which represent a significant source of genetic diversity in mammals, have been shown to be associated with phenotypes of clinical relevance and to be causative of disease. Notwithstanding, little is known about the extent to which CNV contributes to genetic variation in cattle.ResultsWe designed and used a set of NimbleGen CGH arrays that tile across the assayable portion of the cattle genome with approximately 6.3 million probes, at a median probe spacing of 301 bp. This study reports the highest resolution map of copy number variation in the cattle genome, with 304 CNV regions (CNVRs) being identified among the genomes of 20 bovine samples from 4 dairy and beef breeds. The CNVRs identified covered 0.68% (22 Mb) of the genome, and ranged in size from 1.7 to 2,031 kb (median size 16.7 kb). About 20% of the CNVs co-localized with segmental duplications, while 30% encompass genes, of which the majority is involved in environmental response. About 10% of the human orthologous of these genes are associated with human disease susceptibility and, hence, may have important phenotypic consequences.ConclusionsTogether, this analysis provides a useful resource for assessment of the impact of CNVs regarding variation in bovine health and production traits.

A snapshot of CNVs in the pig genome.

Recent studies of mammalian genomes have uncovered the extent of copy number variation (CNV) that contributes to phenotypic diversity, including health and disease status. Here we report a first account of CNVs in the pig genome covering part of the chromosomes 4, 7, 14, and 17 already sequenced and assembled. A custom tiling oligonucleotide array was used with a median probe spacing of 409 bp for screening 12 unrelated Duroc boars that are founders of a large family material. After a strict CNV calling pipeline, 37 copy number variable regions (CNVRs) across all four chromosomes were identified, with five CNVRs overlapping segmental duplications, three overlapping pig unigenes and one overlapping a RefSeq pig mRNA. This CNV snapshot analysis is the first of its kind in the porcine genome and constitutes the basis for a better understanding of porcine phenotypes and genotypes with the prospect of identifying important economic traits.

A Central Role for GRB10 in Regulation of Islet Function in Man

Variants in the growth factor receptor-bound protein 10 (GRB10) gene were in a GWAS meta-analysis associated with reduced glucose-stimulated insulin secretion and increased risk of type 2 diabetes (T2D) if inherited from the father, but inexplicably reduced fasting glucose when inherited from the mother. GRB10 is a negative regulator of insulin signaling and imprinted in a parent-of-origin fashion in different tissues. GRB10 knock-down in human pancreatic islets showed reduced insulin and glucagon secretion, which together with changes in insulin sensitivity may explain the paradoxical reduction of glucose despite a decrease in insulin secretion. Together, these findings suggest that tissue-specific methylation and possibly imprinting of GRB10 can influence glucose metabolism and contribute to T2D pathogenesis. The data also emphasize the need in genetic studies to consider whether risk alleles are inherited from the mother or the father.

TCF7L2 is a master regulator of insulin production and processing

Genome-wide association studies have revealed >60 loci associated with type 2 diabetes (T2D), but the underlying causal variants and functional mechanisms remain largely elusive. Although variants in TCF7L2 confer the strongest risk of T2D among common variants by presumed effects on islet function, the molecular mechanisms are not yet well understood. Using RNA-sequencing, we have identified a TCF7L2-regulated transcriptional network responsible for its effect on insulin secretion in rodent and human pancreatic islets. ISL1 is a primary target of TCF7L2 and regulates proinsulin production and processing via MAFA, PDX1, NKX6.1, PCSK1, PCSK2 and SLC30A8, thereby providing evidence for a coordinated regulation of insulin production and processing. The risk T-allele of rs7903146 was associated with increased TCF7L2 expression, and decreased insulin content and secretion. Using gene expression profiles of 66 human pancreatic islets donors’, we also show that the identified TCF7L2-ISL1 transcriptional network is regulated in a genotype-dependent manner. Taken together, these results demonstrate that not only synthesis of proinsulin is regulated by TCF7L2 but also processing and possibly clearance of proinsulin and insulin. These multiple targets in key pathways may explain why TCF7L2 has emerged as the gene showing one of the strongest associations with T2D.

The diabetes susceptibility gene Clec16a regulates mitophagy

Clec16a has been identified as a disease susceptibility gene for type 1 diabetes, multiple sclerosis, and adrenal dysfunction, but its function is unknown. Here we report that Clec16a is a membrane-associated endosomal protein that interacts with E3 ubiquitin ligase Nrdp1. Loss of Clec16a leads to an increase in the Nrdp1 target Parkin, a master regulator of mitophagy. Islets from mice with pancreas-specific deletion of Clec16a have abnormal mitochondria with reduced oxygen consumption and ATP concentration, both of which are required for normal β cell function. Indeed, pancreatic Clec16a is required for normal glucose-stimulated insulin release. Moreover, patients harboring a diabetogenic SNP in the Clec16a gene have reduced islet Clec16a expression and reduced insulin secretion. Thus, Clec16a controls β cell function and prevents diabetes by controlling mitophagy. This pathway could be targeted for prevention and control of diabetes and may extend to the pathogenesis of other Clec16a- and Parkin-associated diseases.

Global assessment of genomic variation in cattle by genome resequencing and high-throughput genotyping

BackgroundIntegration of genomic variation with phenotypic information is an effective approach for uncovering genotype-phenotype associations. This requires an accurate identification of the different types of variation in individual genomes.ResultsWe report the integration of the whole genome sequence of a single Holstein Friesian bull with data from single nucleotide polymorphism (SNP) and comparative genomic hybridization (CGH) array technologies to determine a comprehensive spectrum of genomic variation. The performance of resequencing SNP detection was assessed by combining SNPs that were identified to be either in identity by descent (IBD) or in copy number variation (CNV) with results from SNP array genotyping. Coding insertions and deletions (indels) were found to be enriched for size in multiples of 3 and were located near the N- and C-termini of proteins. For larger indels, a combination of split-read and read-pair approaches proved to be complementary in finding different signatures. CNVs were identified on the basis of the depth of sequenced reads, and by using SNP and CGH arrays.ConclusionsOur results provide high resolution mapping of diverse classes of genomic variation in an individual bovine genome and demonstrate that structural variation surpasses sequence variation as the main component of genomic variability. Better accuracy of SNP detection was achieved with little loss of sensitivity when algorithms that implemented mapping quality were used. IBD regions were found to be instrumental for calculating resequencing SNP accuracy, while SNP detection within CNVs tended to be less reliable. CNV discovery was affected dramatically by platform resolution and coverage biases. The combined data for this study showed that at a moderate level of sequencing coverage, an ensemble of platforms and tools can be applied together to maximize the accurate detection of sequence and structural variants.

A Syntenic Cross Species Aneuploidy Genetic Screen Links RCAN1 Expression to β-Cell Mitochondrial Dysfunction in Type 2 Diabetes

Type 2 diabetes (T2D) is a complex metabolic disease associated with obesity, insulin resistance and hypoinsulinemia due to pancreatic β-cell dysfunction. Reduced mitochondrial function is thought to be central to β-cell dysfunction. Mitochondrial dysfunction and reduced insulin secretion are also observed in β-cells of humans with the most common human genetic disorder, Down syndrome (DS, Trisomy 21). To identify regions of chromosome 21 that may be associated with perturbed glucose homeostasis we profiled the glycaemic status of different DS mouse models. The Ts65Dn and Dp16 DS mouse lines were hyperglycemic, while Tc1 and Ts1Rhr mice were not, providing us with a region of chromosome 21 containing genes that cause hyperglycemia. We then examined whether any of these genes were upregulated in a set of ~5,000 gene expression changes we had identified in a large gene expression analysis of human T2D β-cells. This approach produced a single gene, RCAN1, as a candidate gene linking hyperglycemia and functional changes in T2D β-cells. Further investigations demonstrated that RCAN1 methylation is reduced in human T2D islets at multiple sites, correlating with increased expression. RCAN1 protein expression was also increased in db/db mouse islets and in human and mouse islets exposed to high glucose. Mice overexpressing RCAN1 had reduced in vivo glucose-stimulated insulin secretion and their β-cells displayed mitochondrial dysfunction including hyperpolarised membrane potential, reduced oxidative phosphorylation and low ATP production. This lack of β-cell ATP had functional consequences by negatively affecting both glucose-stimulated membrane depolarisation and ATP-dependent insulin granule exocytosis. Thus, from amongst the myriad of gene expression changes occurring in T2D β-cells where we had little knowledge of which changes cause β-cell dysfunction, we applied a trisomy 21 screening approach which linked RCAN1 to β-cell mitochondrial dysfunction in T2D.