Gut dysbiosis has been linked to type 1 diabetes (T1D); however, microbial capacity in T1D remains unclear. Here, we integratively profiled gut microbial functional and metabolic alterations in children with new-onset T1D in independent cohorts and investigated the underlying mechanisms. In T1D, the microbiota was characterized by decreased butyrate production and bile acid metabolism and increased lipopolysaccharide biosynthesis at the species, gene, and metabolite levels. The combination of 18 bacterial species and fecal metabolites provided excellently discriminatory power for T1D. Gut microbiota from children with T1D induced elevated fasting glucose levels and declined insulin sensitivity in antibiotic-treated mice. In streptozotocin-induced T1D mice, butyrate and lipopolysaccharide exerted protective and destructive effects on islet structure and function, respectively. Lipopolysaccharide aggravated the pancreatic inflammatory response, while butyrate activated Insulin1 and Insulin2 gene expression. Our study revealed perturbed microbial functional and metabolic traits in T1D, providing potential avenues for microbiome-based prevention and intervention for T1D.
OBJECTIVE Type 1 diabetes (T1D) is a highly heritable disease with much lower incidence but more adult-onset cases in the Chinese population. Although genome-wide association studies (GWAS) have identified >60 T1D loci in Caucasians, less is known in Asians. RESEARCH DESIGN AND METHODS We performed the first two-stage GWAS of T1D using 2,596 autoantibody-positive T1D case subjects and 5,082 control subjects in a Chinese Han population and evaluated the associations between the identified T1D risk loci and age and fasting C-peptide levels at T1D diagnosis. RESULTS We observed a high genetic correlation between children/adolescents and adult T1D case subjects (rg = 0.87), as well as subgroups of autoantibody status (rg ≥ 0.90). We identified four T1D risk loci reaching genome-wide significance in the Chinese Han population, including two novel loci, rs4320356 near BTN3A1 (odds ratio [OR] 1.26, P = 2.70 × 10−8) and rs3802604 in GATA3 (OR 1.24, P = 2.06 × 10−8), and two previously reported loci, rs1770 in MHC (OR 4.28, P = 2.25 × 10−232) and rs705699 in SUOX (OR 1.46, P = 7.48 × 10−20). Further fine mapping in the MHC region revealed five independent variants, including another novel locus, HLA-C position 275 (omnibus P = 9.78 × 10−12), specific to the Chinese population. Based on the identified eight variants, we achieved an area under the curve value of 0.86 (95% CI 0.85–0.88). By building a genetic risk score (GRS) with these variants, we observed that the higher GRS were associated with an earlier age of T1D diagnosis (P = 9.08 × 10−11) and lower fasting C-peptide levels (P = 7.19 × 10−3) in individuals newly diagnosed with T1D. CONCLUSIONS Our results extend current knowledge on genetic contributions to T1D risk. Further investigations in different populations are needed for genetic heterogeneity and subsequent precision medicine.