Quantitative transcriptomic and epigenomic data analysis: a primer.

The advent of microarray and second generation sequencing technology has revolutionized the field of molecular biology, allowing researchers to quantitatively assess transcriptomic and epigenomic features in a comprehensive and cost-efficient manner. Moreover, technical advancements have pushed the resolution of these sequencing techniques to the single cell level. As a result, the bottleneck of molecular biology research has shifted from the bench to the subsequent omics data analysis. Even though most methodologies share the same general strategy, state-of-the-art literature typically focuses on data type specific approaches and already assumes expert knowledge. Here, however, we aim at providing conceptual insight in the principles of genome-wide quantitative transcriptomic and epigenomic (including open chromatin assay) data analysis by describing a generic workflow. By starting from a general framework and its assumptions, the need for alternative or additional data-analytical solutions when working with specific data types becomes clear, and are hence introduced. Thus, we aim to enable readers with basic omics expertise to deepen their conceptual and statistical understanding of general strategies and pitfalls in omics data analysis and to facilitate subsequent progression to more specialized literature.

The phenylalanine ammonia-lyase inhibitor AIP induces rice defence against the root-knot nematode Meloidogyne graminicola.

The phenylalanine ammonia-lyase (PAL) enzyme catalyses the conversion of l-phenylalanine to trans-cinnamic acid. This conversion is the first step in phenylpropanoid biosynthesis in plants. The phenylpropanoid pathway produces diverse plant metabolites that play essential roles in various processes, including structural support and defence. Previous studies have shown that mutation of the PAL genes enhances disease susceptibility. Here, we investigated the functions of the rice PAL genes using 2-aminoindan-2-phosphonic acid (AIP), a strong competitive inhibitor of PAL enzymes. We show that the application of AIP can significantly reduce the PAL activity of rice crude protein extracts in vitro. However, when AIP was applied to intact rice plants, it reduced infection of the root-knot nematode Meloidogyne graminicola. RNA-seq showed that AIP treatment resulted in a rapid but transient upregulation of defence-related genes in roots. Moreover, targeted metabolomics demonstrated higher levels of jasmonates and antimicrobial flavonoids and diterpenoids accumulating after AIP treatment. Furthermore, chemical inhibition of the jasmonate pathway abolished the effect of AIP on nematode infection. Our results show that disturbance of the phenylpropanoid pathway by the PAL inhibitor AIP induces defence in rice against M. graminicola by activating jasmonate-mediated defence.

Promoter hypermethylation of neural-related genes is compatible with stemness in solid cancers.

BACKGROUND: DNA hypermethylation is an epigenetic feature that modulates gene expression, and its deregulation is observed in cancer. Previously, we identified a neural-related DNA hypermethylation fingerprint in colon cancer, where most of the top hypermethylated and downregulated genes have known functions in the nervous system. To evaluate the presence of this signature and its relevance to carcinogenesis in general, we considered 16 solid cancer types available in The Cancer Genome Atlas (TCGA). RESULTS: All tested cancers showed significant enrichment for neural-related genes amongst hypermethylated genes. This signature was already present in two premalignant tissue types and could not be explained by potential confounders such as bivalency status or tumor purity. Further characterization of the neural-related DNA hypermethylation signature in colon cancer showed particular enrichment for genes that are overexpressed during neural differentiation. Lastly, an analysis of upstream regulators identified RE1-Silencing Transcription factor (REST) as a potential mediator of this DNA methylation signature. CONCLUSION: Our study confirms the presence of a neural-related DNA hypermethylation fingerprint in various cancers, of genes linked to neural differentiation, and points to REST as a possible regulator of this mechanism. We propose that this fingerprint indicates an involvement of DNA hypermethylation in the preservation of neural stemness in cancer cells.

Liver transcriptomic and methylomic analyses identify transcriptional mitogen-activated protein kinase regulation in facultative hibernation of Syrian hamster.

Hibernation consists of alternating torpor-arousal phases, during which animals cope with repetitive hypothermia and ischaemia-reperfusion. Due to limited transcriptomic and methylomic information for facultative hibernators, we here conducted RNA and whole-genome bisulfide sequencing in liver of hibernating Syrian hamster (Mesocricetus auratus). Gene ontology analysis was performed on 844 differentially expressed genes and confirmed the shift in metabolic fuel utilization, inhibition of RNA transcription and cell cycle regulation as found in seasonal hibernators. Additionally, we showed a so far unreported suppression of mitogen-activated protein kinase (MAPK) and protein phosphatase 1 pathways during torpor. Notably, hibernating hamsters showed upregulation of MAPK inhibitors (dual-specificity phosphatases and sproutys) and reduced levels of MAPK-induced transcription factors (TFs). Promoter methylation was found to modulate the expression of genes targeted by these TFs. In conclusion, we document gene regulation between hibernation phases, which may aid the identification of pathways and targets to prevent organ damage in transplantation or ischaemia-reperfusion.

Value of ultrasound and computed tomography in the diagnosis of acute appendicitis with histopathology as gold standard.

BACKGROUND: To evaluate the value of ultrasound (US), computed tomography (CT) and laboratory testing in diagnosing acute appendicitis and to evaluate the impact on the treatment of acute appendicitis, keeping histopathology as the gold standard. Retrospective evaluation of all patients treated surgically for appendicitis from 1 January 2016 to 31 December 2018 at Sint-Andries Hospital, Tielt, Belgium. METHODS: For each patient, we recorded: age, sex, duration of hospitalization, type of radiological investigation, fever (>37.3 °C), C-reactive protein (CRP), leukocyte particle count (LPC), Alvarado score and the type of surgery. For US and CT sensitivity, specificity and negative appendectomy rate were evaluated. RESULTS: Over a period of 3 years, 304 appendectomies were performed. The overall prevalence of appendicitis was 95.1%. Mean age was 31 years. All patients underwent radiological examination. US was performed in 35.9% (109), CT in 50.3% (153) and CT after US in 13.8% (42) of all cases. The sensitivity and specificity of CT were 99.4% and 80.0%, respectively. For US, the respective figures were 74.8 and 62.5%. CRP and LPC were significantly higher in the appendicitis group, compared to the non-appendicitis group. The negative appendectomy rate was slightly higher in the CT after US group i.e. 7.1% (3/42) compared to 4.67% (5/107) in the US group and 4.50% (7/155) in the CT group. CONCLUSIONS: The diagnostic value of US and CT are both very high. The diagnostic value of clinical features and biochemistry alone is quite low. By combining radiological examination, clinical examination and laboratory values, low negative appendectomy rates can be achieved.

The genome of the extremophile Artemia provides insight into strategies to cope with extreme environments.

BACKGROUND: Brine shrimp Artemia have an unequalled ability to endure extreme salinity and complete anoxia. This study aims to elucidate its strategies to cope with these stressors. RESULTS AND DISCUSSION: Here, we present the genome of an inbred A. franciscana Kellogg, 1906. We identified 21,828 genes of which, under high salinity, 674 genes and under anoxia, 900 genes were differentially expressed (42%, respectively 30% were annotated). Under high salinity, relevant stress genes and pathways included several Heat Shock Protein and Leaf Embryogenesis Abundant genes, as well as the trehalose metabolism. In addition, based on differential gene expression analysis, it can be hypothesized that a high oxidative stress response and endocytosis/exocytosis are potential salt management strategies, in addition to the expression of major facilitator superfamily genes responsible for transmembrane ion transport. Under anoxia, genes involved in mitochondrial function, mTOR signalling and autophagy were differentially expressed. Both high salt and anoxia enhanced degradation of erroneous proteins and protein chaperoning. Compared with other branchiopod genomes, Artemia had 0.03% contracted and 6% expanded orthogroups, in which 14% of the genes were differentially expressed under high salinity or anoxia. One phospholipase D gene family, shown to be important in plant stress response, was uniquely present in both extremophiles Artemia and the tardigrade Hypsibius dujardini, yet not differentially expressed under the described experimental conditions. CONCLUSIONS: A relatively complete genome of Artemia was assembled, annotated and analysed, facilitating research on its extremophile features, and providing a reference sequence for crustacean research.

Clinically relevant aberrant Filip1l DNA methylation detected in a murine model of cutaneous squamous cell carcinoma.

BACKGROUND: Cutaneous squamous cell carcinomas (cSCC) are among the most common and highly mutated human malignancies. Understanding the impact of DNA methylation in cSCC may provide avenues for new therapeutic strategies. METHODS: We used reduced-representation bisulfite sequencing for DNA methylation analysis of murine cSCC. Differential methylation was assessed at the CpG level using limma. Next, we compared with human cSCC Infinium HumanMethylation BeadArray data. Genes were considered to be of major relevance when they featured at least one significantly differentially methylated CpGs (RRBS) / probes (Infinium) with at least a 30% difference between tumour vs. control in both a murine gene and its human orthologue. The human EPIC Infinium data were used to distinguish two cSCC subtypes, stem-cell-like and keratinocyte-like tumours. FINDINGS: We found increased average methylation in mouse cSCC (by 12.8%, p = 0.0011) as well as in stem-cell like (by 3.1%, p=0.002), but not keratinocyte-like (0.2%, p = 0.98), human cSCC. Comparison of differentially methylated genes revealed striking similarities between human and mouse cSCC. Locus specific methylation changes in mouse cSCC often occurred in regions of potential regulatory function, including enhancers and promoters. A key differentially methylated region was located in a potential enhancer of the tumour suppressor gene Filip1l and its expression was reduced in mouse tumours. Moreover, the FILIP1L locus showed hypermethylation in human cSCC and lower expression in human cSCC cell lines. INTERPRETATION: Deregulation of DNA methylation is an important feature of murine and human cSCC that likely contributes to silencing of tumour suppressor genes, as shown for Filip1l. FUNDING: British Skin Foundation, Cancer Research UK.

DNA Methylation Regulates Transcription Factor-Specific Neurodevelopmental but Not Sexually Dimorphic Gene Expression Dynamics in Zebra Finch Telencephalon.

Song learning in zebra finches (Taeniopygia guttata) is a prototypical example of a complex learned behavior, yet knowledge of the underlying molecular processes is limited. Therefore, we characterized transcriptomic (RNA-sequencing) and epigenomic (RRBS, reduced representation bisulfite sequencing; immunofluorescence) dynamics in matched zebra finch telencephalon samples of both sexes from 1 day post hatching (1 dph) to adulthood, spanning the critical period for song learning (20 and 65 dph). We identified extensive transcriptional neurodevelopmental changes during postnatal telencephalon development. DNA methylation was very low, yet increased over time, particularly in song control nuclei. Only a small fraction of the massive differential expression in the developing zebra finch telencephalon could be explained by differential CpG and CpH DNA methylation. However, a strong association between DNA methylation and age-dependent gene expression was found for various transcription factors (i.e., OTX2, AR, and FOS) involved in neurodevelopment. Incomplete dosage compensation, independent of DNA methylation, was found to be largely responsible for sexually dimorphic gene expression, with dosage compensation increasing throughout life. In conclusion, our results indicate that DNA methylation regulates neurodevelopmental gene expression dynamics through steering transcription factor activity, but does not explain sexually dimorphic gene expression patterns in zebra finch telencephalon.

Molecular correlates of hypothalamic development in songbird ontogeny in comparison with the telencephalon.

Development of the songbird brain provides an excellent experimental model for understanding the regulation of sex differences in ontogeny. Considering the regulatory role of the hypothalamus in endocrine, in particular reproductive, physiology, we measured the structural (volume) and molecular correlates of hypothalamic development during ontogeny of male and female zebra finches. We quantified by relative quantitative polymerase chain reaction (rqPCR) the expression of 14 genes related to thyroid and steroid hormones actions as well as 12 genes related to brain plasticity at four specific time points during ontogeny and compared these expression patterns with the expression of the same genes as detected by transcriptomics in the telencephalon. These two different methodological approaches detected specific changes with age and demonstrated that in a substantial number of cases changes observed in both brain regions are nearly identical. Other genes however had a tissue-specific developmental pattern. Sex differences or interactions of sex by age were detected in the expression of a subset of genes, more in hypothalamus than telencephalon. These results correlate with multiple known aspects of the developmental and reproductive physiology but also raise a number of new functional questions.

Exploratory analysis of the human breast DNA methylation profile upon soymilk exposure.

Upon soy consumption, isoflavone metabolites attain bioactive concentrations in breast tissue possibly affecting health. Though in vitro epigenetic activity of soy metabolites has been described, the in vivo impact on the epigenome is largely unknown. Therefore, in this case-control study, the breast glandular tissue DNA methylome was explored in women undergoing an aesthetic breast reduction. After a run-in phase, 10 generally healthy Belgian or Dutch women received soymilk for 5 days. MethylCap-seq methylation profiles were compared with those of 10 matched controls. Isoflavones and their microbial metabolites were quantified in urine, serum, and glandular breast tissue (liquid chromatography-mass spectrometry) and 17ß-estradiol in glandular breast tissue (immunoassay). Global DNA methylation levels were obtained for 6 cases and 5 controls using liquid chromatography-mass spectrometry. Although lower MethylCap-seq coverages were observed, mass spectrometry results and computational LINE-1 methylation analysis did not provide evidence supporting global methylation alterations upon treatment. At a false discovery rate of 0.05, no differentially methylated loci were identified. Moreover, a set of previously identified loci was specifically tested, but earlier reported results could not be validated. In conclusion, after a 5-day soymilk treatment, no major general epigenetic reprogramming in breast tissue could be found in this exploratory study.