Reversion of islet autoimmunity (IA) may point to mechanisms that prevent IA progression. We followed 199 individuals who developed IA during the Diabetes Autoimmunity Study in the Young. Untargeted metabolomics was performed in serum samples following IA. Cox proportional hazards models were used to test whether the metabolites (2,487) predicted IA reversion: two or more consecutive visits negative for all autoantibodies. We conducted a principal components analysis (PCA) of the top metabolites; |hazard ratio (HR) >1.25| and nominal P < 0.01. Phosphatidylcholine (16:0_18:1(9Z)) was the strongest individual metabolite (HR per 1 SD 2.16, false discovery rate (FDR)-adjusted P = 0.0037). Enrichment analysis identified four clusters (FDR P < 0.10) characterized by an overabundance of sphingomyelin (d40:0), phosphatidylcholine (16:0_18:1(9Z)), phosphatidylcholine (30:0), and l-decanoylcarnitine. Overall, 63 metabolites met the criteria for inclusion in the PCA. PC1 (HR 1.4, P < 0.0001), PC2 (HR 0.85, P = 0.0185), and PC4 (HR 1.28, P = 0.0103) were associated with IA reversion. Given the potential influence of diet on the metabolome, we investigated whether nutrients were correlated with PCs. We identified 20 nutrients that were correlated with the PCs (P < 0.05). Total sugar intake was the top nutrient. Overall, we identified an association between phosphatidylcholine, sphingomyelin, and carnitine levels and reversion of IA.

BACKGROUND: Type 1 diabetes (T1D) is preceded by a heterogenous pre-clinical phase, islet autoimmunity (IA). We aimed to identify pre vs. post-IA seroconversion (SV) changes in DNAm that differed across three IA progression phenotypes, those who lose autoantibodies (reverters), progress to clinical T1D (progressors), or maintain autoantibody levels (maintainers). METHODS: This epigenome-wide association study (EWAS) included longitudinal DNAm measurements in blood (Illumina 450K and EPIC) from participants in Diabetes Autoimmunity Study in the Young (DAISY) who developed IA, one or more islet autoantibodies on at least two consecutive visits. We compared reverters - individuals who sero-reverted, negative for all autoantibodies on at least two consecutive visits and did not develop T1D (n=41); maintainers - continued to test positive for autoantibodies but did not develop T1D (n=60); progressors - developed clinical T1D (n=42). DNAm data were measured before (pre-SV visit) and after IA (post-SV visit). Linear mixed models were used to test for differences in pre- vs post-SV changes in DNAm across the three groups. Linear mixed models were also used to test for group differences in average DNAm. Cell proportions, age, and sex were adjusted for in all models. Median follow-up across all participants was 15.5 yrs. (interquartile range (IQR): 10.8-18.7). RESULTS: The median age at the pre-SV visit was 2.2 yrs. (IQR: 0.8-5.3) in progressors, compared to 6.0 yrs. (IQR: 1.3-8.4) in reverters, and 5.7 yrs. (IQR: 1.4-9.7) in maintainers. Median time between the visits was similar in reverters 1.4 yrs. (IQR: 1-1.9), maintainers 1.3 yrs. (IQR: 1.0-2.0), and progressors 1.8 yrs. (IQR: 1.0-2.0). Changes in DNAm, pre- vs post-SV, differed across the groups at one site (cg16066195) and 11 regions. Average DNAm (mean of pre- and post-SV) differed across 22 regions. CONCLUSION: Differentially changing DNAm regions were located in genomic areas related to beta cell function, immune cell differentiation, and immune cell function.

In our recent article [1], it has come to our attention that the sample labels are not consistent between Table 1, the data labels deposited in the Sequence Read Archive, and Additional file 1: Table S2. We are therefore providing an updated Additional file 1: Table S2 so identical samples now have the same label.

Background

Respiratory illness caused by viral infection is associated with the development and exacerbation of childhood asthma. Little is known about the effects of respiratory viral infections in the absence of illness. Using quantitative PCR (qPCR) for common respiratory viruses and for two genes known to be highly upregulated in viral infections (CCL8/CXCL11), we screened 92 asthmatic and 69 healthy children without illness for respiratory virus infections.

Results

We found 21 viral qPCR-positive and 2 suspected virus-infected subjects with high expression of CCL8/CXCL11. We applied a dual RNA-seq workflow to these subjects, together with 25 viral qPCR-negative subjects, to compare qPCR with sequencing-based virus detection and to generate the airway transcriptome for analysis. RNA-seq virus detection achieved 86% sensitivity when compared to qPCR-based screening. We detected additional respiratory viruses in the two CCL8/CXCL11-high subjects and in two of the qPCR-negative subjects. Viral read counts varied widely and were used to stratify subjects into Virus-High and Virus-Low groups. Examination of the host airway transcriptome found that the Virus-High group was characterized by immune cell airway infiltration, downregulation of cilia genes, and dampening of type 2 inflammation. Even the Virus-Low group was differentiated from the No-Virus group by 100 genes, some involved in eIF2 signaling.

Conclusions

Respiratory virus infection without illness is not innocuous but may determine the airway function of these subjects by driving immune cell airway infiltration, cellular remodeling, and alteration of asthmogenic gene expression.

Type 2 diabetes (T2D) is a heterogeneous disorder with contributions from peripheral insulin resistance and β-cell dysfunction. For minimization of phenotypic heterogeneity, quantitative intermediate phenotypes characterizing basal glucose homeostasis (insulin resistance and HOMA of insulin resistance [HOMAIR] and of β-cell function [HOMAB]) have shown promise in relatively large samples. We investigated the utility of dynamic measures of glucose homeostasis (insulin sensitivity [SI] and acute insulin response [AIRg]) evaluating T2D-susceptibility variants (n = 57) in Hispanic Americans from the GUARDIAN Consortium (n = 2,560). Basal and dynamic measures were genetically correlated (HOMAB-AIRg: ρG = 0.28-0.73; HOMAIR-SI: ρG = -0.73 to -0.83) with increased heritability for the dynamic measure AIRg Significant association of variants with dynamic measures (P < 8.77 × 10(-4)) was observed. A pattern of superior performance of AIRg was observed for well-established loci including MTNR1B (P = 9.46 × 10(-12)), KCNQ1 (P = 1.35 × 10(-4)), and TCF7L2 (P = 5.10 × 10(-4)) with study-wise statistical significance. Notably, significant association of MTNR1B with AIRg (P < 1.38 × 10(-9)) was observed in a population one-fourteenth the size of the initial discovery cohort. These observations suggest that basal and dynamic measures provide different views and levels of sensitivity to discrete elements of glucose homeostasis. Although more costly to obtain, dynamic measures yield significant results that could be considered physiologically "closer" to causal pathways and provide insight into the discrete mechanisms of action.

Unlabelled

Blood transcriptional signatures are promising for tuberculosis (TB) diagnosis but have not been evaluated among U.S.

Patients

To be used clinically, transcriptional classifiers need reproducible accuracy in diverse populations that vary in genetic composition, disease spectrum and severity, and comorbidities. In a prospective case-control study, we identified novel transcriptional classifiers for active TB among U.S. patients and systematically compared their accuracy to classifiers from published studies. Blood samples from HIV-uninfected U.S. adults with active TB, pneumonia, or latent TB infection underwent whole-transcriptome microarray. We used support vector machines to classify disease state based on transcriptional patterns. We externally validated our classifiers using data from sub-Saharan African cohorts and evaluated previously published transcriptional classifiers in our population. Our classifier distinguishing active TB from pneumonia had an area under the concentration-time curve (AUC) of 96.5% (95.4% to 97.6%) among U.S. patients, but the AUC was lower (90.6% [89.6% to 91.7%]) in HIV-uninfected Sub-Saharan Africans. Previously published comparable classifiers had AUC values of 90.0% (87.7% to 92.3%) and 82.9% (80.8% to 85.1%) when tested in U.S.

Patients

Our classifier distinguishing active TB from latent TB had AUC values of 95.9% (95.2% to 96.6%) among U.S. patients and 95.3% (94.7% to 96.0%) among Sub-Saharan Africans. Previously published comparable classifiers had AUC values of 98.0% (97.4% to 98.7%) and 94.8% (92.9% to 96.8%) when tested in U.S.

Patients

Blood transcriptional classifiers accurately detected active TB among U.S. adults. The accuracy of classifiers for active TB versus that of other diseases decreased when tested in new populations with different disease controls, suggesting additional studies are required to enhance generalizability. Classifiers that distinguish active TB from latent TB are accurate and generalizable across populations and can be explored as screening assays.

Background

It is unknown whether immunosuppression influences the physiologic state of Mycobacterium tuberculosis in vivo. We evaluated the impact of host immunity by comparing M. tuberculosis and human gene transcription in sputum between human immunodeficiency virus (HIV)-infected and uninfected patients with tuberculosis.

Methods

We collected sputum specimens before treatment from Gambians and Ugandans with pulmonary tuberculosis, revealed by positive results of acid-fast bacillus smears. We quantified expression of 2179 M. tuberculosis genes and 234 human immune genes via quantitative reverse transcription-polymerase chain reaction. We summarized genes from key functional categories with significantly increased or decreased expression.

Results

A total of 24 of 65 patients with tuberculosis were HIV infected. M. tuberculosis DosR regulon genes were less highly expressed among HIV-infected patients with tuberculosis than among HIV-uninfected patients with tuberculosis (Gambia, P < .0001; Uganda, P = .037). In profiling of human genes from the same sputa, HIV-infected patients had 3.4-fold lower expression of IFNG (P = .005), 4.9-fold higher expression of ARG1 (P = .0006), and 3.4-fold higher expression of IL10 (P = .0002) than in HIV-uninfected patients with tuberculosis.

Conclusions

M. tuberculosis in HIV-infected patients had lower expression of the DosR regulon, a critical metabolic and immunomodulatory switch induced by NO, carbon monoxide, and hypoxia. Our human data suggest that decreased DosR expression may result from alternative pathway activation of macrophages, with consequent decreased NO expression and/or by poor granuloma formation with consequent decreased hypoxic stress.

Objective

Populations of Mexican American ancestry are at an increased risk for nonalcoholic fatty liver disease. The objective of this study was to determine whether loci in known and novel genes were associated with variation in aspartate aminotransferase (AST) (n = 3,644), alanine aminotransferase (ALT) (n = 3,595), and gamma-glutamyl transferase (GGT) (n = 1,577) levels by conducting the first genome-wide association study (GWAS) of liver enzymes, which commonly measure liver function, in individuals of Mexican American ancestry.

Methods

Levels of AST, ALT, and GGT were determined by enzymatic colorimetric assays. A multi-cohort GWAS of individuals of Mexican American ancestry was performed. Single-nucleotide polymorphisms (SNP) were tested for association with liver outcomes by multivariable linear regression using an additive genetic model. Association analyses were conducted separately in each cohort, followed by a nonparametric meta-analysis.

Results

In the PNPLA3 gene, rs4823173 (P = 3.44 × 10^-10 ), rs2896019 (P = 7.29 × 10^-9 ), and rs2281135 (P = 8.73 × 10^-9 ) were significantly associated with AST levels. Although not genome-wide significant, these same SNPs were the top hits for ALT (P = 7.12 × 10^-8 , P = 1.98 × 10^-7 , and P = 1.81 × 10^-7 , respectively). The strong correlation (r²  = 1.0) for these SNPs indicated a single hit in the PNPLA3 gene. No genome-wide significant associations were found for GGT.

Conclusions

PNPLA3, a locus previously identified with ALT, AST, and nonalcoholic fatty liver disease in European and Japanese GWAS, is also associated with liver enzymes in populations of Mexican American ancestry.

Background

Bronchial airway expression profiling has identified inflammatory subphenotypes of asthma, but the invasiveness of this technique has limited its application to childhood asthma.

Objectives

We sought to determine whether the nasal transcriptome can proxy expression changes in the lung airway transcriptome in asthmatic patients. We also sought to determine whether the nasal transcriptome can distinguish subphenotypes of asthma.

Methods

Whole-transcriptome RNA sequencing was performed on nasal airway brushings from 10 control subjects and 10 asthmatic subjects, which were compared with established bronchial and small-airway transcriptomes. Targeted RNA sequencing nasal expression analysis was used to profile 105 genes in 50 asthmatic subjects and 50 control subjects for differential expression and clustering analyses.

Results

We found 90.2% overlap in expressed genes and strong correlation in gene expression (ρ = .87) between the nasal and bronchial transcriptomes. Previously observed asthmatic bronchial differential expression was strongly correlated with asthmatic nasal differential expression (ρ = 0.77, P = 5.6 × 10(-9)). Clustering analysis identified TH2-high and TH2-low subjects differentiated by expression of 70 genes, including IL13, IL5, periostin (POSTN), calcium-activated chloride channel regulator 1 (CLCA1), and serpin peptidase inhibitor, clade B (SERPINB2). TH2-high subjects were more likely to have atopy (odds ratio, 10.3; P = 3.5 × 10(-6)), atopic asthma (odds ratio, 32.6; P = 6.9 × 10(-7)), high blood eosinophil counts (odds ratio, 9.1; P = 2.6 × 10(-6)), and rhinitis (odds ratio, 8.3; P = 4.1 × 10(-6)) compared with TH2-low subjects. Nasal IL13 expression levels were 3.9-fold higher in asthmatic participants who experienced an asthma exacerbation in the past year (P = .01). Several differentially expressed nasal genes were specific to asthma and independent of atopic status.

Conclusion

Nasal airway gene expression profiles largely recapitulate expression profiles in the lung airways. Nasal expression profiling can be used to identify subjects with IL13-driven asthma and a TH2-skewed systemic immune response.

Rationale: Common genetic variants have been associated with idiopathic pulmonary fibrosis (IPF). Objectives: To determine functional relevance of the 10 IPF-associated common genetic variants we previously identified. Methods: We performed expression quantitative trait loci (eQTL) and methylation quantitative trait loci (mQTL) mapping, followed by co-localization of eQTL and mQTL with genetic association signals and functional validation by luciferase reporter assays. Illumina multi-ethnic genotyping arrays, mRNA sequencing, and Illumina 850k methylation arrays were performed on lung tissue of participants with IPF (234 RNA and 345 DNA samples) and non-diseased controls (188 RNA and 202 DNA samples). Measurements and Main Results: Focusing on genetic variants within 10 IPF-associated genetic loci, we identified 27 eQTLs in controls and 24 eQTLs in cases (false-discovery-rate-adjusted P < 0.05). Among these signals, we identified associations of lead variants rs35705950 with expression of MUC5B and rs2076295 with expression of DSP in both cases and controls. mQTL analysis identified CpGs in gene bodies of MUC5B (cg17589883) and DSP (cg08964675) associated with the lead variants in these two loci. We also demonstrated strong co-localization of eQTL/mQTL and genetic signal in MUC5B (rs35705950) and DSP (rs2076295). Functional validation of the mQTL in MUC5B using luciferase reporter assays demonstrates that the CpG resides within a putative internal repressor element. Conclusions: We have established a relationship of the common IPF genetic risk variants rs35705950 and rs2076295 with respective changes in MUC5B and DSP expression and methylation. These results provide additional evidence that both MUC5B and DSP are involved in the etiology of IPF.