Prenatal diagnosis and treatment for fetal angiotensin converting enzyme deficiency.

OBJECTIVE: To elucidate an etiology in a case with persistent oligohydramnios by prenatal diagnosis and actively treat the case to achieve good prognosis. METHODS: We performed whole exome sequencing (WES) of DNA from the fetus and parents. Serial amnioinfusions were conducted until birth. Pressors were required to maintain normal blood pressure. The infant angiotensin-converting enzyme (ACE) activity, angiotensin II (Ang II, a downstream product of ACE), and compensatory enzymes (CEs) activities were measured. Compensatory enzyme activities in plasma from age-matched healthy controls were also detected. RESULTS: We identified a fetus with a severe ACE mutation prenatally. The infant was born prematurely without pulmonary dysplasia. Hypotension and anuria resolved spontaneously. He had almost no ACE activity, but his Ang II level and CE activity exceeded the upper limit of the normal range and the upper limit of the 95% confidence interval of controls, respectively. His renal function also largely recovered. CONCLUSION: Fetuses with ACE mutations can be diagnosed prenatally through WES. Serial amnioinfusion permits the continuation of pregnancy in fetal ACE deficiency. Compensatory enzymes for defective ACE appeared postnatally. Renal function may be spared by preterm delivery; furthermore, for postnatal vasopressor therapy to begin, improving renal perfusion pressure before nephrogenesis has been completed.

Identification of the serum metabolites associated with cow milk consumption in Chinese Peri-/Postmenopausal women.

Cow milk consumption (CMC) and downstream alterations of serum metabolites are commonly considered important factors regulating human health status. Foods may lead to metabolic changes directly or indirectly through remodelling gut microbiota (GM). We sought to identify the metabolic alterations in Chinese Peri-/Postmenopausal women with habitual CMC and explore if the GM mediates the CMC-metabolite associations. 346 Chinese Peri-/Postmenopausal women participants were recruited in this study. Fixed effects regression and partial least squares discriminant analysis (PLS-DA) were applied to reveal alterations of serum metabolic features in different CMC groups. Spearman correlation coefficient was computed to detect metabolome-metagenome association. 36 CMC-associated metabolites including palmitic acid (FA(16:0)), 7alpha-hydroxy-4-cholesterin-3-one (7alphaC4), citrulline were identified by both fixed effects regression (FDR < 0.05) and PLS-DA (VIP score > 2). Some significant metabolite-GM associations were observed, including FA(16:0) with gut species Bacteroides ovatus, Bacteroides sp.D2. These findings would further prompt our understanding of the effect of cow milk on human health.

Gut microbiota accelerates obesity in peri-/post-menopausal women via Bacteroides fragilis and acetic acid.

OBJECTIVE: Many animal experiments and epidemiological studies have shown that the gut microbiota (GM) plays an important role in the development of obesity, but the specific biological mechanism involved in the pathogenesis of disease remain unknown. We aimed to examine the relationships and functional mechanisms of GM on obesity in peri- and post-menopausal women. METHODS: We recruited 499 Chinese peri- and post-menopausal women and performed comprehensive analyses of the gut microbiome, targeted metabolomics for short-chain fatty acids in serum, and host whole-genome sequencing by various association analysis methods. RESULTS: Through constrained linear regression analysis, we found that an elevated abundance of Bacteroides fragilis (B. fragilis) was associated with obesity. We also found that serum levels of acetic acid were negatively associated with obesity, and that B. fragilis was negatively associated with serum acetic acid levels by partial Spearman correlation analysis. Mendelian randomization analysis indicated that B. fragilis increases the risk of obesity and may causally down-regulate acetic acid levels. CONCLUSIONS: We found the gut with B. fragilis may accelerate obesity, in part, by suppressing acetic acid levels. Therefore, B. fragilis and acetic acid may represent important therapeutic targets for obesity intervention in peri- and post-menopausal women.

Combining artificial intelligence: deep learning with Hi-C data to predict the functional effects of non-coding variants.

MOTIVATION: Although genome-wide association studies (GWASs) have identified thousands of variants for various traits, the causal variants and the mechanisms underlying the significant loci are largely unknown. In this study, we aim to predict non-coding variants that may functionally affect translation initiation through long-range chromatin interaction. RESULTS: By incorporating the Hi-C data, we propose a novel and powerful deep learning model of artificial intelligence to classify interacting and non-interacting fragment pairs and predict the functional effects of sequence alteration of single nucleotide on chromatin interaction and thus on gene expression. The changes in chromatin interaction probability between the reference sequence and the altered sequence reflect the degree of functional impact for the variant. The model was effective and efficient with the classification of interacting and non-interacting fragment pairs. The predicted causal SNPs that had a larger impact on chromatin interaction were more likely to be identified by GWAS and eQTL analyses. We demonstrate that an integrative approach combining artificial intelligence-deep learning with high throughput experimental evidence of chromatin interaction leads to prioritizing the functional variants in disease- and phenotype-related loci and thus will greatly expedite uncover of the biological mechanism underlying the association identified in genomic studies. AVAILABILITY AND IMPLEMENTATION: Source code used in data preparing and model training is available at the GitHub website (https://github.com/biocai/DeepHiC). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Epigenomic and Transcriptomic Prioritization of Candidate Obesity-Risk Regulatory GWAS SNPs.

Concern about rising rates of obesity has prompted searches for obesity-related single nucleotide polymorphisms (SNPs) in genome-wide association studies (GWAS). Identifying plausible regulatory SNPs is very difficult partially because of linkage disequilibrium. We used an unusual epigenomic and transcriptomic analysis of obesity GWAS-derived SNPs in adipose versus heterologous tissues. From 50 GWAS and 121,064 expanded SNPs, we prioritized 47 potential causal regulatory SNPs (Tier-1 SNPs) for 14 gene loci. A detailed examination of seven loci revealed that four (CABLES1, PC, PEMT, and FAM13A) had Tier-1 SNPs positioned so that they could regulate use of alternative transcription start sites, resulting in different polypeptides being generated or different amounts of an intronic microRNA gene being expressed. HOXA11 and long noncoding RNA gene RP11-392O17.1 had Tier-1 SNPs in their 3' or promoter region, respectively, and strong preferences for expression in subcutaneous versus visceral adipose tissue. ZBED3-AS1 had two intragenic Tier-1 SNPs, each of which could contribute to mediating obesity risk through modulating long-distance chromatin interactions. Our approach not only revealed especially credible novel regulatory SNPs, but also helped evaluate previously highlighted obesity GWAS SNPs that were candidates for transcription regulation.

Cell-specific network analysis of human folliculogenesis reveals network rewiring in antral stage oocytes.

Although previous studies have explored the gene expression profiles of human oocytes and granulosa cells by single-cell RNA sequencing (scRNA-seq), the dynamic regulatory network at a single-cell resolution during folliculogenesis remains largely unknown. We identified 10 functional modules by WGCNA, four of which were significantly correlated with primary/antral oocyte and antral/pre-ovulatory granulosa cells. Functional enrichment analysis showed that the brown module, which was correlated with antral oocyte, was enriched in oocyte differentiation, and two core subnetworks identified by MCODE were involved in cell cycle (blue subnetwork) and oogenesis (red subnetwork). The cell-specific network (CSN) analysis demonstrated a distinct gene network structure associated with the antral follicular stage, which was notably different from other developmental stages. To our knowledge, this is the first study to explore gene functions during folliculogenesis at single-cell network level. We uncovered two potential gene subnetworks, which may play an important role in oocyte function beginning at the antral stage, and further established their rewiring process at intra-network/whole transcriptome level. The findings provide crucial insights from a novel network perspective to be further explored in functional mechanistic studies.

Mutant Zp1 impedes incorporation of ZP3 and ZP4 in the zona pellucida, resulting in zona absence and female infertility in rats†.

The zona pellucida (ZP) plays vital roles in reproductive processes including oogenesis, fertilization, and preimplantation development. Both human and rat ZP consist of four glycoproteins, called ZP1, ZP2, ZP3, and ZP4. Our previous research reported a novel Zp1 mutation in cases of human infertility, associated with an abnormal phenotype involving the absence of the ZP. Here, we developed a homologous rat strain to investigate the pathogenic effect. The ovaries of homozygous (Zp1MT/MT) females possessed both growing and fully grown oocytes; the oocytes completely lacked a ZP, but ZP1 was detectable inside the cytoplasm. Only 1-2 eggs were recovered from oviducts of superovulated Zp1MT/MT females, while an average of 21 eggs were recovered from superovulated Zp1WT/WT per female. The eggs of Zp1MT/MT females were not surrounded by a ZP and lost their fertilization capacity in vitro. Zp1MT/MT females mated with wild-type males failed to become pregnant. Studies in 293T cells showed that mutant Zp1 resulted in a truncated ZP1 protein, which might be intracellularly sequestered and interacted with wild-type ZP3 or ZP4. Our results suggest that the Zp1 point mutation led to infertility and loss of the ZP in oocytes in rats.

Zp4 is completely dispensable for fertility in female rats†.

Zona pellucida (ZP), which is composed of at most four extracellular glycoproteins (ZP1, ZP2, ZP3, and ZP4) in mammals, shelters the oocytes and is vital in female fertility. Several studies have identified the indispensable roles of ZP1-3 in maintaining normal female fertility. However, the understanding of ZP4 is still very poor because only one study on ZP4-associated infertility performed in rabbits has been reported up to date. Here we investigated the function of mammalian Zp4 by creating a knockout (KO) rat strain (Zp4-/- rat) using CRISPR-Cas9-mediated DNA-editing method. The influence of Zp4 KO on ZP morphology and some pivotal processes of reproduction, including oogenesis, ovulation, fertilization, and pup production, were studied using periodic acid-Schiff's staining, superovulation, in vitro fertilization, and natural mating. The ZP morphology in Zp4-/- rats was normal, and none of these pivotal processes was affected. This study renewed the knowledge of mammalian Zp4 by suggesting that Zp4 was completely dispensable for female fertility.

System network analysis of genomics and transcriptomics data identified type 1 diabetes-associated pathway and genes.

Genome-wide association studies (GWASs) have discovered >50 risk loci for type 1 diabetes (T1D). However, those variations only have modest effects on the genetic risk of T1D. In recent years, accumulated studies have suggested that gene-gene interactions might explain part of the missing heritability. The purpose of our research was to identify potential and novel risk genes for T1D by systematically considering the gene-gene interactions through network analyses. We carried out a novel system network analysis of summary GWAS statistics jointly with transcriptomic gene expression data to identify some of the missing heritability for T1D using weighted gene co-expression network analysis (WGCNA). Using WGCNA, seven modules for 1852 nominally significant (P ≤ 0.05) GWAS genes were identified by analyzing microarray data for gene expression profile. One module (tagged as green module) showed significant association (P ≤ 0.05) between the module eigengenes and the trait. This module also displayed a high correlation (r = 0.45, P ≤ 0.05) between module membership (MM) and gene significant (GS), which indicated that the green module of co-expressed genes is of significant biological importance for T1D status. By further describing the module content and topology, the green module revealed a significant enrichment in the "regulation of immune response" (GO:0050776), which is a crucially important pathway in T1D development. Our findings demonstrated a module and several core genes that act as essential components in the etiology of T1D possibly via the regulation of immune response, which may enhance our fundamental knowledge of the underlying molecular mechanisms for T1D.

Influence of mouse defective zona pellucida in folliculogenesis on apoptosis of granulosa cells and developmental competence of oocytes†.

Zona pellucida (ZP), which enwraps the oocyte during folliculogenesis, initially forms in the primary follicle and plays an important role in female fertility. Here, we investigated a mouse strain ("mutant mice" for short) carrying two types of ZP defects in folliculogenesis, i.e., ZP thinned (but intact) and ZP cracked, caused by targeted mutation in the Zp1 gene. Using this mutant mouse strain and wild-type mouse as control, we studied the effects of the ZP defects on the development of oocytes and granulosa cells during folliculogenesis. For each ZP defect, we examined the morphology of transzonal projections and apoptosis of granulosa cells in the corresponding growing follicles, as well as the morphology of corresponding ovulated eggs and their abilities to develop into viable individuals. Our results suggested that ZP integrity rather than thickness or porosity is crucial for preventing the ectopia of granulosa cells, maintaining adequate routine bilateral signaling between oocyte and surrounding granulosa cells, and thus for ensuring the survival of granulosa cells and the establishment of the full developmental competence of oocytes. This is the first study to elucidate the effects of different degrees of ZP defects caused by the same gene mutation, on the apoptosis of granulosa cells and developmental competence of oocytes, and to explore the potential mechanisms underlying these effects.

Identification of a 1p21 independent functional variant for abdominal obesity.

OBJECTIVES: Aiming to uncover the genetic basis of abdominal obesity, we performed a genome-wide association study (GWAS) meta-analysis of trunk fat mass adjusted by trunk lean mass (TFMadj) and followed by a series of functional investigations. SUBJECTS: A total of 11,569 subjects from six samples were included into the GWAS meta-analysis. METHODS: Meta-analysis was performed by a weighted fixed-effects model. In silico replication analysis was performed in the UK-Biobank (UKB) sample (N = 331,093) and in the GIANT study (N up to 110,204). Cis-expression QTL (cis-eQTL) analysis, dual-luciferase reporter assay and electrophoresis mobility shift assay (EMSA) were conducted to examine the functional relevance of the identified SNPs. At last, differential gene expression analysis (DGEA) was performed. RESULTS: We identified an independent SNP rs12409479 at 1p21 (MAF = 0.07, p = 7.26 × 10-10), whose association was replicated by the analysis of TFM in the UKB sample (one-sided p = 3.39 × 10-3), and was cross-validated by the analyses of BMI (one-sided p = 0.03) and WHRadj (one-sided p = 0.04) in the GIANT study. Cis-eQTL analysis demonstrated that allele A at rs12409479 was positively associated with PTBP2 expression level in subcutaneous adipose tissue (N = 385, p = 4.15 × 10-3). Dual-luciferase reporter assay showed that the region repressed PTBP2 gene expression by downregulating PTBP2 promoter activity (p < 0.001), and allele A at rs12409479 induced higher luciferase activity than allele G did (p = 4.15 × 10-3). EMSA experiment implied that allele A was more capable of binding to unknown transcription factors than allele G. Lastly, DGEA showed that the level of PTBP2 expression was higher in individuals with obesity than in individuals without obesity (N = 20 and 11, p = 0.04 and 9.22 × 10-3), suggesting a regulatory role in obesity development. CONCLUSIONS: Taken together, we hypothesize a regulating path from rs12409479 to trunk fat mass development through its allelic specific regulation of PTBP2 gene expression, thus providing some novel insight into the genetic basis of abdominal obesity.

Inferring causal relationships between phenotypes using summary statistics from genome-wide association studies.

Genome-wide association studies (GWAS) have successfully identified numerous genetic variants associated with diverse complex phenotypes and diseases, and provided tremendous opportunities for further analyses using summary association statistics. Recently, Pickrell et al. developed a robust method for causal inference using independent putative causal SNPs. However, this method may fail to infer the causal relationship between two phenotypes when only a limited number of independent putative causal SNPs identified. Here, we extended Pickrell's method to make it more applicable for the general situations. We extended the causal inference method by replacing the putative causal SNPs with the lead SNPs (the set of the most significant SNPs in each independent locus) and tested the performance of our extended method using both simulation and empirical data. Simulations suggested that when the same number of genetic variants is used, our extended method had similar distribution of test statistic under the null model as well as comparable power under the causal model compared with the original method by Pickrell et al. But in practice, our extended method would generally be more powerful because the number of independent lead SNPs was often larger than the number of independent putative causal SNPs. And including more SNPs, on the other hand, would not cause more false positives. By applying our extended method to summary statistics from GWAS for blood metabolites and femoral neck bone mineral density (FN-BMD), we successfully identified ten blood metabolites that may causally influence FN-BMD. We extended a causal inference method for inferring putative causal relationship between two phenotypes using summary statistics from GWAS, and identified a number of potential causal metabolites for FN-BMD, which may provide novel insights into the pathophysiological mechanisms underlying osteoporosis.

Characterization of TaDREB1 in wheat genotypes with different seed germination under osmotic stress.

BACKGROUND: The cis-acting element DRE/CRT plays an important role in activating gene expression responsive to osmotic stress, low temperature and high-salinity. DREB1/CBF genes encode DRE-binding proteins with the function as transcript activators. TaDREB1 was also found to be induced by osmotic stress. METHODS: The dates of osmotic stress was assessed by seed germination drought resistance index; the full-length cDNA sequences of TaDREB1 gene were downloaded from NCBI datebase; identification of allelic variation and transcript expression were assessed by PCR and semi-quantitive RT-PCR analysis, respectively. RESULTS: Total 13 new allele variations of TaDREB1 were identified in the germplasms tested in the paper, including 5 TaDREB1-A on chromosome 3AL, 4 TaDREB1-B on chromosome 3BL and 4 TaDREB1-D on chromosome 3DL. In each variety, there existed two loci of TaDREB1-D genes, named TaDREB1-D1 and TaDREB1-D2, both of which had the similar nucleotide sequence except an 11 bp insertion in the former. In wheat seeds under osmotic stress, we did not detect the transcript expression level of TaDREB1-A and TaDREB1-B, but that of TaDREB1-D. CONCLUSIONS: The capacity of osmotic stress tolerance was closely correlated with the expression level and tendency of TaDREB1-D.

Mutant ZP1 in familial infertility.

The human zona pellucida is composed of four glycoproteins (ZP1, ZP2, ZP3, and ZP4) and has an important role in reproduction. Here we describe a form of infertility with an autosomal recessive mode of inheritance, characterized by abnormal eggs that lack a zona pellucida. We identified a homozygous frameshift mutation in ZP1 in six family members. In vitro studies showed that defective ZP1 proteins and normal ZP3 proteins colocalized throughout the cells and were not expressed at the cell surface, suggesting that the aberrant ZP1 results in the sequestration of ZP3 in the cytoplasm, thereby preventing the formation of the zona pellucida around the oocyte.

CXADR promote epithelial-mesenchymal transition in endometriosis by modulating AKT/GSK-3ß signaling.

Endometriosis is a benign high prevalent disease exhibiting malignant features. However, the underlying pathogenesis and key molecules of endometriosis remain unclear. By integrating and analysis of existing expression profile datasets, we identified coxsackie and adenovirus receptor (CXADR), as a novel key gene in endometriosis. Based on the results of immunohistochemistry (IHC), we confirmed significant down-regulation of CXADR in ectopic endometrial tissues obtained from women with endometriosis compared with healthy controls. Further in vitro investigation indicated that CXADR regulated the stability and function of the phosphatases and AKT inhibitors PHLPP2 (pleckstrin homology domain and leucine-rich repeat protein phosphatase 2) and PTEN (phosphatase and tensin homolog). Loss of CXADR led to phosphorylation of AKT and glycogen synthase kinase-3ß (GSK-3ß), which resulted in stabilization of an epithelial-mesenchymal transition (EMT) factor, SNAIL1 (snail family transcriptional repressor 1). Therefore, EMT processs was induced, and the proliferation, migration and invasion of Ishikawa cells were enhanced. Over-expression of CXADR showed opposite effects. These findings suggest a previously undefined role of AKT/GSK-3ß signaling axis in regulating EMT and reveal the involvement of a CXADR-induced EMT, in pathogenic progression of endometriosis.

Single-cell RNA-seq identified novel genes involved in primordial follicle formation.

Introduction: The number of primordial follicles (PFs) in mammals determines the ovarian reserve, and impairment of primordial follicle formation (PFF) will cause premature ovarian insufficiency (POI). Methods: By analyzing public single-cell RNA sequencing performed during PFF on mice and human ovaries, we identified novel functional genes and novel ligand-receptor interaction during PFF. Based on immunofluorescence and in vitro ovarian culture, we confirmed mechanisms of genes and ligand-receptor interaction in PFF. We also applied whole exome sequencing (WES) in 93 cases with POI and whole genome sequencing (WGS) in 465 controls. Variants in POI patients were further investigated by in silico analysis and functional verification. Results: We revealed ANXA7 (annexin A7) and GTF2F1 (general transcription factor IIF subunit 1) in germ cells to be novel potentially genes in promoting PFF. Ligand Mdk (midkine) in germ cells and its receptor Sdc1 (syndecan 1) in granulosa cells are novel interaction crucial for PFF. Based on immunofluorescence, we confirmed significant up-regulation of ANXA7 in PFs compared with germline cysts, and uniform expression of GTF2F1, MDK and SDC1 during PFF, in 25 weeks human fetal ovary. In vitro investigation indicated that Anxa7 and Gtf2f1 are vital for mice PFF by regulating Jak/Stat3 and Jnk signaling pathways, respectively. Ligand-receptor (Mdk-Sdc1) are crucial for PFF by regulating Pi3k-akt signaling pathway. Two heterozygous variants in GTF2F1, and one heterozygous variants in SDC1 were identified in cases, but no variant were identified in controls. The protein level of GTF2F1 or SDC1 in POI cases are significantly lower than that of controls, indicating the pathogenic effects of the two genes on ovarian function were dosage dependent. Discussion: Our study identified novel genes and novel ligand-receptor interaction during PFF, and further expanding the genetic architecture of POI.

Novel SETBP1 mutation in a chinese family with intellectual disability.

BACKGROUND: Intellectual disability (ID) is characterized by an IQ < 70, which implies below-average intellectual function and a lack of skills necessary for daily living. ID may occur due to multiple causes, such as metabolic, infectious, and chromosomal causes. ID affects approximately 1-3% of the population; however, the cause can be identified in only 25% of clinical patients. METHODS: To find the cause of genetic ID in a family, we performed whole-exome sequencing and Sanger sequencing to confirm the presence of a SETBP1 variant and real-time quantitative polymerase chain reaction to detect SETBP1 expression in the proband and normal controls. RESULTS: A novel variant, c.942_943insGT (p. Asp316TrpfsTer28), was found in SETBP1. Furthermore, we observed that SETBP1 expression in patients was only 20% that of normal controls (P < 0.05). CONCLUSION: A heterozygous variant in SETBP1 associated with ID was found. This report provides further evidence for its genetic basis and support for clinical genetic diagnosis.

Systematic metabolomic studies identified adult adiposity biomarkers with acetylglycine associated with fat loss in vivo.

Obesity is associated with various adverse health outcomes. Body fat (BF) distribution is recognized as an important factor of negative health consequences of obesity. Although metabolomics studies, mainly focused on body mass index (BMI) and waist circumference, have explored the biological mechanisms involved in the development of obesity, these proxy composite measures are not accurate and cannot reflect BF distribution, and thus may hinder accurate assessment of metabolic alterations and differential risk of metabolic disorders among individuals presenting adiposity differently throughout the body. Thus, the exact relations between metabolites and BF remain to be elucidated. Here, we aim to examine the associations of metabolites and metabolic pathways with BF traits which reflect BF distribution. We performed systematic untargeted serum metabolite profiling and dual-energy X-ray absorptiometry (DXA) whole body fat scan for 517 Chinese women. We jointly analyzed DXA-derived four BF phenotypes to detect cross-phenotype metabolite associations and to prioritize important metabolomic factors. Topology-based pathway analysis was used to identify important BF-related biological processes. Finally, we explored the relationships of the identified BF-related candidate metabolites with BF traits in different sex and ethnicity through two independent cohorts. Acetylglycine, the top distinguished finding, was validated for its obesity resistance effect through in vivo studies of various diet-induced obese (DIO) mice. Eighteen metabolites and fourteen pathways were discovered to be associated with BF phenotypes. Six of the metabolites were validated in varying sex and ethnicity. The obesity-resistant effects of acetylglycine were observed to be highly robust and generalizable in both human and DIO mice. These findings demonstrate the importance of metabolites associated with BF distribution patterns and several biological pathways that may contribute to obesity and obesity-related disease etiology, prevention, and intervention. Acetylglycine is highlighted as a potential therapeutic candidate for preventing excessive adiposity in future studies.

Enteric nervous system damage caused by abnormal intestinal butyrate metabolism may lead to functional constipation.

Functional constipation (FC) is a high morbidity gastrointestinal disease for which dysfunction in the enteric nervous system is a major pathogenesis mechanism. To enhance our understanding of the involvement of intestinal microbiota and its metabolites in the pathogenesis of FC, we conducted a shotgun metagenomic sequencing analysis of gut microbiota and serum short-chain fatty acids (SCFAs) analysis in 460 Chinese women with different defecation frequencies. We observed that the abundance ofFusobacterium_varium, a butyric acid-producing bacterium, was positively correlated (P = 0.0096) with the frequency of defecation; however, the concentrations of serum butyric acid was negatively correlated (P = 3.51E-05) with defecation frequency. These results were verified in an independent cohort (6 patients with FC and 6 controls). To further study the effects of butyric acid on intestinal nerve cells, we treated mouse intestinal neurons in vitro with various concentrations of butyrate (0.1, 0.5, 1, and 2.5 mM). We found that intestinal neurons treated with 0.5 mM butyrate proliferated better than those in the other treatment groups, with significant differences in cell cycle and oxidative phosphorylation signal pathways. We suggest that the decreased butyrate production resulting from the reduced abundance of Fusobacterium in gut microbiota affects the proliferation of intestinal neurons and the energy supply of intestinal cells. However, with FC disease advancing, the consumption and excretion of butyric acid reduce, leading to its accumulation in the intestine. Moreover, the accumulation of an excessively high amount of butyric acid inhibits the proliferation of nerve cells and subsequently exacerbates the disease.