Abstract:Objective: Disuse osteoporosis is known to be primarily caused by a lack of exercise. However, the causal relationships between zinc and immunity and disuse osteoporosis remain unknown. This study investigated these relationships and their potential mechanisms. Methods: This study was an integrative study combining genome-wide association studies and transcriptomics. Two-sample Mendelian randomization analysis (MR) was used to analyze the causal relationships between exposures (zinc, immunity, physical activity) and the outcome (osteoporosis) with the aid of single-nucleotide polymorphisms (SNPs) as instrumental variables (IVs). Four models, MR-Egger, inverse variance weighted, weighted median and MR-Pleiotrophy RESidual Sum and Outlier (MRPRESSO), were used to calculate odds ratio values. Sensitivity and heterogeneity analyses were also performed using MRPRESSO and MR-Egger methods. The mRNA transcriptomic analysis was subsequently conducted. Zinc metabolism scores were acquired through single-sample Gene Set Enrichment Analysis algorithms. Stromal scores were obtained using the R Package "estimate" algorithms. Important Kyoto Encyclopedia of Genes and Genomes and Gene Ontology pathways were also derived through gene set variation analysis. Cytoscape software helped construct the transcription factor (TF)-mRNA-microRNA (miRNA) network. Virtual screening and molecular docking were performed. Polymerase chain reaction validation was also carried out in vivo. Results: Causal relationships were demonstrated between zinc and exercise (95% confidence interval [CI] = 1.30-2.95, p = 0.001), exercise and immunity (95% CI = 0.36-0.80, p = 0.002), exercise and osteoporosis (95% CI = 0.97-0.99, p = 0.0007), and immunity disorder and osteoporosis (95% CI = 1.30-2.03, p = 0.00002). One hundred and seventy-nine mRNAs in important modules were screened. Combining the differential expressional genes (DEGs) and the Boruta selection, six DEGs were screened (AHNAK, CSF2, ADAMTS12, SRA1, RUNX2, and SLC39A14). TF HOXC10 and miRNA hsa-miR-204 were predicted. Then, the TF-mRNA-miRNA network was successfully constructed. RUNX2 and SLC39A14 were identified as hub mRNAs in the TF-mRNA-miRNA network. Eventually, the novel small drug C6O4NH5 was designed according to the pharmacophore structure of SLC39A14. The docking energy for the novel drug was -5.83 kcal/mol. SLC39A14 and RUNX2 were downregulated (of statistical significance p-value < 0.05) in our animal experiment. Conclusion: This study revealed that zinc had a protective causal relationship with disuse osteoporosis by promoting exercise and immunity. SLC39A14 and RUNX2 mRNA participated in this zinc-related mechanism.

Transcription Factor Enrichment Analysis in Enhancers Identifies EZH2 As a Susceptibility Gene for Osteoporosis.

Integrative high-throughput enhancer surveying and functional verification divulges a YY2-condensed regulatory axis conferring risk for osteoporosis

Integrative Genomic Analysis Predicts Novel Functional Enhancer-Snps for Bone Mineral Density

Integrating Epigenomic Elements and GWASs Identifies BDNF Gene Affecting Bone Mineral Density and Osteoporotic Fracture Risk

Transcriptomic Data Identified Key Transcription Factors for Osteoporosis in Caucasian Women

The Osteoporosis Susceptible SNP Rs4325274 Remotely Regulates the SOX6 Gene Through Enhancers.

Gene Expression and RNA Splicing Imputation Identifies Novel Candidate Genes Associated with Osteoporosis.

Author Response for "an Osteoporosis Susceptibility Allele at 11p15 Regulates SOX6 Expression by Modulating TCF4 Chromatin Binding"

Multiple Functional Variants at 13q14 Risk Locus for Osteoporosis Regulate RANKL Expression Through Long-Range Super-Enhancer

Analysis of Recently Identified Osteoporosis Susceptibility Genes in Han Chinese Women.

Integrating genome-wide association and transcriptome prediction model identifies novel target genes for osteoporosis

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

Genetics association study and functional analysis on osteoporosis susceptibility gene BDNF

Integration of Summary Data from GWAS and Eqtl Studies Identified Novel Causal BMD Genes with Functional Predictions.

Integrating Genome- Wide Association Study with Regulatory SNP Annotations Identified Novel Candidate Genes for Osteoporosis

Genetic Association Study Identified a 20 Kb Regulatory Element in WLS Associated with Osteoporosis and Bone Mineral Density in Han Chinese

The osteoporosis risk variant rs9820407 at 3p22.1 acts as an allele-specific enhancer to regulate CTNNB1 expression by long-range chromatin loop formation

LncRNA ZBTB40-IT1 Modulated by Osteoporosis GWAS Risk SNPs Suppresses Osteogenesis

Identification of Key Osteoporosis Genes Through Comparative Analysis of Men's and Women's Osteoblast Transcriptomes

Polymorphisms in the Estrogen Receptor Genes Are Associated with Hip Fractures in Chinese.

Analysis of Zinc and Stromal Immunity in Disuse Osteoporosis: Mendelian Randomization and Transcriptomic Analysis