Genetic regulation of RNA splicing in human pancreatic islets

BACKGROUND

Non-coding genetic variants that influence gene transcription in pancreatic islets play a major role in the susceptibility to type 2 diabetes (T2D), and likely also contribute to type 1 diabetes (T1D) risk. For many loci, however, the mechanisms through which non-coding variants influence diabetes susceptibility are unknown.

RESULTS

We examine splicing QTLs (sQTLs) in pancreatic islets from 399 human donors and observe that common genetic variation has a widespread influence on the splicing of genes with established roles in islet biology and diabetes. In parallel, we profile expression QTLs (eQTLs) and use transcriptome-wide association as well as genetic co-localization studies to assign islet sQTLs or eQTLs to T2D and T1D susceptibility signals, many of which lack candidate effector genes. This analysis reveals biologically plausible mechanisms, including the association of T2D with an sQTL that creates a nonsense isoform in ERO1B, a regulator of ER-stress and proinsulin biosynthesis. The expanded list of T2D risk effector genes reveals overrepresented pathways, including regulators of G-protein-mediated cAMP production. The analysis of sQTLs also reveals candidate effector genes for T1D susceptibility such as DCLRE1B, a senescence regulator, and lncRNA MEG3.

CONCLUSIONS

These data expose widespread effects of common genetic variants on RNA splicing in pancreatic islets. The results support a role for splicing variation in diabetes susceptibility, and offer a new set of genetic targets with potential therapeutic benefit.

Fine-mapping, trans-ancestral and genomic analyses identify causal variants, cells, genes and drug targets for type 1 diabetes

We report the largest and most diverse genetic study of type 1 diabetes (T1D) to date (61,427 participants), yielding 78 genome-wide-significant (P < 5 × 10-8) regions, including 36 that are new. We define credible sets of T1D-associated variants and show that they are enriched in immune-cell accessible chromatin, particularly CD4+ effector T cells. Using chromatin-accessibility profiling of CD4+ T cells from 115 individuals, we map chromatin-accessibility quantitative trait loci and identify five regions where T1D risk variants co-localize with chromatin-accessibility quantitative trait loci. We highlight rs72928038 in BACH2 as a candidate causal T1D variant leading to decreased enhancer accessibility and BACH2 expression in T cells. Finally, we prioritize potential drug targets by integrating genetic evidence, functional genomic maps and immune protein-protein interactions, identifying 12 genes implicated in T1D that have been targeted in clinical trials for autoimmune diseases. These findings provide an expanded genomic landscape for T1D.

Analysis of overlapping genetic association in type 1 and type 2 diabetes

AIMS/HYPOTHESIS

Given the potential shared aetiology between type 1 and type 2 diabetes, we aimed to identify any genetic regions associated with both diseases. For associations where there is a shared signal and the allele that increases risk to one disease also increases risk to the other, inference about shared aetiology could be made, with the potential to develop therapeutic strategies to treat or prevent both diseases simultaneously. Alternatively, if a genetic signal co-localises with divergent effect directions, it could provide valuable biological insight into how the association affects the two diseases differently.

METHODS

Using publicly available type 2 diabetes summary statistics from a genome-wide association study (GWAS) meta-analysis of European ancestry individuals (74,124 cases and 824,006 controls) and type 1 diabetes GWAS summary statistics from a meta-analysis of studies on individuals from the UK and Sardinia (7467 cases and 10,218 controls), we identified all regions of 0.5 Mb that contained variants associated with both diseases (false discovery rate <0.01). In each region, we performed forward stepwise logistic regression to identify independent association signals, then examined co-localisation of each type 1 diabetes signal with each type 2 diabetes signal using coloc. Any association with a co-localisation posterior probability of ≥0.9 was considered a genuine shared association with both diseases.

RESULTS

Of the 81 association signals from 42 genetic regions that showed association with both type 1 and type 2 diabetes, four association signals co-localised between both diseases (posterior probability ≥0.9): (1) chromosome 16q23.1, near CTRB1/BCAR1, which has been previously identified; (2) chromosome 11p15.5, near the INS gene; (3) chromosome 4p16.3, near TMEM129 and (4) chromosome 1p31.3, near PGM1. In each of these regions, the effect of genetic variants on type 1 diabetes was in the opposite direction to the effect on type 2 diabetes. Use of additional datasets also supported the previously identified co-localisation on chromosome 9p24.2, near the GLIS3 gene, in this case with a concordant direction of effect.

CONCLUSIONS/INTERPRETATION

Four of five association signals that co-localise between type 1 diabetes and type 2 diabetes are in opposite directions, suggesting a complex genetic relationship between the two diseases.

Somatic selection of poorly differentiating variant stem cell clones could be a key to human ageing

Any replicating system in which heritable variants with differing replicative potentials can arise is subject to a Darwinian evolutionary process. The continually replicating adult tissue stem cells that control the integrity of many tissues of long-lived, multicellular, complex vertebrate organisms, including humans, constitute such a replicating system. Our suggestion is that somatic selection for mutations (or stable epigenetic changes) that cause an increased rate of adult tissue stem cell proliferation, and their long-term persistence, at the expense of normal differentiation, is a major key to the ageing process. Once an organism has passed the reproductive age, there is no longer any significant counterselection at the organismal level to this inevitable cellular level Darwinian process.

Genetic Variants Predisposing Most Strongly to Type 1 Diabetes Diagnosed Under Age 7 Years Lie Near Candidate Genes That Function in the Immune System and in Pancreatic β-Cells

OBJECTIVE

Immunohistological analyses of pancreata from patients with type 1 diabetes suggest distinct autoimmune islet β-cell pathology between those diagnosed at <7 years (<7 group) and those diagnosed at age ≥13 years (≥13 group), with both B- and T-lymphocyte islet inflammation common in children in the <7 group, whereas B cells are rare in the ≥13 group. Based on these observations, we sought to identify differences in genetic susceptibility between these prespecified age-at-diagnosis groups to inform on the etiology of the most aggressive form of type 1 diabetes that initiates in the first years of life.

RESEARCH DESIGN AND METHODS

Using multinomial logistic regression models, we tested if known type 1 diabetes loci (17 within the HLA and 55 non-HLA loci) had significantly stronger effect sizes in the <7 group compared with the ≥13 group, using genotype data from 27,071 individuals (18,485 control subjects and 3,121 case subjects diagnosed at <7 years, 3,757 at 7-13 years, and 1,708 at ≥13 years).

RESULTS

Six HLA haplotypes/classical alleles and six non-HLA regions, one of which functions specifically in β-cells (GLIS3) and the other five likely affecting key T-cell (IL2RA, IL10, IKZF3, and THEMIS), thymus (THEMIS), and B-cell development/functions (IKZF3 and IL10) or in both immune and β-cells (CTSH), showed evidence for stronger effects in the <7 group.

CONCLUSIONS

A subset of type 1 diabetes-associated variants are more prevalent in children diagnosed under the age of 7 years and are near candidate genes that act in both pancreatic β- and immune cells.

Chronic Immune Activation in Systemic Lupus Erythematosus and the Autoimmune PTPN22 Trp620 Risk Allele Drive the Expansion of FOXP3+ Regulatory T Cells and PD-1 Expression

In systemic lupus erythematosus (SLE), perturbed immunoregulation underpins a pathogenic imbalance between regulatory and effector CD4⁺ T-cell activity. However, to date, the characterization of the CD4⁺ regulatory T cell (Treg) compartment in SLE has yielded conflicting results. Here we show that patients have an increased frequency of CD4⁺FOXP3⁺ cells in circulation owing to a specific expansion of thymically-derived FOXP3⁺HELIOS⁺ Tregs with a demethylated FOXP3 Treg-specific demethylated region. We found that the Treg expansion was strongly associated with markers of recent immune activation, including PD-1, plasma concentrations of IL-2 and the type I interferon biomarker soluble SIGLEC-1. Since the expression of the negative T-cell signaling molecule PTPN22 is increased and a marker of poor prognosis in SLE, we tested the influence of its missense risk allele Trp⁶²⁰ (rs2476601C>T) on Treg frequency. Trp⁶²⁰ was reproducibly associated with increased frequencies of thymically-derived Tregs in blood, and increased PD-1 expression on both Tregs and effector T cells (Teffs). Our results support the hypothesis that FOXP3⁺ Tregs are increased in SLE patients as a consequence of a compensatory mechanism in an attempt to regulate pathogenic autoreactive Teff activity. We suggest that restoration of IL-2-mediated homeostatic regulation of FOXP3⁺ Tregs by IL-2 administration could prevent disease flares rather than treating at the height of a disease flare. Moreover, stimulation of PD-1 with specific agonists, perhaps in combination with low-dose IL-2, could be an effective therapeutic strategy in autoimmune disease and in other immune disorders.

Genetics of the human face: Identification of large-effect single gene variants

To discover specific variants with relatively large effects on the human face, we have devised an approach to identifying facial features with high heritability. This is based on using twin data to estimate the additive genetic value of each point on a face, as provided by a 3D camera system. In addition, we have used the ethnic difference between East Asian and European faces as a further source of face genetic variation. We use principal components (PCs) analysis to provide a fine definition of the surface features of human faces around the eyes and of the profile, and chose upper and lower 10% extremes of the most heritable PCs for looking for genetic associations. Using this strategy for the analysis of 3D images of 1,832 unique volunteers from the well-characterized People of the British Isles study and 1,567 unique twin images from the TwinsUK cohort, together with genetic data for 500,000 SNPs, we have identified three specific genetic variants with notable effects on facial profiles and eyes.

Statistical aspects of evolution under natural selection, with implications for the advantage of sexual reproduction

The prevalence of sexual reproduction remains mysterious, as it poses clear evolutionary drawbacks compared to reproducing asexually. Several possible explanations exist, with one of the most likely being that finite population size causes linkage disequilibria to randomly generate and impede the progress of natural selection, and that these are eroded by recombination via sexual reproduction. Previous investigations have either analysed this phenomenon in detail for small numbers of loci, or performed population simulations for many loci. Here we present a quantitative genetic model for fitness, based on the Price Equation, in order to examine the theoretical consequences of randomly generated linkage disequilibria when there are many loci. In addition, most previous work has been concerned with the long-term consequences of deleterious linkage disequilibria for population fitness. The expected change in mean fitness between consecutive generations, a measure of short-term evolutionary success, is shown under random environmental influences to be related to the autocovariance in mean fitness between the generations, capturing the effects of stochastic forces such as genetic drift. Interaction between genetic drift and natural selection, due to randomly generated linkage disequilibria, is demonstrated to be one possible source of mean fitness autocovariance. This suggests a possible role for sexual reproduction in reducing the negative effects of genetic drift, thereby improving the short-term efficacy of natural selection.

Inferring separate parental admixture components in unknown DNA samples using autosomal SNPs

The identification of ancestral admixture proportions for human DNA samples has recently had success in forensic cases. Current methods infer admixture proportions for the target sample, but not for their parents, which provides an additional layer of information that may aid certain forensic investigations. We describe new maximum likelihood methods (LEAPFrOG and LEAPFrOG Expectation Maximisation), for inferring both an individual's admixture proportions and the admixture proportions possessed by the unobserved parents, with respect to two or more source populations, using single-nucleotide polymorphism data typed only in the target individual. This is achieved by examining the increase in heterozygosity in the offspring of parents who are from different populations or who represent different mixtures from a number of source populations. We validated the methods via simulation; combining chromosomes from different Hapmap Phase III population samples to emulate first-generation admixture. Performance was strong for individuals with mixed African/European (YRI/CEU) ancestry, but poor for mixed Japanese/Chinese (JPT/CHB) ancestry, reflecting the difficulty in distinguishing closely related source populations. A total of 11 African-American trios were used to compare the parental admixture inferred from their own genotypes against that inferred purely from their offspring genotypes. We examined the performance of 34 ancestry informative markers from a multiplex kit for ancestry inference. Simulations showed that estimates were unreliable when parents had similar admixture, suggesting more markers are needed. Our results demonstrate that ancestral backgrounds of case samples and their parents are obtainable to aid in forensic investigations, provided that high-throughput methods are adopted by the forensic community.