Network-Based Meta-Analyses of Associations of Multiple Gene Expression Profiles with Bone Mineral Density Variations in Women

Background

Existing microarray studies of bone mineral density (BMD) have been critical for understanding the pathophysiology of osteoporosis, and have identified a number of candidate genes. However, these studies were limited by their relatively small sample sizes and were usually analyzed individually. Here, we propose a novel network-based meta-analysis approach that combines data across six microarray studies to identify functional modules from human protein-protein interaction (PPI) data, and highlight several differentially expressed genes (DEGs) and a functional module that may play an important role in BMD regulation in women.

Methods

Expression profiling studies were identified by searching PubMed, Gene Expression Omnibus (GEO) and ArrayExpress. Two meta-analysis methods were applied across different gene expression profiling studies. The first, a nonparametric Fisher’s method, combined p-values from individual experiments to identify genes with large effect sizes. The second method combined effect sizes from individual datasets into a meta-effect size to gain a higher precision of effect size estimation across all datasets. Genes with Q test’s p-values < 0.05 or I² values > 50% were assessed by a random effects model and the remainder by a fixed effects model. Using Fisher’s combined p-values, functional modules were identified through an integrated analysis of microarray data in the context of large protein–protein interaction (PPI) networks. Two previously published meta-analysis studies of genome-wide association (GWA) datasets were used to determine whether these module genes were genetically associated with BMD. Pathway enrichment analysis was performed with a hypergeometric test.

Results

Six gene expression datasets were identified, which included a total of 249 (129 high BMD and 120 low BMD) female subjects. Using a network-based meta-analysis, a consensus module containing 58 genes (nodes) and 83 edges was detected. Pathway enrichment analysis of the 58 module genes revealed that these genes were enriched in several important KEGG pathways including Osteoclast differentiation, B cell receptor signaling pathway, MAPK signaling pathway, Chemokine signaling pathway and Insulin signaling pathway. The importance of module genes was replicated by demonstrating that most module genes were genetically associated with BMD in the GWAS data sets. Meta-analyses were performed at the individual gene level by combining p-values and effect sizes. Five candidate genes (ESR1, MAP3K3, PYGM, RAC1 and SYK) were identified based on gene expression meta-analysis, and their associations with BMD were also replicated by two BMD meta-analysis studies.

Conclusions

In summary, our network-based meta-analysis not only identified important differentially expressed genes but also discovered biologically meaningful functional modules for BMD determination. Our study may provide novel therapeutic targets for osteoporosis in women.     

Divergences in gene repertoire among the reference Prevotella genomes derived from distinct body sites of human

Background

The community composition of the human microbiome is known to vary at distinct anatomical niches. But little is known about the nature of variations, if any, at the genome/sub-genome levels of a specific microbial community across different niches. The present report aims to explore, as a case study, the variations in gene repertoire of 28 Prevotella reference genomes derived from different body-sites of human, as reported earlier by the Human Microbiome Consortium.

Results

The pan-genome for Prevotella remains “open”. On an average, 17% of predicted protein-coding genes of any particular Prevotella genome represent the conserved core genes, while the remaining 83% contribute to the flexible and singletons. The study reveals exclusive presence of 11798, 3673, 3348 and 934 gene families and exclusive absence of 17, 221, 115 and 645 gene families in Prevotella genomes derived from human oral cavity, gastro-intestinal tracts (GIT), urogenital tract (UGT) and skin, respectively. Distribution of various functional COG categories differs significantly among the habitat-specific genes. No niche-specific variations could be observed in distribution of KEGG pathways.

Conclusions

Prevotella genomes derived from different body sites differ appreciably in gene repertoire, suggesting that these microbiome components might have developed distinct genetic strategies for niche adaptation within the host. Each individual microbe might also have a component of its own genetic machinery for host adaptation, as appeared from the huge number of singletons.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1350-6) contains supplementary material, which is available to authorized users.     

Prospective Multicenter Study of Community-Associated Skin and Skin Structure Infections due to Methicillin-Resistant Staphylococcus aureus in Buenos Aires,Argentina

Background

Community-associated methicillin-resistant Staphylococcus aureus (CA-MRSA) is now the most common cause of skin and skin structure infections (SSSI) in several world regions. In Argentina prospective, multicenter clinical studies have only been conducted in pediatric populations.

Objective

Primary: describe the prevalence, clinical and demographic characteristics of adult patients with community acquired SSSI due to MRSA; secondary: molecular evaluation of CA-MRSA strains. Patients with MRSA were compared to those without MRSA.

Materials and Methods

Prospective, observational, multicenter, epidemiologic study, with molecular analysis, conducted at 19 sites in Argentina (18 in Buenos Aires) between March 2010 and October 2011. Patients were included if they were ≥14 years, were diagnosed with SSSI, a culture was obtained, and there had no significant healthcare contact identified. A logistic regression model was used to identify factors associated with CA-MRSA. Pulse field types, SCCmec, and PVL status were also determined.

Results

A total of 311 patients were included. CA-MRSA was isolated in 70% (218/311) of patients. Clinical variables independently associated with CA-MRSA were: presence of purulent lesion (OR 3.29; 95%CI 1.67, 6.49) and age <50 years (OR 2.39; 95%CI 1.22, 4.70). The vast majority of CA-MRSA strains causing SSSI carried PVL genes (95%) and were SCCmec type IV. The sequence type CA-MRSA ST30 spa t019 was the predominant clone.

Conclusions

CA-MRSA is now the most common cause of SSSI in our adult patients without healthcare contact. ST30, SCCmec IV, PVL+, spa t019 is the predominant clone in Buenos Aires, Argentina.     

Comparing Cutaneous Research Funded by the National Institute of Arthritis and Musculoskeletal and Skin Diseases with 2010 Global Burden of Disease Results

Importance

Disease burden data helps guide research prioritization.

Objective

To determine the extent to which grants issued by the National Institute of Arthritis and Musculoskeletal and Skin Diseases (NIAMS) reflect disease burden, measured by disability-adjusted life years (DALYs) from Global Burden of Disease (GBD) 2010 project.

Design

Two investigators independently assessed 15 skin conditions studied by GBD 2010 in the NIAMS database for grants issued in 2013. The 15 skin diseases were matched to their respective DALYs from GBD 2010.

Setting

The United States NIAMS database and GBD 2010 skin condition disability data.

Main Outcome(s) and Measure(s)

Relationship of NIAMS grant database topic funding with percent total GBD 2010 DALY and DALY rank for 15 skin conditions.

Results

During fiscal year 2013, 1,443 NIAMS grants were issued at a total value of $424 million. Of these grants, 17.7% covered skin topics. Of the total skin disease funding, 82% (91 grants) were categorized as “general cutaneous research.” Psoriasis, leprosy, and “other skin and subcutaneous diseases” (ie; immunobullous disorders, vitiligo, and hidradenitis suppurativa) were over-represented when funding was compared with disability. Conversely, cellulitis, decubitus ulcer, urticaria, acne vulgaris, viral skin diseases, fungal skin diseases, scabies, and melanoma were under-represented. Conditions for which disability and funding appeared well-matched were dermatitis, squamous and basal cell carcinoma, pruritus, bacterial skin diseases, and alopecia areata.

Conclusions and Relevance

Degree of representation in NIAMS is partly correlated with DALY metrics. Grant funding was well-matched with disability metrics for five of the 15 studied skin diseases, while two skin diseases were over-represented and seven were under-represented. Global burden estimates provide increasingly transparent and important information for investigating and prioritizing national research funding allocations.     

Microarray Analysis of Gene Expression Profiles of Schistosoma japonicum Derived from Less-Susceptible Host Water Buffalo and Susceptible Host Goat

Background

Water buffalo and goats are natural hosts for S. japonicum in endemic areas of China. The susceptibility of these two hosts to schistosome infection is different, as water buffalo are less conducive to S. japonicum growth and development. To identify genes that may affect schistosome development and survival, we compared gene expression profiles of schistosomes derived from these two natural hosts using high-throughput microarray technology.

Results

The worm recovery rate was lower and the length and width of worms from water buffalo were smaller compared to those from goats following S. japonicum infection for 7 weeks. Besides obvious morphological difference between the schistosomes derived from the two hosts, differences were also observed by scanning and transmission electron microscopy. Microarray analysis showed differentially expressed gene patterns for parasites from the two hosts, which revealed that genes related to lipid and nucleotide metabolism, as well as protein folding, sorting, and degradation were upregulated, while others associated with signal transduction, endocrine function, development, immune function, endocytosis, and amino acid/carbohydrate/glycan metabolism were downregulated in schistosomes from water buffalo. KEGG pathway analysis deduced that the differentially expressed genes mainly involved lipid metabolism, the MAPK and ErbB signaling pathways, progesterone-mediated oocyte maturation, dorso-ventral axis formation, reproduction, and endocytosis, etc.

Conclusion

The microarray gene analysis in schistosomes derived from water buffalo and goats provide a useful platform to disclose differences determining S. japonicum host compatibility to better understand the interplay between natural hosts and parasites, and identify schistosome target genes associated with susceptibility to screen vaccine candidates.     

Immunolocalization of the Tumor-Sensitive Calmodulin-Like Protein CALML3 in Normal Human Skin and Hyperproliferative Skin Disorders

Background and Objective

Calmodulin-like protein CALML3 is an epithelial-specific protein regulated during keratinocyte differentiation in vitro. CALML3 expression is downregulated in breast cancers and transformed cell lines making it an attractive marker for tumor formation. The objective of this study was to survey CALML3 localization in normal epidermis and in hyperproliferative skin diseases including actinic keratosis, squamous and basal cell carcinoma as well as verruca and psoriasis and to compare CALML3 immunoreactivity with the proliferation marker Ki-67.

Methods

Paraffin-embedded tissue sections from normal human skin and hyperproliferative skin disorders were examined by immunohistochemistry and analyzed for localization and expression of CALML3 and Ki-67.

Results

CALML3 was strongly expressed in differentiating layers of normal skin, staining the periphery in suprabasal cells and exhibiting nuclear localization in the stratum granulosum. CALML3 nuclear localization was inversely correlated to Ki-67 staining in each disease, indicating that CALML3 nuclear presence is related to terminal cell differentiation and postmitotic state.

Conclusions

Increased CALML3 expression in suprabasal layers is characteristic for differentiating keratinocytes in normal epidermis, and nuclear expression of CALML3 inversely correlates with expression of the proliferation marker Ki-67. This suggests that CALML3 is a useful marker for normal and benign hyperplastic epidermal development, whereas the loss of nuclear CALML3 indicates progression to a proliferative and potentially malignant phenotype.     

Dietary monounsaturated Fatty acids intake and risk of skin photoaging

Background

Intake of monounsaturated fatty acids has been reported to reduce oxidative stress, insulin resistance and related inflammatory processes and may thus protect from skin photoaging. The objective of this study was to investigate the association between the risk of photoaging, monounsaturated fatty acids intake and the sources of monounsaturated fatty acids.

Methodology/Principal Findings

A cross sectional study was conducted within the framework of the SUVIMAX cohort. The survey included 1264 women and 1655 men aged between 45 and 60 years old. Dietary monounsaturated fatty acids intakes were estimated by dietary source through at least ten 24-h diet records completed during the first 2.5 years of the follow-up period. Severity of facial skin photoaging was graded by trained investigators at baseline during a clinical examination using a 6-grade scale illustrated by photographs. A lower risk of severe photoaging was associated with higher intakes of monounsaturated fatty acids from olive oil in both sexes. Strikingly, no association was found with intake of monounsaturated fatty acids from animal sources whether from dairy products, meat or processed meat.

Conclusion/Significance

These findings support the beneficial effect of dietary olive oil or healthy diet habits associated with olive oil consumption on the severity of facial photoaging.     

Intestine-Targeted DGAT1 Inhibition Improves Obesity and Insulin Resistance without Skin Aberrations in Mice

Objective

Diacylglycerol O-acyltransferase 1 (DGAT1) catalyzes the final committed step in triglyceride biosynthesis. DGAT1 null mice are known to be resistant to diet-induced obesity, and more insulin sensitive relative to the wild-type; however, the mice exhibit abnormalities in the skin. This work determined whether the intestine-targeted DGAT1 inhibitor could improve obesity and insulin resistance without skin aberrations in mice.

Design and Methods

We synthesized 2 DGAT1 inhibitors: Compound A, described in the patent application from the Japan Tobacco, and Compound B (A-922500), reported by Abbott Laboratories. Both compounds were evaluated for inhibitory activities against DGAT1 enzymes and effects on the skin in mice in vivo. Compound B was further investigated for effects on obesity and insulin resistance in diet-induced-obese (DIO) mice.

Results

The 2 compounds comparably inhibited the DGAT1 enzyme activity and the cellular triglyceride synthesis in vitro, while they showed different distribution patterns in mice in vivo. Compound A, which distributed systemically, caused skin aberrations, while Compound B, which preferentially distributed to the intestine, improved obesity and insulin resistance without skin aberrations in DIO mice.

Conclusions

Our results suggest that the intestine is the key tissue in which DGAT1 plays a role in promoting obesity and insulin resistance.     

Circadian rhythm and cartilage extracellular matrix genes in osseointegration: a genome-wide screening of implant failure by vitamin D deficiency

Background

Successful dental and orthopedic implants require the establishment of an intimate association with bone tissue; however, the mechanistic explanation of how biological systems accomplish osseointegration is still incomplete. We sought to identify critical gene networks involved in osseointegration by exploring the implant failure model under vitamin D deficiency.

Methodology

Adult male Sprague-Dawley rats were exposed to control or vitamin D-deficient diet prior to the osteotomy surgery in the femur bone and the placement of T-shaped Ti4Al6V implant. Two weeks after the osteotomy and implant placement, tissue formed at the osteotomy site or in the hollow chamber of T-shaped implant was harvested and total RNA was evaluated by whole genome microarray analyses.

Principal Findings

Two-way ANOVA of microarray data identified 103 genes that were significantly (>2 fold) modulated by the implant placement and vitamin D deficiency. Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses assigned the highest z-score to the circadian rhythm pathway including neuronal PAS domain 2 (NPAS2), and period homolog 2 (Per2). NPAS2 and Aryl hydrocarbon receptor nuclear translocator-like (ARNTL/Bmal 1) were upregulated around implant and diminished by vitamin D deficiency, whereas the expression pattern of Per2 was complementary. Hierarchical cluster analysis further revealed that NPAS2 was in a group predominantly composed of cartilage extracellular matrix (ECM) genes. Whereas the expression of bone ECM genes around implant was not significantly affected by vitamin D deficiency, cartilage ECM genes were modulated by the presence of the implant and vitamin D status. In a proof-of-concept in vitro study, the expression of cartilage type II and X collagens was found upregulated when mouse mesenchymal stem cells were cultured on implant disk with 1,25D supplementation.

Conclusions

This study suggests that the circadian rhythm system and cartilage extracellular matrix may be involved in the establishment of osseointegration under vitamin D regulation.     

Differential Gene Repertoire in Mycobacterium ulcerans Identifies Candidate Genes for Patho-Adaptation

Background

Based on large genomic sequence polymorphisms, several haplotypes belonging to two major lineages of the human pathogen Mycobacterium ulcerans could be distinguished among patient isolates from various geographic origins. However, the biological relevance of insertional/deletional diversity is not understood.

Methodology

Using comparative genomics, we have investigated the genes located in regions of difference recently identified by DNA microarray based hybridisation analysis. The analysed regions of difference comprise ∼7% of the entire M. ulcerans genome.

Principal Findings

Several different mechanisms leading to loss of functional genes were identified, ranging from pseudogenization, caused by frame shift mutations or mobile genetic element interspersing, to large sequence polymorphisms. Four hot spot regions for genetic instability were unveiled. Altogether, 229 coding sequences were found to be differentially inactivated, constituting a repertoire of coding sequence variation in the rather monomorphic M. ulcerans.

Conclusions/Significance

The differential gene inactivation patterns associated with the M. ulcerans haplotypes identified candidate genes that may confer enhanced adaptation upon ablation of expression. A number of gene conversions confined to the classical lineage may contribute to particular virulence of this group comprising isolates from Africa and Australia. Identification of this spectrum of anti-virulence gene candidates expands our understanding of the pathogenicity and ecology of the emerging infectious disease Buruli ulcer.     

Cross-platform microarray data normalisation for regulatory network inference

Background

Inferring Gene Regulatory Networks (GRNs) from time course microarray data suffers from the dimensionality problem created by the short length of available time series compared to the large number of genes in the network. To overcome this, data integration from diverse sources is mandatory. Microarray data from different sources and platforms are publicly available, but integration is not straightforward, due to platform and experimental differences.

Methods

We analyse here different normalisation approaches for microarray data integration, in the context of reverse engineering of GRN quantitative models. We introduce two preprocessing approaches based on existing normalisation techniques and provide a comprehensive comparison of normalised datasets.

Conclusions

Results identify a method based on a combination of Loess normalisation and iterative K-means as best for time series normalisation for this problem.