Recent Evolution of a Maternally Acting Sex-Determining Supergene in a Fly with Single-Sex Broods.

Sex determination is a key developmental process, yet it is remarkably variable across the tree of life. The dipteran family Sciaridae exhibits one of the most unusual sex determination systems in which mothers control offspring sex through selective elimination of paternal X chromosomes. Whereas in some members of the family females produce mixed-sex broods, others such as the dark-winged fungus gnat Bradysia coprophila are monogenic, with females producing single-sex broods. Female-producing females were previously found to be heterozygous for a large X-linked paracentric inversion (X'), which is maternally inherited and absent from male-producing females. Here, we assembled and characterized the X' sequence. As close sequence homology between the X and X' made identification of the inversion challenging, we developed a k-mer-based approach to bin genomic reads before assembly. We confirmed that the inversion spans most of the X' chromosome (∼55 Mb) and encodes ∼3,500 genes. Analysis of the divergence between the inversion and the homologous region of the X revealed that it originated very recently (<0.5 Ma). Surprisingly, we found that the X' is more complex than previously thought and is likely to have undergone multiple rearrangements that have produced regions of varying ages, resembling a supergene composed of evolutionary strata. We found functional degradation of ∼7.3% of genes within the region of recombination suppression, but no evidence of accumulation of repetitive elements. Our findings provide an indication that sex-linked inversions are driving turnover of the strange sex determination system in this family of flies.

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.

Characterizing and controlling intrinsic biases of lambda exonuclease in nascent strand sequencing reveals phasing between nucleosomes and G-quadruplex motifs around a subset of human replication origins.

Nascent strand sequencing (NS-seq) is used to discover DNA replication origins genome-wide, allowing identification of features for their specification. NS-seq depends on the ability of lambda exonuclease (λ-exo) to efficiently digest parental DNA while leaving RNA-primer protected nascent strands intact. We used genomics and biochemical approaches to determine if λ-exo digests all parental DNA sequences equally. We report that λ-exo does not efficiently digest G-quadruplex (G4) structures in a plasmid. Moreover, λ-exo digestion of nonreplicating genomic DNA (LexoG0) enriches GC-rich DNA and G4 motifs genome-wide. We used LexoG0 data to control for nascent strand-independent λ-exo biases in NS-seq and validated this approach at the rDNA locus. The λ-exo-controlled NS-seq peaks are not GC-rich, and only 35.5% overlap with 6.8% of all G4s, suggesting that G4s are not general determinants for origin specification but may play a role for a subset. Interestingly, we observed a periodic spacing of G4 motifs and nucleosomes around the peak summits, suggesting that G4s may position nucleosomes at this subset of origins. Finally, we demonstrate that use of Na(+) instead of K(+) in the λ-exo digestion buffer reduced the effect of G4s on λ-exo digestion and discuss ways to increase both the sensitivity and specificity of NS-seq.

Bradysia (Sciara) coprophila larvae up-regulate DNA repair pathways and down-regulate developmental regulators in response to ionizing radiation.

The level of resistance to radiation and the developmental and molecular responses can vary between species, and even between developmental stages of one species. For flies (order: Diptera), prior studies concluded that the fungus gnat Bradysia (Sciara) coprophila (sub-order: Nematocera) is more resistant to irradiation-induced mutations that cause visible phenotypes than the fruit fly Drosophila melanogaster (sub-order: Brachycera). Therefore, we characterized the effects of and level of resistance to ionizing radiation on B. coprophila throughout its life cycle. Our data show that B. coprophila embryos are highly sensitive to even low doses of gamma-irradiation, whereas late-stage larvae can tolerate up to 80 Gy (compared to 40 Gy for D. melanogaster) and still retain their ability to develop to adulthood, though with a developmental delay. To survey the genes involved in the early transcriptional response to irradiation of B. coprophila larvae, we compared larval RNA-seq profiles with and without radiation treatment. The up-regulated genes were enriched for DNA damage response genes, including those involved in DNA repair, cell cycle arrest, and apoptosis, whereas the down-regulated genes were enriched for developmental regulators, consistent with the developmental delay of irradiated larvae. Interestingly, members of the PARP and AGO families were highly up-regulated in the B. coprophila radiation response. We compared the transcriptome responses in B. coprophila to the transcriptome responses in D. melanogaster from 3 previous studies: whereas pathway responses are highly conserved, specific gene responses are less so. Our study lays the groundwork for future work on the radiation responses in Diptera.

Aliens in the CYPome of the black fungus gnat, Bradysia coprophila.

The diverse cytochrome P450 enzymes of insects play essential physiological roles and also play important roles in the metabolism of environmental chemicals such as insecticides. We manually curated the complement of P450 (CYP) genes, or CYPome, of the black fungus gnat, Bradysia (Sciara) coprophila (Diptera, Sciaroidea), a species with a variable number of chromosomes. This CYPome carries two types of "alien" P450 genes. The first type of alien P450s was found among the 163 CYP genes of the core genome (autosomes and X). They consist of 28 sequences resulting from horizontal gene transfer, with closest sequences not found in insects, but in other arthropods, often Collembola. These genes are not contaminants, because they are expressed genes with introns, found in synteny with regular dipteran genes, also found in B. odoriphaga and B. hygida. Two such "alien" genes are representatives of CYP clans not otherwise found in insects, a CYP53 sequence related to fungal CYP53 genes, and a CYP19-like sequence similar to some collembolan sequences but of unclear origin. The second type of alien P450s are represented by 99 sequences from germline-restricted chromosomes (GRC). While most are P450 pseudogenes, 33 are apparently intact, with half being more closely related to P450s from Cecidomyiidae than from Sciaridae, thus supporting the hypothesis of a cross-family hybridization origin of the GRC.

Chromatin and gene expression changes during female Drosophila germline stem cell development illuminate the biology of highly potent stem cells.

Highly potent animal stem cells either self renew or launch complex differentiation programs, using mechanisms that are only partly understood. Drosophila female germline stem cells (GSCs) perpetuate without change over evolutionary time and generate cystoblast daughters that develop into nurse cells and oocytes. Cystoblasts initiate differentiation by generating a transient syncytial state, the germline cyst, and by increasing pericentromeric H3K9me3 modification, actions likely to suppress transposable element activity. Relatively open GSC chromatin is further restricted by Polycomb repression of testis or somatic cell-expressed genes briefly active in early female germ cells. Subsequently, Neijre/CBP and Myc help upregulate growth and reprogram GSC metabolism by altering mitochondrial transmembrane transport, gluconeogenesis, and other processes. In all these respects GSC differentiation resembles development of the totipotent zygote. We propose that the totipotent stem cell state was shaped by the need to resist transposon activity over evolutionary timescales.

Most animals are made up of two cell types: germline stem cells, which give rise to reproductive cells (egg and sperm) and pass their DNA to the next generation, and somatic cells, which make up the rest of the body. Transposable elements ­ fragments of DNA that can copy themselves and integrate into different parts of the genome ­ can greatly disrupt the integrity of the germ cell genome. Systems involving small RNAs and DNA methylation, which respectively modify the sequence and structure of the genome, can protect germ cells from the activity of transposable elements. While these systems have been studied extensively in late germ cells, less is known about how they work in germ cells generated early on in development. To investigate, Pang et al. studied the germline stem cells that give rise to eggs in female fruit flies. Techniques that measure DNA modifications showed that these germline stem cells and the cells they give rise to early on are better protected against transposable elements. This is likely due to the unusual cell cycle of early germ cells, which display a very short initial growth phase and special DNA replication timing during the synthesis phase. Until now, the purpose of these long-known cell cycle differences between early and late germ cells was not understood. Experiments also showed known transposable element defences are upregulated before the cell division that produces reproductive cells. DNA becomes more densely packed and germ cells connect with one another, forming germline 'cysts' that allow them to share small RNAs that can suppress transposable elements. Pang et al. propose that these changes compensate for the loss of enhanced repression that occurs in the earlier stem cell stage. Very similar changes also take place in the cells generated from fertilized eggs and in mammalian reproductive cells. Further experiments investigated how these changes impact the transition from stem cell to egg cell, revealing that germline stem cells express a wide diversity of genes, including most genes whose transcripts will be stored in the mature egg later on. Another type of cell produced by germline stem cells known as nurse cells, which synthesize most of the contents of the egg, dramatically upregulate genes supporting growth. Meanwhile, 25% of genes initially expressed in germline stem cells are switched off during the transition, partly due to a mechanism called Polycomb-mediated repression. The findings advance fundamental knowledge of how germline stem cells become egg cells, and could lead to important findings in developmental biology. Furthermore, understanding that for practical applications germline stem cells do not need to retain transposable element controls designed for evolutionary time scales means that removing them may make it easier to obtain and manipulate new stem cell lines and to develop new medical therapies.

The Drosophila CLAMP protein associates with diverse proteins on chromatin.

Gaining new insights into gene regulation involves an in-depth understanding of protein-protein interactions on chromatin. A powerful model for studying mechanisms of gene regulation is dosage compensation, a process that targets the X-chromosome to equalize gene expression between XY males and XX females. We previously identified a zinc finger protein in Drosophila melanogaster that plays a sex-specific role in targeting the Male-specific lethal (MSL) dosage compensation complex to the male X-chromosome, called the Chromatin-Linked Adapter for MSL Proteins (CLAMP). More recently, we established that CLAMP has non-sex-specific roles as an essential protein that regulates chromatin accessibility at promoters genome-wide. To identify associations between CLAMP and other factors in both male and female cells, we used two complementary mass spectrometry approaches. We demonstrate that CLAMP associates with the transcriptional regulator complex Negative Elongation Factor (NELF) in both sexes and determine that CLAMP reduces NELF recruitment to several target genes. In sum, we have identified many new CLAMP-associated factors and provide a resource for further study of this little understood essential protein.

The hunt for origins of DNA replication in multicellular eukaryotes.

Origins of DNA replication (ORIs) occur at defined regions in the genome. Although DNA sequence defines the position of ORIs in budding yeast, the factors for ORI specification remain elusive in metazoa. Several methods have been used recently to map ORIs in metazoan genomes with the hope that features for ORI specification might emerge. These methods are reviewed here with analysis of their advantages and shortcomings. The various factors that may influence ORI selection for initiation of DNA replication are discussed.

MinION Analysis and Reference Consortium: Phase 1 data release and analysis.

The advent of a miniaturized DNA sequencing device with a high-throughput contextual sequencing capability embodies the next generation of large scale sequencing tools. The MinION™ Access Programme (MAP) was initiated by Oxford Nanopore Technologies™ in April 2014, giving public access to their USB-attached miniature sequencing device. The MinION Analysis and Reference Consortium (MARC) was formed by a subset of MAP participants, with the aim of evaluating and providing standard protocols and reference data to the community. Envisaged as a multi-phased project, this study provides the global community with the Phase 1 data from MARC, where the reproducibility of the performance of the MinION was evaluated at multiple sites. Five laboratories on two continents generated data using a control strain of Escherichia coli K-12, preparing and sequencing samples according to a revised ONT protocol. Here, we provide the details of the protocol used, along with a preliminary analysis of the characteristics of typical runs including the consistency, rate, volume and quality of data produced. Further analysis of the Phase 1 data presented here, and additional experiments in Phase 2 of E. coli from MARC are already underway to identify ways to improve and enhance MinION performance.