Coexpression networks implicate human midfetal deep cortical projection neurons in the pathogenesis of autism.

Autism spectrum disorder (ASD) is a complex developmental syndrome of unknown etiology. Recent studies employing exome- and genome-wide sequencing have identified nine high-confidence ASD (hcASD) genes. Working from the hypothesis that ASD-associated mutations in these biologically pleiotropic genes will disrupt intersecting developmental processes to contribute to a common phenotype, we have attempted to identify time periods, brain regions, and cell types in which these genes converge. We have constructed coexpression networks based on the hcASD "seed" genes, leveraging a rich expression data set encompassing multiple human brain regions across human development and into adulthood. By assessing enrichment of an independent set of probable ASD (pASD) genes, derived from the same sequencing studies, we demonstrate a key point of convergence in midfetal layer 5/6 cortical projection neurons. This approach informs when, where, and in what cell types mutations in these specific genes may be productively studied to clarify ASD pathophysiology.

Haploinsufficiency underlies the neurodevelopmental consequences of SLC6A1 variants.

Heterozygous variants in SLC6A1, encoding the GAT-1 GABA transporter, are associated with seizures, developmental delay, and autism. The majority of affected individuals carry missense variants, many of which are recurrent germline de novo mutations, raising the possibility of gain-of-function or dominant-negative effects. To understand the functional consequences, we performed an in vitro GABA uptake assay for 213 unique variants, including 24 control variants. De novo variants consistently resulted in a decrease in GABA uptake, in keeping with haploinsufficiency underlying all neurodevelopmental phenotypes. Where present, ClinVar pathogenicity reports correlated well with GABA uptake data; the functional data can inform future reports for the remaining 72% of unscored variants. Surface localization was assessed for 86 variants; two-thirds of loss-of-function missense variants prevented GAT-1 from being present on the membrane while GAT-1 was on the surface but with reduced activity for the remaining third. Surprisingly, recurrent de novo missense variants showed moderate loss-of-function effects that reduced GABA uptake with no evidence for dominant-negative or gain-of-function effects. Using linear regression across multiple missense severity scores to extrapolate the functional data to all potential SLC6A1 missense variants, we observe an abundance of GAT-1 residues that are sensitive to substitution. The extent of this missense vulnerability accounts for the clinically observed missense enrichment; overlap with hypermutable CpG sites accounts for the recurrent missense variants. Strategies to increase the expression of the wild-type SLC6A1 allele are likely to be beneficial across neurodevelopmental disorders, though the developmental stage and extent of required rescue remain unknown.

Enhancer variants on chromosome 2p14 regulating SPRED2 and ACTR2 act as a signal amplifier to protect against rheumatoid arthritis.

Genome-wide association studies (GWASs) have repeatedly reported multiple non-coding single-nucleotide polymorphisms (SNPs) at 2p14 associated with rheumatoid arthritis (RA), but their functional roles in the pathological mechanisms of RA remain to be explored. In this study, we integrated a series of bioinformatics and functional experiments and identified three intronic RA SNPs (rs1876518, rs268131, and rs2576923) within active enhancers that can regulate the expression of SPRED2 directly. At the same time, SPRED2 and ACTR2 influence each other as a positive feedback signal amplifier to strengthen the protective role in RA by inhibiting the migration and invasion of rheumatoid fibroblast-like synoviocytes (FLSs). In particular, the transcription factor CEBPB preferentially binds to the rs1876518-T allele to increase the expression of SPRED2 in FLSs. Our findings decipher the molecular mechanisms behind the GWAS signals at 2p14 for RA and emphasize SPRED2 as a potential candidate gene for RA, providing a potential target and direction for precise treatment of RA.

High-throughput functional dissection of noncoding SNPs with biased allelic enhancer activity for insulin resistance-relevant phenotypes.

Most of the single-nucleotide polymorphisms (SNPs) associated with insulin resistance (IR)-relevant phenotypes by genome-wide association studies (GWASs) are located in noncoding regions, complicating their functional interpretation. Here, we utilized an adapted STARR-seq to evaluate the regulatory activities of 5,987 noncoding SNPs associated with IR-relevant phenotypes. We identified 876 SNPs with biased allelic enhancer activity effects (baaSNPs) across 133 loci in three IR-relevant cell lines (HepG2, preadipocyte, and A673), which showed pervasive cell specificity and significant enrichment for cell-specific open chromatin regions or enhancer-indicative markers (H3K4me1, H3K27ac). Further functional characterization suggested several transcription factors (TFs) with preferential allelic binding to baaSNPs. We also incorporated multi-omics data to prioritize 102 candidate regulatory target genes for baaSNPs and revealed prevalent long-range regulatory effects and cell-specific IR-relevant biological functional enrichment on them. Specifically, we experimentally verified the distal regulatory mechanism at IRS1 locus, in which rs952227-A reinforces IRS1 expression by long-range chromatin interaction and preferential binding to the transcription factor HOXC6 to augment the enhancer activity. Finally, based on our STARR-seq screening data, we predicted the enhancer activity of 227,343 noncoding SNPs associated with IR-relevant phenotypes (fasting insulin adjusted for BMI, HDL cholesterol, and triglycerides) from the largest available GWAS summary statistics. We further provided an open resource (http://www.bigc.online/fnSNP-IR) for better understanding genetic regulatory mechanisms of IR-relevant phenotypes.

Systematic evaluation of genome sequencing for the diagnostic assessment of autism spectrum disorder and fetal structural anomalies.

Short-read genome sequencing (GS) holds the promise of becoming the primary diagnostic approach for the assessment of autism spectrum disorder (ASD) and fetal structural anomalies (FSAs). However, few studies have comprehensively evaluated its performance against current standard-of-care diagnostic tests: karyotype, chromosomal microarray (CMA), and exome sequencing (ES). To assess the clinical utility of GS, we compared its diagnostic yield against these three tests in 1,612 quartet families including an individual with ASD and in 295 prenatal families. Our GS analytic framework identified a diagnostic variant in 7.8% of ASD probands, almost 2-fold more than CMA (4.3%) and 3-fold more than ES (2.7%). However, when we systematically captured copy-number variants (CNVs) from the exome data, the diagnostic yield of ES (7.4%) was brought much closer to, but did not surpass, GS. Similarly, we estimated that GS could achieve an overall diagnostic yield of 46.1% in unselected FSAs, representing a 17.2% increased yield over karyotype, 14.1% over CMA, and 4.1% over ES with CNV calling or 36.1% increase without CNV discovery. Overall, GS provided an added diagnostic yield of 0.4% and 0.8% beyond the combination of all three standard-of-care tests in ASD and FSAs, respectively. This corresponded to nine GS unique diagnostic variants, including sequence variants in exons not captured by ES, structural variants (SVs) inaccessible to existing standard-of-care tests, and SVs where the resolution of GS changed variant classification. Overall, this large-scale evaluation demonstrated that GS significantly outperforms each individual standard-of-care test while also outperforming the combination of all three tests, thus warranting consideration as the first-tier diagnostic approach for the assessment of ASD and FSAs.

Transcriptome-wide association study identifies multiple genes and pathways associated with thyroid function.

Thyroid dysfunction is a common endocrine disease measured by thyroid-stimulating hormone (TSH) level. Although >70 genetic loci associated with TSH have been reported through genome-wide association studies (GWASs), the variants can only explain a small fraction of the thyroid function heritability. To identify novel candidate genes for thyroid function, we conducted the first large-scale transcriptome-wide association study (TWAS) for thyroid function using GWAS-summary data for TSH levels in up to 119 715 individuals combined with precomputed gene expression weights of six panels from four tissue types. The candidate genes identified by TWAS were further validated by TWAS replication and gene expression profiles. We identified 74 conditionally independent genes significantly associated with thyroid function, such as PDE8B (P = 1.67 × 10-282), PDE10A (P = 7.61 × 10-119), NR3C2 (P = 1.50 × 10-92) and CAPZB (P = 3.13 × 10-79). After TWAS replication using UKBB datasets, 26 genes were replicated for significant associations with thyroid-relevant diseases/traits. Among them, 16 genes were causal for their associations to thyroid-relevant diseases/traits and further validated in differential expression analyses, including two novel genes (MFSD6 and RBM47) that did not implicate in previous GWASs. Enrichment analyses detected several pathways associated with thyroid function, such as the cAMP signaling pathway (P = 7.27 × 10-4), hemostasis (P = 3.74 × 10-4), and platelet activation, signaling and aggregation (P = 9.98 × 10-4). Our study identified multiple candidate genes and pathways associated with thyroid function, providing novel clues for revealing the genetic mechanisms of thyroid function and disease.

ExsgRNA: reduce off-target efficiency by on-target mismatched sgRNA.

Clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 gene editing technology has been widely used to facilitate efficient genome editing. Current popular sgRNA design tools only consider the sgRNA perfectly matched to the target site and provide the results without any on-target mismatch. We suppose taking on-target gRNA-DNA mismatches into consideration might provide better sgRNA with similar binding activity and reduced off-target sites. Here, we trained a seq2seq-attention model with feedback-loop architecture, to automatically generate sgRNAs with on-target mismatches. Dual-luciferase reporter experiment showed that multiple sgRNAs with three mismatches could achieve the 80% of the relative activity of the perfect matched sgRNA. Meanwhile, it could reduce the number of off-target sites using sgRNAs with on-target mismatches. Finally, we provided a freely accessible web server sgRNA design tool named ExsgRNA. Users could submit their target sequence to this server and get optimal sgRNAs with less off-targets and similar on-target activity compared with the perfect-matched sgRNA.

Recent Advances in Preparation and Structure of Polyurethane Porous Materials for Sound Absorbing Application.

Various acoustic materials are developed to resolve noise pollution problem in many industries. Especially, materials with porous structure are broadly used to absorb sound energy in civil construction and transportation area. Polyurethane (PU) porous materials possess excellent damping properties, good toughness, and well-developed pore structures, which have a broad application prospect in sound absorption field. This work aims to summarize the recent progress of fabrication and structure for PU porous materials in sound absorption application. The sound absorption mechanisms of porous materials are introduced. Different kinds of structure for typical PU porous materials in sound absorption application are covered and highlighted, which include PU foam, modified PU porous materials, aerogel, templated PU, and special PU porous materials. Finally, the development direction and existing problems of PU material in sound absorption application are briefly prospected. It can be expected that porous PU with high sound absorption coefficient can be obtained by using some facile methods. The design and accurate regulation of porous structures or construction of multilayer sound absorption structure is favorably recommended to fulfill the high demand of industrial and commercial applications in the future work.

The value of ictal scalp EEG in focal epilepsies surgery: a retrospective analysis.

OBJECTIVE: To describe the association between preoperative ictal scalp electroencephalogram (EEG) results and surgical outcomes in patients with focal epilepsies. METHODS: The data of consecutive patients with focal epilepsies who received surgical treatments at our center from January 2012 to December 2021 were retrospectively analyzed. RESULTS: Our data showed that 44.2% (322/729) of patients had ictal EEG recorded on video EEG monitoring during preoperative evaluation, of which 60.6% (195/322) had a concordant ictal EEG results. No significant difference of surgery outcomes between patients with and without ictal EEG was discovered. Among MRI-negative patients, those with concordant ictal EEG had a significantly better outcome than those without ictal EEG (75.7% vs. 43.8%, p = 0.024). Further logistic regression analysis showed that concordant ictal EEG was an independent predictor for a favorable outcome (OR = 4.430, 95%CI 1.175-16.694, p = 0.028). Among MRI-positive patients, those with extra-temporal lesions and discordant ictal EEG results had a worse outcome compared to those without an ictal EEG result (44.7% vs. 68.8%, p = 0.005). Further logistic regression analysis showed that discordant ictal EEG was an independent predictor of worse outcome (OR = 0.387, 95%CI 0.186-0.807, p = 0.011) in these patients. Furthermore, our data indicated that the number of seizures was not associated with the concordance rates of the ictal EEG, nor the surgical outcomes. CONCLUSIONS: The value of ictal scalp EEG for epilepsy surgery varies widely among patients. A concordant ictal EEG predicts a good surgical outcome in MRI-negative patients, whereas a discordant ictal EEG predicts a poor postoperative outcome in lesional extratemporal lobe epilepsy.

Assessment of bidirectional relationships between brain imaging-derived phenotypes and stroke: a Mendelian randomization study.

BACKGROUND: Stroke is a major cause of mortality and long-term disability worldwide. Whether the associations between brain imaging-derived phenotypes (IDPs) and stroke are causal is uncertain. METHODS: We performed two-sample bidirectional Mendelian randomization (MR) analyses to explore the causal associations between IDPs and stroke. Summary data of 587 brain IDPs (up to 33,224 individuals) from the UK Biobank and five stroke types (sample size range from 301,663 to 446,696, case number range from 5,386 to 40,585) from the MEGASTROKE consortium were used. RESULTS: Forward MR indicated 14 IDPs belong to projection fibers or association fibers were associated with stroke. For example, higher genetically determined mean diffusivity (MD) in the right external capsule was causally associated with an increased risk of small vessel stroke (IVW OR = 2.76, 95% CI 2.07 to 3.68, P = 5.87 × 10-12). Reverse MR indicated that genetically determined higher risk of any ischemic stroke was associated with increased isotropic or free water volume fraction (ISOVF) in body of corpus callosum (IVW ß = 0.23, 95% CI 0.14 to 0.33, P = 3.22 × 10-7). This IDP is a commissural fiber and it is not included in the IDPs identified by forward MR. CONCLUSIONS: We identified 14 IDPs with statistically significant evidence of causal effects on stroke or stroke subtypes. We also identified potential causal effects of stroke on one IDP of commissural fiber. These findings might guide further work toward identifying preventative strategies at the brain imaging levels.

LDBlockShow: a fast and convenient tool for visualizing linkage disequilibrium and haplotype blocks based on variant call format files.

The triangular correlation heatmap aiming to visualize the linkage disequilibrium (LD) pattern and haplotype block structure of SNPs is ubiquitous component of population-based genetic studies. However, current tools suffered from the problem of time and memory consuming. Here, we developed LDBlockShow, an open source software, for visualizing LD and haplotype blocks from variant call format files. It is time and memory saving. In a test dataset with 100 SNPs from 60 000 subjects, it was at least 10.60 times faster and used only 0.03-13.33% of physical memory as compared with other tools. In addition, it could generate figures that simultaneously display additional statistical context (e.g. association P-values) and genomic region annotations. It can also compress the SVG files with a large number of SNPs and support subgroup analysis. This fast and convenient tool will facilitate the visualization of LD and haplotype blocks for geneticists.

Surgical methods influence on the risk of anastomotic fistula after pancreaticoduodenectomy: a systematic review and network meta-analysis.

BACKGROUND: Pancreaticoduodenectomy is the first choice surgical intervention for the radical treatment of pancreatic tumors. However, an anastomotic fistula is a common complication after pancreaticoduodenectomy with a high mortality rate. With the development of minimally invasive surgery, open pancreaticoduodenectomy (OPD), laparoscopic pancreaticoduodenectomy (LPD), and robotic pancreaticoduodenectomy (RPD) are gaining interest. But the impact of these surgical methods on the risk of anastomosis has not been confirmed. Therefore, we aimed to integrate relevant clinical studies and explore the effects of these three surgical methods on the occurrence of anastomotic fistula after pancreaticoduodenectomy. METHODS: A systematic literature search was conducted for studies reporting the RPD, LPD, and OPD. Network meta-analysis of postoperative anastomotic fistula (Pancreatic fistula, biliary leakage, gastrointestinal fistula) was performed. RESULTS: Sixty-five studies including 10,026 patients were included in the network meta-analysis. The rank of risk probability of pancreatic fistula for RPD (0.00) was better than LPD (0.37) and OPD (0.62). Thus, the analysis suggests the rank of risk of the postoperative pancreatic fistula for RPD, LPD, and OPD. The rank of risk probability for biliary leakage was similar for RPD (0.15) and LPD (0.15), and both were better than OPD (0.68). CONCLUSIONS: This network meta-analysis provided ranking for three different types of pancreaticoduodenectomy. The RPD and LPD can effectively improve the quality of surgery and are safe as well as feasible for OPD.

Subdural osteoma in an adolescent patient with epilepsy: an unusual case report and literature review.

OBJECTIVE: Subdural osteoma (SO) is a rarely reported benign tumor, and there is no report of SO manifested with epileptic seizures. We aim to further the understanding of SO-related epilepsy. METHODS: Here, we report a meaningful case of epilepsy secondary to SO. A systematic review of the literature about SO using the electronic database PubMed and Web of science up to December 2022 was conducted. RESULTS: A 15-year-old girl presented with epileptic seizures for 8 years. Magnetic resonance imaging revealed an irregular lesion with heterogeneous signal in the right frontal convexity. Right frontal craniotomy was performed to remove the lesion. The pathological diagnosis was SO. Histological analysis revealed that the mechanosensitive ion channels Piezo 1/2 were upregulated in the brain tissue compressed by the osteoma, compared with the levels in the osteoma-free region. Seizure freedom was obtained during the 6-month follow-up after the surgery. We identified 24 cases of SO in 23 articles. With our case, a total of 25 cases with 32 SOs was included. Of 25 cases, 24 are adults, and 1 is a child. Seizure has been reported only in our case. Frontal osteoma was found in 76% of the patients. Symptoms were cured in 56% of the patients after surgery. CONCLUSION: Surgery is a safe and effective approach to the treatment of symptomatic osteoma. Mechanical compression on cerebral cortex may be a predisposing factor of the epileptogenesis caused by the SO.

Multivariate analysis of seizure outcomes after resective surgery for focal epilepsy: a single-center study on 833 patients.

The predictors of seizure outcomes after resective surgery for focal epilepsy, for an update on the features of good and poor outcomes, are investigated. A retrospective study of patients with focal epilepsy undergoing resective surgery from March 2011 to April 2019 was performed. There were 3 groups according to the seizure outcomes: seizure freedom, seizure improvement, and no improvement. Predictors of seizure outcomes were identified by multivariate logistic regression analysis. Of all 833 patients, 561 (67.3%) patients remained seizure-free at the last follow-up, 203 (24.4%) patients had seizure improvement, and 69 (8.3%) had no improvement. The mean follow-up duration was 5.2 years (range: 2.7 to 9.6). Predictors of better outcomes included epilepsy duration < 5 years, localized discharge, no. of antiepileptic drugs at surgery < 3, and temporal lobe resection. However, predictors of worse outcomes included intracranial hemorrhage in infancy, interictal abnormal discharge, intracranial electrode monitoring, and acute postoperative seizure. Our study suggests that resective surgery for focal epilepsy has satisfactory outcomes. Short epilepsy duration, localized discharge, and temporal lobe resection are positive predictors of seizure freedom. Patients with these predictors are intensively recommended for surgery.

A novel homozygous C-terminal deletion in BTG4 causes zygotic cleavage failure and female infertility.

PURPOSE: We aimed to identify pathogenic variants in a female patient with primary infertility and recurrent failure of in vitro fertilization with zygotic cleavage failure. METHODS: The genomic DNA from the affected individual was subjected to whole-exome sequencing and the variant was confirmed by Sanger sequencing. The functional effect of the identified variant was further investigated in 293 T cells. RESULTS: We identified a novel homozygous deletion in BTG4 (c.580_616del) in the affected individual. The deletion results in frameshift and replacement of the last 29 residues (aa195-223) with 66 random amino acids. The mutated amino acid residues are highly conserved among mammalian species. Co-immunoprecipitation in 293 T cells showed that the mutation abolished the interaction between BTG4 and PABPN1L. CONCLUSION: This study conforms previous studies and expands the mutational spectrum of BTG4. Our findings prove the functional importance of the C-terminal of BTG4. BTG4 is a potential diagnostic and therapeutic target for patients suffering from zygotic cleavage failure.

The clinical features of patients with seizure freedom and failure after total corpus callosotomy for childhood-onset refractory epilepsy.

BACKGROUND: Corpus callosotomy is a palliative surgery for medically refractory epilepsy. We aim to analyze the clinical features of patients with seizure freedom and failure after total corpus callosotomy for childhood-onset refractory epilepsy. METHODS: We retrospectively reviewed the clinical courses of patients with childhood-onset refractory epilepsy undergoing total corpus callosotomy between May 2009 and March 2019. Seizure outcome at the last follow-up was the primary outcome. The clinical features of patients with seizure freedom and failure after callosotomy were compared. RESULTS: Eighty patients with childhood-onset refractory epilepsy underwent total corpus callosotomy; 15 (18.8%) obtained freedom from all seizures and 19 (23.8%) had unworthwhile improvement and failure. The mean ages at seizure onset in patients with seizure freedom and failure after callosotomy were 5.7 and 5.9 years; and mean seizure durations were 9.4 and 11.5 years, respectively. Univariate analysis found epilepsy syndrome (p = 0.047), mental retardation (p = 0.007), previous medical history (p = 0.004), ≥10 seizures per day (p = 0.024), theta waves in the background electroencephalogram (p = 0.024), and acute postoperative seizure (p = 0.000) were associated with failure after callosotomy. Seizure freedom after callosotomy was more common among patients with less than 10 seizures per day. CONCLUSIONS: Total corpus callosotomy is an effective palliative procedure for childhood-onset refractory epilepsy, particularly for patients with specific clinical characteristics. Callosotomy has a high seizure-free rate in well-selected patients.

Causal Associations Between Blood Lipids and COVID-19 Risk: A Two-Sample Mendelian Randomization Study.

[Figure: see text].

Human brucellosis and fever of unknown origin.

BACKGROUND: Human brucellosis has become one of the major public health problems in China, and increases atypical manifestations, such as fever of unknown origin (FUO), and misdiagnosis rates has complicated the diagnosis of brucellosis. To date, no relevant study on the relationship between brucellosis and FUO has been conducted. METHODS: We retrospectively reviewed the medical charts of 35 patients with confirmed human brucellosis and prospectively recorded their outcomes by telephone interview. The patients were admitted to the Second Affiliated Hospital of Nanchang University between January 01, 2013 and October 31, 2019. Patient data were collected from hospital medical records. RESULTS: The percentage of males was significantly higher than that of female in FUO (78.95% vs. 21.05%, P < 0.05), and 80% of the patients had a clear history of exposure to cattle and sheep. Moreover, 19 (54%) cases were hospitalized with FUO, among which the patients with epidemiological histories were significantly more than those without (P < 0.05). The incidence of toxic hepatitis in FUO patients was higher than that in non-FUO patients (89% vs. 50%, P < 0.05). Meanwhile, the misdiagnosis rate was considerably higher in the FUO group than in the non-FUO group (100% vs. 63%; P < 0.05). CONCLUSION: Brucellosis is predominantly FUO admission in a non-endemic area of China, accompanied by irregular fever and toxic hepatitis. Careful examination of the epidemiological history and timely improvement of blood and bone marrow cultures can facilitate early diagnosis and prevent misdiagnosis.

An Osteoporosis Risk SNP at 1p36.12 Acts as an Allele-Specific Enhancer to Modulate LINC00339 Expression via Long-Range Loop Formation.

Genome-wide association studies (GWASs) have reproducibly associated variants within intergenic regions of 1p36.12 locus with osteoporosis, but the functional roles underlying these noncoding variants are unknown. Through an integrative functional genomic and epigenomic analyses, we prioritized rs6426749 as a potential causal SNP for osteoporosis at 1p36.12. Dual-luciferase assay and CRISPR/Cas9 experiments demonstrate that rs6426749 acts as a distal allele-specific enhancer regulating expression of a lncRNA (LINC00339) (∼360 kb) via long-range chromatin loop formation and that this loop is mediated by CTCF occupied near rs6426749 and LINC00339 promoter region. Specifically, rs6426749-G allele can bind transcription factor TFAP2A, which efficiently elevates the enhancer activity and increases LINC00339 expression. Downregulation of LINC00339 significantly increases the expression of CDC42 in osteoblast cells, which is a pivotal regulator involved in bone metabolism. Our study provides mechanistic insight into how a noncoding SNP affects osteoporosis by long-range interaction, a finding that could indicate promising therapeutic targets for osteoporosis.

Transcriptome-wide association study identifies multiple genes associated with childhood body mass index.

BACKGROUND: Childhood obesity is one of the most common and costly nutritional problems with high heritability. The genetic mechanism of childhood obesity remains unclear. Here, we conducted a transcriptome-wide association study (TWAS) to identify novel genes for childhood obesity. METHODS: By integrating the GWAS summary of childhood body mass index (BMI), we conducted TWAS analyses with pre-computed gene expression weights in 39 obesity priority tissues. The GWAS summary statistics of childhood BMI were derived from the early growth genetics consortium with 35,668 children from 20 studies. RESULTS: We identified 15 candidate genes for childhood BMI after Bonferroni corrections. The most significant gene, ADCY3, was identified in 13 tissues, including adipose, brain, and blood. Interestingly, eight genes were only identified in the specific tissue, such as FAIM2 in the brain (P = 2.04 × 10-7) and fat mass and obesity-associated gene (FTO) in the muscle (P = 1.93 × 10-8). Compared with the TWAS results of adult BMI, we found that one gene TUBA1B with predominant influence only on childhood BMI in the muscle (P = 1.12 × 10-7). We evaluated the candidate genes by querying public databases and identified 12 genes functionally related to obesity phenotypes, including nine differentially expressed genes during the differentiation of human preadipocyte cells. The remaining genes (FAM150B, KNOP1, and LMBR1L) were regarded as novel candidate genes for childhood BMI. CONCLUSIONS: Our study identified multiple candidate genes for childhood BMI, providing novel clues for understanding the genetic mechanism of childhood obesity.