Identification of extracellular vesicle microRNA signatures specifically linked to inflammatory and metabolic mechanisms in obesity-associated low type-2 asthma.

RATIONALE AND OBJECTIVE: Plasma extracellular vesicles (EVs) represent a vital source of molecular information about health and disease states. Due to their heterogenous cellular sources, EVs and their cargo may predict specific pathomechanisms behind disease phenotypes. Here we aimed to utilize EV microRNA (miRNA) signatures to gain new insights into underlying molecular mechanisms of obesity-associated low type-2 asthma. METHODS: Obese low type-2 asthma (OA) and non-obese low type-2 asthma (NOA) patients were selected from an asthma cohort conjointly with healthy controls. Plasma EVs were isolated and characterised by nanoparticle tracking analysis. EV-associated small RNAs were extracted, sequenced and bioinformatically analysed. RESULTS: Based on EV miRNA expression profiles, a clear distinction between the three study groups could be established using a principal component analysis. Integrative pathway analysis of potential target genes of the differentially expressed miRNAs revealed inflammatory cytokines (e.g., interleukin-6, transforming growth factor-beta, interferons) and metabolic factors (e.g., insulin, leptin) signalling pathways to be specifically associated with OA. The miR-17-92 and miR-106a-363 clusters were significantly enriched only in OA. These miRNA clusters exhibited discrete bivariate correlations with several key laboratory (e.g., C-reactive protein) and lung function parameters. Plasma EV miRNA signatures mirrored blood-derived CD4+ T-cell transcriptome data, but achieved an even higher sensitivity in identifying specifically affected biological pathways. CONCLUSION: The identified plasma EV miRNA signatures and particularly the miR-17-92 and -106a-363 clusters were capable to disentangle specific mechanisms of the obesity-associated low type-2 asthma phenotype, which may serve as basis for stratified treatment development.

DNA hypomethylation of the host tree impairs interaction with mutualistic ectomycorrhizal fungus.

Ectomycorrhizas are an intrinsic component of tree nutrition and responses to environmental variations. How epigenetic mechanisms might regulate these mutualistic interactions is unknown. By manipulating the level of expression of the chromatin remodeler DECREASE IN DNA METHYLATION 1 (DDM1) and two demethylases DEMETER-LIKE (DML) in Populus tremula × Populus alba lines, we examined how host DNA methylation modulates multiple parameters of the responses to root colonization with the mutualistic fungus Laccaria bicolor. We compared the ectomycorrhizas formed between transgenic and wild-type (WT) trees and analyzed their methylomes and transcriptomes. The poplar lines displaying lower mycorrhiza formation rate corresponded to hypomethylated overexpressing DML or RNAi-ddm1 lines. We found 86 genes and 288 transposable elements (TEs) differentially methylated between WT and hypomethylated lines (common to both OX-dml and RNAi-ddm1) and 120 genes/1441 TEs in the fungal genome suggesting a host-induced remodeling of the fungal methylome. Hypomethylated poplar lines displayed 205 differentially expressed genes (cis and trans effects) in common with 17 being differentially methylated (cis). Our findings suggest a central role of host and fungal DNA methylation in the ability to form ectomycorrhizas including not only poplar genes involved in root initiation, ethylene and jasmonate-mediated pathways, and immune response but also terpenoid metabolism.

Intranasal administration of Acinetobacter lwoffii in a murine model of asthma induces IL-6-mediated protection associated with cecal microbiota changes.

BACKGROUND: Early-life exposure to certain environmental bacteria including Acinetobacter lwoffii (AL) has been implicated in protection from chronic inflammatory diseases including asthma later in life. However, the underlying mechanisms at the immune-microbe interface remain largely unknown. METHODS: The effects of repeated intranasal AL exposure on local and systemic innate immune responses were investigated in wild-type and Il6-/- , Il10-/- , and Il17-/- mice exposed to ovalbumin-induced allergic airway inflammation. Those investigations were expanded by microbiome analyses. To assess for AL-associated changes in gene expression, the picture arising from animal data was supplemented by in vitro experiments of macrophage and T-cell responses, yielding expression and epigenetic data. RESULTS: The asthma preventive effect of AL was confirmed in the lung. Repeated intranasal AL administration triggered a proinflammatory immune response particularly characterized by elevated levels of IL-6, and consequently, IL-6 induced IL-10 production in CD4+ T-cells. Both IL-6 and IL-10, but not IL-17, were required for asthma protection. AL had a profound impact on the gene regulatory landscape of CD4+ T-cells which could be largely recapitulated by recombinant IL-6. AL administration also induced marked changes in the gastrointestinal microbiome but not in the lung microbiome. By comparing the effects on the microbiota according to mouse genotype and AL-treatment status, we have identified microbial taxa that were associated with either disease protection or activity. CONCLUSION: These experiments provide a novel mechanism of Acinetobacter lwoffii-induced asthma protection operating through IL-6-mediated epigenetic activation of IL-10 production and with associated effects on the intestinal microbiome.

A comparison of high-throughput SARS-CoV-2 sequencing methods from nasopharyngeal samples.

The COVID-19 pandemic caused by the new Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) continues to threaten public health and burden healthcare systems worldwide. Whole SARS-CoV-2 genome sequencing has become essential for epidemiological monitoring and identification of new variants, which could represent a risk of increased transmissibility, virulence, or resistance to vaccines or treatment. Different next-generation sequencing approaches are used in SARS-CoV-2 sequencing, although with different ability to provide whole genome coverage without gaps and to reliably detect new variants. In this study, we compared the performance of three target enrichment methods (two multiplex amplification methods and one hybridization capture) using nasopharyngeal swabs from infected individuals. We applied these target enrichment methods to the same set of nasopharyngeal samples (N = 93) in high-throughput mode. SARS-CoV-2 genome was obtained using short-read next-generation sequencing. We observed that each method has some advantages, such as high mapping rate (CleanPlex and COVIDSeq) or absence of systematic variant calling error (SureSelect) as well as their limitations such as suboptimal uniformity of coverage (CleanPlex), high cost (SureSelect) or supply shortages (COVIDSeq). Nevertheless, each of the three target enrichment kits tested in this study yielded acceptable results of whole SARS-CoV-2 genome sequencing and either of them can therefore be used in prospective programs of genomic surveillance of SARS-CoV-2. Genomic surveillance will be crucial to overcoming the ongoing pandemic of COVID-19, despite its successive waves and continually emerging variants.

From Methylome to Integrative Analysis of Tissue Specificity.

DNA methylation is the most studied epigenetic mark in both plants and animals. The gold standard for assaying genome-wide DNA methylation at single-base resolution is whole-genome bisulfite sequencing (WGBS). Here, we describe an improved procedure for WGBS and original bioinformatic workflows applied to unravel tissue-specific variations of the methylome in relation to gene expression and accumulation of secondary metabolites in the medicinal plant Catharanthus roseus.

RNAi suppression of DNA methylation affects the drought stress response and genome integrity in transgenic poplar.

Trees are long-lived organisms that continuously adapt to their environments, a process in which epigenetic mechanisms are likely to play a key role. Via downregulation of the chromatin remodeler DECREASED IN DNA METHYLATION 1 (DDM1) in poplar (Populus tremula × Populus alba) RNAi lines, we examined how DNA methylation coordinates genomic and physiological responses to moderate water deficit. We compared the growth and drought response of two RNAi-ddm1 lines to wild-type (WT) trees under well-watered and water deficit/rewatering conditions, and analyzed their methylomes, transcriptomes, mobilomes and phytohormone contents in the shoot apical meristem. The RNAi-ddm1 lines were more tolerant to drought-induced cavitation but did not differ in height or stem diameter growth. About 5000 differentially methylated regions were consistently detected in both RNAi-ddm1 lines, colocalizing with 910 genes and 89 active transposable elements. Under water deficit conditions, 136 differentially expressed genes were found, including many involved in phytohormone pathways; changes in phytohormone concentrations were also detected. Finally, the combination of hypomethylation and drought led to the mobility of two transposable elements. Our findings suggest major roles for DNA methylation in regulation of genes involved in hormone-related stress responses, and the maintenance of genome integrity through repression of transposable elements.

Developmental Methylome of the Medicinal Plant Catharanthus roseus Unravels the Tissue-Specific Control of the Monoterpene Indole Alkaloid Pathway by DNA Methylation.

Catharanthus roseus produces a wide spectrum of monoterpene indole alkaloids (MIAs). MIA biosynthesis requires a tightly coordinated pathway involving more than 30 enzymatic steps that are spatio-temporally and environmentally regulated so that some MIAs specifically accumulate in restricted plant parts. The first regulatory layer involves a complex network of transcription factors from the basic Helix Loop Helix (bHLH) or AP2 families. In the present manuscript, we investigated whether an additional epigenetic layer could control the organ-, developmental- and environmental-specificity of MIA accumulation. We used Whole-Genome Bisulfite Sequencing (WGBS) together with RNA-seq to identify differentially methylated and expressed genes among nine samples reflecting different plant organs and experimental conditions. Tissue specific gene expression was associated with specific methylation signatures depending on cytosine contexts and gene parts. Some genes encoding key enzymatic steps from the MIA pathway were found to be simultaneously differentially expressed and methylated in agreement with the corresponding MIA accumulation. In addition, we found that transcription factors were strikingly concerned by DNA methylation variations. Altogether, our integrative analysis supports an epigenetic regulation of specialized metabolisms in plants and more likely targeting transcription factors which in turn may control the expression of enzyme-encoding genes.

Whole-Genome Bisulfite Sequencing for the Analysis of Genome-Wide DNA Methylation and Hydroxymethylation Patterns at Single-Nucleotide Resolution.

The analysis of genome-wide epigenomic alterations including DNA methylation and hydroxymethylation has become a subject of intensive research for many biological and disease-associated investigations. Whole-genome bisulfite sequencing (WGBS) using next-generation sequencing technologies is currently considered as the gold standard for a comprehensive and quantitative analysis of DNA methylation throughout the genome. However, bisulfite conversion does not allow distinguishing between cytosine methylation and hydroxymethylation requiring an additional chemical or enzymatic step to identify hydroxymethylated cytosines. Here we provide two detailed protocols based on commercial kits for the preparation of sequencing libraries for the comprehensive whole-genome analysis of DNA methylation and/or hydroxymethylation. If only DNA methylation is of interest, sequencing libraries can be constructed from limited amounts of input DNA by ligation of methylated adaptors to the fragmented DNA prior to bisulfite conversion. For samples with significant levels of hydroxymethylation such as stem cells or brain tissue, we describe the protocol of oxidative bisulfite sequencing (OxBs-seq), which in its current version uses a post-bisulfite adaptor tagging (PBAT) approach. Two methylomes need to be generated: a classic methylome following bisulfite conversion and analyzing both methylated and hydroxymethylated cytosines and a methylome analyzing only methylated cytosines, respectively. We also provide a step-by-step description of the data analysis using publicly available bioinformatic tools. The described protocols have been successfully applied to different human samples and yield robust and reproducible results.

Deregulation of microRNA expression in purified T and B lymphocytes from patients with primary Sjögren's syndrome.

OBJECTIVE: MicroRNAs (miRNAs) play an important role in the pathogenesis of autoimmune diseases such as primary Sjögren's syndrome (pSS). This study is the first to investigate miRNA expression patterns in purified T and B lymphocytes from patients with pSS using a high-throughput quantitative PCR (qPCR) approach. METHODS: Two independent cohorts of both patients with pSS and controls, one for discovery and one for replication, were included in this study. CD4+ T cells and CD19+ B cells were isolated from peripheral blood mononuclear cells by magnetic microbeads and expression of miRNAs was profiled using the Exiqon Human miRNome panel I analysing 372 miRNAs. A selection of differentially expressed miRNAs was replicated in the second cohort using specific qPCR assays. RESULTS: A major difference in miRNA expression patterns was observed between the lymphocyte populations from patients with pSS and controls. In CD4 T lymphocytes, hsa-let-7d-3p, hsa-miR-155-5 p, hsa-miR-222-3 p, hsa-miR-30c-5p, hsa-miR-146a-5p, hsa-miR-378a-3p and hsa-miR-28-5 p were significantly differentially expressed in both the discovery and the replication cohort. In B lymphocytes, hsa-miR-378a-3p, hsa-miR-222-3 p, hsa-miR-26a-5p, hsa-miR-30b-5p and hsa-miR-19b-3p were significantly differentially expressed. Potential target mRNAs were enriched in disease relevant pathways. Expression of B-cell activating factor (BAFF) mRNA was inversely correlated with the expression of hsa-miR-30b-5p in B lymphocytes from patients with pSS and functional experiments showed increased expression of BAFF after inhibiting hsa-miR-30b-5p. CONCLUSIONS: This study demonstrates major miRNAs deregulation in T and B cells from patients with pSS in two independent cohorts, which might target genes known to be involved in the pathogenesis of pSS.

Whole-Genome Bisulfite Sequencing Using the Ovation® Ultralow Methyl-Seq Protocol.

The analysis of genome-wide epigenomic alterations including DNA methylation has become a subject of intensive research for many complex diseases. Whole-genome bisulfite sequencing (WGBS) using next-generation sequencing technologies can be considered the gold standard for a comprehensive and quantitative analysis of cytosine methylation throughout the genome. Several approaches including tagmentation- and post bisulfite adaptor tagging (PBAT)-based WGBS have been devised. Here, we provide a detailed protocol based on a commercial kit for the preparation of libraries for WGBS from limited amounts of input DNA (50-100 ng) using the classical approach of WGBS by ligation of methylated adaptors to the fragmented DNA prior to bisulfite conversion. The converted library is then amplified with an optimal number of PCR cycles to ensure high sequence diversity and low duplicate rates. Spike-in of unmethylated DNA allows for the precise estimation of bisulfite conversion rates. We also provide a step-by-step description of the data analysis using publicly available bioinformatic tools. The described protocol has been successfully applied to different human samples as well as DNA extracted from plant tissues and yields robust and reproducible results.

Comprehensive evaluation of methods to isolate, quantify, and characterize circulating cell-free DNA from small volumes of plasma.

Circulating cell-free DNA (ccfDNA) has great potential for non-invasive diagnostics, and prediction and monitoring of treatment response, but its amount is usually limited. Therefore, the choice of methods to extract and characterize ccfDNA is crucial. In the current study, we performed the most comprehensive comparison of methods for ccfDNA extraction (11 methods), quantification (3 methods), and estimation of the integrity index (2 methods) from small quantities of different kinds of plasma. The QIAamp® Circulating Nucleic Acid Kit and the Norgen Plasma/Serum Circulating DNA Purification Mini Kit showed the best accuracy and reproducibility, but the Norgen kit allowed to extract a higher amount of ccfDNA. This workflow provides a reliable protocol for the multiple applications of ccfDNA in biomedicine.

Accurate CpG and non-CpG cytosine methylation analysis by high-throughput locus-specific pyrosequencing in plants.

Pyrosequencing permits accurate quantification of DNA methylation of specific regions where the proportions of the C/T polymorphism induced by sodium bisulfite treatment of DNA reflects the DNA methylation level. The commercially available high-throughput locus-specific pyrosequencing instruments allow for the simultaneous analysis of 96 samples, but restrict the DNA methylation analysis to CpG dinucleotide sites, which can be limiting in many biological systems. In contrast to mammals where DNA methylation occurs nearly exclusively on CpG dinucleotides, plants genomes harbor DNA methylation also in other sequence contexts including CHG and CHH motives, which cannot be evaluated by these pyrosequencing instruments due to software limitations. Here, we present a complete pipeline for accurate CpG and non-CpG cytosine methylation analysis at single base-resolution using high-throughput locus-specific pyrosequencing. The devised approach includes the design and validation of PCR amplification on bisulfite-treated DNA and pyrosequencing assays as well as the quantification of the methylation level at every cytosine from the raw peak intensities of the Pyrograms by two newly developed Visual Basic Applications. Our method presents accurate and reproducible results as exemplified by the cytosine methylation analysis of the promoter regions of two Tomato genes (NOR and CNR) encoding transcription regulators of fruit ripening during different stages of fruit development. Our results confirmed a significant and temporally coordinated loss of DNA methylation on specific cytosines during the early stages of fruit development in both promoters as previously shown by WGBS. The manuscript describes thus the first high-throughput locus-specific DNA methylation analysis in plants using pyrosequencing.

Plant genetics. Early allopolyploid evolution in the post-Neolithic Brassica napus oilseed genome.

Oilseed rape (Brassica napus L.) was formed ~7500 years ago by hybridization between B. rapa and B. oleracea, followed by chromosome doubling, a process known as allopolyploidy. Together with more ancient polyploidizations, this conferred an aggregate 72× genome multiplication since the origin of angiosperms and high gene content. We examined the B. napus genome and the consequences of its recent duplication. The constituent An and Cn subgenomes are engaged in subtle structural, functional, and epigenetic cross-talk, with abundant homeologous exchanges. Incipient gene loss and expression divergence have begun. Selection in B. napus oilseed types has accelerated the loss of glucosinolate genes, while preserving expansion of oil biosynthesis genes. These processes provide insights into allopolyploid evolution and its relationship with crop domestication and improvement.

Ultrasensitive detection and identification of BRAF V600 mutations in fresh frozen, FFPE, and plasma samples of melanoma patients by E-ice-COLD-PCR.

A number of molecular diagnostic methods have been developed for the detection and identification of mutations in tumor samples, which are important for the choice of treatment in the context of personalized medicine. For the treatment of metastatic melanoma, Vemurafenib is recommended for patients with BRAF V600 activating mutations. However, the different assays developed to date for the detection of these mutations lack sensitivity or specificity or do not allow a sequencing-based identification or validation of the mutation. Recently, enhanced improved and complete enrichment co-amplification at lower denaturation temperature-polymerase chain reaction (E-ice-COLD-PCR) has been developed as a sensitive method for the detection and identification of mutations in KRAS codons 12/13. Here, we present the first E-ice-COLD-PCR assay for the detection and identification of BRAF codon 600 mutations, which has a large dynamic range, as 25 pg to 25 ng can be used as DNA input without any reduction in mutation enrichment efficiency, and which can detect down to 0.01 % of mutated alleles in a wild-type background. The assay has been validated on fresh frozen, formalin-fixed paraffin-embedded (FFPE), and plasma samples of melanoma patients and has allowed the detection and identification of BRAF mutations present in samples appearing as wild type using standard pyrosequencing, endpoint genotyping, or Sanger sequencing. Thus, the BRAF V600 E-ice-COLD-PCR assay is currently one of the most powerful molecular diagnostic tools for the ultrasensitive detection and identification of BRAF codon 600 mutations.

In utero undernutrition in male mice programs liver lipid metabolism in the second-generation offspring involving altered Lxra DNA methylation.

Obesity and type 2 diabetes have a heritable component that is not attributable to genetic factors. Instead, epigenetic mechanisms may play a role. We have developed a mouse model of intrauterine growth restriction (IUGR) by in utero malnutrition. IUGR mice developed obesity and glucose intolerance with aging. Strikingly, offspring of IUGR male mice also developed glucose intolerance. Here, we show that in utero malnutrition of F1 males influenced the expression of lipogenic genes in livers of F2 mice, partly due to altered expression of Lxra. In turn, Lxra expression is attributed to altered DNA methylation of its 5' UTR region. We found the same epigenetic signature in the sperm of their progenitors, F1 males. Our data indicate that in utero malnutrition results in epigenetic modifications in germ cells (F1) that are subsequently transmitted and maintained in somatic cells of the F2, thereby influencing health and disease risk of the offspring.

Integrated analysis of high-resolution DNA methylation profiles, gene expression, germline genotypes and clinical end points in breast cancer patients.

Breast cancer is a heterogeneous disease for which alterations in DNA methylation patterns have been shown to be of biological and clinical importance. Here we report on the integrated analysis of molecular alterations including the methylation status of 27 gene promoters analyzed by highly quantitative pyrosequencing, and the association to gene expression, germline genotype and clinical parameters including survival. Breast cancer specific deregulation of DNA methylation (both hyper- and hypomethylation) was found in twenty genes including ACVR1, OGG1, IL8 and TFF1. The methylation level in the promoter regions was significantly negatively correlated to gene expression for twelve genes (such as MST1R, ST6GAL1 and TFF1) indicating that a gain of aberrant methylation (hypermethylation) inhibits gene expression. Multiple associations between molecular and clinical parameters were identified, and multivariate statistical analysis demonstrated that methylation was more strongly associated to clinical parameters than gene expression for the investigated genes. The methylation level of BCAP31 and OGG1 showed significant association to survival, and these associations were validated in a larger patient cohort (The Cancer Genome Atlas). Our study provides evidence for the promise of DNA methylation alterations for clinical applications.

Methylation profiling with a panel of cancer related genes: association with estrogen receptor, TP53 mutation status and expression subtypes in sporadic breast cancer.

Breast cancer is a heterogeneous disease that can be divided in subtypes based on histology, gene expression profiles as well as differences in genomic aberrations. Distinct global DNA methylation profiles have been reported in normal breast epithelial cells as well as in breast tumors. However, the influence of the tumor methylome on the previously described subgroups of breast cancer is not fully understood. Here we report the DNA methylation profiles of 80 breast tumors using a panel of 807 cancer related genes interrogating 1505 CpG sites. We identified three major clusters based on the methylation profiles; one consisting of mainly tumors of myoepithelial origin and two other clusters with tumors of predominantly luminal epithelial origin. The clusters were different with respect to estrogen receptor status, TP53 status, ErbB2 status and grade. The most significantly differentially methylated genes including HDAC1, TFF1, OGG1, BMP3, FZD9 and HOXA11 were confirmed by pyrosequencing. Gene Ontology analysis revealed enrichment for genes involved in developmental processes including homeobox domain genes (HOXA9, HOXA11, PAX6, MYBL2, ISL1 and IPF1) and (ETS1, HDAC1, CREBBP, GAS7, SPI1 and TBX1). Extensive correlation to mRNA expression was observed. Pathway analyses identified a significant association with canonical (curated) pathways such as hepatic fibrosis including genes like EGF, NGFR and TNF, dendritic cell maturation and the NF-κB signaling pathway. Our results show that breast tumor expression subtypes harbor major epigenetic differences and tumors with similar gene expression profiles might belong to epigenetically different subtypes. Some of the transcription factors identified, with key roles in differentiation and development might play a role in inducing and maintaining the different phenotypes.