5-methylcytosine in the DNA of the polytene chromosomes of the diptera Sciara coprophila, Drosophila melanogaster and D. persimilis

High contiguity de novo genome assembly and DNA modification analyses for the fungus fly, Sciara coprophila, using single-molecule sequencing

BACKGROUND

The lower Dipteran fungus fly, Sciara coprophila, has many unique biological features that challenge the rule of genome DNA constancy. For example, Sciara undergoes paternal chromosome elimination and maternal X chromosome nondisjunction during spermatogenesis, paternal X elimination during embryogenesis, intrachromosomal DNA amplification of DNA puff loci during larval development, and germline-limited chromosome elimination from all somatic cells. Paternal chromosome elimination in Sciara was the first observation of imprinting, though the mechanism remains a mystery. Here, we present the first draft genome sequence for Sciara coprophila to take a large step forward in addressing these features.

RESULTS

We assembled the Sciara genome using PacBio, Nanopore, and Illumina sequencing. To find an optimal assembly using these datasets, we generated 44 short-read and 50 long-read assemblies. We ranked assemblies using 27 metrics assessing contiguity, gene content, and dataset concordance. The highest-ranking assemblies were scaffolded using BioNano optical maps. RNA-seq datasets from multiple life stages and both sexes facilitated genome annotation. A set of 66 metrics was used to select the first draft assembly for Sciara. Nearly half of the Sciara genome sequence was anchored into chromosomes, and all scaffolds were classified as X-linked or autosomal by coverage.

CONCLUSIONS

We determined that X-linked genes in Sciara males undergo dosage compensation. An entire bacterial genome from the Rickettsia genus, a group known to be endosymbionts in insects, was co-assembled with the Sciara genome, opening the possibility that Rickettsia may function in sex determination in Sciara. Finally, the signal level of the PacBio and Nanopore data support the presence of cytosine and adenine modifications in the Sciara genome, consistent with a possible role in imprinting.

Comparative cytogenetic study on two species of the genus Entedon Dalman, 1820 (Hymenoptera, Eulophidae) using DNA-binding fluorochromes and molecular and immunofluorescent markers

Karyotypes of Entedon cionobius Thomson, 1878 and Entedon cioni Thomson, 1878 (Hymenoptera: Eulophidae) were studied using DNA-binding ligands with different base specificity (propidium iodide, chromomycin A3, methyl green and DAPI; all these ligands, except for the last one, were used for the first time in parasitic wasps), C-banding, fluorescence in situ hybridization (FISH) with a 45S rDNA probe and 5-methylcytosine immunodetection. Female karyotypes of both species contain five pairs of relatively large metacentric chromosomes and a pair of smaller acrocentric chromosomes (2n = 12). As in many other Hymenoptera, males of both Entedon Dalman, 1820 species have haploid chromosome sets (n = 6). Fluorochrome staining revealed chromosome-specific banding patterns that were similar between the different fluorochromes, except for the CMA3- and PI-positive and DAPI-negative band in the pericentromeric regions of the long arms of both acrocentric chromosomes. The obtained banding patterns were virtually identical in both species and allowed for the identification of each individual chromosome. C-banding revealed a pattern similar to DAPI staining, although centromeric and telomeric regions were stained more intensively using the former technique. FISH detected a single rDNA site in the same position on the acrocentric chromosomes as the bright CMA3-positive band. Immunodetection of 5-methylcytosine that was performed for the first time in the order Hymenoptera revealed 5-methylcytosine-rich sites in the telomeric, centromeric and certain interstitial regions of most of the chromosomes.

Distribution of 5-methylcytosine-rich regions in the metaphase chromosomes of Vicia faba

The DNA methylation pattern of Vicia faba metaphase chromosomes was examined with a specific monoclonal antibody. 5-methylcytosine (5-mC) residues are present in different chromosomal sites, and are particularly abundant in telomeric and/or subtelomeric regions and in certain intercalary bands. Chromosomal localization of methylated regions enables a better knowledge of the lengthwise differentiation of this chromosome complement. Our results also indicate that there may be differences in monoclonal antibody binding between corresponding regions of homologous chromosomes in V. faba. This behaviour is detectable in specific regions with different frequencies. The data support results previously obtained for Allium cepa metaphase chromosomes using the same monoclonal antibody.

Changes in the methylation of amplified esterase DNA during loss and reselection of insecticide resistance in peach-potato aphids, Myzus persicae

Insecticide resistance in peach-potato aphids, Myzus persicae, results from the amplification of genes encoding an esterase that hydrolyses and sequesters insecticides. Resistance is normally stable, but highly resistant aphid clones sometimes lose resistance when insecticidal selection pressure is removed. This loss of resistance, termed reversion, arises from a loss of elevated esterase enzyme through transcriptional control, i.e. without loss of the amplified esterase DNA sequences. We have shown that loss of the elevated enzyme occurred simultaneously with loss of methylation at CCGG sites in the amplified DNA sequences. During reselection of resistance in these revertant clones, enzyme levels increased, but there was no corresponding return of methylation to DNA sequences. Thus, although DNA methylation is closely correlated with expression of the amplified esterase genes during reversion, it may not be a factor in the reverse process.

Direct visualization of the sites of DNA methylation in human, and mosquito chromosomes

Human and mosquito fixed chromosomes were digested with restriction endonucleases that are inhibited by the presence of 5-methylcytosine in their restriction sites (Hha I, Hin PI, Hpa II), and with endonucleases for which cleavage is less dependent on the state of methylation (Taq I, Msp I). Methylation-dependent enzymes extracted low DNA amounts from human chromosomes, while methylation-independent enzymes extracted moderate to high amounts of DNA. After DNA demethylation with 5-azacytidine the isoschizomers Hpa II (methylation-dependent) and Msp I (methylation-independent) extracted 12-fold and 1.4-fold amounts of DNA from human chromosomes, respectively. These findings indicate that human DNA has a high concentration of Hpa II and Msp I restriction sites (CCGG), and that the internal C of this sequence is methylated in most cases, while the external cytosine is methylated less often. All the enzymes tested released moderate amounts of DNA from mosquito chromosomes whether or not the DNA was demethylated with 5-azacytidine. Hpa II induced banding in the centromere chromosome regions. After demethylation with 5-azacytidine this banding disappeared. Mosquito DNA has therefore, moderate to high frequencies of nonmethylated CpG duplets. The only exception is the centromeric DNA, in which the high levels of C methylation present produce cleavage by Hpa II and the appearance of banding. Centromere regions of human chromosomes 1 have a moderately low concentration of Hpa II-Msp I restriction sites.

Distribution of 5-methylcytosine-rich regions in the polytene chromosomes of Phaseolus coccineus embryo suspensor as shown by the immunoperoxidase technique

5-Methylcytosine (5-mC) has been visualized in polytene chromosomes of Phaseolus coccineus, scarlet bean using specific antibodies to 5-mC and the immunoperoxidase technique. The results obtained indicate that most heterochromatic regions are methylated, even though the frequency of methylation is highly variable and sometimes low. A preferential binding of anti-5-mC to centromeric heterochromatic blocks was observed. Comparison between anti-5-mC binding and the results of hybridization with highly repetitive DNA and satellite DNA shows, moreover, that centrometric heterochromatic regions hybridize in particular with both DNAs. This finding is consistent with the fact that repetitive DNA and satellite DNA are methylated to a considerably greater extent than main band DNA, in line with many data to be found in the literature. The binding pattern of anti-5-mC that we observed also suggests that methylation does not occur in all classes of repetitive DNA. The high variability of band methylation frequency is discussed in relation to a possible characteristic DNA composition of the band.

Caenorhabditis elegans DNA does not contain 5-methylcytosine at any time during development or aging

DNA, isolated from age-synchronous senescent populations of Caenorhabditis elegans has been quantitatively and qualitatively analyzed for the presence of 5-methylcytosine. High performance liquid chromatography on two wild-type and several mutant strains of C. elegans failed to detect any 5-methylcytosine. The restriction endonuclease isoschizomers, HpaII and MspI, were used to digest genomic DNA after CsCl purification and failed to detect any 5' cytosine methylation at any age. We conclude that C. elegans does not contain detectable (0.01 mole percent) levels of 5-methylcytosine.

Genomic 5-methylcytosine determination by 32P-postlabeling analysis

A simple and sensitive method for the quantitation of 5-methyldeoxycytidine in DNA has been developed by the adaptation of the Randerath 32P-postlabeling technique. Nucleic acids were digested to 3'-monophosphate nucleotides, which were converted to 32P-labeled 3',5'-bisphosphate nucleotides, the 3'-phosphate was cleaved by the action of nuclease P1, and the resultant 5'-[32P]-monophosphate nucleotides were separated by two-dimensional thin-layer chromatography. Less than 1 microgram of DNA was required for the precise quantitation of 5-methyldeoxycytidine content to a detectable limit of 0.01% of the total cytidine residues methylated. The genomic 5-methyldeoxycytidine content may thus be quantitated in tissue samples, small or selective cell populations, senescing or terminally differentiating cells, or DNA from any source. We report here, for the first time, the genomic 5-methyldeoxycytidine content of normal human bronchial epithelial and normal human pulmonary mesothelial cells. The chromatographic separation of all of the normal and some of the rare monophosphate deoxyribonucleotides and ribonucleotides has been characterized. Thus, 5-bromodeoxyuridine and the RNA contamination of DNA or the DNA contamination of RNA can also be quantitated during the same analysis.

Antibodies to DNA

Antibodies that recognize specific conformational variations of DNA structure provide sensitive reagents for testing the extent to which such conformational heterogeneity occurs in nature. A most dramatic recent example has been the development and application of antibodies to left-handed Z-DNA. They provided the first identification of Z-DNA in fixed nuclei and chromosomes, and of DNA sequences that form Z-DNA under the influence of supercoiling. Antibodies have also been induced by chemically modified DNA and by synthetic polydeoxyribonucleotides that differ from the average B-DNA structure. These antibodies recognize only the features that differ from native DNA. In most experiments, native DNA itself is not immunogenic. Antibodies that do react with native DNA occur in sera of patients with autoimmune disease, but even monoclonal anti-DNA autoantibodies usually react with other polynucleotides as well. Anti-DNA antibodies, especially those of monoclonal origin, provide a model for the study of protein-nucleic acid recognition.

Absence of detectable 5-methylcytosine in DNA of embryos of the brine shrimp, Artemia

DNA from three developmental stages of the brine shrimp, Artemia was found to be undermethylated compared to DNA from calf thymus and salmon sperm. Using high-performance liquid chromatography we found that DNA hydrolysates from both dormant cysts and prefeeding nauplii contain less than 1 residue of 5-methylcytosine per 59 kilobases, placing Artemia DNA in the "insect type" category. The absence of detectable 5-methylcytosine in DNA of developing Artemia supports the view that methylation status alone cannot account for regulation of transcription in protostomes, and that DNA methylation may be more common among deuterostomes.

Eukaryotic DNA methylation

Eukaryotic genomes contain 5-methylcytosine (5mC) as a rare base.5mC arises by postsynthetic modification of cytosine and occurs, at least in animals, predominantly in the dinucleotide CpG. The base is not distributed randomly in these genomes but conforms to a pattern. This pattern varies between taxa but appears to be inherited in a semi-conservative fashion. At the level of the genome, gross changes in the level of DNA methylation have been noted. This has encouraged speculation that the modification may play a role in cellular differentiation. Tissue-specific patterns of DNA methylation, predicted by various models of differentiation, have been found for most vertebrate genes so far examined. A correlation has emerged between the undermethylation of these regions and their transcription, but this is not always the case. While data for eukaryotic viral sequences are less equivocal, studies of this kind cannot in isolation distinguish between undermethylation being a cause or a consequence of gene activity. If it were a cause, it is probable that the demethylation of specific CpG sites would be a necessary yet not a sufficient condition for transcription to occur. The introduction of artificially methylated DNA sequences into individual eukaryotic cells by microinjection or transformation may provide the means to elucidate these questions in the future. In the meantime, the study of eukaryotic DNA methylation promises to contribute much to our understanding of the regulation of gene expression in these organisms.

5-Methylcytosine in eukaryotic DNA

A small portion of the cytosine residues in the DNA of higher eukaryotes as well as in that of many lowe eukaryotes if methylated. The resulting 5-methylcytosine residues occur in specific in the DNA, usually adjacent to guanine residues on the 3' side. This methylation of eukaryotic DNA has been proposed to function in many ways, including control of transcription, maintenance of chromosome structure, repair of DNA, establishment of preferred sites for mutation, oncogenic transformation, and, in certain systems, protection of DNA against enzymatic degradation.

The two-dimensional restriction analysis of Drosophila DNAs: males and females

Two-dimensional restriction analysis has been applied to terminally labeled restriction segments from adult male and female flies that were the progeny of ten successive single-pair, brother-sister matings. Many of the 2500-4000 different bands that can be seen on autoradiographs are of widely different density, a fact that suggests the superposition of multiple restriction segments of similar sequence. The number of bands observed is significantly fewer than expected. Nearly all of the bands seen in male-DNA displays appear equivalent to those seen in female-DNA displays, but a few sex differences (male-dense and female-dense bands) can be identified--much fewer than expected. When the HpaII-MspI test was applied to fractionated BamHI fragments, between 900 and 1800 discrete C-C-G-G sites were found to be cleaved in a completely equivalent fashion by each isoschizomer.

The organization of DNA in the mitotic and polytene chromosomes of Sciara corprophila

The organization of DNA in the mitotic metaphase and polytene chromosomes of the fungus gnat, Sciara coprophila, has been studied using base-specific DNA ligands, including anti-nucleoside antibodies. The DNA of metaphase and polytene chromosomes reacts with AT-specific probes (quinacrine, DAPI, Hoechst 33258 and anti-adenosine) and to a somewhat lesser extent with GC-specific probes (mithramycin, chromomycin A3 and anti-cytidine). In virtually every band of the polytene chromosomes chromomycin A3 fluorescence is almost totally quenched by counterstaining with the AT-specific ligand methyl green. This indicates that GC base pairs in most bands are closely interspersed with AT base pairs. The only exceptions are band IV-8A3 and the nucleolus organizer on the X. In contrast, quinacrine and DAPI fluorescence in every band is only slightly quenched by counterstaining with the GC-specific ligand actinomycin D. Thus, each band contains a moderate proportion of AT-rich DNA sequences with few interspersed GC base pairs. - The C-bands in mitotic and polytene chromosomes can be visualized by Giemsa staining after differential extraction of DNA and those in polytene chromosomes by the use of base-specific fluorochromes or antibodies without prior extraction of DNA. C-bands are located in the centromeric region of every chromosome, and the telomeric region of some. The C-bands in the polytene chromosomes contain AT-rich DNA sequences without closely interspered GC base pairs and lack relatively GC-rich sequences. However, one C-band in the centromeric region of chromosome IV contains relatively GC-rich sequences with closely interspersed AT base pairs. - C-bands make up less than 1% of polytene chromosomes compared to nearly 20% of mitotic metaphase chromosomes. The C-bands in polytene chromosomes are detectable with AT-specific or GC-specific probes while those in metaphase chromosomes are not. Thus, during polytenization there is selective replication of highly At-rich and relatively GC-rich sequences and underreplication of the remainder of the DNA sequences in the constitutive heterochromatin.