Developmental progression of DNA puffs in Sciara coprophila: amplification and transcription

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.

Beyond DNA puffs: What can we learn from studying sciarids?

Members of the Sciaridae family attracted the interest of researchers because of the demonstration that the DNA puff regions, which are formed in the salivary gland polytene chromosomes at the end of the fourth larval instar, constitute sites of developmentally regulated gene amplification. Besides contributing to a deeper understanding of the process of gene amplification, the study of sciarids has also provided important insights on other biological processes such as sex determination, programmed cell death, insect immunity, telomere maintenance, and nucleolar organizing regions (NOR) formation. Open questions in sciarids include among others, early development, the role of noncoding RNAs in gene amplification and the relationship between gene amplification and transcription in DNA puff forming regions. These and other questions can now be pursued with next generation sequencing techniques and experiments using RNAi experiments, since this latter technique has been shown to be feasible in sciarids. These new perspectives in the field of sciarid biology open the opportunity to consolidate sciarid species as important emerging models. genesis 54:361-378, 2016. © 2016 Wiley Periodicals, Inc.

The Balbiani Ring Story: Synthesis, Assembly, Processing, and Transport of Specific Messenger RNA-Protein Complexes

Eukaryotic gene expression is the result of the integrated action of multimolecular machineries. These machineries associate with gene transcripts, often already nascent precursor messenger RNAs (pre-mRNAs). They rebuild the transcript and convey properties allowing the processed transcript, the mRNA, to be exported to the cytoplasm, quality controlled, stored, translated, and degraded. To understand these integrated processes, one must understand the temporal and spatial aspects of the fate of the gene transcripts in relation to interacting molecular machineries. Improved methodology is necessary to study gene expression in vivo for endogenous genes. A complementary approach is to study biological systems that provide exceptional experimental possibilities. We describe such a system, the Balbiani ring (BR) genes in polytene cells in the dipteran Chironomus tentans. The BR genes, along with their pre-mRNA-protein complexes (pre-mRNPs) and mRNA-protein complexes (mRNPs), allow the visualization of intact cell nuclei and enable analyses of where and when different molecular machineries associate with and act on the BR pre-mRNAs and mRNAs.

The expanding implications of polyploidy

Polyploid cells, which contain more than two genome copies, occur throughout nature. Beyond well-established roles in increasing cell size/metabolic output, polyploidy can also promote nonuniform genome, transcriptome, and metabolome alterations. Polyploidy also frequently confers resistance to environmental stresses not tolerated by diploid cells. Recent progress has begun to unravel how this fascinating phenomenon contributes to normal physiology and disease.

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.

Developmental control of the DNA replication and transcription programs

Polyploid or polytene cells, which have more than 2C DNA content, are widespread throughout nature and present in most differentiated Drosophila tissues. These cells also can display differential replication, that is, genomic regions of increased or decreased DNA copy number relative to overall genomic ploidy. How frequently differential replication is used as a developmental strategy remains unclear. Here, we use genome-wide array-based comparative genomic hybridization (aCGH) to profile differential DNA replication in isolated and purified larval fat body and midgut tissues of Drosophila, and we compare them with recent aCGH profiles of the larval salivary gland. We identify sites of euchromatic underreplication that are common to all three tissues and others that are tissue specific. We demonstrate that both common and tissue-specific underreplicated sites are dependent on the Suppressor of Underreplication protein, SUUR. mRNA-seq profiling shows that whereas underreplicated regions are generally transcriptionally silent in the larval midgut and salivary gland, transcriptional silencing and underreplication have been uncoupled in the larval fat body. In addition to revealing the prevalence of differential replication, our results show that transcriptional silencing and underreplication can be mechanistically uncoupled.

Indirect immune detection of ecdysone receptor (EcR) during the formation of DNA puffs in Bradysia hygida (Diptera, Sciaridae)

Gene amplification occurs in Bradysia hygida salivary glands, at the end of the fourth larval instar. The hormone 20-hydroxyecdysone (20E) triggers this process, which results in DNA puff formation. Amplified genes are activated in two distinct groups. The activity of the first group is dependent on high levels of 20E, while the second group needs low hormone levels. Consequently, the salivary glands of B. hygida constitute an interesting biological model to study how 20E, and its receptors, affect gene amplification and activity. We produced polyclonal antibodies against B. hygida EcR (BhEcR). In western blots a polypeptide of about 66 kDa was detected in salivary gland extracts. The antibodies were also used for indirect immune-localization of BhEcR in polytene chromosomes. RNA-polymerase II was also immune-detected. We did not detect the receptor in chromosome C where the first and second groups of DNA puffs form during DNA puff anlage formation, but it was present during puff expansion. During the active phase of both groups of DNA puffs, RNA polymerase II co-localized with BhEcR. After puff regression, these antigens were not detected. Apparently, EcR plays a direct role in the transcription of amplified genes, but its role in gene amplification remains enigmatic.

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.

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.

Gene amplification as a developmental strategy: isolation of two developmental amplicons in Drosophila

Gene amplification is known to be critical for upregulating gene expression in a few cases, but the extent to which amplification is utilized in the development of diverse organisms remains unknown. By quantifying genomic DNA hybridization to microarrays to assay gene copy number, we identified two additional developmental amplicons in the follicle cells of the Drosophila ovary. Both amplicons contain genes which, following their amplification, are expressed in the follicle cells, and the expression of three of these genes becomes restricted to specialized follicle cells late in differentiation. Genetic analysis establishes that at least one of these genes, yellow-g, is critical for follicle cell function, because mutations in yellow-g disrupt eggshell integrity. Thus, during follicle cell differentiation the entire genome is overreplicated as the cells become polyploid, and subsequently specific genomic intervals are overreplicated to facilitate gene expression.

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 is differentially expressed in various tissues of Bradysia hygida late larvae and constitutively transcribed in transgenic Drosophila

We extended the characterization of the DNA puff BhB10-1 gene of Bradysia hygida by showing that, although its mRNA is detected only at the end of the fourth larval instar, BhB10-1 expression is not restricted to the salivary gland, the tissue in which this gene is amplified. Different amounts of BhB10-1 mRNA were detected in other larval tissues such as gut, Malpighian tubules, fat body, brain and cuticle, suggesting that this gene is expressed differentially in the various tissues analyzed. Analysis of transgenic Drosophila carrying the BhB10-1 transcription unit and flanking sequences revealed that the tested fragment promotes transcription in a constitutive manner. We suggest that either cis-regulatory elements are missing in the transgene or factors that temporally regulate the BhB10-1 gene in B. hygida are not conserved in Drosophila.

Mapping of intrinsic bent DNA sites in the upstream region of DNA puff BhC4-1 amplified gene

We have identified bent DNA sites in the distal and proximal DNA puff BhC4-1 amplified gene promoter region of Bradysia hygida. The 2D modeling of the 3D DNA path and the ENDS ratio values calculated in this promoter region resulted in the identification of ten pronounced bent sites named BhC4B - 9 to + 1. The 1847 bp fragment (- 3697 to - 1850) in relation to the transcription start site shows multiple bending sites, BhC4B - 9 to BhC4B - 4, with periodicity approximately 300 bp. The analysis of the other identified bent region, starting at position - 957, reveals that the BhC4B + 1 bent site colocalizes with the putative BhC4-1 minimal promoter. The sequence analysis of bent site BhC4B - 4 shows a distribution of dA*dT at approximately 10 bp intervals between the middle of each tract, but intervals with more than one turn, approximately 20 bp, two helix turns, were detected in the other bent sites described here. The bent sites BhC4B - 6 and BhC4B - 4, contain two consensus sequences, with 60 bp each. The apparent molecular weight of fragments in the BhC4-1 promoter region were estimated in agarose gels and compared with the data obtained in polyacrylamide gels without and with ethidium bromide. The mobility reduction ratios (R-values) were determined, and a high R-value, 1.80, for a 1215 bp fragment in the distal promoter region and a 1.23 significant R-value for a 662 bp fragment in the proximal segment were found. To further analyze the predicted bent DNA sites in these fragments, the 2D trajectories of the 3D DNA path and other parameters, AT percentage, roll angle, ENDS ratio and DeltaG, were determined. The role of these bent sites in the BhC4-1 transcription regulation is discussed.

A dual role of 20-hydroxyecdysone in the control of protein synthesis related to DNA puff activity in the anterior region of Bradysia hygida (Diptera, Sciaridae) salivary gland

During the last 30 h of the larval stage, the salivary glands of Bradysia hygida show the amplification of some genes, resulting in the formation of two successive groups of DNA puffs, which direct the synthesis of two different sets of polypeptides. Incubation of anterior (S1) salivary gland regions, at age E7, beginning of first group of DNA puffs activity, in culture medium for 2 to 10 h results in a decrease in the synthesis of the polypeptides characteristic of this period. However, during subsequent incubation (from E7 to E7+12 h-24 h), when the second group of DNA puffs is active, S1 regions were able to synthesize some polypeptides characteristic of this period. The role of 20-OH ecdysone was studied, in vitro and in vivo, during these two periods of protein synthesis in S1 regions. The presence of the hormone was shown to be necessary to maintain, in vitro, the synthesis of the first set of polypeptides and was strongly inhibitory, in vitro and in vivo, to the synthesis of the second set of polypeptides. Thus, it is likely that the activity of the two distinct groups of DNA puffs is under opposite 20-OH-ecdysone control mechanisms.

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.

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.

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.