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

Chromosome-scale scaffolds of the fungus gnat genome reveal multi-Mb-scale chromosome-folding interactions, centromeric enrichments of retrotransposons, and candidate telomere sequences

BACKGROUND

The lower Dipteran fungus gnat, Bradysia (aka Sciara) coprophila, has compelling chromosome biology. Paternal chromosomes are eliminated during male meiosis I and both maternal X sister chromatids are retained in male meiosis II. Embryos start with three copies of the X chromosome, but 1-2 copies are eliminated from somatic cells as part of sex determination, and one is eliminated in the germline to restore diploidy. In addition, there is gene amplification in larval polytene chromosomes, and the X polytene chromosome folds back on itself mediated by extremely long-range interactions between three loci. These developmentally normal events present opportunities to study chromosome behaviors that are unusual in other systems. Moreover, little is known about the centromeric and telomeric sequences of lower Dipterans in general, and there are recent claims of horizontally-transferred genes in fungus gnats. Overall, there is a pressing need to learn more about the fungus gnat chromosome sequences.

RESULTS

We produced the first chromosome-scale models of the X and autosomal chromosomes where each somatic chromosome is represented by a single scaffold. Extensive analysis supports the chromosome identity and structural accuracy of the scaffolds, demonstrating they are co-linear with historical polytene maps, consistent with evolutionary expectations, and have accurate centromere positions, chromosome lengths, and copy numbers. The positions of alleged horizontally-transferred genes in the nuclear chromosomes were broadly confirmed by genomic analyses of the chromosome scaffolds using Hi-C and single-molecule long-read datasets. The chromosomal context of repeats shows family-specific biases, such as retrotransposons correlated with the centromeres. Moreover, scaffold termini were enriched with arrays of retrotransposon-related sequence as well as nucleosome-length (~ 175 bp) satellite repeats. Finally, the Hi-C data captured Mb-scale physical interactions on the X chromosome that are seen in polytene spreads, and we characterize these interesting "fold-back regions" at the sequence level for the first time.

CONCLUSIONS

The chromosome scaffolds were shown to be of exceptional quality, including loci harboring horizontally-transferred genes. Repeat analyses demonstrate family-specific biases and telomere repeat candidates. Hi-C analyses revealed the sequences of ultra-long-range interactions on the X chromosome. The chromosome-scale scaffolds pave the way for further studies of the unusual chromosome movements in Bradysia 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.

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.

Amplification and expression of a salivary gland DNA puff gene in the prothoracic gland of Bradysia hygida (Diptera: Sciaridae)

The DNA puff BhC4-1 gene, located in DNA puff C4 of Bradysiahygida, is amplified and expressed in the salivary gland at the end of the fourth larval instar as a late response to the increase in 20-hydroxyecdysone titer that triggers metamorphosis. Functional studies revealed that the mechanisms that regulate BhC4-1 expression in the salivary gland are conserved in transgenic Drosophila. These studies also led to the identification of a cis-regulatory module that drives developmentally regulated expression of BhC4-1-lacZ in the prothoracic gland cells of the ring gland, a compound organ which in Drosophila results from the fusion of the prothoracic glands, the corpus allatum and the corpus cardiacum. Here we have investigated the occurrence of BhC4-1 expression in B. hygida prothoracic glands. We report the identification of the B. hygida prothoracic gland and demonstrate that it releases ecdysone. Using RT-qPCR, western blots and immunolocalization experiments, we demonstrate that the BhC4-1 mRNA and the BhC4-1 protein are both expressed in the B. hygida prothoracic glands at the same time that DNA puff C4 is formed in the salivary gland. We also show that BhC4-1 is concomitantly amplified 4.8-fold in the prothoracic gland and 23-fold in the salivary gland. Our results reveal the occurrence of stage specific expression of a DNA puff gene in the prothoracic glands of B. hygida, and extend previous studies that have shown that DNA puff genes expression is not restricted to the salivary gland. In addition, the description of stage specific gene amplification in the prothoracic glands of B. hygida constitutes the first demonstration that gene amplification in Diptera might occur concomitantly in two different tissues in the same developmental stage.

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.

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.

Replication in context: dynamic regulation of DNA replication patterns in metazoans

Replication in eukaryotes initiates from discrete genomic regions according to a strict, often tissue-specific temporal programme. However, the locations of initiation events within initiation regions vary, show sequence disparity and are affected by interactions with distal elements. Increasing evidence suggests that specification of replication sites and the timing of replication are dynamic processes that are regulated by tissue-specific and developmental cues, and are responsive to epigenetic modifications. Dynamic specification of replication patterns might serve to prevent or resolve possible spatial and/or temporal conflicts between replication, transcription and chromatin assembly, and facilitate subtle or extensive changes of gene expression during differentiation and development.

Ras signaling modulates activity of the ecdysone receptor EcR during cell migration in the Drosophila ovary

Ecdysone Receptor (EcR) mediates effects of the hormone ecdysone during larval molts, pupal metamorphosis, and adult female oogenesis. In the ovary, egg chamber formation requires interactions between the somatic follicle cell (FC) epithelium and the germ line nurse cell/oocyte cyst. Previous work has shown EcR is required in the germ line for egg chamber maturation, and here we examine EcR requirements in the FC at late stages of oogenesis. EcR protein is ubiquitous in the FC but its activity is restricted, visualized by activity of the "ligand sensor" hs-GAL4-EcR ligand binding domain fusion and EcRE-lacZ reporter gene expression. GAL4-EcR is activated in the FC by an ecdysone agonist and repressed by tissue-specific Ras GTPase signals. To determine the significance of restricted sites of EcR activity in the FC, we used targeted misexpression of the dominant negative EcR (EcR-DN) molecules EcR(F645A) and EcR(W650A). EcR-DN expression at stage 10 reduced EcRE-lacZ expression in the nurse cell FC and resulted in abnormal FC migrations, including aberrant centripetal migration and dorsal appendage tube formation, leading to the formation of cup-shaped eggs with shortened, branched dorsal appendages at stage 14. Clones of FC expressing EcR-DN displayed cell-autonomous increases in DE-cadherin expression and abnormal epithelial junction formation. EcR-DN expression caused thin eggshell phenotypes that correlated with both reduced levels of chorion gene expression and reduction in chorion gene amplification. Our results indicate that tissue-specific modulation of EcR activity by the Ras signaling pathway refines temporal ecdysone signals that regulate FC differentiation and cadherin-mediated epithelial cell shape changes.