Molecular characterization of DNA puff II/9A genes in Sciara coprophila

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.

Anatomy and evolution of a DNA replication origin

DNA amplification occurs at the DNA puff II/9A locus in the fungus fly Sciara coprophila. As a foundation to study the molecular mechanism for the initiating events of II/9A DNA re-replication, we have sequenced 14 kb spanning a DNase hypersensitive site (DHS) upstream of the 1 kb amplification origin and through transcription units II/9-1 and II/9-2 downstream of the origin. These elements are annotated as well as the ORC binding site at the origin and the transition point (TP) between continuous and discontinuous DNA syntheses that marks the origin of bidirectional replication at the nucleotide level. A 9 bp motif found at the TP is repeated near the other end of the 1 kb ORI and may identify a putative second TP. The steroid hormone ecdysone induces DNA amplification as well as transcription and puffing at locus II/9A. Within the 14 kb, several matches to the ecdysone response element (EcRE) consensus sequence were identified, including some in the amplification origin region. EcRE O-P is at a central axis of a remarkable symmetry, equidistant to the TPs that are themselves equidistant to EcRE O-1 and EcRE O-2. DNA sequence alterations have occurred throughout the II/9A region in a newly discovered polymorphism (#2). Polymorphism #2 is not specific to developmental stage, sex, or tissue, and it does not impair DNA amplification. The DHS, both 9 bp TP sequences, and EcREs O-1, O-P, and O-2 are conserved between the polymorphism #1 and #2 sequences, suggesting their functional importance and retention during evolutionary selection. Moreover, a 72 bp sequence in the Sciara DHS at DNA puff II/9A is conserved in DNA puff C-3 of Rhynchosciara americana. Comparisons are discussed between the Sciara II/9A amplicon and the chorion locus amplicon on the third chromosome of Drosophila.

The ecdysone receptor (ScEcR-A) binds DNA puffs at the start of DNA amplification in Sciara coprophila

The steroid hormone ecdysone induces DNA amplification and subsequent DNA puff formation in late fourth larval instar salivary gland polytene chromosomes of the fungus fly, Sciara coprophila. Previous in vitro studies on DNA puff II/9A in Sciara demonstrated that the ecdysone receptor (ScEcR-A) efficiently binds an ecdysone response element adjacent to the origin recognition complex binding site within the II/9A amplification origin, implying a role for ScEcR-A in amplification. Here, we extrapolate the molecular details from locus II/9A to the rest of the genome using immunofluorescence with a ScEcR-A-specific antibody. ScEcR-A binds all DNA puff sites just as amplification begins and persists throughout the processes of amplification, transcription, and puffing. Ecdysone injections into pre-amplification stage larvae prematurely induce both DNA amplification and ScEcR-A binding to DNA puff sites. These data are consistent with a direct role for ScEcR-A in DNA amplification.

Isolation and characterization of the ecdysone receptor and its heterodimeric partner ultraspiracle through development in Sciara coprophila

Regulation of DNA replication is critical, and loss of control can lead to DNA amplification. Naturally occurring, developmentally regulated DNA amplification occurs in the DNA puffs of the late larval salivary gland giant polytene chromosomes in the fungus fly, Sciara coprophila. The steroid hormone ecdysone induces DNA amplification in Sciara, and the amplification origin of DNA puff II/9A contains a putative binding site for the ecdysone receptor (EcR). We report here the isolation, cloning, and characterizing of two ecdysone receptor isoforms in Sciara (ScEcR-A and ScEcR-B) and the heterodimeric partner, ultraspiracle (ScUSP). ScEcR-A is the predominant isoform in larval tissues and ScEcR-B in adult tissues, contrary to the pattern in Drosophila. Moreover, ScEcR-A is produced at amplification but is absent just prior. We discuss these results in relation to the model of ecdysone regulation of DNA amplification.

Functional characterization of the sciarid BhC4-1 core promoter in transgenic Drosophila

Background

Core promoters are cis-regulatory modules to which bind the basal transcriptional machinery and which participate in the regulation of transcription initiation. Although core promoters have not been extensively investigated through functional assays in a chromosomal context, the available data suggested that the response of a given core promoter might vary depending on the promoter context. Previous studies suggest that a (-57/+40) fragment constitutes the core promoter of the BhC4-1 gene which is located in DNA puff C4 of the sciarid fly Bradysia hygida. Here we tested this (-57/+40) fragment in distinct regulatory contexts in order to verify if promoter context affects its core promoter activity.

Results

Consistent with the activity of a core promoter, we showed that in the absence of upstream regulatory sequences the (-57/+40) fragment drives low levels of reporter gene mRNA expression throughout development in transgenic Drosophila. By assaying the (-57/+40) fragment in two distinct regulatory contexts, either downstream of the previously characterized Fbp1 enhancer or downstream of the UAS element, we showed that the BhC4-1 core promoter drives regulated transcription in both the germline and in various tissues throughout development. Furthermore, the use of the BhC4-1 core promoter in a UAS construct significantly reduced salivary gland ectopic expression in third instar larvae, which was previously described to occur in the context of the GAL4/UAS system.

Conclusions

Our results from functional analysis in transgenic Drosophila show that the BhC4-1 core promoter drives gene expression regardless of the promoter context that was assayed. New insights into the functioning of the GAL4/UAS system in Drosophila were obtained, indicating that the presence of the SV40 sequence in the 3' UTR of a UAS construct does not preclude expression in the germline. Furthermore, our analysis indicated that ectopic salivary gland expression in the GAL4/UAS system does not depend only on sequences present in the GAL4 construct, but can also be affected by the core promoter sequences in the UAS construct. In this context, we propose that the sciarid BhC4-1 core promoter constitutes a valuable core promoter which can be employed in functional assays in insects.

Intrinsically bent DNA sites in the Drosophila melanogaster third chromosome amplified domain

Bent DNA sites promote the curvature of DNA in both eukaryotic and prokaryotic chromosomes. Here, we investigate the localization and structure of intrinsically bent DNA sites in the extensively characterized Drosophila melanogaster third chromosome DAFC-66D segment (Drosophila amplicon in the follicle cells). This region contains the amplification control element ACE3, which is a replication enhancer that acts in cis to activate the major replication origin ori-beta. Through both electrophoretic and in silico analysis, we have identified three major bent DNA sites in DAFC-66D. The bent DNA site (b1) is localized in the ACE3 element, whereas the other two bent DNA sites (b2 and b3) are localized in the ori-beta region. Four additional bent DNA sites were identified in the intron of the S18 gene and near the TATA box of the S15, S19, and S16 genes. The identification of DNA bent sites in genomic regions previously characterized as functionally relevant for DNA amplification further supports a function for DNA bent sites in DNA replication in eukaryotes.

The DNA puff 4C expresses a salivary secretion protein in Trichosia pubescens (Diptera; Sciaridae)

DNA puffs are genomic regions of polytene chromosomes that undergo developmentally controlled DNA amplification and transcription in salivary glands of sciarid flies. Here, we tested the hypothesis that DNA puff genes code for salivary proteins in Trichosia pubescens. To do that, we generated antibodies against saliva and immunoscreened a cDNA library made from salivary glands. We isolated clones corresponding to DNA puff regions, including clone D-50 that contained the entire coding sequence of the previously isolated C4B1 gene from puff 4C. Indeed, we showed that puff 4C is a DNA puff region detecting its local transcription and its extra rounds of DNA incorporation compared to neighboring regions. We further confirmed D-50 clone identity in Western blots reacted with the anti-saliva anitiserum. We detected a recombinant protein expressed by this clone that had the expected size for a full-length product of the gene. We end with a discussion of the relationship between DNA puff genes and their products.

Ecdysone induces transcription and amplification in Sciara coprophila DNA puff II/9A

DNA replication is normally tightly regulated to ensure the production of only one copy of the genome per cell cycle. However, DNA puffs of the salivary gland giant polytene chromosomes of Sciara coprophila undergo DNA amplification during the normal course of development, overriding this control. This developmental strategy provides more template for the production of large amounts of protein needed for pupation. We have focused on DNA puff II/9A, which amplifies approximately 17-fold over the rest of the genome. Evidence presented here suggests that DNA amplification at this locus is controlled by the steroid hormone ecdysone, the master regulator of insect development. Explanted, pre-amplification stage salivary glands undergo premature amplification when incubated with ecdysone. Injection of ecdysone into pre-amplification stage larvae induces amplification. Ecdysone also induces transcription of the II/9A genes. We report the presence of a putative ecdysone response element directly adjacent to the origin recognition complex (ORC)-binding site in the II/9A origin and demonstrate that it is efficiently bound by the Sciara ecdysone receptor. These results implicate ecdysone in the regulation of DNA amplification in Sciara and suggest the ecdysone receptor may be the elusive amplification factor. This would be a new role for this transcription factor.

Analysis of expressed sequence tags from Rhynchosciara americana salivary glands

The diptera Rhynchosciara americana (sciaridae) is an important model organism in polyteny and gene amplification research, but up to now a limited amount of data regarding DNA sequences and molecular aspects of this species is available. Considering the importance of going further on the DNA puffs biological meaning, we proposed to generate EST sequences from a DNA library constructed from salivary glands. After their categorization in gene ontology terms, they were used to construct an 'electronic Northern' that represents a general view of the salivary gland metabolic status in an important phase of larval development: the spinning of communal cocoon. In this phase occurs the last polytene DNA replication cycle concomitantly with the specific loci amplification related to protein secretion.

Immunocharacterization of the DNA puff BhC4-1 protein of Bradysia hygida (Diptera: Sciaridae)

The DNA puff BhC4-1 gene is amplified and highly expressed in the salivary gland of Bradysia hygida late larvae. Using affinity-purified polyclonal antibodies we have identified the product of the BhC4-1 gene as a 43 kDa polypeptide which is present in extracts of salivary glands from late fourth instar larvae and in the corresponding gland secretion, but not in glands from earlier stages. We also demonstrate that this protein is produced mainly in the S1 and S3 regions of the salivary gland, where BhC4-1 amplification levels are more pronounced and larger amounts of mRNA are produced. By immunoelectron microscopy the BhC4-1 protein was detected in secretory granules of the S1 and S3 regions, and localized in fibrous structures present in the saliva.

Identification of regulatory regions in the DNA puff BhC4-1 promoter

The mechanisms that control DNA puff BhC4-1 expression in the salivary gland of sciarid late larvae have been shown to be conserved in Drosophila. By analysing Drosophila transformed with constructs carrying progressive deletions of the BhC4-1 promoter fragment (-3314/+40) fused to the lacZ reporter gene we show that the elements required for the correct BhC4-1-lacZ developmental regulation in prepupal salivary glands are contained in a 226 bp fragment (-186/+40). Also, interestingly, this study identified a 67 bp fragment (-253/-187) that activates BhC4-1-lacZ expression specifically in the ring gland.

Developmental changes in the Sciara II/9A initiation zone for DNA replication

Developmentally regulated initiation of DNA synthesis was studied in the fly Sciara at locus II/9A. PCR analysis of nascent strands revealed an initiation zone that spans approximately 8 kb in mitotic embryonic cells and endoreplicating salivary glands but contracts to 1.2 to 2.0 kb during DNA amplification of DNA puff II/9A. Thus, the amplification origin occurs within the initiation zone used for normal replication. The initiation zone left-hand border is constant, but the right-hand border changes during development. Also, there is a shift in the preferred site for initiation of DNA synthesis during DNA amplification compared to that in preamplification stages. This is the first demonstration that once an initiation zone is defined in embryos, its borders and preferred replication start sites can change during development. Chromatin immunoprecipitation showed that the RNA polymerase II 140-kDa subunit occupies the promoter of gene II/9-1 during DNA amplification, even though intense transcription will not start until the next developmental stage. RNA polymerase II is adjacent to the right-hand border of the initiation zone at DNA amplification but not at preamplification, suggesting that it may influence the position of this border. These findings support a relationship between the transcriptional machinery and establishment of the replication initiation zone.

The DNA puff BhB10-1 gene encodes a glycine-rich protein secreted by the late stage larval salivary glands of Bradysia hygida

We present the molecular characterization of a gene of Bradysia hygida DNA puff B10 whose temporal expression in the salivary gland correlates with the puff expansion. The transcription unit of this gene, named BhB10-1, was mapped in a 2-kb EcoRI genomic fragment that is amplified in the salivary gland of late fourth instar larvae. Its 1.3-kb transcript undergoes poly-A tail shortening during development, indicating that post-transcriptional controls as well as transcription activation are involved in the temporal regulation of the BhB10-1 gene. Analysis of the deduced amino acid sequence from the cDNA indicates that the BhB10-1 protein is a glycine-rich secretory protein. A BhB10-1-fusion protein expressed in bacteria was used to raise polyclonal antibodies. Using an immunopurified antibody, we identified the product of the DNA puff BhB10-1 gene as a 23-kDa polypeptide that is produced mainly by the salivary gland regions S1 and S3 and is present in the saliva of late larvae. This is the first direct identification of a protein encoded by a DNA puff amplified gene.

Molecular characterization of the C-3 DNA puff gene of Rhynchosciara americana

We have mapped a region of about 33 kb which includes the transcription unit of the C-3 DNA puff gene of Rhynchosciara americana. The C-3 TU and a region extending approximately 800 bp upstream of the C-3 promoter were characterized. The TU is composed of three exons and produces a 1.1-kb mRNA whose level in salivary glands increases with the expansion of the C-3 puff. The C-3 messenger appears to undergo rapid deadenylation resulting in an RNA of about 0.95 kb which can still be observed in gland cells 15 h after the puff has regressed. The 1.1-kb mRNA codes for a 32.4-kDa, predominantly alpha-helical polypeptide with three conserved parallel coiled-coil stretches. The aa composition and structure of this polypeptide suggests that it is secreted and contributes to the formation of the cocoon in which the larvae pupate. The region upstream of the promoter contains several A-rich sequences with similarity to the ACS of yeast which might have a role in the initiation of replication/amplification.

A 28-fold increase in secretory protein synthesis is associated with DNA puff activity in the salivary gland of Bradysia hygida (Diptera, Sciaridae)

When the first group of DNA puffs is active in the salivary gland regions S1 and S3 of Bradysia hygida larvae, there is a large increase in the production and secretion of new salivary proteins demonstrable by [3H]-Leu incorporation. The present study shows that protein separation by SDS-PAGE and detection by fluorography demonstrated that these polypeptides range in molecular mass from about 23 to 100 kDa. Furthermore, these proteins were synthesized mainly in the S1 and S3 salivary gland regions where the DNA puffs C7, C5, C4 and B10 are conspicuous, while in the S2 region protein synthesis was very low. Others have shown that the extent of amplification for DNA sequences that code for mRNA in the DNA puffs C4 and B10 was about 22 and 10 times, respectively. The present data for this group of DNA puffs are consistent with the proposition that gene amplification is necessary to provide some cells with additional gene copies for the production of massive amounts of proteins within a short period of time.

The complete DNA sequence of lymphocystis disease virus

Lymphocystis disease virus (LCDV) is the causative agent of lymphocystis disease, which has been reported to occur in over 100 different fish species worldwide. LCDV is a member of the family Iridoviridae and the type species of the genus Lymphocystivirus. The virions contain a single linear double-stranded DNA molecule, which is circularly permuted, terminally redundant, and heavily methylated at cytosines in CpG sequences. The complete nucleotide sequence of LCDV-1 (flounder isolate) was determined by automated cycle sequencing and primer walking. The genome of LCDV-1 is 102.653 bp in length and contains 195 open reading frames with coding capacities ranging from 40 to 1199 amino acids. Computer-assisted analyses of the deduced amino acid sequences led to the identification of several putative gene products with significant homologies to entries in protein data banks, such as the two major subunits of the viral DNA-dependent RNA polymerase, DNA polymerase, several protein kinases, two subunits of the ribonucleoside diphosphate reductase, DNA methyltransferase, the viral major capsid protein, insulin-like growth factor, and tumor necrosis factor receptor homolog.

Molecular characterization of an 18 kb segment of DNA puff C4 of Bradysia hygida (Diptera, sciaridae)

The data presented here are an extension of the molecular characterization of DNA puff C4 of Bradysia hygida. A cDNA related to a gene amplified in this puff and expressed when puff C4 expands was cloned and sequenced. Analysis of the amino acid sequence deduced from the open reading frame present in the cDNA indicate that the encoded protein is secreted and comprises mostly alpha-helical coiled-coil. An 18 kb genomic segment containing the transcription unit of this gene was also cloned and the structure and expression of the 1.4 kb mRNA was determined. Quantitative slot blot hybridization of DNA complementary to the transcription unit shows that this gene is amplified about 21 times in the salivary gland, confirming data previously obtained. Fragments upstream of the 5' end, and beyond the 3' end, of the gene transcription unit were also analysed and shown to be amplified at least eight and five times, respectively. Based on these data we discuss how amplification could occur at DNA puffs.

Analysis of an origin of DNA amplification in Sciara coprophila by a novel three-dimensional gel method

The replication origin region for DNA amplification in Sciara coprophila DNA puff II/9A was analyzed with a novel three-dimensional (3D) gel method. Our 3D gel method involves running a neutral/neutral 2D gel and then cutting out vertical gel slices from the area containing replication intermediates, rotating these slices 90 degrees to form the third dimension, and running an alkaline gel for each of the slices. Therefore, replication intermediates are separated into forks and bubbles and then are resolved into parental and nascent strands. We used this technique to determine the size of forks and bubbles and to confirm the location of the major initiation region previously mapped by 2D gels to a 1-kb region. Furthermore, our 3D gel analyses suggest that only one initiation event in the origin region occurs on a single DNA molecule and that the fork arc in the composite fork-plus-bubble pattern in neutral/neutral 2D gels does not result from broken bubbles.

Analysis of an origin of DNA amplification in Sciara coprophila by a novel three-dimensional gel method

The replication origin region for DNA amplification in Sciara coprophila DNA puff II/9A was analyzed with a novel three-dimensional (3D) gel method. Our 3D gel method involves running a neutral/neutral 2D gel and then cutting out vertical gel slices from the area containing replication intermediates, rotating these slices 90 degrees to form the third dimension, and running an alkaline gel for each of the slices. Therefore, replication intermediates are separated into forks and bubbles and then are resolved into parental and nascent strands. We used this technique to determine the size of forks and bubbles and to confirm the location of the major initiation region previously mapped by 2D gels to a 1-kb region. Furthermore, our 3D gel analyses suggest that only one initiation event in the origin region occurs on a single DNA molecule and that the fork arc in the composite fork-plus-bubble pattern in neutral/neutral 2D gels does not result from broken bubbles.

Replication initiates at a confined region during DNA amplification in Sciara DNA puff II/9A

Two independent two-dimensional (2D) gel methods were used to map an origin of replication that is developmentally regulated by the steroid hormone ecdysone, namely an origin for DNA puff amplification in the fungus fly Sciara coprophila. Initiation of replication was found to occur within a small region of no larger than 6 kb by use of the neutral/neutral 2D gel method. Neutral/alkaline 2D gel analyses support the results of the neutral/neutral 2D gels and further define within the origin region an approximately 1-kb area where the majority of replication initiates. This is the first example of an origin of replication in multicellular eukaryotes that has been mapped by 2D gels to such a small defined region. Moreover, replication can be seen by the neutral/alkaline 2D gel method to proceed bidirectionally outward from this replication origin region. These data are consistent with an onion-skin mechanism whereby multiple rounds of DNA replication initiate at a specific origin of replication for Sciara DNA puff amplification.

The promoter of DNA puff gene II/9-1 of Sciara coprophila is inducible by ecdysone in late prepupal salivary glands of Drosophila melanogaster

DNA puffs occur in Sciarid salivary gland chromosomes; they are sites of DNA amplification and intense transcription and they appear to encode secreted structural proteins needed for pupation. In this report we have used P-element transformation of Drosophila to study regulation of a Sciara DNA puff gene. We found that a 718-bp promoter fragment of DNA puff gene II/9-1 from Sciara coprophila directs expression of the bacterial reporter gene CAT in late prepupal salivary glands of transgenic Drosophila melanogaster. The identical tissue and analogous stage specificity indicate that some aspects of the ecdysone response are evolutionarily conserved between Drosophila and Sciara. When transgenic salivary glands are cultured in vitro, CAT activity is rapidly induced by ecdysone, suggesting direct control of gene expression by the ecdysone receptor. Putative stage-specific factors limit expression of the chimeric Sciara-CAT gene in transgenic Drosophila to late prepupae but not to third instar larvae when ecdysone titers are also high.