Search

Epigenetic processes play a crucial role in the regulation of plant stress responses, but their role in plant–pathogen interactions remains poorly understood. Although histone‐modifying enzymes have been observed to be deregulated in galls induced by root‐knot nematodes (RKN, Meloidogyne graminicola) in rice, their influence on plant defence and their genome‐wide impact has not been comprehensively investigated. First, the role of histone modifications in plant–nematode interactions was confirmed by pharmacological inhibition of histone‐modifying enzymes, which all significantly affected rice susceptibility to RKN. For a more specific view, three histone marks, H3K9ac, H3K9me2, and H3K27me3, were subsequently studied by chromatin‐immunoprecipitation‐sequencing on RKN‐induced galls at 3 days postinoculation. While levels of H3K9ac and H3K27me3 were strongly enriched, H3K9me2 was generally depleted in galls versus control root tips. Differential histone peaks were generally associated with plant defence‐related genes. Transcriptome analysis using RNA‐Seq and RT‐qPCR‐based validation revealed that genes marked with H3K9ac or H3K9me2 showed the expected activation or repression gene expression pattern, but this was not the case for H3K27me3 marks. Our results indicate that histone modifications respond dynamically to RKN infection, and that posttranslational modifications mainly at H3K9 specifically target plant defence‐related genes., Histone modifications respond dynamically to root‐knot nematode infection in rice, with levels of H3K9ac and H3K27me3 enriched, and H3K9me2 generally depleted, in galls versus control root tips.

Hirschmanniella oryzae is one of the most devastating nematodes on rice, leading to substantial yield losses. Effector proteins aid the nematode during the infection process by subduing plant defence responses. In this research we characterized two potential H. oryzae effector proteins, chorismate mutase (HoCM) and isochorismatase (HoICM), and investigated their enzymatic activity and their role in plant immunity. Both HoCM and HoICM proved to be enzymatically active in complementation tests in mutant Escherichia coli strains. Infection success by the migratory nematode H. oryzae was significantly higher in transgenic rice lines constitutively expressing HoCM or HoICM. Expression of HoCM, but not HoICM, increased rice susceptibility against the sedentary nematode Meloidogyne graminicola also. Transcriptome and metabolome analyses indicated reductions in secondary metabolites in the transgenic rice plants expressing the potential nematode effectors. The results presented here demonstrate that both HoCM and HoICM suppress the host immune system and that this may be accomplished by lowering secondary metabolite levels in the plant., Chorismate mutase and isochorismatase from the plant parasitic nematode Hirschmanniella oryzae are able to decrease immunity of rice by interfering with the phenylpropanoid pathway.

Summary A role for DNA hypomethylation has recently been suggested in the interaction between bacteria and plants; it is unclear whether this phenomenon reflects a conserved response.Treatment of plants of monocot rice and dicot tomato with nematode‐associated molecular patterns from different nematode species or bacterial pathogen‐associated molecular pattern flg22 revealed global DNA hypomethylation. A similar hypomethylation response was observed during early gall induction by Meloidogyne graminicola in rice. Evidence for the causal impact of hypomethylation on immunity was revealed by a significantly reduced plant susceptibility upon treatment with DNA methylation inhibitor 5‐azacytidine.Whole‐genome bisulphite sequencing of young galls revealed massive hypomethylation in the CHH context, while not for CG or CHG nucleotide contexts. Further, CHH hypomethylated regions were predominantly associated with gene promoter regions, which was not correlated with activated gene expression at the same time point but, rather, was correlated with a delayed transcriptional gene activation. Finally, the relevance of CHH hypomethylation in plant defence was confirmed in rice mutants of the RNA‐directed DNA methylation pathway and DECREASED DNA METHYLATION 1.We demonstrated that DNA hypomethylation is associated with reduced susceptibility in rice towards root‐parasitic nematodes and is likely to be part of the basal pattern‐triggered immunity response in plants.

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.

MicroRNAs miR390-5p, miR7694-3p miR1868 and miR1849 were found to be suitable miRNA reference genes for rice, under either infection with the root-knot nematode Meloidogyne graminicola or treatment with BABA. RT-qPCR is a widely used method to investigate the expression levels of genes under certain conditions. A key step, however, to have reliable results is the normalization of expression. For every experimental condition, suitable reference genes must be chosen. These reference genes must not be affected by differences in experimental conditions. MicroRNAs are regulatory RNA molecules, able to direct the expression levels of protein coding genes. In plants, their attributed functions range from roles in development to immunity. In this work, microRNAs (miRNAs) are evaluated for their suitability as reference genes in rice after infection with root-knot nematode Meloidogyne graminicola or after priming with beta-amino butyric acid. The evaluation was based on their amplification efficiency and their stability estimates according to geNorm, NormFinder and BestKeeper. All tested miRNAs, excluding one, were considered acceptable for normalization. Furthermore, miRNAs were validated using miRNA sequencing data. The set of microRNAs miR390-5p and miR7694-3p was found to be the most stable combination under the tested conditions. Another miRNA set consisting of miR7694-3p, miR1868 and miR1849 also shows potential to be used for miRNA expression normalization under experimental conditions beyond the scope of this study. This work is the first report on reference miRNAs in rice for the purpose of plant defence studies.

A subset of genetic variants found through screening of patients with hereditary breast and ovarian cancer syndrome (HBOC) and Lynch syndrome impact RNA splicing. Through target enrichment of the transcriptome, it is possible to perform deep‐sequencing and to identify the different and even rare mRNA isoforms. A targeted RNA‐seq approach was used to analyse the naturally‐occurring splicing events for a panel of 8 breast and/or ovarian cancer susceptibility genes (BRCA1, BRCA2, RAD51C, RAD51D, PTEN, STK11, CDH1, TP53), 3 Lynch syndrome genes (MLH1, MSH2, MSH6) and the fanconi anaemia SLX4 gene, in which monoallelic mutations were found in non‐BRCA families. For BRCA1, BRCA2, RAD51C and RAD51D the results were validated by capillary electrophoresis and were compared to a non‐targeted RNA‐seq approach. We also compared splicing events from lymphoblastoid cell‐lines with those from breast and ovarian fimbriae tissues. The potential of targeted RNA‐seq to detect pathogenic changes in RNA‐splicing was validated by the inclusion of samples with previously well characterized BRCA1/2 genetic variants. In our study, we update the catalogue of normal splicing events for BRCA1/2, provide an extensive catalogue of normal RAD51C and RAD51D alternative splicing, and list splicing events found for eight other genes. Additionally, we show that our approach allowed the identification of aberrant splicing events due to the presence of BRCA1/2 genetic variants and distinguished between complete and partial splicing events. In conclusion, targeted‐RNA‐seq can be very useful to classify variants based on their putative pathogenic impact on splicing., What's new? Hereditary familial breast/ovarian cancer (HBOC) syndrome involves numerous pathogenic variants, including variants of uncertain clinical significance (VUS). A subset of VUS, however, is suspected to influence RNA splicing, leading to the expression of potentially pathological transcript isoforms. Here, using a targeted RNA‐seq approach, naturally occurring splice isoforms were described for BRCA1/2, RAD51C, RAD51D, and eight additional tumor‐suppressor genes that are associated with HBOC and Lynch syndrome. The targeted RNA‐seq approach also identified aberrant splicing events associated with the presence of BRCA1/2 genetic variants and successfully distinguished complete from incomplete splicing events, which is of major importance in determining pathogenicity.

Existing compendia of non-coding RNA (ncRNA) are incomplete, in part because they are derived almost exclusively from small and polyadenylated RNAs. Here we present a more comprehensive atlas of the human transcriptome, which includes small and polyA RNA as well as total RNA from 300 human tissues and cell lines. We report thousands of previously uncharacterized RNAs, increasing the number of documented ncRNAs by approximately 8%. To infer functional regulation by known and newly characterized ncRNAs, we exploited pre-mRNA abundance estimates from total RNA sequencing, revealing 316 microRNAs and 3,310 long non-coding RNAs with multiple lines of evidence for roles in regulating protein-coding genes and pathways. Our study both refines and expands the current catalog of human ncRNAs and their regulatory interactions. All data, analyses and results are available for download and interrogation in the R2 web portal, serving as a basis for future exploration of RNA biology and function.

Purpose Colonoscopy and the fecal immunochemical test (FIT) are currently the most widely used screening modalities for colorectal cancer (CRC), however, both with their own limitations. Here we aim to identify and validate stool-based DNA methylation markers for the early detection of CRC and investigate the biological pathways prone to DNA methylation. Methods DNA methylation marker discovery was performed using The Cancer Genome Atlas (TCGA) colon adenocarcinoma data set consisting of normal and primary colon adenocarcinoma tissue. The performance of the five best candidate markers and a previously identified marker, NDRG4, was evaluated on tissues and whole stool samples of healthy subjects and CRC patients using quantitative MSP assays. The results were compared and combined with FIT data. Finally, pathway and gene ontology enrichment analyses were performed using ToppFun, GOrilla and clusterProfiler. Results GDNF, HAND2, SLC35F3, SNAP91 and SORCS1 were ranked as the best performing markers. Gene combinations of all five markers, NDRG4 and FIT were evaluated to establish the biomarker panel with the highest diagnostic potential, resulting in the identification of GDNF/SNAP91/NDRG4/FIT as the best performing marker panel. Pathway and gene ontology enrichment analyses revealed that genes associated with the nervous system were enriched in the set of best performing CRC-specific biomarkers. Conclusion In silico discovery analysis using TCGA-derived data yielded a novel DNA-methylation-based assay for the early detection of CRC, potentially improving current screening modalities. Additionally, nervous system-related pathways were enriched in the identified genes, indicating an epigenetic regulation of neuronal genes in CRC.

The European starling, Sturnus vulgaris, is an ecologically significant, globally invasive avian species that is also suffering from a major decline in its native range. Here, we present the genome assembly and long-read transcriptome of an Australian-sourced European starling (S. vulgaris vAU), and a second North American genome (S. vulgaris vNA), as complementary reference genomes for population genetic and evolutionary characterisation. S. vulgaris vAU combined 10x Genomics linked-reads, low-coverage Nanopore sequencing, and PacBio Iso-Seq full-length transcript scaffolding to generate a 1050 Mb assembly on 1,628 scaffolds (72.5 Mb scaffold N50). Species-specific transcript mapping and gene annotation revealed high structural and functional completeness (94.6% BUSCO completeness). Further scaffolding against the high-quality zebra finch (Taeniopygia guttata) genome assigned 98.6% of the assembly to 32 putative nuclear chromosome scaffolds. Rapid, recent advances in sequencing technologies and bioinformatics software have highlighted the need for evidence-based assessment of assembly decisions on a case-by-case basis. Using S. vulgaris vAU, we demonstrate how the multifunctional use of PacBio Iso-Seq transcript data and complementary homology-based annotation of sequential assembly steps (assessed using a new tool, SAAGA) can be used to assess, inform, and validate assembly workflow decisions. We also highlight some counter-intuitive behaviour in traditional BUSCO metrics, and present Buscomp, a complementary tool for assembly comparison designed to be robust to differences in assembly size and base-calling quality. Finally, we present a second starling assembly, S. vulgaris vNA, to facilitate comparative analysis and global genomic research on this ecologically important species.

Meiosis drives reciprocal genetic exchanges and produces gametes with halved chromosome number, which is important for the genetic diversity, plant viability, and ploidy consistency of flowering plants. Alterations in chromosome dynamics and/or cytokinesis during meiosis may lead to meiotic restitution and the formation of unreduced microspores. In this study, we isolated an Arabidopsis mutant male meiotic restitution 1 (mmr1), which produces a small subpopulation of diploid or polyploid pollen grains. Cytological analysis revealed that mmr1 produces dyads, triads, and monads indicative of male meiotic restitution. Both homologous chromosomes and sister chromatids in mmr1 are separated normally, but chromosome condensation at metaphase I is slightly affected. The mmr1 mutant displayed incomplete meiotic cytokinesis. Supportively, immunostaining of the microtubular cytoskeleton showed that the spindle organization at anaphase II and mini-phragmoplast formation at telophase II are aberrant. The causative mutation in mmr1 was mapped to chromosome 1 at the chromatin regulator Male Meiocyte Death 1 (MMD1/DUET) locus. mmr1 contains a C-to-T transition at the third exon of MMD1/DUET at the genomic position 2168 bp from the start codon, which causes an amino acid change G618D that locates in the conserved PHD-finger domain of histone binding proteins. The F1 progenies of mmr1 crossing with knockout mmd1/duet mutant exhibited same meiotic defects and similar meiotic restitution rate as mmr1. Taken together, we here report a hypomorphic mmd1/duet allele that typically shows defects in microtubule organization and cytokinesis.

Background Root knot nematodes (RKN) are plant parasitic nematodes causing major yield losses of widely consumed food crops such as rice (Oryza sativa). Because non-coding RNAs, including small interfering RNAs (siRNA), microRNAs (miRNAs) and long non-coding RNAs (lncRNAs), are key regulators of various plant processes, elucidating their regulation during this interaction may lead to new strategies to improve crop protection. In this study, we aimed to identify and characterize rice siRNAs, miRNAs and lncRNAs responsive to early infection with RKN Meloidogyne graminicola (Mg), based on sequencing of small RNA, degradome and total RNA libraries from rice gall tissues compared with uninfected root tissues. Results We found 425 lncRNAs, 3739 siRNAs and 16 miRNAs to be differentially expressed between both tissues, of which a subset was independently validated with RT-qPCR. Functional prediction of the lncRNAs indicates that a large part of their potential target genes code for serine/threonine protein kinases and transcription factors. Differentially expressed siRNAs have a predominant size of 24 nts, suggesting a role in DNA methylation. Differentially expressed miRNAs are generally downregulated and target transcription factors, which show reduced degradation according to the degradome data. Conclusions To our knowledge, this work is the first to focus on small and long non-coding RNAs in the interaction between rice and Mg, and provides an overview of rice non-coding RNAs with the potential to be used as a resource for the development of new crop protection strategies.

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.

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

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.

Background The HIV-1 proviral genome harbors multiple CpG islands (CpGIs), both in the promoter and intragenic regions. DNA methylation in the promoter region has been shown to be heavily involved in HIV-1 latency regulation in cultured cells. However, its exact role in proviral transcriptional regulation in infected individuals is poorly understood or characterized. Moreover, methylation at intragenic CpGIs has never been studied in depth. Results A large, well-characterized HIV-1 patient cohort (n = 72), consisting of 17 long-term non-progressors and 8 recent seroconverters (SRCV) without combination antiretroviral therapy (cART), 15 early cART-treated, and 32 late cART-treated patients, was analyzed using a next-generation bisulfite sequencing DNA methylation method. In general, we observed low level of promoter methylation and higher levels of intragenic methylation. Additionally, SRCV showed increased promoter methylation and decreased intragenic methylation compared with the other patient groups. This data indicates that increased intragenic methylation could be involved in proviral transcriptional regulation. Conclusions Contrasting in vitro studies, our results indicate that intragenic hypermethylation of HIV-1 proviral DNA is an underestimated factor in viral control in HIV-1-infected individuals, showing the importance of analyzing the complete proviral genome in future DNA methylation studies.

Background Sugar beet (Beta vulgaris subsp. vulgaris) is an economically important crop that provides nearly one third of the global sugar production. The beet cyst nematode (BCN), Heterodera schachtii, causes major yield losses in sugar beet and other crops worldwide. The most effective and economic approach to control this nematode is growing tolerant or resistant cultivars. To identify candidate genes involved in susceptibility and resistance, the transcriptome of sugar beet and BCN in compatible and incompatible interactions at two time points was studied using mRNA-seq. Results In the susceptible cultivar, most defense-related genes were induced at 4 dai while suppressed at 10 dai but in the resistant cultivar Nemakill, induction of genes involved in the plant defense response was observed at both time points. In the compatible interaction, alterations in phytohormone-related genes were detected. The effect of exogenous application of Methyl Jasmonate and ET-generator ethephon on susceptible plants was therefore investigated and the results revealed significant reduction in plant susceptibility. Genes putatively involved in the resistance of Nemakill were identified, such as genes involved in phenylpropanoid pathway and genes encoding CYSTM domain-containing proteins, F-box proteins, chitinase, galactono-1,4-lactone dehydrogenase and CASP-like protein. Also, the transcriptome of the BCN was analyzed in infected root samples and several novel potential nematode effector genes were found. Conclusions Our data provides detailed insights into the plant and nematode transcriptional changes occurring during compatible and incompatible interactions between sugar beet and BCN. Many important genes playing potential roles in susceptibility or resistance of sugar beet against BCN, as well as some BCN effectors with a potential role as avr proteins were identified. In addition, our findings indicate the effective role of jasmonate and ethylene in enhancing sugar beet defense response against BCN. This research provides new molecular insights into the plant-nematode interactions that can be used to design novel management strategies against BCN.

Genomic imprinting plays an important role in growth and development. Loss of imprinting (LOI) has been found in cancer, yet systematic studies are impeded by data-analytical challenges. We developed a methodology to detect monoallelically expressed loci without requiring genotyping data, and applied it on The Cancer Genome Atlas (TCGA, discovery) and Genotype-Tissue expression project (GTEx, validation) breast tissue RNA-seq data. Here, we report the identification of 30 putatively imprinted genes in breast. In breast cancer (TCGA), HM13 is featured by LOI and expression upregulation, which is linked to DNA demethylation. Other imprinted genes typically demonstrate lower expression in cancer, often associated with copy number variation and aberrant DNA methylation. Downregulation in cancer frequently leads to higher relative expression of the (imperfectly) silenced allele, yet this is not considered canonical LOI given the lack of (absolute) re-expression. In summary, our novel methodology highlights the massive deregulation of imprinting in breast cancer., In growth and development genomic imprinting is important in regulating gene expression. Here, the authors study loss of imprinting (LOI) in cancer, developing a mixture model to detect monoallelically expressed loci without genotyping data - applying this novel methodology to TCGA breast cancer data they find massive deregulation of imprinting.

Objective: To investigate the effect of a swimming break on airway hyperresponsiveness (AHR) and to evaluate perception of bronchoconstriction-related symptoms after methacholine challenge testing (MCT). Design: Observational, controlled study. Setting: University department. Participants: Overall, 25 healthy, elite, competitive swimmers aged ≥14 years from the local Olympic training center and 25 healthy control subjects. Independent Variable: Duration of water abstinence. Main Outcome Measures: Primary outcome measures were a positive MCT with a provocation dose (PD) and the perception score of bronchoconstriction-related symptoms at visit 1 (V1, training period) and visit 2 (V2, after swimming break). Results: In the study, 13 swimmers (52%) and 10 control subjects (40%) showed AHR. The PD did not differ significantly between V1 and V2 in the swimmers (P = 0.81) or in the control subjects (P = 0.74). No correlation of fraction of exhaled nitric oxide with the PD could be established in both the groups (swimmers: P = 0.97; controls: P = 0.99). The majority of swimmers with AHR were asymptomatic, and only minimal differences in perception of bronchoconstriction-related symptoms between swimmers and control subjects were observed (P = 0.23). Conclusions: A swimming break of 25 ± 8 days does not seem sufficient to significantly reduce AHR in elite swimmers. Thus, relevant and systematic effects of short-term water abstinence on AHR seem unlikely. Because the majority of swimmers remain asymptomatic, AHR may escape routine screening examinations. The impact of AHR on athletic performance and the long-term clinical consequences remain to be clarified.

The human transcriptome consists of various RNA biotypes including multiple types of non-coding RNAs (ncRNAs). Current ncRNA compendia remain incomplete partially because they are almost exclusively derived from the interrogation of small- and polyadenylated RNAs. Here, we present a more comprehensive atlas of the human transcriptome that is derived from matching polyA-, total-, and small-RNA profiles of a heterogeneous collection of nearly 300 human tissues and cell lines. We report on thousands of novel RNA species across all major RNA biotypes, including a hitherto poorly-cataloged class of non-polyadenylated single-exon long non-coding RNAs. In addition, we exploit intron abundance estimates from total RNA-sequencing to test and verify functional regulation by novel non-coding RNAs. Our study represents a substantial expansion of the current catalogue of human ncRNAs and their regulatory interactions. All data, analyses, and results are available in the R2 web portal and serve as a basis to further explore RNA biology and function.

The recent growth in the number of publicly available cancer omics databases has been accompanied by the development of various tools that allow researchers to visually explore these data. In 2015, we built MEXPRESS, an online tool for the integration and visualization of gene expression, DNA methylation and clinical data from The Cancer Genome Atlas (TCGA), a large collection of publicly available multi-omics cancer data. MEXPRESS addresses the need for an easy-to-use, interactive application that allows researchers to identify dysregulated genes and their clinical relevance in cancer. Furthermore, while other tools typically do not support integrated visualization of expression and DNA methylation data in combination with the precise genomic location of the methylation, MEXPRESS is unique in how it depicts these diverse data types together. Motivated by the large number of users MEXPRESS has managed to attract over the past 3 years and the recent migration of all TCGA data to a new data portal, we developed a new version of MEXPRESS (https://mexpress.be). It contains the latest TCGA data, additional types of omics and clinical data and extra functionality, allowing users to explore mechanisms of gene dysregulation beyond expression and DNA methylation., SCOPUS: ar.j, info:eu-repo/semantics/published

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.

Telomere length is a prognostic biomarker for aging diseases. As it is unknown whether diet plays a role in these associations, we aimed to assess the impact of diet on telomere length. Moreover, given that telomere length is modulated by oxidative stress and inflammation, an additional goal was to evaluate whether the latter may mediate possible telomere – diet associations. Southern blot measured leukocyte telomere length and food frequency questionnaire data were compared for 2509 apparently healthy men and women (~35 to 55 years) from the Asklepios population. No significant associations were found between telomere length and overall dietary characteristics, such as dietary diversity, quality, equilibrium, and the dietary inflammatory index. Exploratory analysis of individual dietary variables revealed that a higher daily intake of deep fried potato products was associated with shorter telomeres (P = 0.002, 151 bp per 100 g/day), also in both sexes separately. Deep fried potato product consumption was also significantly associated with C-reactive protein (P = 0.032) and uric acid (P = 0.042), but not other inflammation and oxidative stress markers. These results suggest an at most limited association between overall dietary patterns and telomere length in the general population. Nevertheless, the association between telomere length and deep fried potato product intake warrants additional research.

Genomic imprinting, the parent-of-origin specific monoallelic expression of genes, plays an important role in growth and development. Loss of imprinting of individual genes has been found in varying cancers, yet data-analytical challenges have impeded systematic studies so far. We developed a mixture distribution model to detect monoallelically expressed loci in a genome-wide manner without the need for genotyping data, and applied the methodology on TCGA breast tissue RNA-seq data. We identified 35 putatively imprinted genes in healthy breast. In breast cancer however, HM13 was featured by significant loss of imprinting and expression upregulation, which could be linked to DNA demethylation. Other imprinted genes (25 out of 35) demonstrated consistent expression downregulation in breast cancer, which often correlated with loss of imprinting. A breast imprinted gene network, deregulated in cancer, might hence be present. In summary, our novel methodology highlights the massive deregulation of imprinting in breast cancer.

Telomere length shortens with aging, and short telomeres have been linked to a wide variety of pathologies. Previous studies suggested a discrepancy in age-associated telomere shortening rate estimated by cross-sectional studies versus the rate measured in longitudinal studies, indicating a potential bias in cross-sectional estimates. Intergenerational changes in initial telomere length, such as that predicted by the previously described effect of a father’s age at birth of his offspring (FAB), could explain the discrepancy in shortening rate measurements. We evaluated whether changes occur in initial telomere length over multiple generations in three large datasets and identified paternal birth year (PBY) as a variable that reconciles the difference between longitudinal and cross-sectional measurements. We also clarify the association between FAB and offspring telomere length, demonstrating that this effect is substantially larger than reported in the past. These results indicate the presence of a downward secular trend in telomere length at birth over generational time with potential public health implications.

This work was supported by the Chief Scientist Office, Scotland (ETM37; GDS, DJH), Cancer Research UK (Experimental Cancer Medicine Centre; TP, London, DJH, Edinburgh), Medical Research Council (AL, DJH), Royal College of Surgeons of Edinburgh Robertson Trust (AL), Melville Trust AL), Renal Cancer Research Fund (GDS), Kidney Cancer Scotland (GDS) and an educational grant from Pfizer (TP). CVN and TDM were funded by Ghent University Multidisciplinary Research Partnership 'Bioinformatics: from nucleotides to networks' Background: Sample collection for clinical epigenetics research often occurs at the time of surgical excision of a diseased organ. However, this approach may compromise the epigenetic profile under study. The use of different sampling procedures during a study can often not be avoided, but may in theory lead to biased results. The effect of tissue sampling approach on DNA methylation is studied here using clear cell renal cell cancer (ccRCC) as a model. A comparison of the DNA methylation profiles between vascularised tumour biopsy samples and subsequent devascularised nephrectomy samples obtained from the same two individuals (total of 6 samples per individual) was undertaken. Validation of the results was performed using biopsy and nephrectomy samples obtained from 14 patients included in a ccRCC clinical trial (SuMR; ClinicalTrials.gov identifier: NCT01024205). Findings: Using MBD2 sequencing, the methylome was analysed for all samples and differential methylome regions were retrieved. The results, from the test set, show six differentially methylated genes, of which four were clearly linked to ischaemia or hypoxia (REXO1L1, TLR4, hsa-mir-1299, and ANKRD2). To validate these findings, it was evaluated whether these loci were also featured by differential methylation in the clinical trial cohort with a similar experimental design to the test set. Three of the six genes are again significantly differentially methylated, showing an overall clear impact of renal artery clamping on DNA methylation. Conclusions: Renal artery ligation modulates the ccRCC methylome, impacting methylation of ischaemia and hypoxia associated genes. Results from devascularised surgical resection specimens do not accurately reflect findings from tumour biopsies. Postprint

Purpose: The currently used prognostic models for patients with nonmetastatic clear cell renal cell carcinoma (ccRCC) are based on clinicopathologic features and might be improved by adding molecular markers. Epigenetic alterations occur frequently in ccRCC and are promising biomarkers. The aim of this study is to identify prognostic promoter methylation markers for ccRCC. Experimental Design: We integrated data generated by massive parallel sequencing of methyl-binding domain enriched DNA and microarray-based RNA expression profiling of 5-aza-2′-deoxycytidine–treated ccRCC cell lines to comprehensively characterize the ccRCC methylome. A selection of the identified methylation markers was evaluated in two independent series of primary ccRCC (n = 150 and n = 185) by methylation-specific PCR. Kaplan–Meier curves and log-rank tests were used to estimate cause-specific survival. HRs and corresponding 95% confidence intervals (CI) were assessed using Cox proportional hazard models. To assess the predictive capacity and fit of models combining several methylation markers, HarrellC statistic and the Akaike Information Criterion were used. Results: We identified four methylation markers, that is, GREM1, NEURL, LAD1, and NEFH, that individually predicted prognosis of patients with ccRCC. The four markers combined were associated with poorer survival in two independent patient series (HR, 3.64; 95% CI, 1.02–13.00 and HR, 7.54; 95% CI, 2.68–21.19). These findings were confirmed in a third series of ccRCC cases from The Cancer Genome Atlas (HR, 3.60; 95% CI, 2.02–6.40). Conclusions: A four-gene promoter methylation marker panel consisting of GREM1, NEURL, LAD1, and NEFH predicts outcome of patients with ccRCC and might be used to improve current prognostic models. Clin Cancer Res; 23(8); 2006–18. ©2016 AACR.

Magnaporthe oryzae (rice blast) and the root-knot nematode Meloidogyne graminicola are causing two of the most important pathogenic diseases jeopardizing rice production. Here, we show that root-knot nematode infestation on rice roots leads to important above-ground changes in plant immunity gene expression, which is correlated with significantly enhanced susceptibility to blast disease. A detailed metabolic analysis of oxidative stress responses and hormonal balances demonstrates that the above-ground tissues have a disturbed oxidative stress level, with accumulation of H2O2, as well as hormonal disturbances. Moreover, double infection experiments on an oxidative stress mutant and an auxin-deficient rice line indicate that the accumulation of auxin in the above-ground tissue is at least partly responsible for the blast-promoting effect of root-knot nematode infection.

Epigenetics refers to the collection of heritable features that modulate the genome-environment interaction without being encoded in the actual DNA sequence. While being mitotically and sometimes even meiotically transmitted, epigenetic traits often demonstrate extensive flexibility. This allows cells to acquire diverse gene expression patterns during differentiation, but also to adapt to a changing environment. However, epigenetic alterations are not always beneficial to the organism, as they are, for example, frequently identified in human diseases such as cancer. Accurate and cost-efficient genome-scale profiling of epigenetic features is thus of major importance to pinpoint these "epimutations," for example, to monitor the epigenetic impact of environmental exposure. Over the last decade, the field of epigenetics has been revolutionized by several innovative "epigenomics" technologies exactly addressing this need. In this review, we discuss and compare widely used next-generation methods to assess DNA methylation and hydroxymethylation, noncoding RNA expression, histone modifications, and nucleosome positioning. Although recent methods are typically based on "second-generation" sequencing, we also pay attention to still commonly used array- and PCR-based methods, and look forward to the additional advantages of single-molecule sequencing. As the current bottleneck in epigenomics research is the analysis rather than generation of data, the basic difficulties and problem-solving strategies regarding data preprocessing and statistical analysis are introduced for the different technologies. Finally, we also consider the complications associated with epigenomic studies of species with yet unsequenced genomes and possible solutions.

† Background and Aims The combination of clonality and a mating system promoting outcrossing is considered advantageous because outcrossing avoids the fitness costs of selfing within clones (geitonogamy) while clonality assures local persistence and increases floral display. The spatial spread of genetically identical plants (ramets) may, however, also decrease paternal diversity (the number of sires fertilizing a given dam) and fertility, particularly towards the centre of large clumped clones. This study aimed to quantify the impact of extensive clonal growth on fine-scale paternity patterns in a population of the allogamous Convallaria majalis. † Methods A full analysis of paternity was performed by genotyping all flowering individuals and all viable seeds produced during a single season using AFLP. Mating patterns were examined and the spatial position of ramets was related to the extent of multiple paternity, fruiting success and seed production. † Key Results The overall outcrossing rate was high (91 %) and pollen flow into the population was considerable (27 %). Despite extensive clonal growth, multiple paternity was relatively common (the fraction of siblings sharing the same father was 0.53 within ramets). The diversity of offspring collected from reproductive ramets surrounded by genetically identical inflorescences was as high as among offspring collected from ramets surrounded by distinct genets. There was no significant relationship between the similarity of the pollen load received by two ramets and the distance between them. Neither the distance of ramets with respect to distinct genets nor the distance to the genet centre significantly affected fruiting success or seed production. † Conclusions Random mating and considerable pollen inflow most probably implied that pollen dispersal distances were sufficiently high to mitigate local mate scarcity despite extensive clonal spread. The data provide no evidence for the intrusion of clonal growth on fine-scale plant mating patterns.

Although the CpG island methylator phenotype (CIMP) was first identified and has been most extensively studied in colorectal cancer, the term “CIMP” has been repeatedly used over the past decade to describe CpG island promoter methylation in other tumor types, including bladder, breast, endometrial, gastric, glioblastoma (gliomas), hepatocellular, lung, ovarian, pancreatic, renal cell, and prostate cancers, as well as for leukemia, melanoma, duodenal adenocarninomas, adrenocortical carcinomas, and neuroblastomas. CIMP has been reported to be useful for predicting prognosis and response to treatment in a variety of tumor types, but it remains unclear whether or not CIMP is a universal phenomenon across human neoplasia or if there should be cancer-specific definitions of the phenotype. Recently, it was shown that somatic isocitrate dehydrogenase-1 (IDH1) mutations, frequently observed in gliomas, establish CIMP in primary human astrocytes by remodeling the methylome. Interestingly, somatic IDH1 and IDH2 mutations, and loss-of-function mutations in ten-eleven translocation (TET) methylcytosine dioxygenase-2 (TET2) associated with a hypermethylation phenotype, are also found in multiple enchondromas of patients with Ollier disease and Mafucci syndrome, and leukemia, respectively. These data provide the first clues for the elucidation of a molecular basis for CIMP. Although CIMP appears as a phenomenon that occurs in various cancer types, the definition is poorly defined and differs for each tumor. The current perspective discusses the use of the term CIMP in cancer, its significance in clinical practice, and future directions that may aid in identifying the true cause and definition of CIMP in different forms of human neoplasia. Cancer Res; 73(19); 5858–68. ©2013 AACR.

We agree with Heidbuchel and colleagues that our findings do not clearly refute the concept of exercise-induced right ventricular (RV) damage because no prospective longitudinal study has been performed so far. However, our finding of marked RV remodeling without functional RV impairment or scarring challenges the notion of an exercise-induced arrhythmogenic right ventricular cardiomyopathy. First, we studied an exclusive cohort of world-class master endurance athletes with a training history of 29±8 years. Such athletes are not as easily available as ambitious recreational athletes, who would not demonstrate an extensive cardiac remodeling as our studied athletes. Our athletes demonstrated clear structural exercise-induced cardiac remodeling with a mild RV chamber size predominance of about 6% (left ventricular end-diastolic volume/right ventricular end-diastolic volume ratio

Background: Ipilimumab (Ipi) improves the survival of advanced melanoma patients with an incremental long-term benefit in 10-15 % of patients. A tumor signature that correlates with this survival benefit could help optimizing individualized treatment strategies. Methods: Freshly frozen melanoma metastases were collected from patients treated with either Ipi alone (n: 7) or Ipi combined with a dendritic cell vaccine (TriMixDC-MEL) (n: 11). Samples were profiled by immunohistochemistry (IHC), whole transcriptome (RNA-seq) and methyl-DNA sequencing (MBD-seq). Results: Patients were divided in two groups according to clinical evolution: durable benefit (DB; 5 patients) and no clinical benefit (NB; 13 patients). 20 metastases were profiled by IHC and 12 were profiled by RNA- and MBD-seq. 325 genes were identified as differentially expressed between DB and NB. Many of these genes reflected a humoral and cellular immune response. MBD-seq revealed differences between DB and NB patients in the methylation of genes linked to nervous system development and neuron differentiation. DB tumors were more infiltrated by CD8+ and PD-L1+ cells than NB tumors. B cells (CD20+) and macrophages (CD163+) co-localized with T cells. Focal loss of HLA class I and TAP-1 expression was observed in several NB samples. Conclusion: Combined analyses of melanoma metastases with IHC, gene expression and methylation profiling can potentially identify durable responders to Ipi-based immunotherapy., SCOPUS: ar.j, info:eu-repo/semantics/published

Learning and memory formation are known to require dynamic CpG (de)methylation and gene expression changes. Here, we aimed at establishing a genome-wide DNA methylation map of the zebra finch genome, a model organism in neuroscience, as well as identifying putatively epigenetically regulated genes. RNA- and MethylCap-seq experiments were performed on two zebra finch cell lines in presence or absence of 5-aza-2′-deoxycytidine induced demethylation. First, the MethylCap-seq methodology was validated in zebra finch by comparison with RRBS-generated data. To assess the influence of (variable) methylation on gene expression, RNA-seq experiments were performed as well. Comparison of RNA-seq and MethylCap-seq results showed that at least 357 of the 3,457 AZA-upregulated genes are putatively regulated by methylation in the promoter region, for which a pathway analysis showed remarkable enrichment for neurological networks. A subset of genes was validated using Exon Arrays, quantitative RT-PCR and CpG pyrosequencing on bisulfite-treated samples. To our knowledge, this study provides the first genome-wide DNA methylation map of the zebra finch genome as well as a comprehensive set of genes of which transcription is under putative methylation control.

Objective— Shorter telomere length is associated with the occurrence of cardiovascular events, but the question of causality is complicated by the intertwined effects of inheritance, aging, and lifestyle factors on both telomere length and cardiovascular disease (CVD). Some studies indicated that healthy offspring of coronary artery disease patients exhibited shorter telomeres than subjects without a family history. Importantly, this result would imply that inheritance of shorter telomeres is a primary abnormality associated with an increased risk of CVD, the so-called Telomere Hypothesis of CVD. Therefore, we aimed at further validating the latter results in the large, population-representative Asklepios Study. Methods and Results— Peripheral blood leukocyte telomere length was measured using telomere restriction fragment analysis in the young to middle-aged (≈35–55 years old) Asklepios study population, free from overt CVD, and could be successfully combined with data from the Asklepios Family History Database for 2136 subjects. No shorter telomere length could be found in healthy subjects with a family history of CVD compared with those without. Conclusion— These findings cast serious doubt on the hypothesis that telomere length is shorter in families with an increased risk of CVD and do not support the Telomere Hypothesis of CVD.

Rice is one of the most important staple crops worldwide, but its yield is compromised by different pathogens, including plant-parasitic nematodes. In this study we have characterized specific and general responses of rice (Oryza sativa) roots challenged with two endoparasitic nematodes with very different modes of action. Local transcriptional changes in rice roots upon root knot (Meloidogyne graminicola) and root rot nematode (RRN, Hirschmanniella oryzae) infection were studied at two time points (3 and 7 d after infection, dai), using mRNA-seq. Our results confirm that root knot nematodes (RKNs), which feed as sedentary endoparasites, stimulate metabolic pathways in the root, and enhance nutrient transport towards the induced root gall. The migratory RRNs, on the other hand, induce programmed cell death and oxidative stress, and obstruct the normal metabolic activity of the root. While RRN infection causes up-regulation of biotic stress-related genes early in the infection, the sedentary RKNs suppress the local defense pathways (e.g. salicylic acid and ethylene pathways). Interestingly, hormone pathways mainly involved in plant development were strongly induced (gibberellin) or repressed (cytokinin) at 3 dai. These results uncover previously unrecognized nematode-induced expression profiles related to their specific infection strategy.

Background— Carotid-femoral pulse wave velocity (CFPWV) is a heritable measure of aortic stiffness that is strongly associated with increased risk for major cardiovascular disease events. Methods and Results— We conducted a meta-analysis of genome-wide association data in 9 community-based European ancestry cohorts consisting of 20 634 participants. Results were replicated in 2 additional European ancestry cohorts involving 5306 participants. Based on a preliminary analysis of 6 cohorts, we identified a locus on chromosome 14 in the 3′- BCL11B gene desert that is associated with CFPWV (rs7152623, minor allele frequency=0.42, β=−0.075±0.012 SD/allele, P =2.8×10 −10 ; replication β=−0.086±0.020 SD/allele, P =1.4×10 −6 ). Combined results for rs7152623 from 11 cohorts gave β=−0.076±0.010 SD/allele, P =3.1×10 −15 . The association persisted when adjusted for mean arterial pressure (β=−0.060±0.009 SD/allele, P =1.0×10 −11 ). Results were consistent in younger (P =2.3×10 −9 ) and older (9 cohorts, n=12 026, β=−0.061±0.014 SD/allele, P =9.4×10 −6 ) participants. In separate meta-analyses, the locus was associated with increased risk for coronary artery disease (hazard ratio=1.05; confidence interval=1.02–1.08; P =0.0013) and heart failure (hazard ratio=1.10, CI=1.03–1.16, P =0.004). Conclusions— Common genetic variation in a locus in the BCL11B gene desert that is thought to harbor 1 or more gene enhancers is associated with higher CFPWV and increased risk for cardiovascular disease. Elucidation of the role this novel locus plays in aortic stiffness may facilitate development of therapeutic interventions that limit aortic stiffening and related cardiovascular disease events.

Despite the major physiological dissimilarities between mature root regions and their tips, differences in their gene expression profiles remain largely unexplored. In this research, the transcriptome of rice (Oryza sativa L. subsp. japonica) mature root tissue and root tips was monitored using mRNA-Seq at two time points. Almost 50 million 76 bp reads were mapped onto the rice genome sequence, expression patterns for different tissues and time points were investigated, and at least 1106 novel transcriptionally active regions (nTARs) expressed in rice root tissue were detected. More than 30 000 genes were found to be expressed in rice roots, among which were 1761 root-enriched and 306 tip-enriched transcripts. Mature root tissue appears to respond more strongly to external stimuli than tips, showing a higher expression of, for instance, auxin-responsive and abscisic acid-responsive genes, as well as the phenylpropanoid pathway and photosynthesis upon light. The root tip-enriched transcripts are mainly involved in mitochondrial electron transport, organelle development, secondary metabolism, DNA replication and metabolism, translation, and cellular component organization. During root maturation, genes involved in cell wall biosynthesis and modification, response to oxidative stress, and secondary metabolism were activated. For some nTARs, a potential role in root development can be put forward based on homology to genes involved in CLAVATA signalling, cell cycle regulators, and hormone signalling. A subset of differentially expressed genes and novel transcripts was confirmed using (quantitative) reverse transcription-PCR. These results uncover previously unrecognized tissue-specific expression profiles and provide an interesting starting point to study the different regulation of transcribed regions of these tissues.

See related article, pp 420–425 The discovery that critical telomere shortening initiates replicative senescence triggered a vast body of epidemiological studies exploring its implications for human aging. These studies demonstrated that telomeres—typically measured in the DNA isolated from blood cells—clearly shorten with age. On average, shorter telomeres were found in subjects with atherosclerosis and independently predicted atherosclerotic cardiovascular disease (ACVD)–associated mortality. Moreover, Mendelian randomization studies point to a causal role for shorter telomere length in ACVD, lending support to the hypothesis that the protective antitumor mechanism elicited by critical telomere attrition may act at the expense of unsuccessful cardiovascular aging.1 Yet, consensus is lacking on the mechanism underlying the presumed associations between shorter telomeres, atherosclerosis, and ACVD. This problem was addressed in the epidemiological research by Toupance et al,2 published in the current issue of Hypertension . In a longitudinal study, both leukocyte telomere length and carotid atherosclerosis were quantified in 154 subjects at baseline and after a 9.5-year follow-up, and their interrelation was evaluated. More specifically, 2 competitive hypotheses on the role of telomeres as primary risk factor for ACVD were scrutinized: (1) shorter (inherited) telomeres versus (2) higher telomere attrition rates (Figure). The study provides evidence supporting the first hypothesis. In this commentary, we focus on embedding these results in the current state of the art and highlight the potential of inferring insight from epidemiological surveys. Figure. Schematic overview depicting either ( A ) shorter inherited telomere length (TL), or ( B ) accelerated telomere attrition (ΔTL), as the …

In small populations of plant species with separate sexes, it can be expected that besides the local environment also stochastic events influence population sex ratios. Biased sex ratios may in turn negatively affect genetic diversity due to increased genetic drift and, in clonal plants, due to reduced sexual reproductive output. Empirical evidence for these processes is scarce, however. We investigated the pattern of sex ratio variation and the distribution of genetic variation of the dioecious clonal forest herb Mercurialis perennis using AFLP markers. Analysis of molecular variance indicated a pronounced genetic structure. Overall within-population genetic diversity was moderate and local sex ratios were slightly male biased. The proportion of male to female plants in large populations slightly increased with increasing light penetration to the herb layer. Small populations, on the contrary, displayed high variability in sex ratios, unrelated to the local light environment. Genotypic diversity decreased with more male-biased sex ratios. We conclude that stochastic events related to small population size and the local forest environment, related to canopy closure, affect the proportion of female plants and indirectly influence local genotypic diversity, likely through the degree of sexual reproduction. This is one of the first studies to report a clear association between gender proportions and genetic diversity of a dioecious plant species in a fairly large survey.

Although telomere biology was revealed to play an important role in several hematopoietic disorders, its impact on the age-dependent dynamics of regular hematopoiesis is poorly understood. In vitro results suggest that particularly the erythropoietic capacity might be limited by critically short telomere length (TL). However, it remains unclear whether TL also affects erythropoiesis in healthy individuals in vivo. Therefore, we analyzed the associations between relevant hematopoietic parameters and peripheral blood leukocyte TL in the apparently healthy Asklepios study population, aged approximately 35-55 years (N2500). Our data indicate a clear positive, age and paternal age at birth adjusted, correlation between TL and red blood cell count, both in men (p0.001) and women (p = 0.011). This association was particularly significant in the older segment of the population (45 years old, both sexes: p = 0.003) and in younger men (p = 0.013), but not in younger women (p = 0.521). Further adjustment for known determinants in a general linear model revealed that peripheral blood leukocyte TL is most probably an independent predictor of red blood cell count (p0.001), suggesting that critical telomere shortening might also limit erythropoiesis in vivo. While negligible in a middle-aged population, the clinical consequences might be important in the elderly (e.g. in anemia of chronic disease). Further studies are required to confirm the impact of our results.

Epigenetics, and more specifically DNA methylation is a fast evolving research area. In almost every cancer type, each month new publications confirm the differentiated regulation of specific genes due to methylation and mention the discovery of novel methylation markers. Therefore, it would be extremely useful to have an annotated, reviewed, sorted and summarized overview of all available data. PubMeth is a cancer methylation database that includes genes that are reported to be methylated in various cancer types. A query can be based either on genes (to check in which cancer types the genes are reported as being methylated) or on cancer types (which genes are reported to be methylated in the cancer (sub) types of interest). The database is freely accessible at http://www.pubmeth.org. PubMeth is based on text-mining of Medline/PubMed abstracts, combined with manual reading and annotation of preselected abstracts. The text-mining approach results in increased speed and selectivity (as for instance many different aliases of a gene are searched at once), while the manual screening significantly raises the specificity and quality of the database. The summarized overview of the results is very useful in case more genes or cancer types are searched at the same time.

Most next generation sequencing experiments generate more data than is usable for the experimental set up. For example, methyl-CpG binding domain (MBD) affinity purification based sequencing is often used for DNA-methylation profiling, but up to 30% of the sequenced fragments cannot be mapped uniquely to the reference genome. Here we present and evaluate a methodology for the identification of viruses in these otherwise unused paired-end MBD-seq data. Viral detection is accomplished by mapping non-reference alignable reads to a comprehensive set of viral genomes. As viruses play an important role in epigenetics and cancer development, 92 (pre)malignant and benign samples, originating from two different collections of cervical samples and related cell lines, were used in this study. These samples include primary carcinomas (n = 22), low- and high-grade cervical intraepithelial neoplasia (CIN1 and CIN2/3 - n = 2/n = 30) and normal tissue (n = 20), as well as control samples (n = 17). Viruses that were detected include phages, adenoviruses, herpesviridae and HPV. HPV, which causes virtually all cervical cancers, was identified in 95% of the carcinomas, 100% of the CIN2/3 samples, both CIN1 samples and in 55% of the normal samples. Comparing the amount of mapped fragments on HPV for each HPV-infected sample yielded a significant difference between normal samples and carcinomas or CIN2/3 samples (adjusted p-values resp.

Two cost-efficient genome-scale methodologies to assess DNA-methylation are MethylCap-seq and Illumina’s Infinium HumanMethylation450 BeadChips (HM450). Objective information regarding the best-suited methodology for a specific research question is scant. Therefore, we performed a large-scale evaluation on a set of 70 brain tissue samples, i.e. 65 glioblastoma and 5 non-tumoral tissues. As MethylCap-seq coverages were limited, we focused on the inherent capacity of the methodology to detect methylated loci rather than a quantitative analysis. MethylCap-seq and HM450 data were dichotomized and performances were compared using a gold standard free Bayesian modelling procedure. While conditional specificity was adequate for both approaches, conditional sensitivity was systematically higher for HM450. In addition, genome-wide characteristics were compared, revealing that HM450 probes identified substantially fewer regions compared to MethylCap-seq. Although results indicated that the latter method can detect more potentially relevant DNA-methylation, this did not translate into the discovery of more differentially methylated loci between tumours and controls compared to HM450. Our results therefore indicate that both methodologies are complementary, with a higher sensitivity for HM450 and a far larger genome-wide coverage for MethylCap-seq, but also that a more comprehensive character does not automatically imply more significant results in biomarker studies.

The natural ends of linear chromosomes, the telomeres, recruit specific proteins in the formation of protective caps that preserve the integrity of the genome. Unprotected chromosomes induce DNA damage checkpoint cascades and ultimately lead to senescence both in mouse and man in a p53 dependent manner and initial telomere length setting therefore determines the proliferative capacity of each cell. Yet, only little information is available on telomere biology during embryonic development. We have previously shown that the p53 gene plays a crucial role in the development of malformations (exencephaly, gastroschisis, polydactyly, cleft palate and dwarfism) in control and irradiated mouse embryos. Here, we investigated telomere biology and the outcome of radiation exposure in wild type (p53+/+) and p53-mutant (p53+/-- and--/--) C57BL mouse foetuses irradiated at three different developmental stages. We show that telomeres are significantly shorter in malformed foetuses as compared to normal counterparts. In addition, our results indicate that the observed telomere attrition is primarily associated with p53-deficiency but is also modulated by irradiation, more specifically during the gastrulation and organogenesis stages. In conclusion, we formulate a hypothesis in which telomere shortening is linked to the absence of p53 in mouse foetuses and that when, in the presence of shorter telomeres, these foetuses are irradiated, the chance for the occurrence of developmental defects increases substantially.

Hirschmanniella oryzae is the most common plant-parasitic nematode in flooded rice cultivation systems. These migratory animals penetrate the plant roots and feed on the root cells, creating large cavities, extensive root necrosis and rotting. The objective of this study was to investigate the systemic response of the rice plant upon root infection by this nematode. RNA sequencing was applied on the above-ground parts of the rice plants at 3 and 7 days post inoculation. The data revealed significant modifications in the primary metabolism of the plant shoot, with a general suppression of for instance chlorophyll biosynthesis, the brassinosteroid pathway, and amino acid production. In the secondary metabolism, we detected a repression of the isoprenoid and shikimate pathways. These molecular changes can have dramatic consequences for the growth and yield of the rice plants, and could potentially change their susceptibility to above-ground pathogens and pests.

Telomere length (TL) has been associated with aging and mortality, but individual differences are also influenced by genetic factors, with previous studies reporting heritability estimates ranging from 34 to 82%. Here we investigate the heritability, mode of inheritance and the influence of parental age at birth on TL in six large, independent cohort studies with a total of 19,713 participants. The meta-analysis estimate of TL heritability was 0.70 (95% CI 0.64-0.76) and is based on a pattern of results that is highly similar for twins and other family members. We observed a stronger mother-offspring (r=0.42; P-value=3.60 × 10(-61)) than father-offspring correlation (r=0.33; P-value=7.01 × 10(-5)), and a significant positive association with paternal age at offspring birth (β=0.005; P-value=7.01 × 10(-5)). Interestingly, a significant and quite substantial correlation in TL between spouses (r=0.25; P-value=2.82 × 10(-30)) was seen, which appeared stronger in older spouse pairs (mean age ≥55 years; r=0.31; P-value=4.27 × 10(-23)) than in younger pairs (mean age55 years; r=0.20; P-value=3.24 × 10(-10)). In summary, we find a high and very consistent heritability estimate for TL, evidence for a maternal inheritance component and a positive association with paternal age.

Monoallelic gene expression is typically initiated early in the development of an organism. Dysregulation of monoallelic gene expression has already been linked to several non-Mendelian inherited genetic disorders. In humans, DNA-methylation is deemed to be an important regulator of monoallelic gene expression, but only few examples are known. One important reason is that current, cost-affordable truly genome-wide methods to assess DNA-methylation are based on sequencing post enrichment. Here, we present a new methodology that combines methylomic data from MethylCap-seq with associated SNP profiles to identify monoallelically methylated loci. Using the Hardy-Weinberg theorem for each SNP locus, it could be established whether the observed frequency of samples featured by biallelic methylation was lower than randomly expected. Applied on 334 MethylCap-seq samples of very diverse origin, this resulted in the identification of 80 genomic regions featured by monoallelic DNA-methylation. Of these 80 loci, 49 are located in genic regions of which 25 have already been linked to imprinting. Further analysis revealed statistically significant enrichment of these loci in promoter regions, further establishing the relevance and usefulness of the method. Additional validation of the found loci was done using 14 whole-genome bisulfite sequencing data sets. Importantly, the developed approach can be easily applied to other enrichment-based sequencing technologies, such as the ChIP-seq-based identification of monoallelic histone modifications.