Meta-Analysis of MicroRNA-146a rs2910164 G>C Polymorphism Association with Autoimmune Diseases Susceptibility,an Update Based on 24 Studies

Background

Published data showed that the susceptibility of autoimmune diseases (ADs) was associated with the polymorphism rs2910164 in microRNA-146a (miR-146a). However, the results remain controversial so far. Two meta-analyses published in 2013 and 2014 came to opposite conclusions. In order to derive a more precise estimation of the relationship, we performed this meta-analysis.

Methods

We searched the PubMed, OvidSP and CNKI databases (published prior to September 8th, 2014) and extracted data from eligible studies. The procedure of meta-analysis was performed by using the Stata 12.0 software. Random effect model or fixed effect model were chosen respectively, according to the between study heterogeneities.

Results

A total of 24 case-control studies, 11 more than previous meta-analysis on this topic, were involved. We took stratified analyses by different ethnicities and different types of diseases in different genetic models. In Caucasian subgroup, significant increased risks of GC genotype and GC+CC genotype with ADs susceptibility were found in heterozygote model (GC vs GG, OR = 1.38, 95% CI 1.04–1.83, p = 0.024) and dominant model (GC+CC vs GG, OR = 1.37, 95% CI 1.01–1.85, p = 0.041), respectively. Meanwhile, in other disease subgroup, significant increased risks of C allele, CC genotype and GC+CC genotype were found in allele model (C vs G, OR = 1.16, 95% CI 1.04–1.31, p = 0.010), homozygote model (CC vs GG, OR = 1.42, 95% CI 1.10–1.84, p = 0.006) and dominant model (GC+CC vs GG, OR = 1.25, 95% CI 1.04–1.51, p = 0.020), respectively.

Conclusions

MiR-146a rs2910164 G>C polymorphism was associated with the susceptibility of ADs.     

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.     

No Specific Gene Expression Signature in Human Granulosa and Cumulus Cells for Prediction of Oocyte Fertilisation and Embryo Implantation

In human IVF procedures objective and reliable biomarkers of oocyte and embryo quality are needed in order to increase the use of single embryo transfer (SET) and thus prevent multiple pregnancies. During folliculogenesis there is an intense bi-directional communication between oocyte and follicular cells. For this reason gene expression profile of follicular cells could be an important indicator and biomarker of oocyte and embryo quality. The objective of this study was to identify gene expression signature(s) in human granulosa (GC) and cumulus (CC) cells predictive of successful embryo implantation and oocyte fertilization. Forty-one patients were included in the study and individual GC and CC samples were collected; oocytes were cultivated separately, allowing a correlation with IVF outcome and elective SET was performed. Gene expression analysis was performed using microarrays, followed by a quantitative real-time PCR validation. After statistical analysis of microarray data, there were no significantly differentially expressed genes (FDR<0,05) between non-fertilized and fertilized oocytes and non-implanted and implanted embryos in either of the cell type. Furthermore, the results of quantitative real-time PCR were in consent with microarray data as there were no significant differences in gene expression of genes selected for validation. In conclusion, we did not find biomarkers for prediction of oocyte fertilization and embryo implantation in IVF procedures in the present study.     

The rs738409 (I148M) variant of the PNPLA3 gene and cirrhosis: a meta-analysis

The human patatin-like phospholipase domain-containing-3 (PNPLA3) gene rs738409 C>G polymorphism is associated with several types of liver disease. The aim of this meta-analysis was to assess the risk of cirrhosis on the basis of rs738409 allele frequency and genotype. Medline, the Cochrane Library, EMBASE, and Google Scholar were searched for prospective and retrospective studies assessing the effect of the rs738409 polymorphism on liver cirrhosis. Seven studies, involving 2,023 patients with cirrhosis, were included. The G allele was associated with a significantly increased risk of cirrhosis versus the C allele [pooled odds ratio (OR) = 1.86, 95% confidence interval (CI) = 1.64–2.12, Z = 9.55, P < 0.001]. Both the GC and GG genotypes were associated with a significantly increased risk of cirrhosis versus the CC genotype (GC vs. CC: pooled OR = 1.73, 95% CI = 1.51–1.98, Z = 7.86, P < 0.001; GG vs. CC: pooled OR = 3.41, 95% CI = 2.77–4.18, Z = 11.65, P < 0.001). There was no evidence of publication bias. Our findings suggest that patients at risk for liver cirrhosis may benefit from PNPLA3 genotyping and thus more intensive monitoring if the rs738409 C>G polymorphism is identified.     

Increased Expression of Pregnancy Up-Regulated Non-Ubiquitous Calmodulin Kinase Is Associated with Poor Prognosis in Clear Cell Renal Cell Carcinoma

Purpose

The aims of this study were to evaluate the clinical significance and potential prognostic value of pregnancy up-regulated non-ubiquitous calmodulin kinase (PNCK) in clear cell renal cell carcinoma (ccRCC) patients.

Materials and Methods

The expression of PNCK mRNA was determined in 24 paired samples of ccRCCs and adjacent normal tissues using real-time RT-PCR. The expression of PNCK was determined in 248 samples of ccRCCs and 92 paired samples of adjacent normal tissues by immunohistochemical analysis. Statistical analysis was performed to define the relationship between PNCK expression and the clinical features of ccRCC.

Results

The mRNA level of PNCK was significantly higher in tumorous tissues than in the adjacent non-tumorous tissues (p<0.001). An immunohistochemical analysis of 92 paired tissue specimens showed that PNCK expression was higher in tumorous tissues than in the adjacent non-tumorous tissues (p<0.001). Moreover, there was a significant correlation between the PNCK expression and various clinicopathological parameters such as Fuhrman grade (p = 0.011), tumor size (p<0.001), T stage (p<0.001) and N stage (p = 0.015). Patients with higher PNCK expression had shorter overall survival time than those with lower PNCK expression (p<0.001). Multivariate analysis indicated that PNCK expression was an independent predictor for poor survival of ccRCC patients.

Conclusions

To our knowledge, this is the first study that determines the relationship between PNCK and prognosis in ccRCC. We found that increased PNCK expression is associated with poor prognosis in ccRCC. PNCK may represent a novel prognostic marker for ccRCC.     

Expression of PROKR1 and PROKR2 in human enteric neural precursor cells and identification of sequence variants suggest a role in HSCR

Background

The enteric nervous system (ENS) is entirely derived from neural crest and its normal development is regulated by specific molecular pathways. Failure in complete ENS formation results in aganglionic gut conditions such as Hirschsprung''s disease (HSCR). Recently, PROKR1 expression has been demonstrated in mouse enteric neural crest derived cells and Prok-1 was shown to work coordinately with GDNF in the development of the ENS.

Principal Findings

In the present report, ENS progenitors were isolated and characterized from the ganglionic gut from children diagnosed with and without HSCR, and the expression of prokineticin receptors was examined. Immunocytochemical analysis of neurosphere-forming cells demonstrated that both PROKR1 and PROKR2 were present in human enteric neural crest cells. In addition, we also performed a mutational analysis of PROKR1, PROKR2, PROK1 and PROK2 genes in a cohort of HSCR patients, evaluating them for the first time as susceptibility genes for the disease. Several missense variants were detected, most of them affecting highly conserved amino acid residues of the protein and located in functional domains of both receptors, which suggests a possible deleterious effect in their biological function.

Conclusions

Our results suggest that not only PROKR1, but also PROKR2 might mediate a complementary signalling to the RET/GFRα1/GDNF pathway supporting proliferation/survival and differentiation of precursor cells during ENS development. These findings, together with the detection of sequence variants in PROKR1, PROK1 and PROKR2 genes associated to HSCR and, in some cases in combination with RET or GDNF mutations, provide the first evidence to consider them as susceptibility genes for HSCR.     

Five genes encoding surface-exposed LPXTG proteins are enriched in hospital-adapted Enterococcus faecium clonal complex 17 isolates

Most Enterococcus faecium isolates associated with hospital outbreaks and invasive infections belong to a distinct genetic subpopulation called clonal complex 17 (CC17). It has been postulated that the genetic evolution of CC17 involves the acquisition of various genes involved in antibiotic resistance, metabolic pathways, and virulence. To gain insight into additional genes that may have favored the rapid emergence of this nosocomial pathogen, we aimed to identify surface-exposed LPXTG cell wall-anchored proteins (CWAPs) specifically enriched in CC17 E. faecium. Using PCR and Southern and dot blot hybridizations, 131 E. faecium isolates (40 CC17 and 91 non-CC17) were screened for the presence of 22 putative CWAP genes identified from the E. faecium TX0016 genome. Five genes encoding LPXTG surface proteins were specifically enriched in E. faecium CC17 isolates. These five LPXTG surface protein genes were found in 28 to 40 (70 to 100%) of CC17 and in only 7 to 24 (8 to 26%) of non-CC17 isolates (P < 0.05). Three of these CWAP genes clustered together on the E. faecium TX0016 genome, which may comprise a novel enterococcal pathogenicity island covering E. faecium contig 609. Expression at the mRNA level was demonstrated, and immunotransmission electron microscopy revealed an association of the five LPXTG surface proteins with the cell wall. Minimal spanning tree analysis based on the presence and absence of 22 CWAP genes revealed grouping of all 40 CC17 strains together with 18 hospital-derived but evolutionary unrelated non-CC17 isolates in a distinct CWAP-enriched cluster, suggesting horizontal transfer of CWAP genes and a role of these CWAPs in hospital adaptation.     

A Transcriptome-Wide Screen for mRNAs Enriched in Fetal Leydig Cells: CRHR1 Agonism Stimulates Rat and Mouse Fetal Testis Steroidogenesis

Fetal testis steroidogenesis plays an important role in the reproductive development of the male fetus. While regulators of certain aspects of steroidogenesis are known, the initial driver of steroidogenesis in the human and rodent fetal testis is unclear. Through comparative analysis of rodent fetal testis microarray datasets, 54 candidate fetal Leydig cell-specific genes were identified. Fetal mouse testis interstitial expression of a subset of these genes with unknown expression (Crhr1, Gramd1b, Itih5, Vgll3, and Vsnl1) was verified by whole-mount in situ hybridization. Among the candidate fetal Leydig cell-specific factors, three receptors (CRHR1, PRLR, and PROKR2) were tested for a steroidogenic function using ex vivo fetal testes treated with receptor agonists (CRH, PRL, and PROK2). While PRL and PROK2 had no effect, CRH, at low (approximately 1 to 10) nM concentration, increased expression of the steroidogenic genes Cyp11a1, Cyp17a1, Scarb1, and Star in GD15 mouse and GD17 rat testes, and in conjunction, testosterone production was increased. Exposure of GD15 fetal mouse testis to a specific CRHR1 antagonist blunted the CRH-induced steroidogenic gene expression and testosterone responses. Similar to ex vivo rodent fetal testes, ≥10 nM CRH exposure of MA-10 Leydig cells increased steroidogenic pathway mRNA and progesterone levels, showing CRH can enhance steroidogenesis by directly targeting Leydig cells. Crh mRNA expression was observed in rodent fetal hypothalamus, and CRH peptide was detected in rodent amniotic fluid. Together, these data provide a resource for discovering factors controlling fetal Leydig cell biology and suggest that CRHR1 activation by CRH stimulates rat and mouse fetal Leydig cell steroidogenesis in vivo.     

Identification of Highly Methylated Genes across Various Types of B-Cell Non-Hodgkin Lymphoma

Epigenetic alterations of gene expression are important in the development of cancer. In this study, we identified genes which are epigenetically altered in major lymphoma types. We used DNA microarray technology to assess changes in gene expression after treatment of 11 lymphoma cell lines with epigenetic drugs. We identified 233 genes with upregulated expression in treated cell lines and with downregulated expression in B-cell lymphoma patient samples (n = 480) when compared to normal B cells (n = 5). The top 30 genes were further analyzed by methylation specific PCR (MSP) in 18 lymphoma cell lines. Seven of the genes were methylated in more than 70% of the cell lines and were further subjected to quantitative MSP in 37 B-cell lymphoma patient samples (diffuse large B-cell lymphoma (activated B-cell like and germinal center B-cell like subtypes), follicular lymphoma and Burkitt`s lymphoma) and normal B lymphocytes from 10 healthy donors. The promoters of DSP, FZD8, KCNH2, and PPP1R14A were methylated in 28%, 67%, 22%, and 78% of the 36 tumor samples, respectively, but not in control samples. Validation using a second series of healthy donor controls (n = 42; normal B cells, peripheral blood mononuclear cells, bone marrow, tonsils and follicular hyperplasia) and fresh-frozen lymphoma biopsies (n = 25), confirmed the results. The DNA methylation biomarker panel consisting of DSP, FZD8, KCNH2, and PPP1R14A was positive in 89% (54/61) of all lymphomas. Receiver operating characteristic analysis to determine the discriminative power between lymphoma and healthy control samples showed a c-statistic of 0.96, indicating a possible role for the biomarker panel in monitoring of lymphoma patients.     

Poor Prognosis of Phosphatase of Regenerating Liver 3 Expression in Gastric Cancer: A Meta-Analysis

Background

Overexpression of phosphatase of regenerating liver 3 (PRL-3) has been implicated in gastric cancer (GC) metastasis. Epidemiological studies have evaluated the relationship between PRL-3 expression and prognosis in GC. However, results still remains controversial. In this study, a meta-analysis was performed to evaluate the association of PRL-3 expression with overall survival (OS) and clinicopathological characteristics.

Methods

Literature databases were searched to identify eligible studies dated until April 2013. Summary hazard ratios (HRs) or odds ratios (ORs) with 95% confidence interval (95% CI) were calculated to estimate the association.

Results

A total of 1380 GC patients from six studies were included in the meta-analysis. Overall, the combined HR estimate for OS in a random-effect model was 1.89 (95% CI = 1.38–2.60; P<0.001). Results showed that PRL-3 overexpression was significantly associated with OS, indicating that it may be a biomarker for poor prognosis of GC. Both subgroup and sensitivity analyses further identified the prognostic role of PRL-3 expression in GC patients. Moreover, PRL-3 overexpression was significantly associated with tumor stage (OR = 2.25; 95% CI = 1.63–3.12; P<0.001), depth of invasion (OR = 2.03; 95% CI = 1.38–2.98; P<0.001), vascular invasion (OR = 2.52; 95% CI = 1.79–3.56; P<0.001), lymphatic invasion (OR = 3.74; 95% CI = 2.49–5.63; P<0.001), and lymph node metastasis (OR = 4.56; 95% CI = 2.37–8.76; P<0.001). However, when age, sex, tumor size, and tumor differentiation were considered, no obvious association was observed.

Conclusions

This meta-analysis reveals significant association of PRL-3 overexpression with OS and some clinicopathological features in GC. PRL-3 may be a predicative factor of poor prognosis and aggressive tumor behavior in GC patients.     

Gene Polymorphisms of Micrornas in Helicobacter pylori-Induced High Risk Atrophic Gastritis and Gastric Cancer

Background and aims

MicroRNAs (miRNAs) are known for their function as translational regulators of tumor suppressor or oncogenes. Single nucleotide polymorphisms (SNPs) in miRNAs related genes have been shown to affect the regulatory capacity of miRNAs and were linked with gastric cancer (GC) and premalignant gastric conditions. The purpose of this study was to evaluate potential associations between miRNA-related gene polymorphisms (miR-27a, miR-146a, miR-196a-2, miR-492 and miR-608) and the presence of GC or high risk atrophic gastritis (HRAG) in European population.

Methods

Gene polymorphisms were analyzed in 995 subjects (controls: n = 351; GC: n = 363; HRAG: n = 281) of European descent. MiR-27a T>C (rs895819), miR-146a G>C (rs2910164), miR-196a-2 C>T (rs11614913), miR-492 G>C (rs2289030) and miR-608 C>G (rs4919510) SNPs were genotyped by RT-PCR.

Results

Overall, SNPs of miRNAs were not associated with the presence of GC or HRAG. We observed a tendency for miR-196a-2 CT genotype to be associated with higher risk of GC when compared to CC genotype, however, the difference did not reach the adjusted P-value (odds ratio (OR) - 1.46, 95% confidence interval (CI) 1.03-2.07, P = 0.032). MiR-608 GG genotype was more frequent in GC when compared to controls (OR −2.34, 95% CI 1.08–5.04), but significance remained marginal (P = 0.029). A similar tendency was observed in a recessive model for miR-608, where CC + CG vs GG genotype comparison showed a tendency for increased risk of GC with OR of 2.44 (95% CI 1.14–5.22, P = 0.021). The genotypes and alleles of miR-27a, miR-146a, miR-196a-2, miR-492 and miR-608 SNPs had similar distribution between histological subtypes of GC and were not linked with the presence of diffuse or intestinal-type GC.

Conclusions

Gene polymorphisms of miR-27a, miR-146a, miR-196a-2, miR-492, miR-492a and miR-608 were not associated with the presence of HRAG, GC or different histological subtypes of GC in European subjects.     

Restoration of Dlk1 and Rtl1 Is Necessary but Insufficient to Rescue Lethality in Intergenic Differentially Methylated Region (IG-DMR)-deficient Mice

In the Dlk1-Dio3 imprinted domain, an intergenic differentially methylated region (IG-DMR) regulates the parental allele-specific expression of imprinted genes. The maternally inherited deletion of IG-DMR (IG-DMR^(−/+)) results in perinatal lethality because of the overexpression of paternally expressed genes and repression of maternally expressed noncoding RNAs (ncRNAs), including Gtl2. To better understand the possible contribution of paternally expressed genes to the lethality, we attempted to rescue the lethality of IG-DMR^(−/+) mutants by restoring the paternally expressed genes. Because the paternally inherited Gtl2 deletion (Gtl2^(+/−)) induced a decrease in the expression of paternally expressed genes, we crossed female IG-DMR heterozygous mice and male Gtl2 heterozygous mutant mice. The resultant IG-DMR^(−/+)/Gtl2^(+/−) double mutant mice had normal expression levels of paternally expressed genes, and none of them showed perinatal lethality; however, most mice showed postnatal lethality with decreased expression of the maternally expressed ncRNAs. Thus, we inferred that paternally expressed genes are necessary for perinatal survivability and that maternally expressed ncRNAs are involved in postnatal lethality.     

Proteome Analysis of Human Sebaceous Follicle Infundibula Extracted from Healthy and Acne-Affected Skin

Acne vulgaris is a very common disease of the pilosebaceous unit of the human skin. The pathological processes of acne are not fully understood. To gain further insight sebaceous follicular casts were extracted from 18 healthy and 20 acne-affected individuals by cyanoacrylate-gel biopsies and further processed for mass spectrometry analysis, aiming at a proteomic analysis of the sebaceous follicular casts. Human as well as bacterial proteins were identified. Human proteins enriched in acne and normal samples were detected, respectively. Normal follicular casts are enriched in proteins such as prohibitins and peroxiredoxins which are involved in the protection from various stresses, including reactive oxygen species. By contrast, follicular casts extracted from acne-affected skin contained proteins involved in inflammation, wound healing and tissue remodeling. Among the most distinguishing proteins were myeloperoxidase, lactotransferrin, neutrophil elastase inhibitor and surprisingly, vimentin. The most significant biological process among all acne-enriched proteins was ‘response to a bacterium’. Identified bacterial proteins were exclusively from Propionibacterium acnes. The most abundant P. acnes proteins were surface-exposed dermatan sulphate adhesins, CAMP factors, and a so far uncharacterized lipase in follicular casts extracted from normal as well as acne-affected skin. This is a first proteomic study that identified human proteins together with proteins of the skin microbiota in sebaceous follicular casts.