Single-embryo RNA sequencing for continuous and sex-specific gene expression analysis on Drosophila

Exploring early embryonic gene expression is challenging due to the rate of development and the limited material available. Here, we present a protocol for ordering Drosophila embryos along a developmental pseudo-time trajectory and determining the sex of the embryos using RNA-seq data. We describe steps for sample collection, RNA isolation, RNA-seq, and RNA-seq data processing. We then detail the establishment of a continuous transcriptome dataset for assessing gene expression throughout early development and in a sex-specific manner. For complete details on the use and execution of this protocol, please refer to Pérez-Mojica et al.1.

Continuous transcriptome analysis reveals novel patterns of early gene expression in Drosophila embryos

The transformative events during early organismal development lay the foundation for body formation and long-term phenotype. The rapid progression of events and the limited material available present major barriers to studying these earliest stages of development. Herein, we report an operationally simple RNA sequencing approach for high-resolution, time-sensitive transcriptome analysis in early (≤3 h) Drosophila embryos. This method does not require embryo staging but relies on single-embryo RNA sequencing and transcriptome ordering along a developmental trajectory (pseudo-time). The resulting high-resolution, time-sensitive mRNA expression profiles reveal the exact onset of transcription and degradation for thousands of transcripts. Further, using sex-specific transcription signatures, embryos can be sexed directly, eliminating the need for Y chromosome genotyping and revealing patterns of sex-biased transcription from the beginning of zygotic transcription. Our data provide an unparalleled resolution of gene expression during early development and enhance the current understanding of early transcriptional processes.

Docosahexaenoic acid and oleic acid induce altered DNA methylation of individual CpG loci in Jurkat T cells

Docosahexaenoic acid (DHA, 22:6n-3) and oleic acid (18:1n-9) can alter the DNA methylation of individual CpG loci in vivo and in vitro, although the targeting mechanism is unknown. We tested the hypothesis that the targeting of altered methylation is associated with putative transcription factor response elements (pTREs) proximal to modified loci. Jurkat cells were treated with 22:6n-3 or 18:1n-9 (both 15 μM) for eight days and DNA methylation measured using the MethylationEPIC 850K array. 1596 CpG loci were altered significantly (508 hypermethylated) by 22:6n-3 and 563 CpG loci (294 hypermethylated) by 18:1n-9. 78 loci were modified by both fatty acids. Induced differential methylation was not modified by the PPARα antagonist GW6471. DNA sequences proximal to differentially methylated CpG loci were enriched in zinc-finger pTREs. These findings suggest that zinc-finger-containing transcription factors may be involved in targeting altered DNA methylation modifying processes induced by fatty acids to individual CpG loci.

Polyunsaturated Fatty Acid Biosynthesis Involving Δ8 Desaturation and Differential DNA Methylation of FADS2 Regulates Proliferation of Human Peripheral Blood Mononuclear Cells

Polyunsaturated fatty acids (PUFAs) are important for immune function. Limited evidence indicates that immune cell activation involves endogenous PUFA synthesis, but this has not been characterised. To address this, we measured metabolism of 18:3n-3 in quiescent and activated peripheral blood mononuclear cells (PBMCs), and in Jurkat T cell leukaemia. PBMCs from men and women (n = 34) were incubated with [1-¹³C]18:3n-3 with or without Concanavalin A (Con. A). 18:3n-3 conversion was undetectable in unstimulated PBMCs, but up-regulated when stimulated. The main products were 20:3n-3 and 20:4n-3, while 18:4n-3 was undetectable, suggesting initial elongation and Δ8 desaturation. PUFA synthesis was 17.4-fold greater in Jurkat cells than PBMCs. The major products of 18:3n-3 conversion in Jurkat cells were 20:4n-3, 20:5n-3, and 22:5n-3. ¹³C Enrichment of 18:4n-3 and 20:3n-3 suggests parallel initial elongation and Δ6 desaturation. The FADS2 inhibitor SC26196 reduced PBMC, but not Jurkat cell, proliferation suggesting PUFA synthesis is involved in regulating mitosis in PBMCs. Con. A stimulation increased FADS2, FADS1, ELOVL5 and ELOVL4 mRNA expression in PBMCs. A single transcript corresponding to the major isoform of FADS2, FADS20001, was detected in PBMCs and Jurkat cells. PBMC activation induced hypermethylation of a 470bp region in the FADS2 5'-regulatory sequence. This region was hypomethylated in Jurkat cells compared to quiescent PBMCs. These findings show that PUFA synthesis involving initial elongation and Δ8 desaturation is involved in regulating PBMC proliferation and is regulated via transcription possibly by altered DNA methylation. These processes were dysregulated in Jurkat cells. This has implications for understanding the regulation of mitosis in normal and transformed lymphocytes.

Associations between whole peripheral blood fatty acids and DNA methylation in humans

Fatty acids (FA) modify DNA methylation in vitro, but limited information is available on whether corresponding associations exist in vivo and reflect any short-term effect of the diet. Associations between global DNA methylation and FAs were sought in blood from lactating infants (LI; n = 49) and adult males (AMM; n = 12) equally distributed across the three conventional BMI classes. AMM provided multiple samples at 2-hour intervals during 8 hours after either a single Western diet-representative meal (post-prandial samples) or no meal (fasting samples). Lipid/glucose profile, HDAC4 promoter and PDK4 5’UTR methylation were determined in AMM. Multiple regression analysis revealed that global (in LI) and both global and PDK4-specific DNA methylation (in AMM) were positively associated with eicosapentaenoic and arachidonic acid. HDAC4 methylation was inversely associated with arachidonic acid post-prandially in AMM. Global DNA methylation did not show any defined within-day pattern that would suggest a short-term response to the diet. Nonetheless, global DNA methylation was higher in normal weight subjects both post-prandially and in fasting and coincided with higher polyunsaturated relative to monounsaturated and saturated FAs. We show for the first time strong associations of DNA methylation with specific FAs in two human cohorts of distinct age, diet and postnatal development stage.