Genome-wide methylome stability and parental effects in the worldwide distributed Lombardy poplar.

BACKGROUND: Despite the increasing number of epigenomic studies in plants, little is known about the forces that shape the methylome in long-lived woody perennials. The Lombardy poplar offers an ideal opportunity to investigate the impact of the individual environmental history of trees on the methylome. RESULTS: We present the results of three interconnected experiments on Lombardy poplar. In the first experiment, we investigated methylome variability during a growing season and across vegetatively reproduced generations. We found that ramets collected over Europe and raised in common conditions have stable methylomes in symmetrical CG-contexts. In contrast, seasonal dynamics occurred in methylation patterns in CHH context. In the second experiment, we investigated whether methylome patterns of plants grown in a non-parental environment correlate with the parental climate. We did not observe a biological relevant pattern that significantly correlates with the parental climate. Finally, we investigated whether the parental environment has persistent carry-over effects on the vegetative offspring's phenotype. We combined new bud set observations of three consecutive growing seasons with former published bud set data. Using a linear mixed effects analysis, we found a statistically significant but weak short-term, parental carry-over effect on the timing of bud set. However, this effect was negligible compared to the direct effects of the offspring environment. CONCLUSIONS: Genome-wide cytosine methylation patterns in symmetrical CG-context are stable in Lombardy poplar and appear to be mainly the result of random processes. In this widespread poplar clone, methylation patterns in CG-context can be used as biomarkers to infer a common ancestor and thus to investigate the recent environmental history of a specific Lombardy poplar. The Lombardy poplar shows high phenotypic plasticity in a novel environment which enabled this clonal tree to adapt and survive all over the temperate regions of the world.

The multi-omics single-cell landscape of sinus mucosa in uncontrolled severe chronic rhinosinusitis with nasal polyps.

Uncontrolled severe chronic rhinosinusitis with nasal polyps (CRSwNP) is associated with elevated levels of type 2 inflammatory cytokines and raised immunoglobulin concentrations in nasal polyp tissue. By using single-cell RNA sequencing, transcriptomics, surface proteomics, and T cell and B cell receptor sequencing, we found the predominant cell types in nasal polyps were shifted from epithelial and mesenchymal cells to inflammatory cells compared to nasal mucosa from healthy controls. Broad expansions of CD4 T effector memory cells, CD4 tissue-resident memory T cells, CD8 T effector memory cells and all subtypes of B cells in nasal polyp tissues. The T and B cell receptor repertoires were skewed in NP. This study highlights the deviated immune response and remodeling mechanisms that contribute to the pathogenesis of uncontrolled severe CRSwNP. CLINICAL IMPLICATIONS: We identified differences in the cellular compositions, transcriptomes, proteomes, and deviations in the immune profiles of T cell and B cell receptors as well as alterations in the intercellular communications in uncontrolled severe CRSwNP patients versus healthy controls, which might help to define potential therapeutic targets in the future.

ADAR1 interaction with Z-RNA promotes editing of endogenous double-stranded RNA and prevents MDA5-dependent immune activation.

Loss of function of adenosine deaminase acting on double-stranded RNA (dsRNA)-1 (ADAR1) causes the severe autoinflammatory disease Aicardi-Goutières syndrome (AGS). ADAR1 converts adenosines into inosines within dsRNA. This process called A-to-I editing masks self-dsRNA from detection by the antiviral dsRNA sensor MDA5. ADAR1 binds to dsRNA in both the canonical A-form and the poorly defined Z conformation (Z-RNA). Mutations in the Z-RNA-binding Zα domain of ADAR1 are common in patients with AGS. How loss of ADAR1/Z-RNA interaction contributes to disease development is unknown. We demonstrate that abrogated binding of ADAR1 to Z-RNA leads to reduced A-to-I editing of dsRNA structures formed by base pairing of inversely oriented short interspersed nuclear elements. Preventing ADAR1 binding to Z-RNA triggers an MDA5/MAVS-mediated type I interferon response and leads to the development of lethal autoinflammation in mice. This shows that the interaction between ADAR1 and Z-RNA restricts sensing of self-dsRNA and prevents AGS development.

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.

Distinct transcriptome signatures of Helicobacter suis and Helicobacter heilmannii strains upon adherence to human gastric epithelial cells.

The porcine Helicobacter suis and canine-feline H. heilmannii are gastric Helicobacter species with zoonotic potential. However, little is known about the pathogenesis of human infections with these Helicobacter species. To gain more insight into the interactions of both zoonotic Helicobacter species with human gastric epithelial cells, we investigated bacterial genes that are differentially expressed in a H. suis and H. heilmannii strain after adhesion to the human gastric epithelial cell line MKN7. In vitro Helicobacter-MKN7 binding assays were performed to obtain bacterial RNA for sequencing analysis. H. suis and H. heilmannii bacteria attached to the gastric epithelial cells (i.e. cases) as well as unbound bacteria (i.e. controls) were isolated, after which prokaryotic RNA was purified and sequenced. Differentially expressed genes were identified using the DESeq2 package and SARTools pipeline in R. A list of 134 (83 up-regulated and 51 down-regulated) and 143 (60 up-regulated and 83 down-regulated) differentially expressed genes (padj ≤ 0.01; fold change ≥ 2) were identified for the adherent H. suis and H. heilmannii strains, respectively. According to BLASTp analyses, only 2 genes were commonly up-regulated and 4 genes commonly down-regulated in both pathogens. Differentially expressed genes of the H. suis and H. heilmannii strains belonged to multiple functional classes, indicating that adhesion of both strains to human gastric epithelial cells evokes pleiotropic adaptive responses. Our results suggest that distinct pathways are involved in human gastric colonization of H. suis and H. heilmannii. Further research is needed to elucidate the clinical significance of these findings.

Author Correction: Development and performance of a targeted whole exome sequencing enrichment kit for the dog (Canis Familiaris Build 3.1).

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Photothermally Triggered Endosomal Escape and Its Influence on Transfection Efficiency of Gold-Functionalized JetPEI/pDNA Nanoparticles.

Plasmonic nanoparticles for drug delivery have attracted increasing interest over the last few years. Their localized surface plasmon resonance causes photothermal effects on laser irradiation, which allows for delivering drugs in a spatio-temporally controlled manner. Here, we explore the use of gold nanoparticles (AuNP) as carriers for pDNA in combination with pulsed laser irradiation to induce endosomal escape, which is currently considered to be one of the major bottlenecks in macromolecular drug delivery on the intracellular level. In particular, we evaluate nanocomplexes composed of JetPEI (polyethylenimine)pDNA and 10 nm AuNP, which do not exhibit endosomal escape by themselves. After incubating HeLa cells with these complexes, we evaluated endosomal escape and transfection efficiency using low- and high-energy laser pulses. At low laser energy heat is produced by the nanocomplexes, while, at higher laser energy, explosive vapour nanobubbles (VNB) are formed. We investigated the ability of heat transfer and VNB formation to induce endosomal escape and we examine the integrity of pDNA cargo after inducing both photothermal effects. We conclude that JetPEI/pDNA/AuNP complexes are unable to induce meaningful transfection efficiencies because laser treatment causes either dysfunctionality of the cargo when VNB are formed or forms too small pores in the endosomal membrane to allow pDNA to escape in case of heating. We conclude that laser-induced VNB is the most suitable to induce effective pDNA endosomal escape, but a different nanocomplex structure will be required to keep the pDNA intact.

A genome-wide search for epigenetically [corrected] regulated genes in zebra finch using MethylCap-seq and RNA-seq.

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.

Improved canine exome designs, featuring ncRNAs and increased coverage of protein coding genes.

By limiting sequencing to those sequences transcribed as mRNA, whole exome sequencing is a cost-efficient technique often used in disease-association studies. We developed two target enrichment designs based on the recently released annotation of the canine genome: the exome-plus design and the exome-CDS design. The exome-plus design combines the exons of the CanFam 3.1 Ensembl annotation, more recently discovered protein-coding exons and a variety of non-coding RNA regions (microRNAs, long non-coding RNAs and antisense transcripts), leading to a total size of ≈ 152 Mb. The exome-CDS was designed as a subset of the exome-plus by omitting all 3' and 5' untranslated regions. This reduced the size of the exome-CDS to ≈ 71 Mb. To test the capturing performance, four exome-plus captures were sequenced on a NextSeq 500 with each capture containing four pre-capture pooled, barcoded samples. At an average sequencing depth of 68.3x, 80% of the regions and well over 90% of the targeted base pairs were completely covered at least 5 times with high reproducibility. Based on the performance of the exome-plus, we estimated the performance of the exome-CDS. Overall, these designs provide flexible solutions for a variety of research questions and are likely to be reliable tools in disease studies.

PROTEOFORMER: deep proteome coverage through ribosome profiling and MS integration.

An increasing amount of studies integrate mRNA sequencing data into MS-based proteomics to complement the translation product search space. However, several factors, including extensive regulation of mRNA translation and the need for three- or six-frame-translation, impede the use of mRNA-seq data for the construction of a protein sequence search database. With that in mind, we developed the PROTEOFORMER tool that automatically processes data of the recently developed ribosome profiling method (sequencing of ribosome-protected mRNA fragments), resulting in genome-wide visualization of ribosome occupancy. Our tool also includes a translation initiation site calling algorithm allowing the delineation of the open reading frames (ORFs) of all translation products. A complete protein synthesis-based sequence database can thus be compiled for mass spectrometry-based identification. This approach increases the overall protein identification rates with 3% and 11% (improved and new identifications) for human and mouse, respectively, and enables proteome-wide detection of 5'-extended proteoforms, upstream ORF translation and near-cognate translation start sites. The PROTEOFORMER tool is available as a stand-alone pipeline and has been implemented in the galaxy framework for ease of use.

A proteogenomics approach integrating proteomics and ribosome profiling increases the efficiency of protein identification and enables the discovery of alternative translation start sites.

Next-generation transcriptome sequencing is increasingly integrated with MS to enhance MS-based protein and peptide identification. Recently, a breakthrough in transcriptome analysis was achieved with the development of ribosome profiling (ribo-seq). This technology is based on the deep sequencing of ribosome-protected mRNA fragments, thereby enabling the direct observation of in vivo protein synthesis at the transcript level. In order to explore the impact of a ribo-seq-derived protein sequence search space on MS/MS spectrum identification, we performed a comprehensive proteome study on a human cancer cell line, using both shotgun and N-terminal proteomics, next to ribosome profiling, which was used to delineate (alternative) translational reading frames. By including protein-level evidence of sample-specific genetic variation and alternative translation, this strategy improved the identification score of 69 proteins and identified 22 new proteins in the shotgun experiment. Furthermore, we discovered 18 new alternative translation start sites in the N-terminal proteomics data and observed a correlation between the quantitative measures of ribo-seq and shotgun proteomics with a Pearson correlation coefficient ranging from 0.483 to 0.664. Overall, this study demonstrated the benefits of ribosome profiling for MS-based protein and peptide identification and we believe this approach could develop into a common practice for next-generation proteomics.

Development and performance of a targeted whole exome sequencing enrichment kit for the dog (Canis Familiaris Build 3.1).

Whole exome sequencing is a technique that aims to selectively sequence all exons of protein-coding genes. A canine whole exome sequencing enrichment kit was designed based on the latest canine reference genome (build 3.1.72). Its performance was tested by sequencing 2 exome captures, each consisting of 4 pre-capture pooled, barcoded Illumina libraries on an Illumina HiSeq 2500. At an average sequencing depth of 102x, 83 to 86% of the target regions were completely sequenced with a minimum coverage of five and 90% of the reads mapped on the target regions. Additionally, it is shown that the reproducibility within and between captures is high and that pooling four samples per capture is a valid option. Overall, we have demonstrated the strong performance of this WES enrichment kit and are confident it will be a valuable tool in future disease association studies.

Massively parallel sequencing for early molecular diagnosis in Leber congenital amaurosis.

PURPOSE: Leber congenital amaurosis (LCA) is a rare congenital retinal dystrophy associated with 16 genes. Recent breakthroughs in LCA gene therapy offer the first prospect of treating inherited blindness, which requires an unequivocal and early molecular diagnosis. While present genetic tests do not address this due to a tremendous genetic heterogeneity, massively parallel sequencing (MPS) strategies might bring a solution. Here, we developed a comprehensive molecular test for LCA based on targeted MPS of all exons of 16 known LCA genes. METHODS: We designed a unique and flexible workflow for targeted resequencing of all 236 exons from 16 LCA genes based on quantitative PCR (qPCR) amplicon ligation, shearing, and parallel sequencing of multiple patients on a single lane of a short-read sequencer. Twenty-two prescreened LCA patients were included, five of whom had a known molecular cause. RESULTS: Validation of 107 variations was performed as proof of concept. In addition, the causal genetic defect and a single heterozygous mutation were identified in 3 and 5, respectively, of 17 patients without previously identified mutations. CONCLUSION: We propose a novel targeted MPS-based approach that is suitable for accurate, fast, and cost-effective early molecular testing in LCA, and easily applicable in other genetic disorders.