Serum Hydroxysteroid (17beta) dehydrogenase 1 concentration in pregnant women correlates with pregnancy-associated plasma protein a but does not serve as an independent marker for preeclampsia

Hydroxysteroid (17beta) dehydrogenase 1 (HSD17B1) is a steroid synthetic enzyme expressed in ovarian granulosa cells and placental syncytiotrophoblasts. Here, HSD17B1 serum concentration was measured with a validated immuno assay during pregnancy at three time points (12-14, 18-20 and 26-28 weeks of gestation). The concentration increased 2.5-fold (p < 0.0001) and 1.7-fold (p = 0.0019) during the follow-up period for control women and women who later developed preeclampsia (PE), respectively, and a significant difference was observed at weeks 26-28 (p = 0.0266). HSD17B1 concentration at all the three time points positively correlated with serum PAPPA measured at the first time point (first time point r = 0.38, p = 1.1x10-10; second time point r = 0.27, p = 5.9x10-6 and third timepoint r = 0.26, p = 2.3x10-5). No correlation was observed between HSD17B1 and placental growth factor (PLGF). Serum HSD17B1, furthermore, negatively correlated with the mother's weight and body mass index (BMI), mirroring the pattern observed for PAPPA. The univariable logistic regression identified a weak association between HSD17B1 at 26-28 weeks and later development of PE (P = 0.04). Also, the best multivariable model obtained using penalized logistic regression with stable iterative variable selection at 26-28 weeks included HSD17B1, together with PLGF, PAPPA and the mother's BMI. While the area under the ROC curve of the model was higher than that of the adjusted PLGF, the difference was not statistically significant. In summary, the serum concentration of HSD17B1 correlated with PAPPA, another protein expressed in syncytiotrophoblasts, and with mother's weight and BMI but could not be considered as an independent marker for PE.

A proximal enhancer regulates RORA expression during early human Th17 cell differentiation

BACKGROUND

Gene regulatory elements, such as enhancers, greatly influence cell identity by tuning the transcriptional activity of specific cell types. Dynamics of enhancer landscape during early human Th17 cell differentiation remains incompletely understood. Leveraging ATAC-seq-based profiling of chromatin accessibility and comprehensive analysis of key histone marks, we identified a repertoire of potential enhancers that potentially exert control over the fate specification of Th17 cells. We found 23 SNPs associated with autoimmune diseases within Th17-enhancers that precisely overlapped with the binding sites of transcription factors actively engaged in T-cell functions. Among the Th17-specific enhancers, we identified an enhancer in the intron of RORA and demonstrated that this enhancer positively regulates RORA transcription. Moreover, CRISPR-Cas9-mediated deletion of a transcription factor binding site-rich region within the identified RORA enhancer confirmed its role in regulating RORA transcription. These findings provide profound insights into the potential mechanism by which the RORA enhancer orchestrates Th17 differentiation.

RUNX transcription factors are essential in maintaining epididymal epithelial differentiation

Apart from the androgen receptor, transcription factors (TFs) that are required for the development and formation of the different segments of the epididymis have remained unknown. We identified TF families expressed in the developing epididymides, of which many showed segment specificity. From these TFs, down-regulation of runt related transcription factors (RUNXs) 1 and 2 expression coincides with epithelial regression in Dicer1 cKO mice. Concomitant deletion of both Runx1 and Runx2 in a mouse epididymal epithelial cell line affected cell morphology, adhesion and mobility in vitro. Furthermore, lack of functional RUNXs severely disturbed the formation of 3D epididymal organoid-like structures. Transcriptomic analysis of the epididymal cell organoid-like structures indicated that RUNX1 and RUNX2 are involved in the regulation of MAPK signaling, NOTCH pathway activity, and EMT-related gene expression. This suggests that RUNXs are master regulators of several essential signaling pathways, and necessary for the maintenance of proper differentiation of the epididymal epithelium.

PIM kinases regulate early human Th17 cell differentiation

The serine/threonine-specific Moloney murine leukemia virus (PIM) kinase family (i.e., PIM1, PIM2, and PIM3) has been extensively studied in tumorigenesis. PIM kinases are downstream of several cytokine signaling pathways that drive immune-mediated diseases. Uncontrolled T helper 17 (Th17) cell activation has been associated with the pathogenesis of autoimmunity. However, the detailed molecular function of PIMs in human Th17 cell regulation has yet to be studied. In the present study, we comprehensively investigated how the three PIMs simultaneously alter transcriptional gene regulation during early human Th17 cell differentiation. By combining PIM triple knockdown with bulk and scRNA-seq approaches, we found that PIM deficiency promotes the early expression of key Th17-related genes while suppressing Th1-lineage genes. Further, PIMs modulate Th cell signaling, potentially via STAT1 and STAT3. Overall, our study highlights the inhibitory role of PIMs in human Th17 cell differentiation, thereby suggesting their association with autoimmune phenotypes.

Profiling of peripheral blood B-cell transcriptome in children who developed coeliac disease in a prospective study

Background

In coeliac disease (CoD), the role of B-cells has mainly been considered to be production of antibodies. The functional role of B-cells has not been analysed extensively in CoD.

Methods

We conducted a study to characterize gene expression in B-cells from children developing CoD early in life using samples collected before and at the diagnosis of the disease. Blood samples were collected from children at risk at 12, 18, 24 and 36 months of age. RNA from peripheral blood CD19⁺ cells was sequenced and differential gene expression was analysed using R package Limma.

Findings

Overall, we found one gene, HNRNPL, modestly downregulated in all patients (logFC -0·7; q = 0·09), and several others downregulated in those diagnosed with CoD already by the age of 2 years.

Interpretation

The data highlight the role of B-cells in CoD development. The role of HNRPL in suppressing enteroviral replication suggests that the predisposing factor for both CoD and enteroviral infections is the low level of HNRNPL expression.

Funding

EU FP7 grant no. 202063, EU Regional Developmental Fund and research grant PRG712, The Academy of Finland Centre of Excellence in Molecular Systems Immunology and Physiology Research (SyMMyS) 2012-2017, grant no. 250114) and, AoF Personalized Medicine Program (grant no. 292482), AoF grants 292335, 294337, 319280, 31444, 319280, 329277, 331790) and grants from the Sigrid Jusélius Foundation (SJF).

SMG6 localizes to the chromatoid body and shapes the male germ cell transcriptome to drive spermatogenesis

Nonsense-mediated RNA decay (NMD) is a highly conserved and selective RNA turnover pathway that depends on the endonuclease SMG6. Here, we show that SMG6 is essential for male germ cell differentiation in mice. Germ-cell conditional knockout (cKO) of Smg6 induces extensive transcriptome misregulation, including a failure to eliminate meiotically expressed transcripts in early haploid cells, and accumulation of NMD target mRNAs with long 3' untranslated regions (UTRs). Loss of SMG6 in the male germline results in complete arrest of spermatogenesis at the early haploid cell stage. We find that SMG6 is strikingly enriched in the chromatoid body (CB), a specialized cytoplasmic granule in male germ cells also harboring PIWI-interacting RNAs (piRNAs) and the piRNA-binding protein PIWIL1. This raises the possibility that SMG6 and the piRNA pathway function together, which is supported by several findings, including that Piwil1-KO mice phenocopy Smg6-cKO mice and that SMG6 and PIWIL1 co-regulate many genes in round spermatids. Together, our results demonstrate that SMG6 is an essential regulator of the male germline transcriptome, and highlight the CB as a molecular platform coordinating RNA regulatory pathways to control sperm production and fertility.

Umbilical cord blood DNA methylation in children who later develop type 1 diabetes

AIMS/HYPOTHESIS

Distinct DNA methylation patterns have recently been observed to precede type 1 diabetes in whole blood collected from young children. Our aim was to determine whether perinatal DNA methylation is associated with later progression to type 1 diabetes.

METHODS

Reduced representation bisulphite sequencing (RRBS) analysis was performed on umbilical cord blood samples collected within the Finnish Type 1 Diabetes Prediction and Prevention (DIPP) Study. Children later diagnosed with type 1 diabetes and/or who tested positive for multiple islet autoantibodies (n = 43) were compared with control individuals (n = 79) who remained autoantibody-negative throughout the DIPP follow-up until 15 years of age. Potential confounding factors related to the pregnancy and the mother were included in the analysis.

RESULTS

No differences in the umbilical cord blood methylation patterns were observed between the cases and controls at a false discovery rate <0.05.

CONCLUSIONS/INTERPRETATION

Based on our results, differences between children who progress to type 1 diabetes and those who remain healthy throughout childhood are not yet present in the perinatal DNA methylome. However, we cannot exclude the possibility that such differences would be found in a larger dataset.

High intratumoral dihydrotestosterone is associated with antiandrogen resistance in VCaP prostate cancer xenografts in castrated mice

Antiandrogen treatment resistance is a major clinical concern in castration-resistant prostate cancer (CRPC) treatment. Using xenografts of VCaP cells we showed that growth of antiandrogen resistant CRPC tumors were characterized by a higher intratumor dihydrotestosterone (DHT) concentration than that of treatment responsive tumors. Furthermore, the slow tumor growth after adrenalectomy was associated with a low intratumor DHT concentration. Reactivation of androgen signaling in enzalutamide-resistant tumors was further shown by the expression of several androgen-dependent genes. The data indicate that intratumor DHT concentration and expression of several androgen-dependent genes in CRPC lesions is an indication of enzalutamide treatment resistance and an indication of the need for further androgen blockade. The presence of an androgen synthesis, independent of CYP17A1 activity, has been shown to exist in prostate cancer cells, and thus, novel androgen synthesis inhibitors are needed for the treatment of enzalutamide-resistant CRPC tumors that do not respond to abiraterone.

A systematic comparison of FOSL1, FOSL2 and BATF-mediated transcriptional regulation during early human Th17 differentiation

Th17 cells are essential for protection against extracellular pathogens, but their aberrant activity can cause autoimmunity. Molecular mechanisms that dictate Th17 cell-differentiation have been extensively studied using mouse models. However, species-specific differences underscore the need to validate these findings in human. Here, we characterized the human-specific roles of three AP-1 transcription factors, FOSL1, FOSL2 and BATF, during early stages of Th17 differentiation. Our results demonstrate that FOSL1 and FOSL2 co-repress Th17 fate-specification, whereas BATF promotes the Th17 lineage. Strikingly, FOSL1 was found to play different roles in human and mouse. Genome-wide binding analysis indicated that FOSL1, FOSL2 and BATF share occupancy over regulatory regions of genes involved in Th17 lineage commitment. These AP-1 factors also share their protein interacting partners, which suggests mechanisms for their functional interplay. Our study further reveals that the genomic binding sites of FOSL1, FOSL2 and BATF harbour hundreds of autoimmune disease-linked SNPs. We show that many of these SNPs alter the ability of these transcription factors to bind DNA. Our findings thus provide critical insights into AP-1-mediated regulation of human Th17-fate and associated pathologies.

Early DNA methylation changes in children developing beta cell autoimmunity at a young age

AIMS/HYPOTHESIS

Type 1 diabetes is a chronic autoimmune disease of complex aetiology, including a potential role for epigenetic regulation. Previous epigenomic studies focused mainly on clinically diagnosed individuals. The aim of the study was to assess early DNA methylation changes associated with type 1 diabetes already before the diagnosis or even before the appearance of autoantibodies.

METHODS

Reduced representation bisulphite sequencing (RRBS) was applied to study DNA methylation in purified CD4⁺ T cell, CD8⁺ T cell and CD4^-CD8^- cell fractions of 226 peripheral blood mononuclear cell samples longitudinally collected from seven type 1 diabetes-specific autoantibody-positive individuals and control individuals matched for age, sex, HLA risk and place of birth. We also explored correlations between DNA methylation and gene expression using RNA sequencing data from the same samples. Technical validation of RRBS results was performed using pyrosequencing.

RESULTS

We identified 79, 56 and 45 differentially methylated regions in CD4⁺ T cells, CD8⁺ T cells and CD4^-CD8^- cell fractions, respectively, between type 1 diabetes-specific autoantibody-positive individuals and control participants. The analysis of pre-seroconversion samples identified DNA methylation signatures at the very early stage of disease, including differential methylation at the promoter of IRF5 in CD4⁺ T cells. Further, we validated RRBS results using pyrosequencing at the following CpG sites: chr19:18118304 in the promoter of ARRDC2; chr21:47307815 in the intron of PCBP3; and chr14:81128398 in the intergenic region near TRAF3 in CD4⁺ T cells.

CONCLUSIONS/INTERPRETATION

These preliminary results provide novel insights into cell type-specific differential epigenetic regulation of genes, which may contribute to type 1 diabetes pathogenesis at the very early stage of disease development. Should these findings be validated, they may serve as a potential signature useful for disease prediction and management.

Differential ATAC-seq and ChIP-seq peak detection using ROTS

Changes in cellular chromatin states fine-tune transcriptional output and ultimately lead to phenotypic changes. Here we propose a novel application of our reproducibility-optimized test statistics (ROTS) to detect differential chromatin states (ATAC-seq) or differential chromatin modification states (ChIP-seq) between conditions. We compare the performance of ROTS to existing and widely used methods for ATAC-seq and ChIP-seq data using both synthetic and real datasets. Our results show that ROTS outperformed other commonly used methods when analyzing ATAC-seq data. ROTS also displayed the most accurate detection of small differences when modeling with synthetic data. We observed that two-step methods that require the use of a separate peak caller often more accurately called enrichment borders, whereas one-step methods without a separate peak calling step were more versatile in calling sub-peaks. The top ranked differential regions detected by the methods had marked correlation with transcriptional differences of the closest genes. Overall, our study provides evidence that ROTS is a useful addition to the available differential peak detection methods to study chromatin and performs especially well when applied to study differential chromatin states in ATAC-seq data.

Evaluation of tools for identifying large copy number variations from ultra-low-coverage whole-genome sequencing data

BACKGROUND

Detection of copy number variations (CNVs) from high-throughput next-generation whole-genome sequencing (WGS) data has become a widely used research method during the recent years. However, only a little is known about the applicability of the developed algorithms to ultra-low-coverage (0.0005-0.8×) data that is used in various research and clinical applications, such as digital karyotyping and single-cell CNV detection.

RESULT

Here, the performance of six popular read-depth based CNV detection algorithms (BIC-seq2, Canvas, CNVnator, FREEC, HMMcopy, and QDNAseq) was studied using ultra-low-coverage WGS data. Real-world array- and karyotyping kit-based validation were used as a benchmark in the evaluation. Additionally, ultra-low-coverage WGS data was simulated to investigate the ability of the algorithms to identify CNVs in the sex chromosomes and the theoretical minimum coverage at which these tools can accurately function. Our results suggest that while all the methods were able to detect large CNVs, many methods were susceptible to producing false positives when smaller CNVs (< 2 Mbp) were detected. There was also significant variability in their ability to identify CNVs in the sex chromosomes. Overall, BIC-seq2 was found to be the best method in terms of statistical performance. However, its significant drawback was by far the slowest runtime among the methods (> 3 h) compared with FREEC (~ 3 min), which we considered the second-best method.

CONCLUSIONS

Our comparative analysis demonstrates that CNV detection from ultra-low-coverage WGS data can be a highly accurate method for the detection of large copy number variations when their length is in millions of base pairs. These findings facilitate applications that utilize ultra-low-coverage CNV detection.

Complex Interplay Between MAZR and Runx3 Regulates the Generation of Cytotoxic T Lymphocyte and Memory T Cells

The BTB zinc finger transcription factor MAZR (also known as PATZ1) controls, partially in synergy with the transcription factor Runx3, the development of CD8 lineage T cells. Here we explored the role of MAZR as well as combined activities of MAZR/Runx3 during cytotoxic T lymphocyte (CTL) and memory CD8+ T cell differentiation. In contrast to the essential role of Runx3 for CTL effector function, the deletion of MAZR had a mild effect on the generation of CTLs in vitro. However, a transcriptome analysis demonstrated that the combined deletion of MAZR and Runx3 resulted in much more widespread downregulation of CTL signature genes compared to single Runx3 deletion, indicating that MAZR partially compensates for loss of Runx3 in CTLs. Moreover, in line with the findings made in vitro, the analysis of CTL responses to LCMV infection revealed that MAZR and Runx3 cooperatively regulate the expression of CD8α, Granzyme B and perforin in vivo. Interestingly, while memory T cell differentiation is severely impaired in Runx3-deficient mice, the deletion of MAZR leads to an enlargement of the long-lived memory subset and also partially restored the differentiation defect caused by loss of Runx3. This indicates distinct functions of MAZR and Runx3 in the generation of memory T cell subsets, which is in contrast to their cooperative roles in CTLs. Together, our study demonstrates complex interplay between MAZR and Runx3 during CTL and memory T cell differentiation, and provides further insight into the molecular mechanisms underlying the establishment of CTL and memory T cell pools.

Exon-level estimates improve the detection of differentially expressed genes in RNA-seq studies

Detection of differentially expressed genes (DEGs) between different biological conditions is a key data analysis step of most RNA-sequencing studies. Conventionally, computational tools have used gene-level read counts as input to test for differential gene expression between sample condition groups. Recently, it has been suggested that statistical testing could be performed with increased power at a lower feature level prior to aggregating the results to the gene level. In this study, we systematically compared the performance of calling the DEGs when using read count data at different levels (gene, transcript, and exon) as input, in the context of two publicly available data sets. Additionally, we tested two different methods for aggregating the lower feature-level p-values to gene-level: Lancaster and empirical Brown’s method. Our results show that detection of DEGs is improved compared to the conventional gene-level approach regardless of the lower feature-level used for statistical testing. The overall best balance between accuracy and false discovery rate was obtained using the exon-level approach with empirical Brown’s aggregation method, which we provide as a freely available Bioconductor package EBSEA (https://bioconductor.org/packages/release/bioc/html/EBSEA.html).

The Transcription Factor MAZR/PATZ1 Regulates the Development of FOXP3+ Regulatory T Cells

Forkhead box protein P3⁺ (FOXP3⁺) regulatory T cells (T_reg cells) play a key role in maintaining tolerance and immune homeostasis. Here, we report that a T cell-specific deletion of the transcription factor MAZR (also known as PATZ1) leads to an increased frequency of T_reg cells, while enforced MAZR expression impairs T_reg cell differentiation. Further, MAZR expression levels are progressively downregulated during thymic T_reg cell development and during in-vitro-induced human T_reg cell differentiation, suggesting that MAZR protein levels are critical for controlling T_reg cell development. However, MAZR-deficient T_reg cells show only minor transcriptional changes ex vivo, indicating that MAZR is not essential for establishing the transcriptional program of peripheral T_reg cells. Finally, the loss of MAZR reduces the clinical score in dextran-sodium sulfate (DSS)-induced colitis, suggesting that MAZR activity in T cells controls the extent of intestinal inflammation. Together, these data indicate that MAZR is part of a T_reg cell-intrinsic transcriptional network that modulates T_reg cell development.

Reproducibility-optimized detection of differential DNA methylation

Aim: DNA methylation is a key epigenetic mechanism regulating gene expression. Identifying differentially methylated regions is integral to DNA methylation analysis and there is a need for robust tools reliably detecting regions with significant differences in their methylation status. Materials & methods: We present here a reproducibility-optimized test statistic (ROTS) for detection of differential DNA methylation from high-throughput sequencing or array-based data. Results: Using both simulated and real data, we demonstrate the ability of ROTS to identify differential methylation between sample groups. Conclusion: Compared with state-of-the-art methods, ROTS shows competitive sensitivity and specificity in detecting consistently differentially methylated regions.

Histone deacetylases 1 and 2 restrain CD4+ cytotoxic T lymphocyte differentiation

Some effector CD4+ T cell subsets display cytotoxic activity, thus breaking the functional dichotomy of CD4+ helper and CD8+ cytotoxic T lymphocytes. However, molecular mechanisms regulating CD4+ cytotoxic T lymphocyte (CD4+ CTL) differentiation are poorly understood. Here we show that levels of histone deacetylases 1 and 2 (HDAC1-HDAC2) are key determinants of CD4+ CTL differentiation. Deletions of both Hdac1 and 1 Hdac2 alleles (HDAC1cKO-HDAC2HET) in CD4+ T cells induced a T helper cytotoxic program that was controlled by IFN-γ-JAK1/2-STAT1 signaling. In vitro, activated HDAC1cKO-HDAC2HET CD4+ T cells acquired cytolytic activity and displayed enrichment of gene signatures characteristic of effector CD8+ T cells and human CD4+ CTLs. In vivo, murine cytomegalovirus-infected HDAC1cKO-HDAC2HET mice displayed a stronger induction of CD4+ CTL features compared with infected WT mice. Finally, murine and human CD4+ T cells treated with short-chain fatty acids, which are commensal-produced metabolites acting as HDAC inhibitors, upregulated CTL genes. Our data demonstrate that HDAC1-HDAC2 restrain CD4+ CTL differentiation. Thus, HDAC1-HDAC2 might be targets for the therapeutic induction of CD4+ CTLs.

Systematic evaluation of differential splicing tools for RNA-seq studies

Differential splicing (DS) is a post-transcriptional biological process with critical, wide-ranging effects on a plethora of cellular activities and disease processes. To date, a number of computational approaches have been developed to identify and quantify differentially spliced genes from RNA-seq data, but a comprehensive intercomparison and appraisal of these approaches is currently lacking. In this study, we systematically evaluated 10 DS analysis tools for consistency and reproducibility, precision, recall and false discovery rate, agreement upon reported differentially spliced genes and functional enrichment. The tools were selected to represent the three different methodological categories: exon-based (DEXSeq, edgeR, JunctionSeq, limma), isoform-based (cuffdiff2, DiffSplice) and event-based methods (dSpliceType, MAJIQ, rMATS, SUPPA). Overall, all the exon-based methods and two event-based methods (MAJIQ and rMATS) scored well on the selected measures. Of the 10 tools tested, the exon-based methods performed generally better than the isoform-based and event-based methods. However, overall, the different data analysis tools performed strikingly differently across different data sets or numbers of samples.

Peripheral blood DNA methylation differences in twin pairs discordant for Alzheimer's disease

BACKGROUND

Alzheimer's disease results from a neurodegenerative process that starts well before the diagnosis can be made. New prognostic or diagnostic markers enabling early intervention into the disease process would be highly valuable. Environmental and lifestyle factors largely modulate the disease risk and may influence the pathogenesis through epigenetic mechanisms, such as DNA methylation. As environmental and lifestyle factors may affect multiple tissues of the body, we hypothesized that the disease-associated DNA methylation signatures are detectable in the peripheral blood of discordant twin pairs.

RESULTS

Comparison of 23 disease discordant Finnish twin pairs with reduced representation bisulfite sequencing revealed peripheral blood DNA methylation differences in 11 genomic regions with at least 15.0% median methylation difference and FDR adjusted p value ≤ 0.05. Several of the affected genes are primarily associated with neuronal functions and pathologies and do not display disease-associated differences in gene expression in blood. The DNA methylation mark in ADARB2 gene was found to be differentially methylated also in the anterior hippocampus, including entorhinal cortex, of non-twin cases and controls. Targeted bisulfite pyrosequencing of the DNA methylation mark in ADARB2 gene in 62 Finnish and Swedish twin pairs revealed that, in addition to the disease status, DNA methylation of this region is influenced by gender, age, zygosity, APOE genotype, and smoking. Further analysis of 120 Swedish twin pairs indicated that this specific DNA methylation mark is not predictive for Alzheimer's disease and becomes differentially methylated after disease onset.

CONCLUSIONS

DNA methylation differences can be detected in the peripheral blood of twin pairs discordant for Alzheimer's disease. These DNA methylation signatures may have value as disease markers and provide insights into the molecular mechanisms of pathogenesis. We found no evidence that the DNA methylation marks would be associated with gene expression in blood. Further studies are needed to elucidate the potential importance of the associated genes in neuronal functions and to validate the prognostic or diagnostic value of the individual marks or marker panels.

RepViz: a replicate-driven R tool for visualizing genomic regions

OBJECTIVE

Visualization of sequencing data is an integral part of genomic data analysis. Although there are several tools to visualize sequencing data on genomic regions, they do not offer user-friendly ways to view simultaneously different groups of replicates. To address this need, we developed a tool that allows efficient viewing of both intra- and intergroup variation of sequencing counts on a genomic region, as well as their comparison to the output of user selected analysis methods, such as peak calling.

RESULTS

We present an R package RepViz for replicate-driven visualization of genomic regions. With ChIP-seq and ATAC-seq data we demonstrate its potential to aid visual inspection involved in the evaluation of normalization, outlier behavior, detected features from differential peak calling analysis, and combined analysis of multiple data types. RepViz is readily available on Bioconductor ( https://www.bioconductor.org/packages/devel/bioc/html/RepViz.html ) and on Github ( https://github.com/elolab/RepViz ).

Intratumoral androgen levels are linked to TMPRSS2-ERG fusion in prostate cancer

Intratumoral androgen biosynthesis is one of the mechanisms involved in the progression of prostate cancer, and an important target for novel prostate cancer therapies. Using gas chromatography-tandem mass spectrometry and genome-wide RNA sequencing, we have analyzed androgen concentrations and androgen-regulated gene expression in cancerous and morphologically benign prostate tissue specimens and serum samples obtained from 48 primary prostate cancer patients. Intratumoral dihydrotestosterone (DHT) concentrations were significantly higher in the cancerous tissues compared to benign prostate (P < 0.001). The tissue/serum ratios of androgens were highly variable between the patients, indicating individual patterns of androgen metabolism and/or uptake of androgens within the prostate tissue. An unsupervised hierarchical clustering analysis of intratissue androgen concentrations indicated that transmembrane protease, serine 2/ETS-related gene (TMPRSS2-ERG)-positive patients have different androgen profiles compared to TMPRSS2-ERG-negative patients. TMPRSS2-ERG gene fusion status was also associated with an enhanced androgen-regulated gene expression, along with altered intratumoral androgen metabolism, demonstrated by reduced testosterone concentrations and increased DHT/testosterone ratios in TMPRSS2-ERG-positive tumors. TMPRSS2-ERG-positive and -negative prostate cancer specimens have distinct intratumoral androgen profiles, possibly due to activation of testosterone-independent DHT biosynthesis via the alternative pathway in TMPRSS2-ERG-positive tumors. Thus, patients with TMPRSS2-ERG-positive prostate cancer may benefit from novel inhibitors targeting the alternative DHT biosynthesis.

RNA Polymerase III Subunit POLR3G Regulates Specific Subsets of PolyA+ and SmallRNA Transcriptomes and Splicing in Human Pluripotent Stem Cells

POLR3G is expressed at high levels in human pluripotent stem cells (hPSCs) and is required for maintenance of stem cell state through mechanisms not known in detail. To explore how POLR3G regulates stem cell state, we carried out deep-sequencing analysis of polyA⁺ and smallRNA transcriptomes present in hPSCs and regulated in POLR3G-dependent manner. Our data reveal that POLR3G regulates a specific subset of the hPSC transcriptome, including multiple transcript types, such as protein-coding genes, long intervening non-coding RNAs, microRNAs and small nucleolar RNAs, and affects RNA splicing. The primary function of POLR3G is in the maintenance rather than repression of transcription. The majority of POLR3G polyA⁺ transcriptome is regulated during differentiation, and the key pluripotency factors bind to the promoters of at least 30% of the POLR3G-regulated transcripts. Among the direct targets of POLR3G, POLG is potentially important in sustaining stem cell status in a POLR3G-dependent manner.

Antiandrogens Reduce Intratumoral Androgen Concentrations and Induce Androgen Receptor Expression in Castration-Resistant Prostate Cancer Xenografts

The development of castration-resistant prostate cancer (CRPC) is associated with the activation of intratumoral androgen biosynthesis and an increase in androgen receptor (AR) expression. We recently demonstrated that, similarly to the clinical CRPC, orthotopically grown castration-resistant VCaP (CR-VCaP) xenografts express high levels of AR and retain intratumoral androgen concentrations similar to tumors grown in intact mice. Herein, we show that antiandrogen treatment (enzalutamide or ARN-509) significantly reduced (10-fold, P < 0.01) intratumoral testosterone and dihydrotestosterone concentrations in the CR-VCaP tumors, indicating that the reduction in intratumoral androgens is a novel mechanism by which antiandrogens mediate their effects in CRPC. Antiandrogen treatment also altered the expression of multiple enzymes potentially involved in steroid metabolism. Identical to clinical CRPC, the expression levels of the full-length AR (twofold, P < 0.05) and the AR splice variants 1 (threefold, P < 0.05) and 7 (threefold, P < 0.01) were further increased in the antiandrogen-treated tumors. Nonsignificant effects were observed in the expression of certain classic androgen-regulated genes, such as TMPRSS2 and KLK3, despite the low levels of testosterone and dihydrotestosterone. However, other genes recently identified to be highly sensitive to androgen-regulated AR action, such as NOV and ST6GalNAc1, were markedly altered, which indicated reduced androgen action. Taken together, the data indicate that, besides blocking AR, antiandrogens modify androgen signaling in CR-VCaP xenografts at multiple levels.

The metabolic syndrome in mice overexpressing neuropeptide Y in noradrenergic neurons

A gain-of-function polymorphism in human neuropeptide Y (NPY) gene (rs16139) associates with metabolic disorders and earlier onset of type 2 diabetes (T2D). Similarly, mice overexpressing NPY in noradrenergic neurons (OE-NPY^DBH) display obesity and impaired glucose metabolism. In this study, the metabolic syndrome-like phenotype was characterized and mechanisms of impaired hepatic fatty acid, cholesterol and glucose metabolism in pre-obese (2-month-old) and obese (4-7-month-old) OE-NPY^DBH mice were elucidated. Susceptibility to T2D was assessed by subjecting mice to high caloric diet combined with low-dose streptozotocin. Contribution of hepatic Y1-receptor to the phenotype was studied using chronic treatment with an Y1-receptor antagonist, BIBO3304. Obese OE-NPY^DBH mice displayed hepatosteatosis and hypercholesterolemia preceded by decreased fatty acid oxidation and accelerated cholesterol synthesis. Hyperinsulinemia in early obese state inhibited pyruvate- and glucose-induced hyperglycemia, and deterioration of glucose metabolism of OE-NPY^DBH mice developed with aging. Furthermore, streptozotocin induced T2D only in OE-NPY^DBH mice. Hepatic inflammation was not morphologically visible, but upregulated hepatic anti-inflammatory pathways and increased 8-isoprostane combined with increased serum resistin and decreased interleukin 10 pointed to increased NPY-induced oxidative stress that may predispose OE-NPY^DBH mice to insulin resistance. Chronic treatment with BIBO3304 did not improve the metabolic status of OE-NPY^DBH mice. Instead, downregulation of beta-1-adrenoceptors suggests indirect actions of NPY via inhibition of sympathetic nervous system. In conclusion, changes in hepatic fatty acid, cholesterol and glucose metabolism favoring energy storage contribute to the development of NPY-induced metabolic syndrome, and the effect is likely mediated by changes in sympathetic nervous system activity.

Gene expression profiles and signaling mechanisms in α2B-adrenoceptor-evoked proliferation of vascular smooth muscle cells

BACKGROUND

α2-adrenoceptors are important regulators of vascular tone and blood pressure. Regulation of cell proliferation is a less well investigated consequence of α2-adrenoceptor activation. We have previously shown that α2B-adrenoceptor activation stimulates proliferation of vascular smooth muscle cells (VSMCs). This may be important for blood vessel development and plasticity and for the pathology and therapeutics of cardiovascular disorders. The underlying cellular mechanisms have remained mostly unknown. This study explored pathways of regulation of gene expression and intracellular signaling related to α2B-adrenoceptor-evoked VSMC proliferation.

RESULTS

The cellular mechanisms and signaling pathways of α2B-adrenoceptor-evoked proliferation of VSMCs are complex and include redundancy. Functional enrichment analysis and pathway analysis identified differentially expressed genes associated with α2B-adrenoceptor-regulated VSMC proliferation. They included the upregulated genes Egr1, F3, Ptgs2 and Serpine1 and the downregulated genes Cx3cl1, Cav1, Rhoa, Nppb and Prrx1. The most highly upregulated gene, Lypd8, represents a novel finding in the VSMC context. Inhibitor library screening and kinase activity profiling were applied to identify kinases in the involved signaling pathways. Putative upstream kinases identified by two different screens included PKC, Raf-1, Src, the MAP kinases p38 and JNK and the receptor tyrosine kinases EGFR and HGF/HGFR. As a novel finding, the Src family kinase Lyn was also identified as a putative upstream kinase.

CONCLUSIONS

α2B-adrenoceptors may mediate their pro-proliferative effects in VSMCs by promoting the activity of bFGF and PDGF and the growth factor receptors EGFR, HGFR and VEGFR-1/2. The Src family kinase Lyn was also identified as a putative upstream kinase. Lyn is known to be expressed in VSMCs and has been identified as an important regulator of GPCR trafficking and GPCR effects on cell proliferation. Identified Ser/Thr kinases included several PKC isoforms and the β-adrenoceptor kinases 1 and 2. Cross-talk between the signaling mechanisms involved in α2B-adrenoceptor-evoked VSMC proliferation thus appears to involve PKC activation, subsequent changes in gene expression, transactivation of EGFR, and modulation of kinase activities and growth factor-mediated signaling. While many of the identified individual signals were relatively small in terms of effect size, many of them were validated by combining pathway analysis and our integrated screening approach.

Gut Microbiota Analysis Results Are Highly Dependent on the 16S rRNA Gene Target Region, Whereas the Impact of DNA Extraction Is Minor

Next-generation sequencing (NGS) is currently the method of choice for analyzing gut microbiota composition. As gut microbiota composition is a potential future target for clinical diagnostics, it is of utmost importance to enhance and optimize the NGS analysis procedures. Here, we have analyzed the impact of DNA extraction and selected 16S rDNA primers on the gut microbiota NGS results. Bacterial DNA from frozen stool specimens was extracted with 5 commercially available DNA extraction kits. Special attention was paid to the semiautomated DNA extraction methods that could expedite the analysis procedure, thus being especially suitable for clinical settings. The microbial composition was analyzed with 2 distinct protocols: 1 targeting the V3-V4 and the other targeting the V4-V5 area of the bacterial 16S rRNA gene. The overall effect of DNA extraction on the gut microbiota 16S rDNA profile was relatively small, whereas the 16S rRNA gene target region had an immense impact on the results. Furthermore, semiautomated DNA extraction methods clearly appeared suitable for NGS procedures, proposing that application of these methods could importantly reduce hands-on time and human errors without compromising the validity of results.

Neuropeptide Y Overexpressing Female and Male Mice Show Divergent Metabolic but Not Gut Microbial Responses to Prenatal Metformin Exposure

Background

Prenatal metformin exposure has been shown to improve the metabolic outcome in the offspring of high fat diet fed dams. However, if this is evident also in a genetic model of obesity and whether gut microbiota has a role, is not known.

Methods

The metabolic effects of prenatal metformin exposure were investigated in a genetic model of obesity, mice overexpressing neuropeptide Y in the sympathetic nervous system and in brain noradrenergic neurons (OE-NPY^DβH). Metformin was given for 18 days to the mated female mice. Body weight, body composition, glucose tolerance and serum parameters of the offspring were investigated on regular diet from weaning and sequentially on western diet (at the age of 5–7 months). Gut microbiota composition was analysed by 16S rRNA sequencing at 10–11 weeks.

Results

In the male offspring, metformin exposure inhibited weight gain. Moreover, weight of white fat depots and serum insulin and lipids tended to be lower at 7 months. In contrast, in the female offspring, metformin exposure impaired glucose tolerance at 3 months, and subsequently increased body weight gain, fat mass and serum cholesterol. In the gut microbiota, a decline in Erysipelotrichaceae and Odoribacter was detected in the metformin exposed offspring. Furthermore, the abundance of Sutterella tended to be decreased and Parabacteroides increased. Gut microbiota composition of the metformin exposed male offspring correlated to their metabolic phenotype.

Conclusion

Prenatal metformin exposure caused divergent metabolic phenotypes in the female and male offspring. Nevertheless, gut microbiota of metformin exposed offspring was similarly modified in both genders.

Gut Microbiota Richness and Composition and Dietary Intake of Overweight Pregnant Women Are Related to Serum Zonulin Concentration, a Marker for Intestinal Permeability

BACKGROUND

Increased intestinal permeability may precede adverse metabolic conditions. The extent to which the composition of the gut microbiota and diet contribute to intestinal permeability during pregnancy is unknown.

OBJECTIVE

The aim was to investigate whether the gut microbiota and diet differ according to serum zonulin concentration, a marker of intestinal permeability, in overweight pregnant women.

METHODS

This cross-sectional study included 100 overweight women [mean age: 29 y; median body mass index (in kg/m(2)): 30] in early pregnancy (<17 wk of gestation; median: 13 wk). Serum zonulin (primary outcome) was determined by using ELISA, gut microbiota by 16S ribosomal RNA sequencing, and dietary intake of macro- and micronutrients from 3-d food diaries. The Mann-Whitney U test was used for pairwise comparisons and linear regression and Spearman's nonparametric correlations for relations between serum zonulin and other outcome variables.

RESULTS

Women were divided into "low" (<46.4 ng/mL) and "high" (≥46.4 ng/mL) serum zonulin groups on the basis of the median concentration of zonulin (46.4 ng/mL). The richness of the gut microbiota (Chao 1, observed species and phylogenetic diversity) was higher in the low zonulin group than in the high zonulin group (P = 0.01). The abundances of Bacteroidaceae and Veillonellaceae, Bacteroides and Blautia, and Blautia sp. were lower and of Faecalibacterium and Faecalibacterium prausnitzii higher (P < 0.05) in the low zonulin group than in the high zonulin group. Dietary quantitative intakes of n-3 (ω-3) polyunsaturated fatty acids (PUFAs), fiber, and a range of vitamins and minerals were higher (P < 0.05) in women in the low zonulin group than those in the high zonulin group.

CONCLUSIONS

The richness and composition of the gut microbiota and the intake of n-3 PUFAs, fiber, and a range of vitamins and minerals in overweight pregnant women are associated with serum zonulin concentration. Modification of the gut microbiota and diet may beneficially affect intestinal permeability, leading to improved metabolic health of both the mother and fetus. This trial was registered at clinicaltrials.gov as NCT01922791.

IRF4 Deficiency Leads to Altered BCR Signalling Revealed by Enhanced PI3K Pathway, Decreased SHIP Expression and Defected Cytoskeletal Responses

The graded expression of transcription factor interferon regulatory factor 4 (IRF4) regulates B cell development and is critical for plasma cell differentiation. However, the mechanisms, by which IRF4 elicits its crucial tasks, are largely unknown. To characterize the molecular targets of IRF4 in B cells, we established an IRF4-deficient DT40 B cell line. We found that in the absence of IRF4, the expression of several molecules involved in BCR signalling was altered. For example, the expression of B cell adaptor for PI3K (BCAP) was upregulated, whereas the SHIP (SH2-containing Inositol 5?-Phosphatase) expression was downregulated. These molecular unbalances were accompanied by increased BCR-induced calcium signalling, attenuated B cell linker protein (BLNK) and ERK activity and enhanced activity of PI3K/protein kinase B (Akt) pathway. Further, the IRF4-deficient cells showed dramatically diminished cytoskeletal responses to anti-IgM cross-linking. Our results show that IRF4 has an important role in the regulation of BCR signalling and help to shed light on the molecular mechanisms of B cell development and germinal centre response.

Retene causes multifunctional transcriptomic changes in the heart of rainbow trout (Oncorhynchus mykiss) embryos

Fish are particularly sensitive to aryl hydrocarbon receptor (AhR)-mediated developmental toxicity. The molecular mechanisms behind these adverse effects have remained largely unresolved in salmonids, and for AhR-agonistic polycyclic aromatic hydrocarbons (PAHs). This study explored the cardiac transcriptome of rainbow trout (Oncorhynchus mykiss) eleuteroembryos exposed to retene, an AhR-agonistic PAH. The embryos were exposed to retene (nominal concentration 32 μg/L) and control, their hearts were collected before, at and after the onset of the visible signs of developmental toxicity, and transcriptomic changes were studied by microarray analysis. Retene up- or down-regulated 122 genes. The largest Gene Ontology groups were signal transduction, transcription, apoptosis, cell growth, cytoskeleton, cell adhesion/mobility, cardiovascular development, xenobiotic metabolism, protein metabolism, lipid metabolism and transport, and amino acid metabolism. Together these findings suggest that retene affects multiple signaling cascades in the heart of rainbow trout embryos, and potentially disturbs processes related to cardiovascular development and function.

Differential Promoter Methylation of Macrophage Genes Is Associated With Impaired Vascular Growth in Ischemic Muscles of Hyperlipidemic and Type 2 Diabetic Mice: Genome-Wide Promoter Methylation Study

RATIONALE

Hyperlipidemia and type 2 diabetes mellitus (T2DM) severely impair adaptive vascular growth responses in ischemic muscles. This is largely attributed to dysregulated gene expression, although details of the changes are unknown.

OBJECTIVE

To define the role of promoter methylation in adaptive vascular growth in hyperlipidemia (LDLR(-/-)ApoB(100/100)) and T2DM (IGF-II/LDLR(-/-)ApoB(100/100)) mouse models of hindlimb ischemia.

METHODS AND RESULTS

Unilateral hindlimb ischemia was induced by ligating femoral artery. Perfusion was assessed using ultrasound, and capillary and arteriole parameters were assessed using immunohistochemistry. Genome-wide methylated DNA sequencing was performed with DNA isolated from ischemic muscle, tissue macrophages (Mϕs), and endothelial cells. Compared with the controls, hyperlipidemia and T2DM mice showed impaired perfusion recovery, which was associated with impaired angiogenesis and arteriogenesis. Genome-wide proximal promoter DNA methylation analysis suggested differential patterns of methylation in Mϕ genes in ischemic muscles. Classically activated M1-Mϕ gene promoters, including Cfb, Serping1, and Tnfsf15, were significantly hypomethylated, whereas alternatively activated M2-Mϕ gene promoters, including Nrp1, Cxcr4, Plxnd1, Arg1, Cdk18, and Fes, were significantly hypermethylated in Mϕs isolated from hyperlipidemia and T2DM ischemic muscles compared with controls. These results combined with mRNA expression and immunohistochemistry showed the predominance of proinflammatory M1-Mϕs, compared with anti-inflammatory and proangiogenic M2-Mϕs in hyperlipidemia and T2DM ischemic muscles.

CONCLUSIONS

We found significant promoter hypomethylation of genes typical for proinflammatory M1-Mϕs and hypermethylation of anti-inflammatory, proangiogenic M2-Mϕ genes in hyperlipidemia and T2DM ischemic muscles. Epigenetic alterations modify Mϕ phenotype toward proinflammatory M1 as opposed to anti-inflammatory, proangiogenic, and tissue repair M2 phenotype, which may contribute to the impaired adaptive vascular growth under these pathological conditions.

Seasonal gene expression kinetics between diapause phases in Drosophila virilis group species and overwintering differences between diapausing and non-diapausing females

Most northern insect species experience a period of developmental arrest, diapause, which enables them to survive over the winter and postpone reproduction until favorable conditions. We studied the timing of reproductive diapause and its long-term effects on the cold tolerance of Drosophila montana, D. littoralis and D. ezoana females in seasonally varying environmental conditions. At the same time we traced expression levels of 219 genes in D. montana using a custom-made microarray. We show that the seasonal switch to reproductive diapause occurs over a short time period, and that overwintering in reproductive diapause has long-lasting effects on cold tolerance. Some genes, such as Hsc70, Jon25Bi and period, were upregulated throughout the diapause, while others, including regucalcin, couch potato and Thor, were upregulated only at its specific phases. Some of the expression patterns induced during the sensitive stage, when the females either enter diapause or not, remained induced regardless of the later conditions. qPCR analyses confirmed the findings of the microarray analysis in D. montana and revealed similar gene expression changes in D. littoralis and D. ezoana. The present study helps to achieve a better understanding of the genetic regulation of diapause and of the plasticity of seasonal responses in general.

Maternal diabetes induces changes in the umbilical cord gene expression

INTRODUCTION

Since maternal diabetes may affect fetal development and the umbilical cord provides an extension of the fetal vasculature, we decided to investigate cords' biological responses to maternal diabetic milieu.

METHODS

Using microarray analysis, we determined the gene expression profiles in the umbilical cords of six neonates born to type 1 diabetic mothers and in six control cords. Umbilical cord tissue was collected immediately after elective cesarean section. Expression data were confirmed by real-time polymerase chain reaction (11 genes). Additionally, the same umbilical cords were analyzed histologically.

RESULTS

Two hundred eighty six genes were differentially expressed in the umbilical cords from diabetic pregnancies compared to the controls (fold change ±1.5 and P < 0.01). Maternal diabetes had a major effect on the expression of genes involved in vascular development (Bone morphogenetic protein 4, Delta-like 1, and Notch homolog 4), vessel wall integrity (Collagen type VIII alpha 1, Myocyte enhancer factor 2C, and Matrix metalloproteinase 2), and vascular function (Natriuretic peptide precursor B, Endothelin 1, Endothelin receptor B, Cyclooxygenase 1, and Phosphodiesterase 5A). Maternal diabetes was associated with thicker umbilical vein intima-media layers and larger umbilical vein and artery intima-media areas compared to the controls.

DISCUSSION

Maternal diabetic environment seems to alter umbilical cord expression of genes involved in the regulation of vascular development and function with simultaneous umbilical vessel muscle layer thickening. These alterations suggest vascular phenotypic modifications, which in turn may lead to long-term vascular consequences in various tissues in infants of diabetic mothers.

Early dental epithelial transcription factors distinguish ameloblastoma from keratocystic odontogenic tumor

The aim of the study was to characterize the molecular relationship between ameloblastoma and keratocystic odontogenic tumor (KCOT) by means of a genome-wide expression analysis. Total RNA from 27 fresh tumor samples of 15 solid/multicystic intraosseous ameloblastomas and 12 sporadic KCOTs was hybridized on Affymetrix whole genome arrays. Hierarchical clustering separated ameloblastomas and KCOTs into 2 distinct groups. The gene set enrichment analysis based on 303 dental genes showed a similar separation of ameloblastomas and KCOTs. Early dental epithelial markers PITX2, MSX2, DLX2, RUNX1, and ISL1 were differentially overexpressed in ameloblastoma, indicating its dental identity. Also, PTHLH, a hormone involved in tooth eruption and invasive growth, was one of the most differentially upregulated genes in ameloblastoma. The most differentially overexpressed genes in KCOT were squamous epithelial differentiation markers SPRR1A, KRTDAP, and KRT4, as well as DSG1, a component of desmosomal cell-cell junctions. Additonally, the epithelial stem cell marker SOX2 was significantly upregulated in KCOT when compared with ameloblastoma. Taken together, the gene expression profile of ameloblastoma reflects differentiation from dental lamina toward the cap/bell stage of tooth development, as indicated by dental epithelium-specific transcription factors. In contrast, gene expression of KCOT indicates differentiation toward keratinocytes.

A note on an exon-based strategy to identify differentially expressed genes in RNA-seq experiments

RNA-sequencing (RNA-seq) has rapidly become the method of choice in many genome-wide transcriptomic studies. To meet the high expectations posed by this technology, powerful computational techniques are needed to translate the measurements into biological and biomedical understanding. A number of statistical procedures have already been developed to identify differentially expressed genes between distinct sample groups. With these methods statistical testing is typically performed after the data has been summarized at the gene level. As an alternative strategy, developed with the aim to improve the results, we demonstrate a method in which statistical testing at the exon level is performed prior to the summary of the results at the gene level. Using publicly available RNA-seq datasets as case studies, we illustrate how this exon-based strategy can improve the performance of the widely used differential expression software packages as compared to the conventional gene-based strategy. In particular, we show how it enables robust detection of moderate but systematic changes that are missed when relying on single gene-level summary counts only.

Fermentable fibres condition colon microbiota and promote diabetogenesis in NOD mice

AIMS/HYPOTHESIS

Gut microbiota (GM) and diet both appear to be important in the pathogenesis of type 1 diabetes. Fermentable fibres (FFs), of which there is an ample supply in natural, diabetes-promoting diets, are used by GM as a source of energy. Our aim was to determine whether FFs modify GM and diabetes incidence in the NOD mouse.

METHODS

Female NOD mice were weaned to a semisynthetic diet and the effects of FF supplementation on diabetes incidence and insulitis were evaluated. Real-time quantitative PCR was employed to determine the effects imposed to gene transcripts in the colon and lymph nodes. Changes to GM were analysed by next-generation sequencing.

RESULTS

NOD mice fed semisynthetic diets free from FFs were largely protected from diabetes while semisynthetic diets supplemented with the FFs pectin and xylan (PX) resulted in higher diabetes incidence. Semisynthetic diet free from FFs altered GM composition significantly; addition of PX changed the composition of the GM towards that found in natural-diet-fed mice and increased production of FF-derived short-chain fatty acid metabolites in the colon. The highly diabetogenic natural diet was associated with expression of proinflammatory and stress-related genes in the colon, while the semisynthetic diet free from FFs promoted Il4, Il22, Tgfβ and Foxp3 transcripts in the colon and/or pancreatic lymph node. PX in the same diet counteracted these effects and promoted stress-related IL-18 activation in gut epithelial cells. 16S RNA sequencing revealed each diet to give rise to its particular GM composition, with different Firmicutes to Bacteroidetes ratios, and enrichment of mucin-degrading Ruminococcaceae following diabetes-protective FF-free diet.

CONCLUSIONS/INTERPRETATION

FFs condition microbiota, affect colon homeostasis and are important components of natural, diabetes-promoting diets in NOD mice.

Promoter-specific alterations of APC are a rare cause for mutation-negative familial adenomatous polyposis

n familial adenomatous polyposis (FAP), 20% of classical and 70% of attenuated/atypical (AFAP) cases remain mutation-negative after routine testing; yet, allelic expression imbalance may suggest an APC alteration. Our aim was to determine the proportion of families attributable to genetic or epigenetic changes in the APC promoter region. We studied 51 unrelated families/cases (26 with classical FAP and 25 with AFAP) with no point mutations in the exons and exon/intron borders and no rearrangements by multiplex ligation-dependent probe amplification (MLPA, P043-B1). Promoter-specific events of APC were addressed by targeted resequencing, MLPA (P043-C1), methylation-specific MLPA, and Sanger sequencing of promoter regions. A novel 132-kb deletion encompassing the APC promoter 1B and upstream sequence occurred in a classical FAP family with allele-specific APC expression. No promoter-specific point mutations or hypermethylation were present in any family. In conclusion, promoter-specific alterations are a rare cause for mutation-negative FAP (1/51, 2%). The frequency and clinical correlations of promoter 1B deletions are poorly defined. This investigation provides frequencies of 1/26 (4%) for classical FAP, 0/25 (0%) for AFAP, and 1/7 (14%) for families with allele-specific expression of APC. Clinically, promoter 1B deletions may associate with classical FAP without extracolonic manifestations.

Gene expression differences between Noccaea caerulescens ecotypes help to identify candidate genes for metal phytoremediation

Populations of Noccaea caerulescens show tremendous differences in their capacity to hyperaccumulate and hypertolerate metals. To explore the differences that could contribute to these traits, we undertook SOLiD high-throughput sequencing of the root transcriptomes of three phenotypically well-characterized N. caerulescens accessions, i.e., Ganges, La Calamine, and Monte Prinzera. Genes with possible contribution to zinc, cadmium, and nickel hyperaccumulation and hypertolerance were predicted. The most significant differences between the accessions were related to metal ion (di-, trivalent inorganic cation) transmembrane transporter activity, iron and calcium ion binding, (inorganic) anion transmembrane transporter activity, and antioxidant activity. Analysis of correlation between the expression profile of each gene and the metal-related characteristics of the accessions disclosed both previously characterized (HMA4, HMA3) and new candidate genes (e.g., for nickel IRT1, ZIP10, and PDF2.3) as possible contributors to the hyperaccumulation/tolerance phenotype. A number of unknown Noccaea-specific transcripts also showed correlation with Zn(2+), Cd(2+), or Ni(2+) hyperaccumulation/tolerance. This study shows that N. caerulescens populations have evolved great diversity in the expression of metal-related genes, facilitating adaptation to various metalliferous soils. The information will be helpful in the development of improved plants for metal phytoremediation.

Whole transcriptome characterization of the effects of dehydration and rehydration on Cladonia rangiferina, the grey reindeer lichen

BACKGROUND

Lichens are symbiotic organisms with a fungal and an algal or a cyanobacterial partner. Lichens inhabit some of the harshest climates on earth and most lichen species are desiccation-tolerant. Lichen desiccation-tolerance has been studied at the biochemical level and through proteomics, but the underlying molecular genetic mechanisms remain largely unexplored. The objective of our study was to examine the effects of dehydration and rehydration on the gene expression of Cladonia rangiferina.

RESULTS

Samples of C. rangiferina were collected at several time points during both the dehydration and rehydration process and the gene expression intensities were measured using a custom DNA microarray. Several genes, which were differentially expressed in one or more time points, were identified. The microarray results were validated using qRT-PCR analysis. Enrichment analysis of differentially expressed transcripts was also performed to identify the Gene Ontology terms most associated with the rehydration and dehydration process.

CONCLUSIONS

Our data identify differential expression patterns for hundreds of genes that are modulated during dehydration and rehydration in Cladonia rangiferina. These dehydration and rehydration events clearly differ from each other at the molecular level and the largest changes to gene expression are observed within minutes following rehydration. Distinct changes are observed during the earliest stage of rehydration and the mechanisms not appear to be shared with the later stages of wetting or with drying. Several of the most differentially expressed genes are similar to genes identified in previous studies that have investigated the molecular mechanisms of other desiccation-tolerant organisms. We present here the first microarray experiment for any lichen species and have for the first time studied the genetic mechanisms behind lichen desiccation-tolerance at the whole transcriptome level.

Comparison of software packages for detecting differential expression in RNA-seq studies

RNA-sequencing (RNA-seq) has rapidly become a popular tool to characterize transcriptomes. A fundamental research problem in many RNA-seq studies is the identification of reliable molecular markers that show differential expression between distinct sample groups. Together with the growing popularity of RNA-seq, a number of data analysis methods and pipelines have already been developed for this task. Currently, however, there is no clear consensus about the best practices yet, which makes the choice of an appropriate method a daunting task especially for a basic user without a strong statistical or computational background. To assist the choice, we perform here a systematic comparison of eight widely used software packages and pipelines for detecting differential expression between sample groups in a practical research setting and provide general guidelines for choosing a robust pipeline. In general, our results demonstrate how the data analysis tool utilized can markedly affect the outcome of the data analysis, highlighting the importance of this choice.

Continuous hypoxic culturing of human embryonic stem cells enhances SSEA-3 and MYC levels

Low oxygen tension (hypoxia) contributes critically to pluripotency of human embryonic stem cells (hESCs) by preventing spontaneous differentiation and supporting self-renewal. However, it is not well understood how hESCs respond to reduced oxygen availability and what are the molecular mechanisms maintaining pluripotency in these conditions. In this study we characterized the transcriptional and molecular responses of three hESC lines (H9, HS401 and HS360) on short (2 hours), intermediate (24 hours) and prolonged (7 days) exposure to low oxygen conditions (4% O₂). In response to prolonged hypoxia the expression of pluripotency surface marker SSEA-3 was increased. Furthermore, the genome wide gene-expression analysis revealed that a substantial proportion (12%) of all hypoxia-regulated genes in hESCs, were directly linked to the mechanisms controlling pluripotency or differentiation. Moreover, transcription of MYC oncogene was induced in response to continuous hypoxia. At the protein level MYC was stabilized through phosphorylation already in response to a short hypoxic exposure. Total MYC protein levels remained elevated throughout all the time points studied. Further, MYC protein expression in hypoxia was affected by silencing HIF2α, but not HIF1α. Since MYC has a crucial role in regulating pluripotency we propose that induction of sustained MYC expression in hypoxia contributes to activation of transcriptional programs critical for hESC self-renewal and maintenance of enhanced pluripotent state.

Transcriptome profiling of the murine testis during the first wave of spermatogenesis

Correct gene expression patterns form the basis for male germ cell differentiation and male fertility. Although previous studies have elucidated the importance of testis specific gene expression, the exact transcripts and comprehensive gene expression patterns remain unknown. Large scale sequencing techniques have enabled cost effective analysis of gene expression and isoform studies. Using the SOLiD 4 next-generation sequencing platform we have investigated the gene expression patterns at five different time points during the first wave on murine spermatogenesis. Our results highlight the upregulation of spermatogenesis related biological processes and associated cellular components. Elucidation of differential gene expression at important time points during the sperm development emphasizes the importance of correct timing of gene expression within biological processes. Differential gene level expression was analyzed with R/Bioconductor’s Limma package and isoform analysis was conducted with the Cufflinks pipeline. At gene level total of 2494 differentially expressed genes were identified and Cufflinks characterized over 160 000 gene isoforms, of which 29% were novel transcripts assigned to known genes. Isoforms were detected for 57% of expressed genes and in a total over 26 000 genes were expressed in the testis. Differential promoter and transcription start site usage appears also to play a role in regulation of gene expression during spermatogenesis. Furthermore, we identified 947 upregulated long non-coding RNAs during the first wave of spermatogenesis. These RNAs appeared to be highly specific to different time points. Transcriptomic analysis of testis tissue samples is highly informative due to the large number of expressed genes and identified isoforms. Our study provides a very valuable basis for investigation of gene isoforms and regulation and factors contributing to male fertility.

SAP30L (Sin3A-associated protein 30-like) is involved in regulation of cardiac development and hematopoiesis in zebrafish embryos

The Sin3A-associated proteins SAP30 and SAP30L share 70% sequence identity and are part of the multiprotein Sin3A corepressor complex. They participate in gene repression events by linking members of the complex and stabilizing interactions among the protein members as well as between proteins and DNA. While most organisms have both SAP30 and SAP30L, the zebrafish is exceptional because it only has SAP30L. Here we demonstrate that SAP30L is expressed ubiquitously in embryonic and adult zebrafish tissues. Knockdown of SAP30L using morpholino-mediated technology resulted in a morphant phenotype manifesting as cardiac insufficiency and defective hemoglobinization of red blood cells. A microarray analysis of gene expression in SAP30L morphant embryos revealed changes in the expression of genes involved in regulation of transcription, TGF-beta signaling, Wnt-family transcription factors, and nuclear genes encoding mitochondrial proteins. The expression of the heart-specific nkx2.5 gene was markedly down-regulated in SAP30L morphants, and the cardiac phenotype could be partially rescued by nkx2.5 mRNA. In addition, changes were detected in the expression of genes known to be important in hemoglobin synthesis and erythropoiesis. Our results demonstrate that SAP30L regulates several transcriptional pathways in zebrafish embryos and is involved in the development of cardiac and hematopoietic systems.

Expression of small nucleolar RNAs in leukemic cells

PURPOSE

Small nucleolar RNAs (snoRNAs) direct sequence-specific modifications to ribosomal RNA. We hypothesized that the expression of snoRNAs may be altered in leukemic cells.

METHODS

The expression of snoRNAs was analyzed in various leukemic cell lines by massive parallel sequencing (SOLiD). Quantitative real-time PCR (RT-qPCR) was used to validate the expression profiles.

RESULTS

Our results show characteristic differences in the expression patterns of snoRNAs between cell lines representing the main subgroups of leukemia, AML, pre-B-ALL and T-ALL, respectively. In RT-qPCR analyses, several snoRNAs were found to be differentially expressed in T-ALL as compared to pre-B-ALL cell lines.

CONCLUSIONS

snoRNAs are differentially expressed in various leukemic cell lines and could, therefore, be potentially useful in the classification of leukemia subgroups.