Genetic Associations with C-peptide Levels before Type 1 Diabetes Diagnosis in At-Risk Relatives.

OBJECTIVE: We sought to determine whether the type 1 diabetes genetic risk score-2 (T1D-GRS2) and single nucleotide polymorphisms (SNPs) are associated with C-peptide preservation before type 1 diabetes diagnosis. METHODS: We conducted a retrospective analysis of 713 autoantibody-positive participants who developed type 1 diabetes in the TrialNet Pathway to Prevention Study who had T1DExomeChip data. We evaluated the relationships of 16 known SNPs and T1D-GRS2 with area under the curve (AUC) C-peptide levels during oral glucose tolerance tests conducted in the 9 months before diagnosis. RESULTS: Higher T1D-GRS2 was associated with lower C-peptide AUC in the 9 months before diagnosis in univariate (ß=-0.06, P<0.0001) and multivariate (ß=-0.03, P=0.005) analyses. Participants with the JAZF1 rs864745 T allele had lower C-peptide AUC in both univariate (ß=-0.11, P=0.002) and multivariate (ß=-0.06, P=0.018) analyses. CONCLUSIONS: The type 2 diabetes-associated JAZF1 rs864745 T allele and higher T1D-GRS2 are associated with lower C-peptide AUC prior to diagnosis of type 1 diabetes, with implications for the design of prevention trials.

HLA Class II (DR, DQ, DP) Genes Were Separately Associated With the Progression From Seroconversion to Onset of Type 1 Diabetes Among Participants in Two Diabetes Prevention Trials (DPT-1 and TN07).

OBJECTIVE: To explore associations of HLA class II genes (HLAII) with the progression of islet autoimmunity from asymptomatic to symptomatic type 1 diabetes (T1D). RESEARCH DESIGN AND METHODS: Next-generation targeted sequencing was used to genotype eight HLAII genes (DQA1, DQB1, DRB1, DRB3, DRB4, DRB5, DPA1, DPB1) in 1,216 participants from the Diabetes Prevention Trial-1 and Randomized Diabetes Prevention Trial with Oral Insulin sponsored by TrialNet. By the linkage disequilibrium, DQA1 and DQB1 are haplotyped to form DQ haplotypes; DP and DR haplotypes are similarly constructed. Together with available clinical covariables, we applied the Cox regression model to assess HLAII immunogenic associations with the disease progression. RESULTS: First, the current investigation updated the previously reported genetic associations of DQA1*03:01-DQB1*03:02 (hazard ratio [HR] = 1.25, P = 3.50*10-3) and DQA1*03:03-DQB1*03:01 (HR = 0.56, P = 1.16*10-3), and also uncovered a risk association with DQA1*05:01-DQB1*02:01 (HR = 1.19, P = 0.041). Second, after adjusting for DQ, DPA1*02:01-DPB1*11:01 and DPA1*01:03-DPB1*03:01 were found to have opposite associations with progression (HR = 1.98 and 0.70, P = 0.021 and 6.16*10-3, respectively). Third, DRB1*03:01-DRB3*01:01 and DRB1*03:01-DRB3*02:02, sharing the DRB1*03:01, had opposite associations (HR = 0.73 and 1.44, P = 0.04 and 0.019, respectively), indicating a role of DRB3. Meanwhile, DRB1*12:01-DRB3*02:02 and DRB1*01:03 alone were found to associate with progression (HR = 2.6 and 2.32, P = 0.018 and 0.039, respectively). Fourth, through enumerating all heterodimers, it was found that both DQ and DP could exhibit associations with disease progression. CONCLUSIONS: These results suggest that HLAII polymorphisms influence progression from islet autoimmunity to T1D among at-risk subjects with islet autoantibodies.

Distinct transcriptomic profiles in children prior to the appearance of type 1 diabetes-linked islet autoantibodies and following enterovirus infection.

Although the genetic basis and pathogenesis of type 1 diabetes have been studied extensively, how host responses to environmental factors might contribute to autoantibody development remains largely unknown. Here, we use longitudinal blood transcriptome sequencing data to characterize host responses in children within 12 months prior to the appearance of type 1 diabetes-linked islet autoantibodies, as well as matched control children. We report that children who present with insulin-specific autoantibodies first have distinct transcriptional profiles from those who develop GADA autoantibodies first. In particular, gene dosage-driven expression of GSTM1 is associated with GADA autoantibody positivity. Moreover, compared with controls, we observe increased monocyte and decreased B cell proportions 9-12 months prior to autoantibody positivity, especially in children who developed antibodies against insulin first. Lastly, we show that control children present transcriptional signatures consistent with robust immune responses to enterovirus infection, whereas children who later developed islet autoimmunity do not. These findings highlight distinct immune-related transcriptomic differences between case and control children prior to case progression to islet autoimmunity and uncover deficient antiviral response in children who later develop islet autoimmunity.

Type 1 Diabetes Risk Phenotypes Using Cluster Analysis.

Background: Although statistical models for predicting type 1 diabetes risk have been developed, approaches that reveal clinically meaningful clusters in the at-risk population and allow for non-linear relationships between predictors are lacking. We aimed to identify and characterize clusters of islet autoantibody-positive individuals that share similar characteristics and type 1 diabetes risk. Methods: We tested a novel outcome-guided clustering method in initially non-diabetic autoantibody-positive relatives of individuals with type 1 diabetes, using the TrialNet Pathway to Prevention (PTP) study data (n=1127). The outcome of the analysis was time to type 1 diabetes and variables in the model included demographics, genetics, metabolic factors and islet autoantibodies. An independent dataset (Diabetes Prevention Trial of Type 1 Diabetes, DPT-1 study) (n=704) was used for validation. Findings: The analysis revealed 8 clusters with varying type 1 diabetes risks, categorized into three groups. Group A had three clusters with high glucose levels and high risk. Group B included four clusters with elevated autoantibody titers. Group C had three lower-risk clusters with lower autoantibody titers and glucose levels. Within the groups, the clusters exhibit variations in characteristics such as glucose levels, C-peptide levels, age, and genetic risk. A decision rule for assigning individuals to clusters was developed. The validation dataset confirms that the clusters can identify individuals with similar characteristics. Interpretation: Demographic, metabolic, immunological, and genetic markers can be used to identify clusters of distinctive characteristics and different risks of progression to type 1 diabetes among autoantibody-positive individuals with a family history of type 1 diabetes. The results also revealed the heterogeneity in the population and complex interactions between variables.

ISLET: individual-specific reference panel recovery improves cell-type-specific inference.

We propose a statistical framework ISLET to infer individual-specific and cell-type-specific transcriptome reference panels. ISLET models the repeatedly measured bulk gene expression data, to optimize the usage of shared information within each subject. ISLET is the first available method to achieve individual-specific reference estimation in repeated samples. Using simulation studies, we show outstanding performance of ISLET in the reference estimation and downstream cell-type-specific differentially expressed genes testing. We apply ISLET to longitudinal transcriptomes profiled from blood samples in a large observational study of young children and confirm the cell-type-specific gene signatures for pancreatic islet autoantibody. ISLET is available at https://bioconductor.org/packages/ISLET .

Epigenetic Age in Peripheral Blood Among Children, Adolescent, and Adult Survivors of Childhood Cancer.

Importance: Certain cancer therapies are risk factors for epigenetic age acceleration (EAA) among survivors of childhood cancer, and EAA is associated with chronic health conditions (CHCs). However, small numbers of younger survivors (aged <20 years) previously evaluated have limited the ability to calculate EAA among this age group. Objective: To evaluate the change rate of epigenetic age (EA) and EAA in younger compared with older survivors and the possible association of EAA with early-onset obesity (aged <20 years), severity/burden of CHCs, and late mortality (>5 years from cancer diagnosis). Design, Setting, and Participants: Study participants were from the St Jude Lifetime Cohort, initiated in 2007 with ongoing follow-up. The present study was conducted from April 17, 2022, to March 23, 2023. Survivors in this cohort of European ancestry with DNA methylation data were included. Cross-sectional annual changes in EA and EAA were compared across 5 different chronologic age groups: age 0 to 9 (children), 10 to 19 (adolescents), 20 to 34 (younger adults), 35 to 49 (middle-aged adults), and greater than or equal to 50 (older adults) years. Logistic regression evaluated the association between EAA and early-onset obesity or severity/burden of CHCs. Cox proportional hazards regression assessed the association between EAA and late mortality. Main Outcomes and Measures: Early-onset obesity, severity/burden of CHCs (graded using the Common Terminology Criteria for Adverse Events (grade 1, mild; 2, moderate; 3, severe/disabling; 4, life-threatening) and were combined into high vs low severity/burden based on frequency and grade), and late mortality were the outcomes based on follow-up until April 2020. Expanded DNA methylation profiling increased the number of survivors younger than 20 years (n = 690). Epigenetic age was calculated primarily using the Levine clock, and EAA was derived from least squares regression of EA against chronologic age and was standardized to a z score (Levine EEA). Results: Among 2846 participants (median age, 30.3 [IQR, 9.3-41.5] years; 53% males), the cross-sectional annual change in EA_Levine was higher in children (1.63 years) and adolescents (1.14 years), and the adjusted least-squares mean of Levine EEA was lower in children (-0.22 years) and older adults (-1.70 years). Each 1-SD increase in Levine EEA was associated with increased risk of developing early-onset obesity (odds ratio [OR], 1.46; 95% CI, 1.19-1.78), high severity/burden of CHCs (OR, 1.13; 95% CI, 1.03-1.24), and late mortality (hazard ratio, 1.75; 95% CI, 1.35-2.26). Conclusions and Relevance: The findings of this study suggest that EAA measured in children and adolescent survivors of childhood cancer is associated with early-onset obesity, severity/burden of all CHCs, and late mortality. Evaluating EAA may help identify survivors of childhood cancer at increased risk for early-onset obesity, morbidity in general, and mortality.

Interaction Between Dietary Iron Intake and Genetically Determined Iron Overload: Risk of Islet Autoimmunity and Progression to Type 1 Diabetes in the TEDDY Study.

OBJECTIVE: To examine whether iron intake and genetically determined iron overload interact in predisposing to the development of childhood islet autoimmunity (IA) and type 1 diabetes (T1D). RESEARCH DESIGN AND METHODS: In The Environmental Determinants of Diabetes in the Young (TEDDY) study, 7,770 genetically high-risk children were followed from birth until the development of IA and progression to T1D. Exposures included energy-adjusted iron intake in the first 3 years of life and a genetic risk score (GRS) for increased circulating iron. RESULTS: We found a U-shaped association between iron intake and risk of GAD antibody as the first autoantibody. In children with GRS ≥2 iron risk alleles, high iron intake was associated with an increased risk of IA, with insulin as first autoantibody (adjusted hazard ratio 1.71 [95% CI 1.14; 2.58]) compared with moderate iron intake. CONCLUSIONS: Iron intake may alter the risk of IA in children with high-risk HLA haplogenotypes.

Data Mining Framework for Discovering and Clustering Phenotypes of Atypical Diabetes.

CONTEXT: Some individuals present with forms of diabetes that are "atypical" (AD), which do not conform to typical features of either type 1 diabetes (T1D) or type 2 diabetes (T2D). These forms of AD display a range of phenotypic characteristics that likely reflect different endotypes based on unique etiologies or pathogenic processes. OBJECTIVE: To develop an analytical approach to identify and cluster phenotypes of AD. METHODS: We developed Discover Atypical Diabetes (DiscoverAD), a data mining framework, to identify and cluster phenotypes of AD. DiscoverAD was trained against characteristics of manually classified patients with AD among 278 adults with diabetes within the Cameron County Hispanic Cohort (CCHC) (Study A). We then tested DiscoverAD in a separate population of 758 multiethnic children with T1D within the Texas Children's Hospital Registry for New-Onset Type 1 Diabetes (TCHRNO-1) (Study B). RESULTS: We identified an AD frequency of 11.5% in the CCHC (Study A) and 5.3% in the pediatric TCHRNO-1 (Study B). Cluster analysis identified 4 distinct groups of AD in Study A: cluster 1, positive for the 65 kDa glutamate decarboxylase autoantibody (GAD65Ab), adult-onset, long disease duration, preserved beta-cell function, no insulin treatment; cluster 2, GAD65Ab negative, diagnosed at age ≤21 years; cluster 3, GAD65Ab negative, adult-onset, poor beta-cell function, lacking central obesity; cluster 4, diabetic ketoacidosis (DKA)-prone participants lacking a typical T1D phenotype. Applying DiscoverAD to the pediatric patients with T1D in Study B revealed 2 distinct groups of AD: cluster 1, autoantibody negative, poor beta-cell function, lower body mass index (BMI); cluster 2, autoantibody positive, higher BMI, higher incidence of DKA. CONCLUSION: DiscoverAD can be adapted to different datasets to identify and define phenotypes of participants with AD based on available clinical variables.

High-throughput muscle fiber typing from RNA sequencing data.

BACKGROUND: Skeletal muscle fiber type distribution has implications for human health, muscle function, and performance. This knowledge has been gathered using labor-intensive and costly methodology that limited these studies. Here, we present a method based on muscle tissue RNA sequencing data (totRNAseq) to estimate the distribution of skeletal muscle fiber types from frozen human samples, allowing for a larger number of individuals to be tested. METHODS: By using single-nuclei RNA sequencing (snRNAseq) data as a reference, cluster expression signatures were produced by averaging gene expression of cluster gene markers and then applying these to totRNAseq data and inferring muscle fiber nuclei type via linear matrix decomposition. This estimate was then compared with fiber type distribution measured by ATPase staining or myosin heavy chain protein isoform distribution of 62 muscle samples in two independent cohorts (n = 39 and 22). RESULTS: The correlation between the sequencing-based method and the other two were rATPas = 0.44 [0.13-0.67], [95% CI], and rmyosin = 0.83 [0.61-0.93], with p = 5.70 × 10-3 and 2.00 × 10-6, respectively. The deconvolution inference of fiber type composition was accurate even for very low totRNAseq sequencing depths, i.e., down to an average of ~ 10,000 paired-end reads. CONCLUSIONS: This new method ( https://github.com/OlaHanssonLab/PredictFiberType ) consequently allows for measurement of fiber type distribution of a larger number of samples using totRNAseq in a cost and labor-efficient way. It is now feasible to study the association between fiber type distribution and e.g. health outcomes in large well-powered studies.

Telomere length is not a main factor for the development of islet autoimmunity and type 1 diabetes in the TEDDY study.

The Environmental Determinants of Diabetes in the Young (TEDDY) study enrolled 8676 children, 3-4 months of age, born with HLA-susceptibility genotypes for islet autoimmunity (IA) and type 1 diabetes (T1D). Whole-genome sequencing (WGS) was performed in 1119 children in a nested case-control study design. Telomere length was estimated from WGS data using five tools: Computel, Telseq, Telomerecat, qMotif and Motif_counter. The estimated median telomere length was 5.10 kb (IQR 4.52-5.68 kb) using Computel. The age when the blood sample was drawn had a significant negative correlation with telomere length (P = 0.003). European children, particularly those from Finland (P = 0.041) and from Sweden (P = 0.001), had shorter telomeres than children from the U.S.A. Paternal age (P = 0.019) was positively associated with telomere length. First-degree relative status, presence of gestational diabetes in the mother, and maternal age did not have a significant impact on estimated telomere length. HLA-DR4/4 or HLA-DR4/X children had significantly longer telomeres compared to children with HLA-DR3/3 or HLA-DR3/9 haplogenotypes (P = 0.008). Estimated telomere length was not significantly different with respect to any IA (P = 0.377), IAA-first (P = 0.248), GADA-first (P = 0.248) or T1D (P = 0.861). These results suggest that telomere length has no major impact on the risk for IA, the first step to develop T1D. Nevertheless, telomere length was shorter in the T1D high prevalence populations, Finland and Sweden.

Transcriptional networks in at-risk individuals identify signatures of type 1 diabetes progression.

Type 1 diabetes (T1D) is a disease of insulin deficiency that results from autoimmune destruction of pancreatic islet ß cells. The exact cause of T1D remains unknown, although asymptomatic islet autoimmunity lasting from weeks to years before diagnosis raises the possibility of intervention before the onset of clinical disease. The number, type, and titer of islet autoantibodies are associated with long-term disease risk but do not cause disease, and robust early predictors of individual progression to T1D onset remain elusive. The Environmental Determinants of Diabetes in the Young (TEDDY) consortium is a prospective cohort study aiming to determine genetic and environmental interactions causing T1D. Here, we analyzed longitudinal blood transcriptomes of 2013 samples from 400 individuals in the TEDDY study before both T1D and islet autoimmunity. We identified and interpreted age-associated gene expression changes in healthy infancy and age-independent changes tracking with progression to both T1D and islet autoimmunity, beginning before other evidence of islet autoimmunity was present. We combined multivariate longitudinal data in a Bayesian joint model to predict individual risk of T1D onset and validated the association of a natural killer cell signature with progression and the model's predictive performance on an additional 356 samples from 56 individuals in the independent Type 1 Diabetes Prediction and Prevention study. Together, our results indicate that T1D is characterized by early and longitudinal changes in gene expression, informing the immunopathology of disease progression and facilitating prediction of its course.

Children's erythrocyte fatty acids are associated with the risk of islet autoimmunity.

Our aim was to investigate the associations between erythrocyte fatty acids and the risk of islet autoimmunity in children. The Environmental Determinants of Diabetes in the Young Study (TEDDY) is a longitudinal cohort study of children at high genetic risk for type 1 diabetes (n = 8676) born between 2004 and 2010 in the U.S., Finland, Sweden, and Germany. A nested case-control design comprised 398 cases with islet autoimmunity and 1178 sero-negative controls matched for clinical site, family history, and gender. Fatty acids composition was measured in erythrocytes collected at the age of 3, 6, and 12 months and then annually up to 6 years of age. Conditional logistic regression models were adjusted for HLA risk genotype, ancestry, and weight z-score. Higher eicosapentaenoic and docosapentaenoic acid (n - 3 polyunsaturated fatty acids) levels during infancy and conjugated linoleic acid after infancy were associated with a lower risk of islet autoimmunity. Furthermore, higher levels of some even-chain saturated (SFA) and monounsaturated fatty acids (MUFA) were associated with increased risk. Fatty acid status in early life may signal the risk for islet autoimmunity, especially n - 3 fatty acids may be protective, while increased levels of some SFAs and MUFAs may precede islet autoimmunity.

Relationship between insulin sensitivity and gene expression in human skeletal muscle.

BACKGROUND: Insulin resistance (IR) in skeletal muscle is a key feature of the pre-diabetic state, hypertension, dyslipidemia, cardiovascular diseases and also predicts type 2 diabetes. However, the underlying molecular mechanisms are still poorly understood. METHODS: To explore these mechanisms, we related global skeletal muscle gene expression profiling of 38 non-diabetic men to a surrogate measure of insulin sensitivity, i.e. homeostatic model assessment of insulin resistance (HOMA-IR). RESULTS: We identified 70 genes positively and 110 genes inversely correlated with insulin sensitivity in human skeletal muscle, identifying autophagy-related genes as positively correlated with insulin sensitivity. Replication in an independent study of 9 non-diabetic men resulted in 10 overlapping genes that strongly correlated with insulin sensitivity, including SIRT2, involved in lipid metabolism, and FBXW5 that regulates mammalian target-of-rapamycin (mTOR) and autophagy. The expressions of SIRT2 and FBXW5 were also positively correlated with the expression of key genes promoting the phenotype of an insulin sensitive myocyte e.g. PPARGC1A. CONCLUSIONS: The muscle expression of 180 genes were correlated with insulin sensitivity. These data suggest that activation of genes involved in lipid metabolism, e.g. SIRT2, and genes regulating autophagy and mTOR signaling, e.g. FBXW5, are associated with increased insulin sensitivity in human skeletal muscle, reflecting a highly flexible nutrient sensing.

First-appearing islet autoantibodies for type 1 diabetes in young children: maternal life events during pregnancy and the child's genetic risk.

AIMS/HYPOTHESIS: Psychological stress has long been considered a possible trigger of type 1 diabetes, although prospective studies examining the link between psychological stress or life events during pregnancy and the child's type 1 diabetes risk are rare. The objective of this study was to examine the association between life events during pregnancy and first-appearing islet autoantibodies (IA) in young children, conditioned by the child's type 1 diabetes-related genetic risk. METHODS: The IA status of 7317 genetically at-risk The Environmental Determinants of Diabetes in the Young (TEDDY) participants was assessed every 3 months from 3 months to 4 years, and bi-annually thereafter. Reports of major life events during pregnancy were collected at study inception when the child was 3 months of age and placed into one of six categories. Life events during pregnancy were examined for association with first-appearing insulin (IAA) (N = 222) or GAD (GADA) (N = 209) autoantibodies in the child until 6 years of age using proportional hazard models. Relative excess risk due to interaction (RERI) by the child's HLA-DR and SNP profile was estimated. RESULTS: Overall, 65% of mothers reported a life event during pregnancy; disease/injury (25%), serious interpersonal (28%) and job-related (25%) life events were most common. The association of life events during pregnancy differed between IAA and GADA as the first-appearing autoantibody. Serious interpersonal life events correlated with increased risk of GADA-first only in HLA-DR3 children with the BACH2-T allele (HR 2.28, p < 0.0001), an additive interaction (RERI 1.87, p = 0.0004). Job-related life events were also associated with increased risk of GADA-first among HLA-DR3/4 children (HR 1.53, p = 0.04) independent of serious interpersonal life events (HR 1.90, p = 0.002), an additive interaction (RERI 1.19, p = 0.004). Job-related life events correlated with reduced risk of IAA-first (HR 0.55, p = 0.004), particularly in children with the BTNL2-GG allele (HR 0.48; 95% CI 0.31, 0.76). CONCLUSIONS/INTERPRETATION: Specific life events during pregnancy are differentially related to IAA vs GADA as first-appearing IA and interact with different HLA and non-HLA genetic factors, supporting the concept of different endotypes underlying type 1 diabetes. However, the mechanisms underlying these associations remain to be discovered. Life events may be markers for other yet-to-be-identified factors important to the development of first-appearing IA.

Longitudinal Metabolome-Wide Signals Prior to the Appearance of a First Islet Autoantibody in Children Participating in the TEDDY Study.

Children at increased genetic risk for type 1 diabetes (T1D) after environmental exposures may develop pancreatic islet autoantibodies (IA) at a very young age. Metabolic profile changes over time may imply responses to exposures and signal development of the first IA. Our present research in The Environmental Determinants of Diabetes in the Young (TEDDY) study aimed to identify metabolome-wide signals preceding the first IA against GAD (GADA-first) or against insulin (IAA-first). We profiled metabolomes by mass spectrometry from children's plasma at 3-month intervals after birth until appearance of the first IA. A trajectory analysis discovered each first IA preceded by reduced amino acid proline and branched-chain amino acids (BCAAs), respectively. With independent time point analysis following birth, we discovered dehydroascorbic acid (DHAA) contributing to the risk of each first IA, and γ-aminobutyric acid (GABAs) associated with the first autoantibody against insulin (IAA-first). Methionine and alanine, compounds produced in BCAA metabolism and fatty acids, also preceded IA at different time points. Unsaturated triglycerides and phosphatidylethanolamines decreased in abundance before appearance of either autoantibody. Our findings suggest that IAA-first and GADA-first are heralded by different patterns of DHAA, GABA, multiple amino acids, and fatty acids, which may be important to primary prevention of T1D.

Plasma ascorbic acid and the risk of islet autoimmunity and type 1 diabetes: the TEDDY study.

AIMS/HYPOTHESIS: We studied the association of plasma ascorbic acid with the risk of developing islet autoimmunity and type 1 diabetes and examined whether SNPs in vitamin C transport genes modify these associations. Furthermore, we aimed to determine whether the SNPs themselves are associated with the risk of islet autoimmunity or type 1 diabetes. METHODS: We used a risk set sampled nested case-control design within an ongoing international multicentre observational study: The Environmental Determinants of Diabetes in the Young (TEDDY). The TEDDY study followed children with increased genetic risk from birth to endpoints of islet autoantibodies (350 cases, 974 controls) and type 1 diabetes (102 cases, 282 controls) in six clinical centres. Control participants were matched for family history of type 1 diabetes, clinical centre and sex. Plasma ascorbic acid concentration was measured at ages 6 and 12 months and then annually up to age 6 years. SNPs in vitamin C transport genes were genotyped using the ImmunoChip custom microarray. Comparisons were adjusted for HLA genotypes and for background population stratification. RESULTS: Childhood plasma ascorbic acid (mean ± SD 10.76 ± 3.54 mg/l in controls) was inversely associated with islet autoimmunity risk (adjusted OR 0.96 [95% CI 0.92, 0.99] per +1 mg/l), particularly islet autoimmunity, starting with insulin autoantibodies (OR 0.94 [95% CI 0.88, 0.99]), but not with type 1 diabetes risk (OR 0.93 [95% Cl 0.86, 1.02]). The SLC2A2 rs5400 SNP was associated with increased risk of type 1 diabetes (OR 1.77 [95% CI 1.12, 2.80]), independent of plasma ascorbic acid (OR 0.92 [95% CI 0.84, 1.00]). CONCLUSIONS/INTERPRETATION: Higher plasma ascorbic acid levels may protect against islet autoimmunity in children genetically at risk for type 1 diabetes. Further studies are warranted to confirm these findings. DATA AVAILABILITY: The datasets generated and analysed during the current study will be made available in the NIDDK Central Repository at https://www.niddkrepository.org/studies/teddy.

An open resource for accurately benchmarking small variant and reference calls.

Benchmark small variant calls are required for developing, optimizing and assessing the performance of sequencing and bioinformatics methods. Here, as part of the Genome in a Bottle (GIAB) Consortium, we apply a reproducible, cloud-based pipeline to integrate multiple short- and linked-read sequencing datasets and provide benchmark calls for human genomes. We generate benchmark calls for one previously analyzed GIAB sample, as well as six genomes from the Personal Genome Project. These new genomes have broad, open consent, making this a 'first of its kind' resource that is available to the community for multiple downstream applications. We produce 17% more benchmark single nucleotide variations, 176% more indels and 12% larger benchmark regions than previously published GIAB benchmarks. We demonstrate that this benchmark reliably identifies errors in existing callsets and highlight challenges in interpreting performance metrics when using benchmarks that are not perfect or comprehensive. Finally, we identify strengths and weaknesses of callsets by stratifying performance according to variant type and genome context.

Serum biomarkers of glucocorticoid response and safety in anti-neutrophil cytoplasmic antibody-associated vasculitis and juvenile dermatomyositis.

Glucocorticoids are standard of care for many chronic inflammatory conditions, including juvenile dermatomyositis (JDM) and anti-neutrophil cytoplasmic antibody (ANCA)-associated vasculitis (AAV). We sought to define pharmacodynamic biomarkers of therapeutic efficacy and safety concerns of glucocorticoid treatment for these two disorders. Previous proteomic profiling of patients with Duchenne muscular dystrophy (DMD) and inflammatory bowel disease (IBD) treated with glucocorticoids identified candidate biomarkers for efficacy and safety concerns of glucocorticoids. Serial serum samples from patients with AAV (n = 30) and JDM (n = 12) were obtained during active disease, and after treatment with glucocorticoids. For AAV, 8 of 11 biomarkers of the anti-inflammatory response to glucocorticoids were validated (P-value ≤0.05; CD23, macrophage-derived cytokine, interleukin-22 binding protein, matrix metalloproteinase-12, T lymphocyte surface antigen Ly9, fibrinogen gamma chain, angiopoietin-2 [all decreased], and protein C [increased]), as were 5 of 7 safety biomarkers (P-value ≤0.05; afamin, matrix metalloproteinase-3, insulin growth factor binding protein-5, angiotensinogen, leptin [all increased]). For JDM, 10 of 11 efficacy biomarkers were validated (P-value ≤0.05; all proteins except fibrinogen gamma chain) and 6 of 7 safety biomarkers (P-value ≤0.05; AAV proteins plus growth hormone binding protein). The identified efficacy biomarkers may be useful as objective outcome measures for early phase proof-of-concept studies when assessing novel anti-inflammatory drugs in JDM and AAV, and likely in other inflammatory disorders. Similarly, safety biomarkers may also be helpful assessing toxicity of alternatives to glucocorticoids.

Metabolic pathways and immunometabolism in rare kidney diseases.

OBJECTIVES: To characterise renal tissue metabolic pathway gene expression in different forms of glomerulonephritis. METHODS: Patients with nephrotic syndrome (NS), antineutrophil cytoplasmic antibody-associated vasculitis (AAV), systemic lupus erythematosus (SLE) and healthy living donors (LD) were studied. Clinically indicated renal biopsies were obtained at time of diagnosis and microdissected into glomerular and tubulointerstitial compartments. Microarray-derived differential gene expression of 88 genes representing critical enzymes of metabolic pathways and 25 genes related to immune cell markers was compared between disease groups. Correlation analyses measured relationships between metabolic pathways, kidney function and cytokine production. RESULTS: Reduced steady state levels of mRNA species were enriched in pathways of oxidative phosphorylation and increased in the pentose phosphate pathway (PPP) with maximal perturbation in AAV and SLE followed by NS, and least in LD. Transcript regulation was isozymes specific with robust regulation in hexokinases, enolases and glucose transporters. Intercorrelation networks were observed between enzymes of the PPP (eg, transketolase) and macrophage markers (eg, CD68) (r=0.49, p<0.01). Increased PPP transcript levels were associated with reduced glomerular filtration rate in the glomerular (r=-0.49, p<0.01) and tubulointerstitial (r=-0.41, p<0.01) compartments. PPP expression and tumour necrosis factor activation were tightly co-expressed (r=0.70, p<0.01). CONCLUSION: This study demonstrated concordant alterations of the renal transcriptome consistent with metabolic reprogramming across different forms of glomerulonephritis. Activation of the PPP was tightly linked with intrarenal macrophage marker expression, reduced kidney function and increased production of cytokines. Modulation of glucose metabolism may offer novel immune-modulatory therapeutic approaches in rare kidney diseases.

Characterising cis-regulatory variation in the transcriptome of histologically normal and tumour-derived pancreatic tissues.

OBJECTIVE: To elucidate the genetic architecture of gene expression in pancreatic tissues. DESIGN: We performed expression quantitative trait locus (eQTL) analysis in histologically normal pancreatic tissue samples (n=95) using RNA sequencing and the corresponding 1000 genomes imputed germline genotypes. Data from pancreatic tumour-derived tissue samples (n=115) from The Cancer Genome Atlas were included for comparison. RESULTS: We identified 38 615 cis-eQTLs (in 484 genes) in histologically normal tissues and 39 713 cis-eQTL (in 237 genes) in tumour-derived tissues (false discovery rate <0.1), with the strongest effects seen near transcriptional start sites. Approximately 23% and 42% of genes with significant cis-eQTLs appeared to be specific for tumour-derived and normal-derived tissues, respectively. Significant enrichment of cis-eQTL variants was noted in non-coding regulatory regions, in particular for pancreatic tissues (1.53-fold to 3.12-fold, p≤0.0001), indicating tissue-specific functional relevance. A common pancreatic cancer risk locus on 9q34.2 (rs687289) was associated with ABO expression in histologically normal (p=5.8×10-8) and tumour-derived (p=8.3×10-5) tissues. The high linkage disequilibrium between this variant and the O blood group generating deletion variant in ABO (exon 6) suggested that nonsense-mediated decay (NMD) of the 'O' mRNA might explain this finding. However, knockdown of crucial NMD regulators did not influence decay of the ABO 'O' mRNA, indicating that a gene regulatory element influenced by pancreatic cancer risk alleles may underlie the eQTL. CONCLUSIONS: We have identified cis-eQTLs representing potential functional regulatory variants in the pancreas and generated a rich data set for further studies on gene expression and its regulation in pancreatic tissues.