Gene-gene interaction between CD40 and CD40L reduces bone mineral density and increases osteoporosis risk in women

Bone Metabolite Profile Differs between Normal and Femur Head Necrosis (FHN/BCO)-Affected Broilers: Implications for Dysregulated Metabolic Cascades in FHN Pathophysiology

Femur head necrosis (FHN), also known as bacterial chondronecrosis with osteomyelitis (BCO), has remained an animal welfare and production concern for modern broilers regardless of efforts to select against it in primary breeder flocks. Characterized by the bacterial infection of weak bone, FHN has been found in birds without clinical lameness and remains only detectable via necropsy. This presents an opportunity to utilize untargeted metabolomics to elucidate potential non-invasive biomarkers and key causative pathways involved in FHN pathology. The current study used ultra-performance liquid chromatography coupled with high-resolution mass spectrometry (UPLC-HRMS) and identified a total of 152 metabolites. Mean intensity differences at p < 0.05 were found in 44 metabolites, with 3 significantly down-regulated and 41 up-regulated in FHN-affected bone. Multivariate analysis and a partial least squares discriminant analysis (PLS-DA) scores plot showed the distinct clustering of metabolite profiles from FHN-affected vs. normal bone. Biologically related molecular networks were predicted using an ingenuity pathway analysis (IPA) knowledge base. Using a fold-change cut off of -1.5 and 1.5, top canonical pathways, networks, diseases, molecular functions, and upstream regulators were generated using the 44 differentially abundant metabolites. The results showed the metabolites NAD+, NADP+, and NADH to be downregulated, while 5-Aminoimidazole-4-carboxamide ribonucleotide (AICAR) and histamine were significantly increased in FHN. Ascorbate recycling and purine nucleotides degradation were the top canonical pathways, indicating the potential dysregulation of redox homeostasis and osteogenesis. Lipid metabolism and cellular growth and proliferation were some of the top molecular functions predicted based on the metabolite profile in FHN-affected bone. Network analysis showed significant overlap across metabolites and predicted upstream and downstream complexes, including AMP-activated protein kinase (AMPK), insulin, collagen type IV, mitochondrial complex, c-Jun N-terminal kinase (Jnk), extracellular signal-regulated kinase (ERK), and 3β-hydroxysteroid dehydrogenase (3β HSD). The qPCR analysis of relevant factors showed a significant decrease in AMPKα2 mRNA expression in FHN-affected bone, supporting the predicted downregulation found in the IPA network analysis. Taken as a whole, these results demonstrate a shift in energy production, bone homeostasis, and bone cell differentiation that is distinct in FHN-affected bone, with implications for how metabolites drive the pathology of FHN.

Analysis of Communal Molecular Mechanism Between Chronic Obstructive Pulmonary Disease and Osteoporosis

Background

Chronic obstructive pulmonary disease (COPD) patients with osteoporosis (OP) usually experience more frequent exacerbations, worse quality of life, and heavier economic burden, however, few studies have investigated common molecular mechanisms of COPD and OP.

Objective

To explore the relationship between COPD and OP through bioinformatics analysis.

Methods

The miRNA microarray data of COPD and OP were retrieved from the Gene Expression Database (GEO), and the differentially expressed microRNAs (DEmiRNAs) were screened and the intersection was obtained. The Targetscan, miRDB, and miRWalk databases were used to predict the target genes of DEmiRNA, and the gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed using the R package clusterProfiler, the STRING database was used to analyze the target protein-protein interaction network (PPI) and screens to determine the core modules and core genes.

Results

Two DEmiRNAs (miR-23a-5p, miR-194-3p) have been found in COPD and OP, which have predicted 76 and 114 target genes, respectively. GO functional annotations of miR-23a-5p were significantly enriched in CD40 signaling pathway, ubiquitin-conjugating enzyme activity, etc; KEGG pathways of miR-23a-5p were significantly enriched in ubiquitin-mediated proteolysis, folate biosynthesis, and regulation of actin cytoskeleton. GO function annotations of miR-194-3p were significantly enriched in T cell activation regulation, ubiquitin protein ligase activity, and DNA transcription factor binding; KEGG pathways of miR-194-3p were significantly enriched in cell adhesion molecules, intercellular tight junctions, and lysosomal pathway. PPI analysis found target coding proteins formed complex regulatory networks. Ten core genes (TP53, SRC, PXN, CHD4, SYK, TNRC6B, PML, KAT5, BRD1 and IGF2) were picked out among them, then we used the MCODE plugin found three core subnetworks.

Conclusion

Two identical DEmiRNAs (miR-23a-5p, miR-194-3p) exist in the peripheral blood of COPD and OP patients, which are important biomarkers for COPD patients with OP and may represent novel targets for diagnosis and treatment of COPD patients with OP.

The rs1883832 Polymorphism (CD40-1C>T) Affects the Intensity of IgA Responses after BNT162b2 Vaccination

The effectiveness of coronavirus disease 2019 (COVID-19) vaccination strategies is affected by several factors, including the genetic background of the host. In our study, we evaluated the contribution of the functional polymorphism rs1883832 affecting the Kozak sequence of the TNFSF5 gene (c.-1C>T), encoding CD40, to humoral immune responses after vaccination with the spike protein of SARS-CoV-2. The rs1883832 polymorphism was analyzed by PCR-RFLP in 476 individuals (male/female: 216/260, median age: 55.0 years, range: 20−105) of whom 342 received the BNT162b2 mRNA vaccine and 134 received the adenovirus-based vector vaccines (67 on ChAdOx1-nCoV-19 vaccine, 67 on Ad.26.COV2.S vaccine). The IgG and IgA responses were evaluated with chemiluminescent microparticle and ELISA assays on days 21, 42, and 90 after the first dose. The T allele of the rs1883832 polymorphism (allele frequency: 32.8%) was significantly associated with lower IgA levels and represented, as revealed by multivariable analysis, an independent risk factor for reduced anti-spike protein IgA levels on days 42 and 90 following BNT162b2 mRNA vaccination. Similar to serum anti-spike IgA levels, a trend of lower anti-spike IgA concentrations in saliva was found in individuals with the T allele of rs1883832. Finally, the intensity of IgA and IgG responses on day 42 significantly affected the prevalence of COVID-19 after vaccination. The rs1883832 polymorphism may be used as a molecular predictor of the intensity of anti-spike IgA responses after BNT162b2 mRNA vaccination.

Association of a single nucleotide polymorphism of RANK gene with blood pressure in Spanish women

In addition to governing key functions in bone metabolism and the immune system, the RANK/RANKL/OPG system plays a role in the vascular system, particularly in vascular calcification and atherosclerosis.Given that these 2 phenotypes are considered a major cause of high blood pressure (BP), in this study we analyzed the association of SNPs in RANK and OPG genes with blood pressure. An observational study was conducted of 2 SNPs in the RANK gene (rs884205 and rs78326403) and 1 in the OPG gene (rs4876869) with systolic (SBP) and diastolic blood pressure (DBP) in a cohort of 695 women.Data analysis revealed a statistically significant association between the SNP rs884205 and BP pressure (SBP and DBP). Analyzing this relationship by the dominant inheritance model for this SNP (allele risk: A), women of the AA/AC genotype showed higher BP than women of the CC genotype, both for SBP (P = .001) and for DBP (P = .003), and these associations both surpassed the Bonferroni threshold for multiple comparisons. Multivariate regression analysis including known predictors of BP as independent variables was performed to evaluate the strength of this association, which in the case of the SNP rs884205 of the RANK gene remained statistically significant after adjustment for both SBP (P = .0006) and DBP (P = .005), demonstrating the key role of this SNP in BP.We report a robust association between the SNP rs884205 in RANK gene and BP in women, and this SNP is validated as a candidate in cardiovascular risk studies.

Regulation of Skeletal Homeostasis

Landmark advances in skeletal biology have arisen mainly from the identification of disease-causing mutations and the advent of rapid and selective gene-targeting technologies to phenocopy human disease in mice. Here, we discuss work on newly identified mechanisms controlling the remodeling of bone, communication of bone cells with cells of other lineages, and crosstalk between bone and vital organs as these relate to the therapeutic targeting of the skeleton.

Heterozygosity for the rs696217 SNP in the Preproghrelin Gene Predicts Weight Loss After Bariatric Surgery in Severely Obese Individuals

BACKGROUND

Several patients encompass a scarce weight loss after Roux-en-Y gastric bypass (RYGB). As such event is not related to surgical complications, finding markers able to identify "well responders" and to predict weight loss outcome is clinically relevant. Ghrelin regulates appetite and energy balance. Common single nucleotide polymorphisms (SNPs) in its encoding genes have been associated with body weight regulation. Other peptides involved in satiety modulation, like the CD40/CD40L complex, are less explored.

METHODS

One hundred, otherwise healthy, obese subjects (aged 45 ± 11 years, 65 females, BMI 48.0 ± 0.7 kg/m²) were sequentially enrolled in years 2014-2015. SNPs rs2241766 for adiponectin gene, rs490683 for ghrelin receptor, rs696217 and rs27647 for the preproghrelin/ghrelin gene, and rs1126535 for the CD40L gene were determined on DNA extracted from circulating lymphomonocytes. Patients were reevaluated at 6 (n = 100), 26 (n = 91), and 52 weeks (n = 79) after RYGB.

RESULTS

Subjects carrying the rs696217 T allele encompassed a significantly greater reduction in BMI 52 weeks after surgery (GG vs GT 30.5 ± 1.1 vs 38.1 ± 2.1 %; p < 0.001). Carrying the rs1126535 C allele in the CD40L gene was associated with a significantly lower BMI reduction at week 52 (TT vs CT 33.2 ± 1.1 vs 28.1 ± 2.3 %, p = 0.049). rs490683 and rs27647 SNPs of ghrelin and rs2241766 for adiponectin gene did not show any difference between carriers and non-carriers of the mutant allele.

CONCLUSION

Carrying a G to T substitution in rs696217 (preproghrelin gene) seems to mark a successful weight loss outcome; we also report for the first time that the rs1126535 C allele (CD40L gene) may predict a worse response to bariatric surgery.

The role of CD40 and CD40L in bone mineral density and in osteoporosis risk: A genetic and functional study

Compelling data are revealing that the CD40/CD40L system is involved in bone metabolism. Furthermore, we have previously demonstrated that polymorphisms in both genes are associated with bone phenotypes. The aim of this study is to further characterize this association and to identify the causal functional mechanism. We conducted an association study of BMD with 15 SNPs in CD40/CD40L genes in a population of 779 women. In addition, we assessed the functionality of this association through the study of the allele-dependent expression of CD40 and CD40L in peripheral blood leukocytes (PBLs) and in human osteoblasts (OBs) obtained from bone explants by qPCR and by sequencing. When an allelic imbalance (AI) was detected, studies on allele-dependent in vitro transcription rate and on CpG methylation in the gene promoter were also performed. Our results confirm the genetic association between SNP rs116535 (T>C) of CD40L gene with LS-BMD. Regarding CD40 gene, two SNPs showed nominal P-values<0.05 for FN- and LS-BMD (Z-scores), although the association was not significant after correcting for multiple testing. Homozygous TT women for SNP rs1883832 (C>T) of CD40 gene showed a trend to have lower levels of OPG (Q-value=0.059), especially when women of BMD-quartile ends were selected (P<0.05). Regarding functionality, we detected an AI for rs1883832 with the C allele the most expressed in OBs and in PBLs. Since the rs116535 of CD40L gene did not show AI, it was not further analyzed. Finally, we described a differential methylation of CpGs in the CD40 promoter among women of high in comparison to low BMD. Our results suggest that the CD40/CD40L system plays a role in regulating BMD. Effectively, our data suggest that a decreased production of OPG could be the cause of the lower BMD observed in TT women for rs1883832 of the CD40 gene and that the degree of methylation of CpGs in the CD40 promoter could contribute to the acquisition of BMD. One possibility that deserves further study is whether the degree of methylation of the CD40 gene affects the level of CD40 expression and, consequently, the level of OPG.

Are polymorphisms of the immunoregulatory factor CD40LG implicated in acute transfusion reactions?

The CD40 ligand (CD40L/CD154), a member of TNF superfamily, is notably expressed on activated CD4+ T-cells and stimulated platelets. CD40L is linked to a variety of pathologies and to acute transfusion reactions (ATR). Mutations in this gene (CD40LG) lead to X-linked hyper-IgM syndrome. Some CD40LG polymorphisms are associated with variable protein expression. The rationale behind this study is that CD40L protein has been observed to be involved in ATR. We wondered whether genetic polymorphisms are implicated. We investigated genetic diversity in the CD40LG using DHPLC and capillary electrophoresis for screening and genotyping (n = 485 French and Tunisian blood donors). We identified significant difference in the CD40LG linkage pattern between the two populations. Variant minor alleles were significantly over-represented in Tunisian donors (P<0.0001 to 0.0270). We found higher heterogeneity in the Tunisian, including three novel low frequency variants. As there was not a particular pattern of CD40LG in single apheresis donors whose platelet components induced an ATR, we discuss how this information may be useful for future disease association studies on CD40LG.

Replication study of three functional polymorphisms associated with bone mineral density in a cohort of Spanish women

Gene candidate and genome-wide association studies have revealed tens of loci of susceptibility for osteoporosis. Some limitations such as sample size, use of confounding variables, and control for multiple testing and for population stratification, however, represent common problems in these studies that make replication in independent cohorts desirable and even necessary. The main objective of the present study is to replicate previous data on three functional polymorphisms in a cohort of Spanish women. To that end, we performed an association study of three functional polymorphisms previously associated with bone phenotypes in the LRP5, TNFRSF11B, and FGFBP1 genes with low bone mineral density (BMD) in a cohort of 721 Spanish women, most of them postmenopausal. We detected a strong significant association, even when correcting for multiple comparisons, for polymorphism rs312009 in the LRP5 gene with low BMD at the lumbar-spine site. These were women with the CC genotype, which showed the worst bone parameters. Moreover, these women had a higher risk of osteoporosis (adjusted odds ratio 2.82, P = 0.001) than women with the TT/TC genotype. This association seems to be caused because the rs312009 single nucleotide polymorphism (SNP) is located at a binding site for the transcription factor RUNX2 at the 5' region of the LRP5 gene, and the T allele seems to be a better transcriber than the C allele. Regarding the other two SNPs, only the rs4876869 SNP in the TNFRSF11B gene showed a suggestive trend for both skeletal sites. These results underscore the significance of the LRP5 gene in bone metabolism and emphasize the significance of the replication of previous results in independent cohorts.

Lack of association between CD40 polymorphisms and acute rejection in German liver transplant recipients

CD40 and its ligand, CD154, are major costimulatory molecules whose interactions are important in alloreactive transplant rejection. The aim of this study was to examine the association of CD40 polymorphisms with the susceptibility to acute rejection episodes in liver transplantation. In total, 112 liver transplant recipients with biopsy proven acute rejections (BPAR), 97 without BPAR (WBPAR), and 112 healthy control individuals were enrolled in the study. Two single nucleotide polymorphisms (SNPs) of CD40 gene (rs1883832 and rs4810485) were genotyped by polymerase chain reaction-allele specific restriction enzyme analysis (PCR-ASRA). Both SNPs has been tested for a recessive and a dominant model. No significant differences were found in the genotype and allele frequencies of the SNPs rs1883832 and rs4810485 between BPAR liver recipients and WBPAR recipients. Our results do not suggest an important role of tested CD40 SNPs in the susceptibility to acute liver transplant rejection in a Caucasian population.

Gene expression profile induced by ovariectomy in bone marrow of mice: a functional approach to identify new candidate genes associated to osteoporosis risk in women

Osteoporosis is a multifactorial skeletal pathology with a main genetic component. To date, however, the majority of genes associated with this pathology remain unknown since genes cataloged to date only explain a part of the heritability of bone phenotypes. In the present study, we have used a genome-wide gene expression approach by means of microarrays to identify new candidate genes involved in the physiopathology of osteoporosis, using as a model the ovariectomized (OVX) mice by comparing global bone marrow gene expression of the OVX mice with those of SHAM operated mice. One hundred and eighty transcripts were found to be differentially expressed between groups. The analysis showed 23 significant regulatory networks, of which the top five canonical pathways included B-cell development, primary immunodeficiency signaling, PI3K signaling in B-cells, phospholipase C signaling, and FcgRIIB signaling in B-cells. Twelve differentially expressed genes were validated by MALDI-TOF mass spectrometry with good reproducibility. Finally, the association to bone phenotypes of SNPs in genes whose expression was increased (IL7R and CD79A) or decreased (GPX3 and IRAK3) by OVX in mice was analyzed in a cohort of 706 postmenopausal women. We detected an association of a SNP in a gene involved in the detoxification of free radicals like glutathione peroxidase 3 (GPX3) with femoral neck BMD (rs8177447, P=0.043) and two SNPs in the Ig-alpha protein of the B-cell antigen component gene (CD79A) with lumbar spine BMD (rs3810153 and rs1428922, P=0.016 and P=0.001, respectively). These results reinforce the role of antioxidant pathways and of B-cells in bone metabolism. Furthermore, it shows that a genome-wide gene expression approach in animal models is a useful method for detecting genes associated to BMD and osteoporosis risk in humans.

Investigation of gene-gene interactions between CD40 and CD40L in Polish multiple sclerosis patients

CD40-CD40L interaction is necessary for the activation of both humoral and cellular immune response and has been suggested to play a role in the pathogenesis of multiple sclerosis (MS). Therefore, we analyzed the combined influence of the CD40 and CD40L variants on MS susceptibility and progression on well-defined Polish population. Our investigation revealed that CT individuals in rs1883832 locus of CD40 possessed almost 1.5-fold higher risk for MS than CC individuals (OR = 1.44; 95%CI = 1.03-2.1; p = 0.032), while this risk for TT individuals was almost 2.5-fold higher (OR = 2.36; 95%CI = 1.19-4.78; p = 0.014). Moreover, for the first time, we observed the association of CD40 gene with MS development and progression. We observed that for the rs1883832CC individuals the age at diagnosis was on average 2 years lower than for the rs1883832CT and rs1883832TT individuals (CI95% = -3.69-(-0.29); p = 0.023). Additionally, we detected that individuals with TT and CT genotypes showed lower risk of developing secondary progressive course in comparison to those with CC genotype. For rs1883832TT individuals this risk was 4-fold lower (HR = 0.24; CI95% = 0.10-0.53; p = 0.00062). Despite the fact that CD40-CD40L pathway plays a key role in development of autoimmune diseases, we were not able to detect gene-gene interactions between CD40 and CD40L polymorphisms associated with multiple sclerosis.

Genome-wide approaches for identifying genetic risk factors for osteoporosis

Osteoporosis, the most common type of bone disease worldwide, is clinically characterized by low bone mineral density (BMD) and increased susceptibility to fracture. Multiple genetic and environmental factors and gene-environment interactions have been implicated in its pathogenesis. Osteoporosis has strong genetic determination, with the heritability of BMD estimated to be as high as 60%. More than 80 genes or genetic variants have been implicated in risk of osteoporosis by hypothesis-free genome-wide studies. However, these genes or genetic variants can only explain a small portion of BMD variation, suggesting that many other genes or genetic variants underlying osteoporosis risk await discovery. Here, we review recent progress in genome-wide studies of osteoporosis and discuss their implications for medicine and the major challenges in the field.

Association of SNPs of CD40 gene with multiple sclerosis in Russians

Multiple sclerosis (MS) is a serious, incurable neurological disease. In 2009, the ANZgene studies detected the suggestive association of located upstream of CD40 gene in chromosome 20q13 (p = 1.3×10(-7)). Identification of the causal variant(s) in the CD40 locus leads to a better understanding of the mechanism underlying the development of autoimmune pathologies. We determined the genotypes of rs6074022, rs1883832, rs1535045, and rs11086996 in patients with MS (n = 1684) and in the control group (n = 879). Two SNPs were significantly associated with MS: rs6074022 (additive model C allele OR = 1.27, 95% CI = [1.12-1.45], p = 3×10(-4)) and rs1883832 (additive model T allele OR = 1.20, 95% CI = [1.05-1.38], p = 7×10(-3)). In the meta-analysis of our results and the results of four previous studies, we obtain the association p-value of 2.34×10(-12), which confirmed the association between MS and rs6074022 at a genome-wide significant level. Next, we demonstrated that the model including rs6074022 only sufficiently described the association. From our analysis, we can speculate that the association between rs1883832 and MS was induced by LD, whereas rs6074022 was a marker in stronger LD with the functional variant or was the functional variant itself. Our results indicated that the functional variants were located in the upstream region of the gene CD40 and were in higher LD with rs6074022 than LD with rs1883832.

Gene-gene interaction between RBMS3 and ZNF516 influences bone mineral density

Osteoporosis is characterized by low bone mineral density (BMD), a highly heritable trait that is determined, in part, by the actions and interactions of multiple genes. Although an increasing number of genes have been identified to have independent effects on BMD, few studies have been performed to identify genes that interact with one another to affect BMD. In this study, we performed gene-gene interaction analyses in selected candidate genes in individuals with extremely high versus low hip BMD (20% tails of the distributions), in two independent U.S. Caucasian samples. The first sample contained 916 unrelated subjects with extreme hip BMD Z-scores selected from a population composed of 2286 subjects. The second sample consisted of 400 unrelated subjects with extreme hip BMD Z-scores selected from a population composed of 1000 subjects. Combining results from these two samples, we found one interacting gene pair (RBMS3 versus ZNF516) which, even after Bonferroni correction for multiple testing, showed consistently significant effects on hip BMD. RMBS3 harbored two single-nucleotide polymorphisms (SNPs), rs6549904 and rs7640046, both of which had significant interactions with an SNP, rs4891159, located on ZNF516 (p = 7.04 × 10(-11) and 1.03 × 10(-10) ). We further validated these results in two additional samples of Caucasian and African descent. The gene pair, RBMS3 versus ZNF516, was successfully replicated in the Caucasian sample (p = 8.07 × 10(-3) and 2.91 × 10(-3) ). For the African sample, a significant interaction was also detected (p = 0.031 and 0.043), but the direction of the effect was opposite to that observed in the three Caucasian samples. By providing evidence for genetic interactions underlying BMD, this study further delineates the genetic architecture of osteoporosis.

Effects of administration of hormone therapy or raloxifene on the immune system and on biochemical markers of bone remodeling

OBJECTIVE

Over the last few years, conclusive evidence of the involvement of the immune system in the regulation of bone metabolism has been identified. Consequently, one question that should be formulated concerns the possible effects of antiresorptive therapies on the immune system. Therefore, the purpose of the present work was to evaluate both the functionality of the immune system and bone turnover in women receiving antiresorptive therapies, such as hormone therapy (HT; n = 33) and raloxifene (RLX; n = 66), acting through estrogen receptors.

METHODS

To that end, this study analyzed bone turnover markers in a population of postmenopausal women before and after beginning therapy and compared these with data of women not treated (NT; n = 102). In a subgroup of participants (NT = 33, RLX = 24, and HT = 26), we analyzed the effects of treatments on immune system parameters such as serum levels of interleukin (IL)-6, tumor necrosis factor α, and IL-1β; lymphocyte subpopulations; cell proliferation by peripheral blood mononuclear cells (PBMCs); in vitro production of IL-1β by PBMCs; and the expression of receptor activator of nuclear factor-κB ligand, transforming growth factor β, and IL-4 genes by PBMCs.

RESULTS

The results showed that bone resorption was inhibited strongly in women in the RLX and HT groups when compared with women in the NT group. Interestingly, the administration of RLX inhibited the production of the Wnt/β-catenin signaling pathway inhibitor Dickkopf Homolog-1 (P < 0.05) and tended to increase the levels of the osteoclastogenesis inhibitor osteoprotegerin at month 6 (P = 0.059). With regard to the immune system, the different treatments did not markedly perturb the parameters analyzed, with the exception of the increase in serum IL-1β detected in the HT group at month 6 (P < 0.05).

CONCLUSIONS

The main conclusions of the present work were that HT or RLX do not disturb the immune system and that both treatments have a similar antiresorptive power.