Detection of significant pathways in osteoporosis based on graph clustering

Identification of common genetic characteristics of rheumatoid arthritis and major depressive disorder by bioinformatics analysis and machine learning

Introduction

Depression is the most common comorbidity of rheumatoid arthritis (RA). In particular, major depressive disorder (MDD) and rheumatoid arthritis share highly overlapping mental and physical manifestations, such as depressed mood, sleep disturbance, fatigue, pain, and worthlessness. This overlap and indistinguishability often lead to the misattribution of physical and mental symptoms of RA patients to depression, and even, the depressive symptoms of MDD patients are ignored when receiving RA treatment. This has serious consequences, since the development of objective diagnostic tools to distinguish psychiatric symptoms from similar symptoms caused by physical diseases is urgent.

Methods

Bioinformatics analysis and machine learning.

Results

The common genetic characteristics of rheumatoid arthritis and major depressive disorder are EAF1, SDCBP and RNF19B.

Discussion

We discovered a connection between RA and MDD through immune infiltration studies: monocyte infiltration. Futhermore, we explored the correlation between the expression of the 3 marker genes and immune cell infiltration using the TIMER 2.0 database. This may help to explain the potential molecular mechanism by which RA and MDD increase the morbidity of each other.

The selected genes NR6A1, RSAD2-CMPK2, and COL3A1 contribute to body size variation in Meishan pigs through different patterns

The high-fertility Meishan pig is currently categorized into medium sized (MMS) and small sized (SMS) based on body size. To identify causal genes responsible for the variation in body size within the two categories, we sequenced individuals representing the entire consanguinity of the existing Meishan pig. This enabled us to conduct genome selective signal analysis. Our findings revealed the genomes of MMS and SMS are stratified, with selective sweep regions formed by differential genomic intervals between the two categories enriched in multiple pig body size related quantitative trait loci (QTLs). Furthermore, the missense mutation c.575T > C of candidate causal gene NR6A1, accounting for the variation in lumbar vertebrae number in pigs, was positively selected in MMS only, leading to an increase in body length of MMS at 6 months of age. To precisely identify causal genes accounting for body size variation through multi-omics, we collected femoral cartilage and liver transcription data from MMS and SMS respectively, and re-sequencing data from pig breeds exhibiting varying body sizes. We found that two selected regions where the RSAD2-CMPK2 and COL3A1 genes are located, respectively, showed different haplotypes in pig breeds of varying body size, and was associated with body or carcass length in hybridized Suhuai pig. Additionally, the above three hub genes, were significantly greater expressed in SMS femoral cartilage and liver tissues compared to MMS. These three genes could strengthen the pathways related to bone resorption and metabolism in SMS, potentially hindering bone and skeletal development and resulting in a smaller body size in SMS. These findings provide valuable insights into the genetic mechanism of body size variation in Meishan pig population.

RNAseq and RNA molecular barcoding reveal differential gene expression in cortical bone following hindlimb unloading in female mice

Disuse-induced bone loss is seen following spinal cord injury, prolonged bed rest, and exposure to microgravity. We performed whole transcriptomic profiling of cortical bone using RNA sequencing (RNAseq) and RNA molecular barcoding (NanoString) on a hindlimb unloading (HLU) mouse model to identify genes whose mRNA transcript abundances change in response to disuse. Eleven-week old female C57BL/6 mice were exposed to ambulatory loading or HLU for 7 days (n = 8/group). Total RNA from marrow-flushed femoral cortical bone was analyzed on HiSeq and NanoString platforms. The expression of several previously reported genes associated with Wnt signaling and metabolism was altered by HLU. Furthermore, the increased abundance of transcripts, such as Pfkfb3 and Mss51, after HLU imply these genes also have roles in the cortical bone’s response to altered mechanical loading. Our study demonstrates that an unbiased approach to assess the whole transcriptomic profile of cortical bone can reveal previously unidentified mechanosensitive genes and may eventually lead to novel targets to prevent disuse-induced osteoporosis.

Natural Killer Lytic-Associated Molecule (NKLAM): An E3 Ubiquitin Ligase With an Integral Role in Innate Immunity

Natural Killer Lytic-Associated Molecule (NKLAM), also designated RNF19B, is a unique member of a small family of E3 ubiquitin ligases. This 14-member group of ligases has a characteristic cysteine-rich RING-IBR-RING (RBR) domain that mediates the ubiquitination of multiple substrates. The consequence of substrate ubiquitination varies, depending on the type of ubiquitin linkages formed. The most widely studied effect of ubiquitination of proteins is proteasome-mediated substrate degradation; however, ubiquitination can also alter protein localization and function. Since its discovery in 1999, much has been deciphered about the role of NKLAM in innate immune responses. We have discerned that NKLAM has an integral function in both natural killer (NK) cells and macrophages in vitro and in vivo. NKLAM expression is required for each of these cell types to mediate maximal killing activity and cytokine production. However, much remains to be determined. In this review, we summarize what has been learned about NKLAM expression, structure and function, and discuss new directions for investigation. We hope that this will stimulate interest in further exploration of NKLAM.

An eQTL in the cystathionine beta synthase gene is linked to osteoporosis in laying hens

BACKGROUND

Skeletal damage is a challenge for laying hens because the physiological adaptations required for egg laying make them susceptible to osteoporosis. Previously, we showed that genetic factors explain 40% of the variation in end of lay bone quality and we detected a quantitative trait locus (QTL) of large effect on chicken chromosome 1. The aim of this study was to combine data from the commercial founder White Leghorn population and the F2 mapping population to fine-map this QTL and understand its function in terms of gene expression and physiology.

RESULTS

Several single nucleotide polymorphisms on chromosome 1 between 104 and 110 Mb (galGal6) had highly significant associations with tibial breaking strength. The alternative genotypes of markers of large effect that flanked the region had tibial breaking strengths of 200.4 vs. 218.1 Newton (P < 0.002) and, in a subsequent founder generation, the higher breaking strength genotype was again associated with higher breaking strength. In a subsequent generation, cortical bone density and volume were increased in individuals with the better bone genotype but with significantly reduced medullary bone quality. The effects on cortical bone density were confirmed in a further generation and was accompanied by increased mineral maturity of the cortical bone as measured by infrared spectrometry and there was evidence of better collagen cross-linking in the cortical bone. Comparing the transcriptome of the tibia from individuals with good or poor bone quality genotypes indicated four differentially-expressed genes at the locus, one gene, cystathionine beta synthase (CBS), having a nine-fold higher expression in the genotype for low bone quality. The mechanism was cis-acting and although there was an amino-acid difference in the CBS protein between the genotypes, there was no difference in the activity of the enzyme. Plasma homocysteine concentration, the substrate of CBS, was higher in the poor bone quality genotype.

CONCLUSIONS

Validated markers that predict bone strength have been defined for selective breeding and a gene was identified that may suggest alternative ways to improve bone health in addition to genetic selection. The identification of how genetic variants affect different aspects of bone turnover shows potential for translational medicine.

Identification of a key gene module associated with glucocorticoid- induced derangement in bone mineral density in patients with asthma

Derangement in bone mineral density (BMD) caused by glucocorticoid is well-known. The present study aimed to find key biological pathways associated with low BMD after glucocorticoid treatment in asthmatics using gene expression profiles of peripheral blood cells. We utilized immortalized B cells (IBCs) from 32 childhood asthmatics after multiple oral glucocorticoid bursts and peripheral blood mononuclear cells (PBMCs) from 17 adult asthmatics after a long-term use of oral glucocorticoid. We searched co-expressed gene modules significantly related with the BMD Z score in childhood asthmatics and tested if these gene modules were preserved and significantly associated with the BMD Z score in adult asthmatics as well. We identified a gene module composed of 199 genes significantly associated with low BMD in both childhood and adult asthmatics. The structure of this module was preserved across gene expression profiles. We found that the cellular metabolic pathway was significantly enriched in this module. Among 18 hub genes in this module, we postulated that 2 genes, CREBBP and EP300, contributed to low BMD following a literature review. A novel biologic pathway identified in this study highlighted a gene module and several genes as playing possible roles in the pathogenesis of glucocorticoid- induced derangement in BMD.

"Omics" Signatures in Peripheral Monocytes from Women with Low BMD Condition

Postmenopausal osteoporosis (PMO) is a result of increased bone resorption compared to formation. Osteoclasts are responsible for bone resorption, which are derived from circulating monocytes that undertake a journey from the blood to the bone for the process of osteoclastogenesis. In recent times, the use of high throughput technologies to explore monocytes from women with low versus high bone density has led to the identification of candidate molecules that may be deregulated in PMO. This review provides a list of molecules in monocytes relevant to bone density which have been identified by "omics" studies in the last decade or so. The molecules in monocytes that are deregulated in low BMD condition may contribute to processes such as monocyte survival, migration/chemotaxis, adhesion, transendothelial migration, and differentiation into the osteoclast lineage. Each of these processes may be crucial to the overall route of osteoclastogenesis and an increase in any/all of these processes can lead to increased bone resorption and subsequently low bone density. Whether these molecules are indeed the cause or effect is an arena currently unexplored.

Transcriptional profile in afferent lymph cells following vaccination with liposomes incorporating CpG

Vaccine formulations incorporating innate immune stimulants are highly immunogenic; however, the biological signals that originate in the peripheral tissues at the site of injection and are transmitted to the local lymph node to induce immunity remain unclear. By directly cannulating the ovine afferent lymphatic vessels, we have previously shown that it takes 72 hr for mature antigen-loaded dendritic cells and monocytes to appear within afferent lymph following injection of a liposomal formulation containing the Toll-like receptor ligand CpG. In this present study, we characterize the global transcriptional signatures at this time-point in ovine afferent lymph cells as they migrate from the injection site into the lymphatics following vaccination with a liposome antigen formulation incorporating CpG. We show that at 72 hr post vaccination, liposomes alone induce no changes in gene expression and inflammatory profiles within afferent lymph; however, the incorporation of CpG drives interferon, antiviral and cytotoxic gene programmes. This study also measures the expression of key genes within individual cell types in afferent lymph. Antiviral gene signatures are most prominent in lymphocytes, which may play a significant and unexpected role in sustaining the immune response to vaccination at the site of injection. These findings provide a comprehensive analysis of the in vivo immunological pathways that connect the injection site with the local draining lymph node following vaccination.

Manganese superoxide dismutase is required to maintain osteoclast differentiation and function under static force

Bone homeostasis is maintained by the balance of osteoblasts (OBs) and osteoclasts (OCs). Increased activity of OCs not only contributes to pathological bone resorption, such as osteoporosis and periodontitis, but also is responsible for physiological conditions like orthodontic tooth movement (OTM). However, the detailed mechanism by which orthodontic force promotes the formation of OCs is still poorly understood. In this study, we confirmed that static force promoted the differentiation of human cord monocytes (HMNCs) into OCs depending on loading time and magnitude. Protein expression profiles among HMNCs, HMNCs subjected to static force and mature OCs were established via 2-DE and MALDI-TOF-MS analyses. Total respective protein spot numbers of 549 ± 13, 612 ± 19 and 634 ± 16 were detected in each of the gels by image analysis. The five proteins identified were plasminogen activator inhibitor 2 (PAI-2, Spot 1), peroxiredoxin-6 (PRD-6, Spot 3), manganese superoxide dismutase (SOD2, Spot 6), Rho GDP-dissociation inhibitor 2 (Rho-GDI2, Spot 11) and L-lactate dehydrogenase B chain (L-LDH, Spot 15). More importantly, we revealed that SOD2 was required to maintain monocyte differentiation into functional OCs and may become a potential target for regulating the efficiency of OTM in the future.

Meta-analysis of genome-wide association studies identifies six new Loci for serum calcium concentrations

Calcium is vital to the normal functioning of multiple organ systems and its serum concentration is tightly regulated. Apart from CASR, the genes associated with serum calcium are largely unknown. We conducted a genome-wide association meta-analysis of 39,400 individuals from 17 population-based cohorts and investigated the 14 most strongly associated loci in ≤ 21,679 additional individuals. Seven loci (six new regions) in association with serum calcium were identified and replicated. Rs1570669 near CYP24A1 (P = 9.1E-12), rs10491003 upstream of GATA3 (P = 4.8E-09) and rs7481584 in CARS (P = 1.2E-10) implicate regions involved in Mendelian calcemic disorders: Rs1550532 in DGKD (P = 8.2E-11), also associated with bone density, and rs7336933 near DGKH/KIAA0564 (P = 9.1E-10) are near genes that encode distinct isoforms of diacylglycerol kinase. Rs780094 is in GCKR. We characterized the expression of these genes in gut, kidney, and bone, and demonstrate modulation of gene expression in bone in response to dietary calcium in mice. Our results shed new light on the genetics of calcium homeostasis.