Dnmt1 is essential for egg production and embryo viability in the large milkweed bug, Oncopeltus fasciatus

Background The function of cytosine (DNA) methylation in insects remains inconclusive due to a lack of mutant and/or genetic studies. Results Here, we provide evidence for the functional role of the maintenance DNA methyltransferase 1 (Dnmt1) in an insect using experimental manipulation. Through RNA interference (RNAi), we successfully posttranscriptionally knocked down Dnmt1 in ovarian tissue of the hemipteran Oncopeltus fasciatus (the large milkweed bug). Individuals depleted for Dnmt1, and subsequently DNA methylation, failed to reproduce. Eggs were inviable and declined in number, and nuclei structure of follicular epithelium was aberrant. Erasure of DNA methylation from gene or transposon element bodies did not reveal a direct causal link to steady-state mRNA levels in somatic cells. These results reveal an important function of Dnmt1 seemingly not contingent on directly controlling gene expression. Conclusions This study provides direct experimental evidence for a functional role of Dnmt1 in egg production and embryo viability and uncovers a trivial role, if any, for DNA methylation in control of gene expression in O. fasciatus. Electronic supplementary material The online version of this article (10.1186/s13072-018-0246-5) contains supplementary material, which is available to authorized users.


Background
DNA methylation in insects has been hypothesized to play numerous functional roles including polyphenism, diapause, longevity, and social behavior and caste differentiation . However, these hypotheses are based primarily on correlational studies; for example, there is evidence that the dependence of social behavior and caste differentiation on DNA methylation is not absolute [11][12][13]. Comparative epigenomic studies have provided valuable insights into patterns of DNA methylation across most insects, but have offered limited insights into its functional significance.
In mammals, DNA methylation typically occurs at CG sites, which is established by DNA methyltransferase 3 (Dnmt3) and subsequently maintained by DNA methyltransferase 1 (Dnmt1). Expression of Dnmt1 and Dnmt3 is spatially and temporally dynamic in Mus musculus (house mouse) and more notably during development where DNA methylation undergoes two rounds of genome-wide erasure and reestablishment [22]. Proper regulation of Dnmt1 during development is necessary to maintain DNA methylation at imprinted genes [23]. Furthermore, mutations causing partial or severe loss of function of Dnmt1 cause abnormal development and embryonic lethality [24]. As in mammals, many insects possess Dnmt1 and Dnmt3, the expression of which are spatially and temporally dynamic with elevated transcript levels in both male and female gonads and throughout the course of development [12,13,[25][26][27]. Dnmt1 might be necessary for development and embryonic viability in insects, but the functional contribution of DNA methylation is unknown [25].
In many insects, DNA methylation is localized to moderately and constitutively expressed genes that are highly conserved between species. This has led some to the hypothesis that DNA methylation functions in transcriptional regulation [6,15,28,29]. However, functional tests of DNA methyltransferases and the consequences to DNA methylation in insects are limited to the posttranscriptional knockdown of Dnmt3 in Apis mellifera, which is associated with modest reductions in DNA methylation, differential gene expression, and alternative splicing [30]. Given the limited genetic and mutant studies, the function of insect DNA methyltransferases, and DNA methylation itself, is obscure.
Here, we show that DNA methylation maintained by Dnmt1 is essential for egg production and embryo viability and is not associated with transcription of the cells profiled in the hemipteran Oncopeltus fasciatus (the large milkweed bug) [31,32]. Oncopeltus fasciatus has been an established laboratory model organism in fields ranging from physiology and development to evolutionary ecology for > 60 years [33]. Its use as a model in functional genomics is apparent: RNA interference (RNAi) is highly effective across different life history stages (reviewed in [31]), and the recent release of the genome allows for "-omic" level exploration and hypothesis testing [32]. Posttranscriptional knockdown of Dnmt1 in O. fasciatus led to reduced egg viability, fecundity, and aberrant follicular epithelium, and thus failure to produce a successive generation. Despite finding levels of methylated CG (mCG) within coding regions reduced by 83.55%, we found no evidence for DNA methylation directly affecting transcription in the cells examined. Our results suggest Dnmt1 is required for reproduction in O. fasciatus that is likely mediated by a gene-regulatory independent function. Oncopeltus fasciatus represents a fruitful model species for functional studies of DNA methylation, and continuation of studies in this system will unravel the insect epigenome and its functional consequences.

Dnmt1 is necessary for reproduction
To assess the function of Dnmt1, double-stranded RNA (dsRNA) targeting either the cytosine-specific DNA methyltransferase replication foci domain (ds-dnmt1-1) or the DNA methylase domain (AdoMet) (ds-dnmt1-2) of Dnmt1 was injected between the abdominal sternites of virgin O. fasciatus females (Additional file 1: Fig. S1 and Fig. S2, and Tables S1 and S2). For the remainder of the text, we report findings from only ds-dnmt1-1 as resulting phenotypes between the two experiments are essentially the same (Additional file 1: Fig. S2). Maternal dsRNA injection has been used in O. fasciatus to posttranscriptionally knockdown expression in embryos developing from eggs laid by injected females [34] and is known to reduce expression of embryonically expressed genes. A control injection of double-stranded Red (ds-red) or buffer was used to confirm that there was no technical effect associated with RNAi treatment. The effects of RNAi treatment on gene expression were assessed in gut, head, thorax, and ovary tissues 10 days after mating to untreated O. fasciatus males. The posttranscriptional knockdown reduced dnmt1 mRNA expression in all examined tissues (Fig. 1a, Additional file 1: Fig. S3, Additional file 1: Table S2). Moreover, there was little variation in expression among biological replicates of ds-dnmt1-1injected and control individuals, respectively (Fig. 1a). Taken together, treatment with ds-dnmt1-1 specifically and reliably knocked down Dnmt1 transcripts.
Posttranscriptional knockdown of Dnmt1 affects egg development and viability. In the early stages following injection, females injected with ds-dnmt1-1 did not lay fewer eggs (F = 2.91, df = 1; p = 1.1e−01) (Additional file 1: Fig. S4). Although eggs laid by females injected with ds-dnmt1-1 within the first 8 days postinjection looked typical, they were significantly less likely to develop than the eggs laid by control females (χ 2 = 8.470; df = 1; p = 4.0e−03). Furthermore, a mean of 93% of the eggs laid by control females initiated development whereas only a mean of 4% of the eggs laid by ds-dnmt1-1 females initiated development (Fig. 1b). Although eggs laid by the control females that initiated development were viable and hatched, the few eggs laid by ds-dnmt1-1-injected females that initiated development were not viable and failed to hatch.
Posttranscriptional knockdown of Dnmt1 affects egg production following the first 10 days post-injection. By 10-12 days, females injected with ds-dnmt1-1 have mainly stopped laying eggs. While ovaries dissected from control females at this stage have intact eggs in the oviduct, the oviducts of ds-dnmt1-1-injected females are either empty or filled with a mass of what appears to be yolk (Fig. 1c), indicating a fault in production of a functional chorion. Analysis of ovarian cell structure indicates that knockdown of Dnmt1 transcripts affects the cell structure of the follicular epithelium, which is responsible for production of the chorion and vitelline envelope. Nuclei of the follicular epithelium in Dnmt1 posttranscriptional knock-downed females are aberrant and fewer in number than that of the control females (Fig. 1c). Therefore, both through disruption of embryonic development from eggs produced within a week of injection and cessation of egg production, reproduction is compromised following knockdown of Dnmt1 transcripts, preventing a successive generation.

Posttranscriptional knockdown of Dnmt1 successfully and severely reduces mCG in ovaries
The successful knockdown of Dnmt1 transcripts in gut, head, thorax, and ovaries prompted the evaluation of the consequences on DNA methylation using wholegenome bisulfite sequencing (WGBS). We used a low coverage sequencing approach, which is sufficient for the detection of changes in bulk levels of DNA methylation [35]. Although Dnmt1 mRNA expression was reduced in all tissues, DNA methylation was only reduced in ovary tissue ( Fig. 2a-c, Additional file 1: Fig. S3). The reduction specifically in ovaries is likely due to the higher rate of cell division in comparison with cells in other tissues surveyed, which would facilitate the passive of loss of DNA b Parental RNAi injection with ds-dnmt1-1 significantly affected the development of eggs laid by injected females compared to eggs laid by control females. Dots indicate mean number of eggs developed, and error bars indicate standard error of the mean. c Whole ovaries from females removed 12-14 days post-injection (i and iv). In control females (i), mature oocytes can be seen collecting in the lateral oviduct (Od). In ds-dnmt1-1 females (iv), no mature oocytes are apparent, but the oviduct has filled with a yolk-like substance. Follicular epithelium surrounding a maturing oocyte at 20-(ii and v) and 40-times (iii and vi) magnification. Nuclei are stained with Hoechst 33258. Nuclei from control follicular epithelium (ii and iii) are round and regular in shape. Nuclei from ds-dnmt1-1 females (v and vi) are highly irregular in shape methylation in the absence of Dnmt1. Very low expression (< 1 fragments per kilobase of transcript per million [FPKM]) of the de novo DNA methyltransferase Dnmt3 was detected in ovaries from RNA-seq data. Additionally, low levels of expression of Dnmt3 were detected in ovaries using qRT-PCR and did not differ between ds-dnmt1-1 and control individuals (Additional file 1: Fig.  S5). High coverage single-base resolution DNA methylomes from ds-dnmt1-1 and control ovaries were generated to understand the impact of the loss of Dnmt1. Greater than 75% reductions of mCG were observed across the genome, consistent with the decrease observed in the low coverage experiments (Additional file 1: Fig. S3). The presence of symmetrical mCG-DNA methylation occurring on both DNA strands at a CpG site-is indicative of the presence of a functional Dnmt1. In control individuals, methylation at CpGs is highly symmetrical and highly methylated. However, even though knockdown of Dnmt1 significantly reduces methylation in ovaries, a minority of CpG sites remain symmetrically methylated (Fig. 2c).
This suggests that there are some cells within ovaries that have wild-type methylomes. CG methylation is passively lost in the absence of Dnmt1 upon DNA replication. The lack of a complete loss of mCG is expected due to the injection of fully developed individuals with some cells not undergoing further replication and hence the presence of fully methylated genomes that existed prior to the injection of dsRNA.

Decreased levels of mCG are not associated with transcriptome-wide changes in gene or transposable element expression in ovaries
We performed RNA-seq analysis from ovary and other tissues in ds-dnmt1-1 and control individuals to better understand the relationship between gene expression and DNA methylation. No relationship between discrete mCG levels and gene expression was observed, and this was consistent in other insects that we investigated (Fig. 3a, Additional file 1: Fig. S6 and Tables S3-S7). Similarly, no correlation between continuous mCG levels and gene expression was observed in ds-dnmt1-1 and control  Tables  S3-S7). Finally, RNAi treatment resulted in no changes to gene expression despite appreciable reductions to genic DNA methylation levels in O. fasciatus ds-dnmt1-1 (Fig. 3c). A maximum of 264 differentially expressed genes (DEG) were observed between O. fasciatus ds-dnmt1-1 and control ovaries (Additional file 1: Table S4), whereas no DEGs are observed between O. fasciatus ds-dnmt1-1 and control for gut, head, and thorax (Additional file 1: Table S8). Of genes with differences in DNA methylation between ds-dnmt1-1 and control ovaries (N = 6590), the majority (N = 6484; 98.39%) had no changes to gene expression (Fig. 3d, Additional file 1: Tables S3 and S4). Furthermore, mCG was always reduced in ds-dnmt1-1 compared to control ovaries in the 6484 gene set. Despite genes being unmethylated (UM) in ds-dnmt1-1 ovaries, changes to gene expression occurred in both directions for the 106 differentially DNA methylated and expressed genes (Additional file 1: Fig. S7). Even for genes that by definition are unmethylated (N = 5982) or similarly CGmethylated (N = 21) in ds-dnmt1-1 and control, a similar proportion (N = 6003; 97.53%) was observed to have no difference in gene expression. The lack of an association between DNA methylation and gene expression is further supported when applying more stringent thresholds to the definition of CG-methylated and unmethylated genes (Additional file 1: Fig. S8). Our observations support a trivial role of DNA methylation in gene expression if any at all. The function of DNA methylation in O. fasciatus might lie within genome defense-the transcriptional regulation of repetitive DNA and transposons-as its genome is composed of a non-trivial amount of repetitive DNA and transposons (6.21%) (Additional file 1: Table S9). However, as in genes, reduction of mCG in ds-dnmt1-1 compared to control ovaries was found across the bodies of the TEs and was not associated with a genome-wide increase in expression ( Fig. 4a-c, Additional file 1: Fig.  S9 and Table S10). This further supports a trivial role of DNA methylation in transcriptional regulation of loci.

Comparative epigenomic analysis of the O. fasciatus methylome with other insects
The O. fasciatus methylome revealed much higher levels of mCG throughout the genome, similar to other hemimetabolous insects, when compared to holometabolous insects (Fig. 5a, b). DNA methylation of gene bodies is found in O. fasciatus similar to other insects that possess DNA methylation; however, the pattern of mCG within exonic regions of hemi-and holometabolous insects is distinct (Fig. 5b). Higher levels of mCG are toward the 5′ end compared to the 3′ end of coding regions in holometabolous insects. This distribution resembles the most recent common ancestor of all insects (crustaceans) represented by Daphnia pulex. Hemimetabolous insects, which include O. fasciatus, have a more uniform distribution of mCG across coding regions, and higher levels toward the 3′ end of coding regions. Previous studies have demonstrated that CG-methylated genes are often conserved within insects. We identified 39.31% (N = 7561) CG-methylated genes in O. fasciatus, and 85.99% of these are CG-methylated in at least one other insect species investigated (Additional file 1: Table S9). Therefore, even though the pattern of DNA methylation within gene bodies of hemimetabolous insects is distinct compared to holometabolous insects, the targeting of specific genes is conserved. Thus, the reductions of mCG across gene bodies in ds-dnmt1-1 individuals provide an opportunity to study its potential function.

Discussion
DNA methylation has been implicated to play numerous roles in insects ; however, its exact function remains uncertain. In this study, we observed that posttranscriptional knockdown of Dnmt1 negatively impacts reproduction through a decrease in egg development and egg production in O. fasciatus.
Posttranscriptional knockdown of Dnmt1 transcripts impacted the maternal somatic gonad and egg maturation. The follicular epithelial cells undergo multiple rounds of mitosis during oogenesis [36]. Oocytes that were in the early stages of maturation would be the most affected, as their follicular epithelium would passively lose DNA methylation with each successive round of mitosis. Loss was not compensated by the de novo DNA methyltransferase Dnmt3 as expression was low in ovaries. Dnmt3 expression during development and across tissues has not been as well explored in hemimetabolous insects compared to holometabolous insects [26]. However, limited Dnmt3 expression in ovaries, which contain developing oocytes and embryos, possibly reflects differences in epigenomic reprogramming and imprinting during oogenesis and early embryo development compared to mammals [37][38][39][40]. The loss of function of follicular epithelial cells, which includes production of the chorion and vitelline membrane, could explain the loss of oocyte integrity.
The knockdown of Dnmt1 transcripts in O. fasciatus ovaries subsequently did not lead to a causal link between DNA methylation of genes and their mRNA levels in the cells we examined. No correlation between mCG and gene expression in O. fasciatus ovaries and other insects was observed (Fig. 3, Additional file 1: Fig. S6) [6,13]. Additionally, knockdown of Dnmt3 in Apis mellifera, although resulting in minor reductions to genic DNA methylation, resulted in no global changes to gene expression (Additional file 1: Fig. S6). A similar lack of relationship between loss of genic CG methylation levels and gene expression has been reported in the plants Arabidopsis thaliana through met1/dnmt1 mutants and Eutrema salsugineum and Conringia planisiliqua through natural, independent losses of gene body DNA methylation [41,42]. Furthermore, in a recent study, maternal transcriptional knockdown of Dnmt1 in the red flour beetle (Tribolium castaneum)-an insect with Dnmt1, but without DNA methylation-resulted in developmental arrest in offspring embryos [43]. However, the effect on gene expression following dnmt1 knockdown in T. castaneum is unknown [43]. Hence, the importance of mCG in O. fasciatus might lie within the orchestration of genome structure during DNA replication, gametogenesis, or cell types not examined in our study and/or Dnmt1 has an alternative function unrelated to the maintenance of CG methylation. We contribute to a growing body of evidence for a potentially alternative, non-catalytic function of Dnmt1 [43,44].
DNA methylation might be important for the transcriptional regulation of 106 genes that are differentially The subset presents the RPKM distribution when all TEs are considered within ds-dnmt1-1 and control ovaries expressed when their DNA methylation is reduced (Fig. 3). Around 70% of O. fasciatus genes have DNA methylation above the non-conversion rate (Table S3) with the corresponding 106 genes representing < 1% of genes with DNA methylation. Thus, it seems unlikely that the function of DNA methylation is to regulate gene expression in O. fasciatus and possibly other insects. One of these genes (OFAS014839) has lower levels of DNA methylation and expression in ds-dnmt1-1 compared to control and is orthologous to Drosophila melanogaster secretory 23 (Sec23) (NP_730978.1). In D. melanogaster, P-element insertion in the 5′ UTR of sec23 eliminates or severely attenuates gene function and is lethal [45]. Thus, it is possible that one of these 106 genes is responsible for the phenotype, although no obvious candidates stood out on this list, based on the developmental literature around insect oogenesis and early development. Furthermore, the similar phenotype arising in T. castaneum [43], an insect without DNA methylation, suggests that the phenotype is likely not due to the loss of DNA methylation. DNA methylation might be required for downstream functions such as proper mitosis, and the segregation of sister chromatids into their respective daughter cells, as there does appear to be aberrant nuclei structure in ds-dnmt1-1 compared to control ovarian cells (Fig. 1c).
This phenotype could be the result of epigenomic defects in scaffold/matrix attachment regions (S/MAR), as some correlate with origins of replication [46]. This phenotype could be exacerbated by multiple mitotic divisions of follicular epithelial cells during oogenesis, and/or the presence of holocentric chromosomes [47]. Future work describing the epigenomic contributions to cellular defects are now possible to study in this newly emerging and tractable model species, O. fasciatus.
We suggest that DNA methylation in O. fasciatus and likely other insects is more important for genome structure, integrity or other cellular processes than it is for the regulation of somatic gene expression. Although it is possible that DNA methylation is important for expression control in a rare cell type that was not examined in our study, we instead propose that regulation of expression is not likely the primary role in insects. Investigating the importance of DNA methylation prior to the onset of meiosis and throughout embryo development is a fundamental next step. Oncopeltus fasciatus represents a fruitful model species for functional studies of DNA methylation, and continuation of studies in this system will unravel the insect epigenome and its functional consequences.

PCR confirmation for the presence of a single Dnmt1 ortholog in O. fasciatus
To determine whether OFAS015351 and OFAS018396 were two parts of a single Dnmt1 ortholog, we designed one sense primer at the 3′ end of OFAS015351 (Of_ DMNT1-1_3603S; Additional file 1: Table S1) and two antisense primers at the 3′ end of OFAS018396 (Of_ DNMT1-2_424A and Of_DNMT1-2_465A) (Additional file 1: Table S1). A fourth sense primer (Of_DNMT1-2_1S) (Additional file 1: Table S1) was designed at the 5′ end of OFAS018396 to confirm the size and sequence of this possibly truncated gene annotation. Polymerase chain reaction (PCR) with primer combinations Of_DMNT1-1_3603S-Of_DNMT1-2_424A, Of_DMNT1-1_3603S-Of_DNMT1-2_465A, and Of_ DNMT1-2_1S-Of_DNMT1-2_424A was performed using Q5 Polymerase (New England BioLabs, Ipswich, MA) per manufactures instructions. Thermacycler conditions were 98 °C for 15 s (s) (denaturing), 60 °C for 30 s (annealing), and 72 °C for 30 s (extension) and repeated for 40 cycles. The PCR products were then purified using QIAquick PCR Purification Kit (Qiagen, Venlo, The Netherlands) and sequenced at the Georgia Genomics and Bioinformatics Core (Athens, GA).

Animal culture
Oncopeltus fasciatus cultures were originally purchased from Carolina Biologicals (Burlington, NC). Mass colonies were maintained in incubators under a 12 h/12 h light/dark cycle at 27 °C. Colonies and individual experimental animals were fed organic raw sunflower seeds and provided with ad libitum deionized water. Late instar nymphs were separated from the mass colonies and housed under the same conditions. Nymph colonies were checked daily for newly emerged adults. Adults were separated by sex and kept with food and water for 7-10 days until females reached sexual maturity.

Parental RNAi
Template for the in vitro transcription of reactions was prepared from a PCR reaction in which T7 phage promoter sequences were added to the gene-specific Dnmt1 primers [52]. For our control sequence, we used the red fluorescent protein (Red) sequence used in previous parental RNAi experiments in O. fasciatus [52] or buffer. Primer sequences can be found in Table S2. Sense and antisense RNA was synthesized in a single reaction using the Ambion MEGAscript kit (ThermoFisher Sci, Waltham, MA). After purification, the double-stranded RNA (dsRNA) concentration was adjusted to 2 μg/μL in injection buffer (5 mM KCl, 0.1 mM NaH 2 PO 4 ) [52]. Females were injected with 5 μL of dsRNA between the abdominal sternites using an insulin syringe. Following injection, females were paired with an uninjected male to stimulate oviposition and fertilize eggs. Individual females with their mate were housed in petri dishes with sunflower seeds, water and cotton wool as an oviposition site. The parental RNAi protocol has been reported to result in 100% penetrance by the third clutch of eggs [52], and this was also our experience.
To test for any potential off-target effect of our RNAi, we replicated the experiment using a second Dnmt1specific dsRNA. The original dsRNA was targeted to the cytosine-specific DNA methyltransferase replication foci domain (RFD) of Dnmt1 (ds-dnmt1-1). We developed a second dsRNA that was homologous to the DNA methylase domain (AdoMet) of Dnmt1 (ds-dnmt1-2).

Reproductive phenotype screening and analysis
Eggs were collected between days 4-10 post-injection and assessed for development. Oncopeltus fasciatus embryos change from a creamy white color to orange as they develop, which indicates viability. Thus, color change is a useful tool for assessing healthy development. We examined the number of developing eggs at 5 days post-oviposition, at which point viable eggs are clearly distinguishable from inviable eggs, as well as hatching rate of eggs from females injected with double-stranded Dnmt1 (ds-dnmt1-1) and double-stranded Red (ds-red). The data for the number of eggs laid were normally distributed, so differences among RNAi-treated groups were tested using analysis of variance (ANOVA). The development data, however, were not normally distributed and consisted of binary states (developed and not developed), and so differences among treatment groups were analyzed with a generalized linear model (GLM) using a Poisson distribution.
A second set of O. fasciatus females were injected in the same manner as described in "Parental RNAi" to assess posttranscriptional knockdown of Dnmt1 on ovarian structure. Oncopeltus fasciatus females were dissected 10 days after injection. By 10 days post-injection, O. fasciatus females are beginning to stop laying recognizable eggs. Ovaries were removed from O. fasciatus females and placed in 1 × PBS. Whole ovaries were imaged with a Leica DFC295 stereomicroscope using Leica Application Suite morphometric software (LAS V4.1; Leica, Wetzlar, Germany). Dissected ovaries were fixed within 15 min of dissection in 4% formaldehyde in 1 × PBS for 25 min. Fixed ovaries were stained with Hoechst 33342 (Sigma-Aldrich) at 0.5 µg/ml. The stained ovarioles were imaged using a Zeiss LSM 710 Confocal Microscope (Zeiss) at the University of Georgia Biomedical Microscopy Core.

Quantitative RT-PCR
To assess the effectiveness of posttranscriptional knockdown of Dnmt1, females were dissected 11 days postinjection. Ovaries were removed from each female, flash frozen in liquid nitrogen, and stored at − 80 °C until processing. Total RNA (and DNA) was extracted from a single ovary per female using a Qiagen Allprep DNA/RNA Mini Kit (Qiagen, Venlo, The Netherlands) per manufacturer's instructions. Complementary DNA (cDNA) was synthesized from 500 ng RNA with qScript cDNA Super-Mix (Quanta Biosciences, Gaithersburg, MD).
Expression level of Dnmt1 was quantified by quantitative real-time PCR (qRT-PCR). Primers were designed for Dnmt1 using the O. fasciatus genome as a Ref. [52]. Actin and GAPDH were used as endogenous reference genes. Primer sequences can be found in Additional file 1: Table S2. We used Roche LightCycler 480 SYBR Green Master Mix with a Roche LightCycler 480 (Roche Applied Science, Indianapolis, IN) for qRT-PCR. All samples were run with 3 technical replicates using 10 μL reactions using the manufacturer's recommended protocol. Primer efficiency calculations, genomic contamination testing, and endogenous control gene selection were performed as described by [53]. We used the ΔΔCT method [54] to examine differences in expression between control and ds-dnmt1-1 injected females.

Whole-genome bisulfite sequencing (WGBS) and analysis of cytosine (DNA) methylation
MethylC-seq libraries for an O. fasciatus ds-dnmt1-1 and control individual were prepared according to the protocol described in [55] using genomic DNA extracted from ovaries (see "Materials and methods" section "Quantitative RT-PCR"). Libraries were single-end 75 bp sequenced on an Illumina NextSeq 500 machine [56,57]. Unmethylated lambda phage DNA was used to as a control for sodium bisulfite conversion, and an error rate of ~ 0.05% was estimated. Oncopeltus fasciatus ds-dnmt1-1 and control were sequenced to a depth of ~ 18 × and ~ 21 × , which corresponded to an actual mapped coverage of ~ 9 × and ~ 11 × , respectively. Additionally, low pass (< 1 ×) WGBS from gut, head, ovary, and thorax was performed for three control and ds-dnmt1-1 biological replicates. Blattella germanica was additionally sequenced to generate equal numbers of hemi-and holometabolous insects investigated in this study. However, DNA was extracted from whole-body minus gastrointestinal tract. MethylC-seq libraries were prepared and sequenced identically to O. fasciatus. An error rate of ~ 0.14% was estimated from unmethylated lambda phage DNA. Blattella germanica was sequenced to a depth of ~ 8 × , which corresponded to an actual mapped coverage of ~ 5 × . WGBS data for O. fasciatus and B. germanica can be found on Gene Expression Omnibus (GEO) under accession GSE109199. Previously published WGBS data for A. mellifera [30,58], Bo. mori [17], Daphnia pulex [59], N. vespilloides [19], and Z. nevadensis [15] were downloaded from the Short Read Archive (SRA) using accessions SRR1270128, SRR1270129, SRR445803-4, SRR027157-9, SRR1552830, SRR2017555, and SRR3139749. Thus, DNA methylation was investigated for six insects from six different orders spread evenly across developmental groups, and a crustacean outgroup. Additionally, DNA methylation (and gene expression, see "Materials and methods" section "RNA-seq and differential expression analysis") changes following posttranscriptional knockdown of Dnmt3 in A. mellifera were investigated [30]. WGBS data were aligned to each species respective genome assembly using the methylpy pipeline [60]. In brief, reads were trimmed of sequencing adapters using Cutadapt v1.9 [61] and then mapped to both a converted forward strand (cytosines to thymines) and converted reverse strand (guanines to adenines) using Bowtie v1.1.1 [62]. Reads that mapped to multiple locations, and clonal reads were removed.
Weighted DNA methylation was calculated for CG sites by dividing the total number of aligned methylated reads by the total number of methylated plus unmethylated reads [63]. For genic metaplots, the gene body (start to stop codon), 1000 base pairs (bp) upstream, and 1000 bp downstream were divided into 20 windows proportional windows based on sequence length (bp). Weighted DNA methylation was calculated for each window and then plotted in R v3.2.4 (https ://www.rproje ct.org/). CG sequence context enrichment for each gene was determined through a binomial test followed by Benjamini-Hochberg false discovery rate [64,65]. A background mCG level was determined from all coding sequences, which was used as a threshold in determining significance with a false discovery rate (FDR) correction. Genes were classified as CG-methylated if they had reads mapping to at least 20 reads mapping to 20 CG sites and a q < 0.05. Using a binomial test can lead to false negatives-highly CG-methylated genes that are classified as unmethylated (UM)-due to a low number of statistically CG-methylated sites (Additional file 1: Table S4). Genes classified as unmethylated, but had a mCG level greater than the lowest CG-methylated gene were dropped from future analyses.

RNA-seq and differential expression analysis
RNA-seq libraries for RNA extracted from ovaries of three biological O. fasciatus ds-dnmt1-1 and control replicates at 11 days post-injection were constructed using Illumina TruSeq Stranded RNA LT Kit (Illumina, San Diego, CA) following the manufacturer's instructions with limited modifications. RNA from ovaries of an additional three biological O. fasciatus ds-dnmt1-1 and control replicates, and three biological O. fasciatus ds-dnmt1-1 and control replicates from gut, head, and thorax were extracted. The starting quantity of total RNA was adjusted to 1.3 µg, and all volumes were reduced to a third of the described quantity. Libraries were single-end 75 bp sequenced on an Illumina NextSeq 500 machine. RNA-seq data for O. fasciatus ds-dnmt1-1 and control can be found on GEO under accession GSE109199. Previously published RNA-seq data for A. mellifera [66] and Z. nevadensis [15] were downloaded from the SRA using accessions SRR893029, SRR893030, SRR893035, SRR893036, SRR2954345, and SRR3139740.
Differentially expressed genes (DEGs) between ds-dnmt1-1 and control libraries were determined using edgeR v3.20.1 [69] implemented in R v3.2.4 (https :// www.r-proje ct.org/). Genes were retained for DEG analysis if they possessed a counts per million (CPM) ≥ 1 in at least ≥ 2 libraries. Significance was determined using the glmQLFTest function, which uses empirical Bayes quasilikelihood F tests. Parameter settings were determined following best practices for DEG analysis as described by [70]. Gene expression metrics for A. mellifera, Z. nevadensis, and O. fasciatus ds-dnmt1-1 and control are located in Additional file 1: Tables S4-S8, respectively.

Transposable element (TE) annotation and expression
TEs were identified using RepeatMasker v4.0.5 (http:// www.repea tmask er.org) provided with the invertebrate repeat library from Repbase (http://www.girin st.org/repba se/) (arguments: -lib < Repbase invertebrate library > -no_is -engine wublast -a -inv -x -gff. Following RepeatMasker, neighboring TEs of the same type were collapsed into a single locus within the outputted GFF. The unmodified GFF is located in Additional file 1: Table S10. To quantify expression from TEs, RNA-seq libraries from O. fasciatus ds-dnmt1-1 and control were independently combined and mapped to the O. fasciatus v1.1 reference genome assembly [32] using Bowtie2 v2.2.9 [73] with the following arguments: -sensitive. Mapped reads overlapping with the top ten most abundant TEs of ≥ 500 bp in length were identified using the intersect command in BEDTools suite v 2.26.0 [74]. TE expression is quantified as reads per kilobase per million mapped reads (RPKM) for each intersected TE type by counting the number reads and dividing by the mapped library read number in millions. Significance in expression of TEs between ds-dnmt1-1 and control tissues was assessed using the Mann-Whitney test with the alternative hypothesis set to "greater" in R v3.2.4 (https ://www.rproje ct.org/).

Additional file
Additional file 1: Fig. S1 Fig. S2. Assessment of RNAi treatment targeting S-adenosyl-L-methionine (AdoMet) region (ds-dnmt1-2) of Dnmt1 produces similar results as ds-dnmt1-1. a Assessment of RNAi treatment targeting Dnmt1 using qRT-PCR demonstrates successful reduction in transcription in ovaries compared to control. b Whole ovaries from ds-dnmt1-2 females removed 12-14 days postinjection. c Control follicular epithelium nuclei. d ds-dnmt1-1 follicular epithelium nuclei. e ds-dnmt1-2 follicular epithelium nuclei. For c-e scale bar corresponds to 100 µm. Fig. S3. DNA methylation consequences following posttranscriptional knockdown of Dnmt1 are restricted to ovaries. a Assessment of RNAi treatment targeting Dnmt1 using qRT-PCR demonstrates successful reduction in transcription compared to control across all tissues sampled. Colored dots indicate independent biological replicates. b Genome-wide CG methylation level across tissues sampled. Numbers at the top of each bar correspond to a unique individual identifier.  Fig. S5. Expression of Dnmt3 is not affected by dsRNA targeting Dnmt1 (ds-dnmt1-1). Dots indicate mean expression level, and error bars indicate standard error of the mean. Fig. S6. mCG in A. mellifera and Z. nevadensis is not associated with transcription. a Gene expression level for deciles of increasing mCG (1-10) and unmethylated genes (UM). Error bars represent 95% confidence interval of the mean. b Regression of gene expression against a continuous measure of mCG among all genes with > 0 FPKM and weighted mCG. Raw p values are provided for each regression, and significance (S) or non-significance (NS) is indicated in brackets following Bonferroni correction. c. The relationship between change in gene expression measured as the log 2 fold-change of FPKM between A. mellifera dnmt3 knockdown and GFP control, and DNA methylation measured as the difference between A. mellifera dnmt3 knockdown and GFP control mCG levels. For gene expression, FPKM values were averaged across libraries. Density of bins corresponds to the number of genes with similar changes to expression and DNA methylation. Fig.  S7. Expression changes of differentially methylated genes. A heatmap showing gene expression changes for genes that are differentially CGmethylated between O. fasciatus ds-dnmt1-1 and control. Expression was standardized by the highest value per gene per biological replicate to produce a relative fragments per kilobase of transcript per million (RFPKM) value. RFPKM were clustered using a hierarchical clustering method. Fig.  S8. Loss of mCG in O. fasciatus ovaries has a limited effect on transcription. Combinational overlap of genes that are differential CG-methylated and expressed, and similarly CG-methylated and expressed between O. fasciatus ds-dnmt1-1 and control. Gene groups: differentially methylated gene (DMG), differentially expressed gene (DEG), similarly methylated gene (SMG)/unmethylated genes (UMG), and non-differentially expressed gene (non-DEG). More stringent thresholds were used to group genes as CGmethylated or unmethylated (see "Materials and methods").  Table S1. PCR primers used to validate the presence of a single Dnmt1 ortholog in the O. fasciatus genome. Table S2. Primer sequences for use in producing template DNA for use in the MegaScript transcription kit to generate double-stranded RNAs for injection and for quantitative real-time PCR to assess expression levels of Dnmt1. Table S3. DNA methylation summary statistics for all species investigated in this study. Table S4. Output from edgeR v3.20.1 [69] with gene expression (fragments per kilobase of transcript per million mapped reads [FPKM]) for O. fasciatus ovaries ds-dnmt1-1 and control biological and technical replicates. Table S5. Gene expression (FPKM) for A. mellifera queen and drone brains. Table S6. Gene expression (FPKM) for A. mellifera dnmt3 knockdown and GFP control. Table S7. Gene expression (FPKM) for Z. nevadensis female worker at the final instar larva. Table S8. Output from edgeR v3.20.1 [69] with gene expression (FPKM) for O. fasciatus gut, head, and thorax ds-dnmt1-1 and control biological replicates. Table S9. Overlap between O. fasciatus (control) CG-methylated and unmethylated genes in none and ≥ 1 other insect species used in this study. Table S10. Output General Features File (GFF) from RepeatMasker v4.0.5 (http://www.repea tmask er.org). Table S11. Significance (p value) of TE expression between ds-dnmt1-1 and control tissues. Table S12. Gene Ontology (GO) terms for O. fasciatus v1.1 reference genome assembly [32] annotated genes. Table S13. Significantly enriched GO terms for the intersections between differentially methylated genes (DMG), similarly methylated genes (SMG), differentially expressed genes (DEG), and nondifferentially expressed genes (non-DEG). Data S1. Unaligned and aligned DNA methyltransferase protein sequences in FASTA format, and a parenthetical phylogenetic tree in nexus format estimated from the aligned DNA methyltransferase protein sequences using BEAST v2.3.2 [51].
Authors' contributions AJB, AJM, PJM, and RJS designed the study. PJM, ZS, and ECM maintained animals and performed RNAi experiments and microscopy analyses. WGBS and RNA-seq data were generated by RJS. AJB analyzed and interpreted the WGBS and RNA-seq data and wrote the manuscript. All authors read and approved the final manuscript.