Meningeal lymphatics affect microglia responses and anti-Aß immunotherapy.

Alzheimer's disease (AD) is the most prevalent cause of dementia1. Although there is no effective treatment for AD, passive immunotherapy with monoclonal antibodies against amyloid beta (Aß) is a promising therapeutic strategy2,3. Meningeal lymphatic drainage has an important role in the accumulation of Aß in the brain4, but it is not known whether modulation of meningeal lymphatic function can influence the outcome of immunotherapy in AD. Here we show that ablation of meningeal lymphatic vessels in 5xFAD mice (a mouse model of amyloid deposition that expresses five mutations found in familial AD) worsened the outcome of mice treated with anti-Aß passive immunotherapy by exacerbating the deposition of Aß, microgliosis, neurovascular dysfunction, and behavioural deficits. By contrast, therapeutic delivery of vascular endothelial growth factor C improved clearance of Aß by monoclonal antibodies. Notably, there was a substantial overlap between the gene signature of microglia from 5xFAD mice with impaired meningeal lymphatic function and the transcriptional profile of activated microglia from the brains of individuals with AD. Overall, our data demonstrate that impaired meningeal lymphatic drainage exacerbates the microglial inflammatory response in AD and that enhancement of meningeal lymphatic function combined with immunotherapies could lead to better clinical outcomes.

A multi-ancestry genome-wide association study in type 1 diabetes.

Type 1 diabetes (T1D) is an autoimmune disease caused by destruction of the pancreatic ß-cells. Genome-wide association (GWAS) and fine mapping studies have been conducted mainly in European ancestry (EUR) populations. We performed a multi-ancestry GWAS to identify SNPs and HLA alleles associated with T1D risk and age at onset. EUR families (N = 3223), and unrelated individuals of African (AFR, N = 891) and admixed (Hispanic/Latino) ancestry (AMR, N = 308) were genotyped using the Illumina HumanCoreExome BeadArray, with imputation to the TOPMed reference panel. The Multi-Ethnic HLA reference panel was utilized to impute HLA alleles and amino acid residues. Logistic mixed models (T1D risk) and frailty models (age at onset) were used for analysis. In GWAS meta-analysis, seven loci were associated with T1D risk at genome-wide significance: PTPN22, HLA-DQA1, IL2RA, RNLS, INS, IKZF4-RPS26-ERBB3, and SH2B3, with four associated with T1D age at onset (PTPN22, HLA-DQB1, INS, and ERBB3). AFR and AMR meta-analysis revealed NRP1 as associated with T1D risk and age at onset, although NRP1 variants were not associated in EUR ancestry. In contrast, the PTPN22 variant was significantly associated with risk only in EUR ancestry. HLA alleles and haplotypes most significantly associated with T1D risk in AFR and AMR ancestry differed from that seen in EUR ancestry; in addition, the HLA-DRB1*08:02-DQA1*04:01-DQB1*04:02 haplotype was 'protective' in AMR while HLA-DRB1*08:01-DQA1*04:01-DQB1*04:02 haplotype was 'risk' in EUR ancestry, differing only at HLA-DRB1*08. These results suggest that much larger sample sizes in non-EUR populations are required to capture novel loci associated with T1D risk.

Leveraging type 1 diabetes human genetic and genomic data in the T1D knowledge portal.

To address the challenge of translating genetic discoveries for type 1 diabetes (T1D) into mechanistic insight, we have developed the T1D Knowledge Portal (T1DKP), an open-access resource for hypothesis development and target discovery in T1D.

Genetic variants in the complement system and their potential link in the aetiology of type 1 diabetes.

Type 1 diabetes is an autoimmune disease in which one's own immune system destroys insulin-secreting beta cells in the pancreas. This process results in life-long dependence on exogenous insulin for survival. Both genetic and environmental factors play a role in disease initiation, progression, and ultimate clinical diagnosis of type 1 diabetes. This review will provide background on the natural history of type 1 diabetes and the role of genetic factors involved in the complement system, as several recent studies have identified changes in levels of these proteins as the disease evolves from pre-clinical through to clinically apparent disease.

Efferocytosis induces a novel SLC program to promote glucose uptake and lactate release.

Development and routine tissue homeostasis require a high turnover of apoptotic cells. These cells are removed by professional and non-professional phagocytes via efferocytosis1. How a phagocyte maintains its homeostasis while coordinating corpse uptake, processing ingested materials and secreting anti-inflammatory mediators is incompletely understood1,2. Here, using RNA sequencing to characterize the transcriptional program of phagocytes actively engulfing apoptotic cells, we identify a genetic signature involving 33 members of the solute carrier (SLC) family of membrane transport proteins, in which expression is specifically modulated during efferocytosis, but not during antibody-mediated phagocytosis. We assessed the functional relevance of these SLCs in efferocytic phagocytes and observed a robust induction of an aerobic glycolysis program, initiated by SLC2A1-mediated glucose uptake, with concurrent suppression of the oxidative phosphorylation program. The different steps of phagocytosis2-that is, 'smell' ('find-me' signals or sensing factors released by apoptotic cells), 'taste' (phagocyte-apoptotic cell contact) and 'ingestion' (corpse internalization)-activated distinct and overlapping sets of genes, including several SLC genes, to promote glycolysis. SLC16A1 was upregulated after corpse uptake, increasing the release of lactate, a natural by-product of aerobic glycolysis3. Whereas glycolysis within phagocytes contributed to actin polymerization and the continued uptake of corpses, lactate released via SLC16A1 promoted the establishment of an anti-inflammatory tissue environment. Collectively, these data reveal a SLC program that is activated during efferocytosis, identify a previously unknown reliance on aerobic glycolysis during apoptotic cell uptake and show that glycolytic by-products of efferocytosis can influence surrounding cells.

Publisher Correction: Functional aspects of meningeal lymphatics in ageing and Alzheimer's disease.

Change history: In this Article, Extended Data Fig. 9 was appearing as Fig. 2 in the HTML, and in Fig. 2, the panel labels 'n' and 'o' overlapped the figure; these errors have been corrected online.

Functional aspects of meningeal lymphatics in ageing and Alzheimer's disease.

Ageing is a major risk factor for many neurological pathologies, but its mechanisms remain unclear. Unlike other tissues, the parenchyma of the central nervous system (CNS) lacks lymphatic vasculature and waste products are removed partly through a paravascular route. (Re)discovery and characterization of meningeal lymphatic vessels has prompted an assessment of their role in waste clearance from the CNS. Here we show that meningeal lymphatic vessels drain macromolecules from the CNS (cerebrospinal and interstitial fluids) into the cervical lymph nodes in mice. Impairment of meningeal lymphatic function slows paravascular influx of macromolecules into the brain and efflux of macromolecules from the interstitial fluid, and induces cognitive impairment in mice. Treatment of aged mice with vascular endothelial growth factor C enhances meningeal lymphatic drainage of macromolecules from the cerebrospinal fluid, improving brain perfusion and learning and memory performance. Disruption of meningeal lymphatic vessels in transgenic mouse models of Alzheimer's disease promotes amyloid-ß deposition in the meninges, which resembles human meningeal pathology, and aggravates parenchymal amyloid-ß accumulation. Meningeal lymphatic dysfunction may be an aggravating factor in Alzheimer's disease pathology and in age-associated cognitive decline. Thus, augmentation of meningeal lymphatic function might be a promising therapeutic target for preventing or delaying age-associated neurological diseases.

Frequent somatic TET2 mutations in chronic NK-LGL leukemia with distinct patterns of cytopenias.

Chronic natural killer large granular lymphocyte (NK-LGL) leukemia, also referred to as chronic lymphoproliferative disorder of NK cells, is a rare disorder defined by prolonged expansion of clonal NK cells. Similar prevalence of STAT3 mutations in chronic T-LGL and NK-LGL leukemia is suggestive of common pathogenesis. We undertook whole-genome sequencing to identify mutations unique to NK-LGL leukemia. The results were analyzed to develop a resequencing panel that was applied to 58 patients. Phosphatidylinositol 3-kinase pathway gene mutations (PIK3CD/PIK3AP1) and TNFAIP3 mutations were seen in 5% and 10% of patients, respectively. TET2 was exceptional in that mutations were present in 16 (28%) of 58 patient samples, with evidence that TET2 mutations can be dominant and exclusive to the NK compartment. Reduced-representation bisulfite sequencing revealed that methylation patterns were significantly altered in TET2 mutant samples. The promoter of TET2 and that of PTPRD, a negative regulator of STAT3, were found to be methylated in additional cohort samples, largely confined to the TET2 mutant group. Mutations in STAT3 were observed in 19 (33%) of 58 patient samples, 7 of which had concurrent TET2 mutations. Thrombocytopenia and resistance to immunosuppressive agents were uniquely observed in those patients with only TET2 mutation (Games-Howell post hoc test, P = .0074; Fisher's exact test, P = .00466). Patients with STAT3 mutation, inclusive of those with TET2 comutation, had lower hematocrit, hemoglobin, and absolute neutrophil count compared with STAT3 wild-type patients (Welch's t test, P ≤ .015). We present the discovery of TET2 mutations in chronic NK-LGL leukemia and evidence that it identifies a unique molecular subtype.

Genetic Regulation of Atherosclerosis-Relevant Phenotypes in Human Vascular Smooth Muscle Cells.

RATIONALE: Coronary artery disease (CAD) is a major cause of morbidity and mortality worldwide. Recent genome-wide association studies revealed 163 loci associated with CAD. However, the precise molecular mechanisms by which the majority of these loci increase CAD risk are not known. Vascular smooth muscle cells (VSMCs) are critical in the development of CAD. They can play either beneficial or detrimental roles in lesion pathogenesis, depending on the nature of their phenotypic changes. OBJECTIVE: To identify genetic variants associated with atherosclerosis-relevant phenotypes in VSMCs. METHODS AND RESULTS: We quantified 12 atherosclerosis-relevant phenotypes related to calcification, proliferation, and migration in VSMCs isolated from 151 multiethnic heart transplant donors. After genotyping and imputation, we performed association mapping using 6.3 million genetic variants. We demonstrated significant variations in calcification, proliferation, and migration. These phenotypes were not correlated with each other. We performed genome-wide association studies for 12 atherosclerosis-relevant phenotypes and identified 4 genome-wide significant loci associated with at least one VSMC phenotype. We overlapped the previously identified CAD loci with our data set and found nominally significant associations at 79 loci. One of them was the chromosome 1q41 locus, which harbors MIA3. The G allele of the lead risk single nucleotide polymorphism (SNP) rs67180937 was associated with lower VSMC MIA3 expression and lower proliferation. Lentivirus-mediated silencing of MIA3 (melanoma inhibitory activity protein 3) in VSMCs resulted in lower proliferation, consistent with human genetics findings. Furthermore, we observed a significant reduction of MIA3 protein in VSMCs in thin fibrous caps of late-stage atherosclerotic plaques compared to early fibroatheroma with thick and protective fibrous caps in mice and humans. CONCLUSIONS: Our data demonstrate that genetic variants have significant influences on VSMC function relevant to the development of atherosclerosis. Furthermore, high MIA3 expression may promote atheroprotective VSMC phenotypic transitions, including increased proliferation, which is essential in the formation or maintenance of a protective fibrous cap.

Macrophages redirect phagocytosis by non-professional phagocytes and influence inflammation.

Professional phagocytes (such as macrophages) and non-professional phagocytes (such as epithelial cells) clear billions of apoptotic cells and particles on a daily basis. Although professional and non-professional macrophages reside in proximity in most tissues, whether they communicate with each other during cell clearance, and how this might affect inflammation, is not known. Here we show that macrophages, through the release of a soluble growth factor and microvesicles, alter the type of particles engulfed by non-professional phagocytes and influence their inflammatory response. During phagocytosis of apoptotic cells or in response to inflammation-associated cytokines, macrophages released insulin-like growth factor 1 (IGF-1). The binding of IGF-1 to its receptor on non-professional phagocytes redirected their phagocytosis, such that uptake of larger apoptotic cells was reduced whereas engulfment of microvesicles was increased. IGF-1 did not alter engulfment by macrophages. Macrophages also released microvesicles, whose uptake by epithelial cells was enhanced by IGF-1 and led to decreased inflammatory responses by epithelial cells. Consistent with these observations, deletion of IGF-1 receptor in airway epithelial cells led to exacerbated lung inflammation after allergen exposure. These genetic and functional studies reveal that IGF-1- and microvesicle-dependent communication between macrophages and epithelial cells can critically influence the magnitude of tissue inflammation in vivo.

Association of Coding Variants in Hydroxysteroid 17-beta Dehydrogenase 14 (HSD17B14) with Reduced Progression to End Stage Kidney Disease in Type 1 Diabetes.

BACKGROUND: Rare variants in gene coding regions likely have a greater impact on disease-related phenotypes than common variants through disruption of their encoded protein. We searched for rare variants associated with onset of ESKD in individuals with type 1 diabetes at advanced kidney disease stage. METHODS: Gene-based exome array analyses of 15,449 genes in five large incidence cohorts of individuals with type 1 diabetes and proteinuria were analyzed for survival time to ESKD, testing the top gene in a sixth cohort (n=2372/1115 events all cohorts) and replicating in two retrospective case-control studies (n=1072 cases, 752 controls). Deep resequencing of the top associated gene in five cohorts confirmed the findings. We performed immunohistochemistry and gene expression experiments in human control and diseased cells, and in mouse ischemia reperfusion and aristolochic acid nephropathy models. RESULTS: Protein coding variants in the hydroxysteroid 17-ß dehydrogenase 14 gene (HSD17B14), predicted to affect protein structure, had a net protective effect against development of ESKD at exome-wide significance (n=4196; P value=3.3 × 10-7). The HSD17B14 gene and encoded enzyme were robustly expressed in healthy human kidney, maximally in proximal tubular cells. Paradoxically, gene and protein expression were attenuated in human diabetic proximal tubules and in mouse kidney injury models. Expressed HSD17B14 gene and protein levels remained low without recovery after 21 days in a murine ischemic reperfusion injury model. Decreased gene expression was found in other CKD-associated renal pathologies. CONCLUSIONS: HSD17B14 gene is mechanistically involved in diabetic kidney disease. The encoded sex steroid enzyme is a druggable target, potentially opening a new avenue for therapeutic development.

Heterogeneous long-term trajectories of glycaemic control in type 1 diabetes.

AIMS: We aimed to identify long-term HbA1c trajectories and examine associated characteristics in an observational, childhood-onset (<17 years) type 1 diabetes cohort. METHODS: Data are from the Epidemiology of Diabetes Complications study, comprising 405 participants with ≥2 of seven possible HbA1c measurements over follow-up (1988-2013) and available DNA (baseline mean diabetes duration 21 years, 53% men). HbA1c trajectories were estimated using latent class growth models. Baseline and change in participant characteristics were compared across trajectories. RESULTS: Five HbA1c trajectories were identified: low (51%), intermediate stable (22%), improved (19%), high stable (6%), and worsened (2%; not included in analyses). Age, diabetes duration, diabetes onset age, and sex did not differ across trajectories. Characteristics did not differ significantly between intermediate stable and low trajectories at baseline, though albumin excretion rate (AER, p = 0.0002) and estimated glomerular filtration rate (eGFR, p = 0.001) worsened slightly more in intermediate stable over time. Improved and high stable trajectories had higher baseline LDL-c (p = 0.002 and 0.003, respectively). Improved trajectory increased median self-monitoring of blood glucose from <1 to 3.5 times/day (p < 0.0001) and had larger LDL-c improvement (p = 0.004) but greater worsening of AER (p < 0.0001) and eGFR (p < 0.0001) than low. The A allele of rs12970134 (near MC4R) was associated with improved (p = 0.0003) or high stable (p = 0.001) HbA1c trajectory, both patterns with high baseline HbA1c. CONCLUSIONS: Long-term HbA1c trajectories were primarily associated with modifiable factors in this type 1 diabetes cohort. The intermediate stable pattern had a risk factor profile that suggests some protection against adverse metabolic effects of chronic hyperglycaemia, warranting further study.

PRF1 mutation alters immune system activation, inflammation, and risk of autoimmunity.

BACKGROUND: Defective alleles within the PRF1 gene, encoding the pore-forming protein perforin, in combination with environmental factors, cause familial type 2 hemophagocytic lymphohistiocytosis (FHL2), a rare, severe autosomal recessive childhood disorder characterized by massive release of cytokines-cytokine storm. OBJECTIVE: The aim of this study was to determine the function of hypomorph PRF1:p.A91V g.72360387 G > A on multiple sclerosis (MS) and type 1 diabetes (T1D). METHODS: We cross-compare the association data for PRF1:p.A91V mutation derived from GWAS on adult MS and pediatric T1D in Sardinians. The novel association with T1D was replicated in metanalysis in 12,584 cases and 17,692 controls from Sardinia, the United Kingdom, and Scotland. To dissect this mutation function, we searched through the coincident association immunophenotypes in additional set of general population Sardinians. RESULTS: We report that PRF1:p.A91V, is associated with increase of lymphocyte levels, especially within the cytotoxic memory T-cells, at general population level with reduced interleukin 7 receptor expression on these cells. The minor allele increased risk of MS, in 2903 cases and 2880 controls from Sardinia p = 2.06 × 10-4, odds ratio OR = 1.29, replicating a previous finding, whereas it protects from T1D p = 1.04 × 10-5, OR = 0.82. CONCLUSION: Our results indicate opposing contributions of the cytotoxic T-cell compartment to MS and T1D pathogenesis.

Disease-specific biases in alternative splicing and tissue-specific dysregulation revealed by multitissue profiling of lymphocyte gene expression in type 1 diabetes.

Genome-wide association studies (GWAS) have identified multiple, shared allelic associations with many autoimmune diseases. However, the pathogenic contributions of variants residing in risk loci remain unresolved. The location of the majority of shared disease-associated variants in noncoding regions suggests they contribute to risk of autoimmunity through effects on gene expression in the immune system. In the current study, we test this hypothesis by applying RNA sequencing to CD4+, CD8+, and CD19+ lymphocyte populations isolated from 81 subjects with type 1 diabetes (T1D). We characterize and compare the expression patterns across these cell types for three gene sets: all genes, the set of genes implicated in autoimmune disease risk by GWAS, and the subset of these genes specifically implicated in T1D. We performed RNA sequencing and aligned the reads to both the human reference genome and a catalog of all possible splicing events developed from the genome, thereby providing a comprehensive evaluation of the roles of gene expression and alternative splicing (AS) in autoimmunity. Autoimmune candidate genes displayed greater expression specificity in the three lymphocyte populations relative to other genes, with significantly increased levels of splicing events, particularly those predicted to have substantial effects on protein isoform structure and function (e.g., intron retention, exon skipping). The majority of single-nucleotide polymorphisms within T1D-associated loci were also associated with one or more cis-expression quantitative trait loci (cis-eQTLs) and/or splicing eQTLs. Our findings highlight a substantial, and previously underrecognized, role for AS in the pathogenesis of autoimmune disorders and particularly for T1D.

The active contribution of OPCs to neuroinflammation is mediated by LRP1.

Oligodendrocyte progenitor cells (OPCs) account for about 5% of total brain and spinal cord cells, giving rise to myelinating oligodendrocytes that provide electrical insulation to neurons of the CNS. OPCs have also recently been shown to regulate inflammatory responses and glial scar formation, suggesting functions that extend beyond myelination. Low-density lipoprotein receptor-related protein 1 (LRP1) is a multifaceted phagocytic receptor that is highly expressed in several CNS cell types, including OPCs. Here, we have generated an oligodendroglia-specific knockout of LRP1, which presents with normal myelin development, but is associated with better outcomes in two animal models of demyelination (EAE and cuprizone). At a mechanistic level, LRP1 did not directly affect OPC differentiation into mature oligodendrocytes. Instead, animals lacking LRP1 in OPCs in the demyelinating CNS were characterized by a robust dampening of inflammation. In particular, LRP1-deficient OPCs presented with impaired antigen cross-presentation machinery, suggesting a failure to propagate the inflammatory response and thus promoting faster myelin repair and neuroprotection. Our study places OPCs as major regulators of neuroinflammation in an LRP1-dependent fashion.

Genome-Wide Association Study of Diabetic Kidney Disease Highlights Biology Involved in Glomerular Basement Membrane Collagen.

BACKGROUND: Although diabetic kidney disease demonstrates both familial clustering and single nucleotide polymorphism heritability, the specific genetic factors influencing risk remain largely unknown. METHODS: To identify genetic variants predisposing to diabetic kidney disease, we performed genome-wide association study (GWAS) analyses. Through collaboration with the Diabetes Nephropathy Collaborative Research Initiative, we assembled a large collection of type 1 diabetes cohorts with harmonized diabetic kidney disease phenotypes. We used a spectrum of ten diabetic kidney disease definitions based on albuminuria and renal function. RESULTS: Our GWAS meta-analysis included association results for up to 19,406 individuals of European descent with type 1 diabetes. We identified 16 genome-wide significant risk loci. The variant with the strongest association (rs55703767) is a common missense mutation in the collagen type IV alpha 3 chain (COL4A3) gene, which encodes a major structural component of the glomerular basement membrane (GBM). Mutations in COL4A3 are implicated in heritable nephropathies, including the progressive inherited nephropathy Alport syndrome. The rs55703767 minor allele (Asp326Tyr) is protective against several definitions of diabetic kidney disease, including albuminuria and ESKD, and demonstrated a significant association with GBM width; protective allele carriers had thinner GBM before any signs of kidney disease, and its effect was dependent on glycemia. Three other loci are in or near genes with known or suggestive involvement in this condition (BMP7) or renal biology (COLEC11 and DDR1). CONCLUSIONS: The 16 diabetic kidney disease-associated loci may provide novel insights into the pathogenesis of this condition and help identify potential biologic targets for prevention and treatment.

Negligible impact of rare autoimmune-locus coding-region variants on missing heritability.

Genome-wide association studies (GWAS) have identified common variants of modest-effect size at hundreds of loci for common autoimmune diseases; however, a substantial fraction of heritability remains unexplained, to which rare variants may contribute. To discover rare variants and test them for association with a phenotype, most studies re-sequence a small initial sample size and then genotype the discovered variants in a larger sample set. This approach fails to analyse a large fraction of the rare variants present in the entire sample set. Here we perform simultaneous amplicon-sequencing-based variant discovery and genotyping for coding exons of 25 GWAS risk genes in 41,911 UK residents of white European origin, comprising 24,892 subjects with six autoimmune disease phenotypes and 17,019 controls, and show that rare coding-region variants at known loci have a negligible role in common autoimmune disease susceptibility. These results do not support the rare-variant synthetic genome-wide-association hypothesis (in which unobserved rare causal variants lead to association detected at common tag variants). Many known autoimmune disease risk loci contain multiple, independently associated, common and low-frequency variants, and so genes at these loci are a priori stronger candidates for harbouring rare coding-region variants than other genes. Our data indicate that the missing heritability for common autoimmune diseases may not be attributable to the rare coding-region variant portion of the allelic spectrum, but perhaps, as others have proposed, may be a result of many common-variant loci of weak effect.

Neuromedin B Expression Defines the Mouse Retrotrapezoid Nucleus.

The retrotrapezoid nucleus (RTN) consists, by definition, of Phox2b-expressing, glutamatergic, non-catecholaminergic, noncholinergic neurons located in the parafacial region of the medulla oblongata. An unknown proportion of RTN neurons are central respiratory chemoreceptors and there is mounting evidence for biochemical diversity among these cells. Here, we used multiplexed in situ hybridization and single-cell RNA-Seq in male and female mice to provide a more comprehensive view of the phenotypic diversity of RTN neurons. We now demonstrate that the RTN of mice can be identified with a single and specific marker, Neuromedin B mRNA (Nmb). Most (∼75%) RTN neurons express low-to-moderate levels of Nmb and display chemoreceptor properties. Namely they are activated by hypercapnia, but not by hypoxia, and express proton sensors, TASK-2 and Gpr4. These Nmb-low RTN neurons also express varying levels of transcripts for Gal, Penk, and Adcyap1, and receptors for substance P, orexin, serotonin, and ATP. A subset of RTN neurons (∼20-25%), typically larger than average, express very high levels of Nmb mRNA. These Nmb-high RTN neurons do not express Fos after hypercapnia and have low-to-undetectable levels of Kcnk5 or Gpr4 transcripts; they also express Adcyap1, but are essentially devoid of Penk and Gal transcripts. In male rats, Nmb is also a marker of the RTN but, unlike in mice, this gene is expressed by other types of nearby neurons located within the ventromedial medulla. In sum, Nmb is a selective marker of the RTN in rodents; Nmb-low neurons, the vast majority, are central respiratory chemoreceptors, whereas Nmb-high neurons likely have other functions.SIGNIFICANCE STATEMENT Central respiratory chemoreceptors regulate arterial PCO2 by adjusting lung ventilation. Such cells have recently been identified within the retrotrapezoid nucleus (RTN), a brainstem nucleus defined by genetic lineage and a cumbersome combination of markers. Using single-cell RNA-Seq and multiplexed in situ hybridization, we show here that a single marker, Neuromedin B mRNA (Nmb), identifies RTN neurons in rodents. We also suggest that >75% of these Nmb neurons are chemoreceptors because they are strongly activated by hypercapnia and express high levels of proton sensors (Kcnk5 and Gpr4). The other RTN neurons express very high levels of Nmb, but low levels of Kcnk5/Gpr4/pre-pro-galanin/pre-pro-enkephalin, and do not respond to hypercapnia. Their function is unknown.

Identification of non-HLA genes associated with development of islet autoimmunity and type 1 diabetes in the prospective TEDDY cohort.

Traditional linkage analysis and genome-wide association studies have identified HLA and a number of non-HLA genes as genetic factors for islet autoimmunity (IA) and type 1 diabetes (T1D). However, the relative risk associated with previously identified non-HLA genes is usually very small as measured in cases/controls from mixed populations. Genetic associations for IA and T1D may be more accurately assessed in prospective cohorts. In this study, 5806 subjects from the TEDDY (The Environmental Determinants of Diabetes in the Young) study, an international prospective cohort study, were genotyped for 176,586 SNPs on the ImmunoChip. Cox proportional hazards analyses were performed to discover the SNPs associated with the risk for IA, T1D, or both. Three regions were associated with the risk of developing any persistent confirmed islet autoantibody: one known region near SH2B3 (HR = 1.35, p = 3.58 × 10-7) with Bonferroni-corrected significance and another known region near PTPN22 (HR = 1.46, p = 2.17 × 10-6) and one novel region near PPIL2 (HR = 2.47, p = 9.64 × 10-7) with suggestive evidence (p < 10-5). Two known regions (PTPN22: p = 2.25 × 10-6, INS; p = 1.32 × 10-7) and one novel region (PXK/PDHB: p = 8.99 × 10-6) were associated with the risk for multiple islet autoantibodies. First appearing islet autoantibodies differ with respect to association. Two regions (INS: p = 5.67 × 10-6 and TTC34/PRDM16: 6.45 × 10-6) were associated if the fist appearing autoantibody was IAA and one region (RBFOX1: p = 8.02 × 10-6) was associated if the first appearing autoantibody was GADA. The analysis of T1D identified one region already known to be associated with T1D (INS: p = 3.13 × 10-7) and three novel regions (RNASET2, PLEKHA1, and PPIL2; 5.42 × 10-6 > p > 2.31 × 10-6). These results suggest that a number of low frequency variants influence the risk of developing IA and/or T1D and these variants can be identified by large prospective cohort studies using a survival analysis approach.

Dissection of a Complex Disease Susceptibility Region Using a Bayesian Stochastic Search Approach to Fine Mapping.

Identification of candidate causal variants in regions associated with risk of common diseases is complicated by linkage disequilibrium (LD) and multiple association signals. Nonetheless, accurate maps of these variants are needed, both to fully exploit detailed cell specific chromatin annotation data to highlight disease causal mechanisms and cells, and for design of the functional studies that will ultimately be required to confirm causal mechanisms. We adapted a Bayesian evolutionary stochastic search algorithm to the fine mapping problem, and demonstrated its improved performance over conventional stepwise and regularised regression through simulation studies. We then applied it to fine map the established multiple sclerosis (MS) and type 1 diabetes (T1D) associations in the IL-2RA (CD25) gene region. For T1D, both stepwise and stochastic search approaches identified four T1D association signals, with the major effect tagged by the single nucleotide polymorphism, rs12722496. In contrast, for MS, the stochastic search found two distinct competing models: a single candidate causal variant, tagged by rs2104286 and reported previously using stepwise analysis; and a more complex model with two association signals, one of which was tagged by the major T1D associated rs12722496 and the other by rs56382813. There is low to moderate LD between rs2104286 and both rs12722496 and rs56382813 (r2 ≃ 0:3) and our two SNP model could not be recovered through a forward stepwise search after conditioning on rs2104286. Both signals in the two variant model for MS affect CD25 expression on distinct subpopulations of CD4+ T cells, which are key cells in the autoimmune process. The results support a shared causal variant for T1D and MS. Our study illustrates the benefit of using a purposely designed model search strategy for fine mapping and the advantage of combining disease and protein expression data.