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.

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.

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.

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 (approximately 55Mb) and encodes around 3500 genes. Analysis of the divergence between the inversion and the homologous region of the X revealed that it originated very recently (<0.5 mya). 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 around 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.

Observation of A Raman mode splitting in few layer black phosphorus encapsulated with hexagonal boron nitride

We investigate the impact of encapsulation with hexagonal boron nitride (h-BN) on the Raman spectrum of few layer black phosphorus. The encapsulation results in a significant reduction of the line width of the Raman modes of black phosphorus, due to a reduced phonon scattering rate. We observe a so far elusive peak in the Raman spectra ∼4 cm^-1 above the A mode in trilayer and thicker flakes, which had not been observed experimentally. The newly observed mode originates from the strong black phosphorus inter-layer interaction, which induces a hardening of the surface atom vibration with respect to the corresponding modes of the inner layers. The observation of this mode suggests a significant impact of h-BN encapsulation on the properties of black phosphorus and can serve as an indicator of the quality of its surface.

Probing the Interlayer Exciton Physics in a MoS2/MoSe2/MoS2 van der Waals Heterostructure

Stacking atomic monolayers of semiconducting transition metal dichalcogenides (TMDs) has emerged as an effective way to engineer their properties. In principle, the staggered band alignment of TMD heterostructures should result in the formation of interlayer excitons with long lifetimes and robust valley polarization. However, these features have been observed simultaneously only in MoSe₂/WSe₂ heterostructures. Here we report on the observation of long-lived interlayer exciton emission in a MoS₂/MoSe₂/MoS₂ trilayer van der Waals heterostructure. The interlayer nature of the observed transition is confirmed by photoluminescence spectroscopy, as well as by analyzing the temporal, excitation power, and temperature dependence of the interlayer emission peak. The observed complex photoluminescence dynamics suggests the presence of quasi-degenerate momentum-direct and momentum-indirect bandgaps. We show that circularly polarized optical pumping results in long-lived valley polarization of interlayer exciton. Intriguingly, the interlayer exciton photoluminescence has helicity opposite to the excitation. Our results show that through a careful choice of the TMDs forming the van der Waals heterostructure it is possible to control the circular polarization of the interlayer exciton emission.

In situ Raman spectroscopy of the graphene/water interface of a solution-gated field-effect transistor: electron-phonon coupling and spectroelectrochemistry

We present a novel measurement approach which combines the electrical characterization of solution-gated field-effect transistors based on epitaxial bilayer graphene on 4H-SiC (0001) with simultaneous Raman spectroscopy. By changing the gate voltage, we observed Raman signatures related to the resonant electron-phonon coupling. An analysis of these Raman bands enabled the extraction of the geometrical capacitance of the system and an accurate calculation of the Fermi levels for bilayer graphene. An intentional application of higher gate voltages allowed us to trigger electrochemical reactions, which we followed in situ by Raman spectroscopy. The reactions showed a partially reversible character, as indicated by an emergence/disappearance of peaks assigned to C-H and Si-H vibration modes as well as an increase/decrease of the defect-related Raman D band intensity. Our setup provides a highly interesting platform for future spectroelectrochemical research on electrically-induced sorption processes of graphene on the micrometer scale.

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.

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.

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.

Novel p21-activated kinase-dependent protrusions characteristically formed at the edge of transformed cells

During long-term culture, certain lines become neoplastic while accumulating changes in cell shape. Early and late cell populations have characteristic shape phenotypes that have been quantified by computerized assay. Phenotypes are determined from variables describing three-dimensional aspects of the subcellular distribution of mass. The features of cells can be recognized by use of latent factors, which are theoretical variables based on the covariance of the primary variables. Factor #7 represented a cell edge feature different from filopodia. We studied the morphological characteristics and morphogenesis of the feature. Brief exposure of cells from rat tracheal epithelium to phorbol 12-myristate 13-acetate (PMA) enhanced #7 values. The time to reach maximal #7 values was prolonged if PMA was administered with calcium ionophore or lysophosphatidic acid (LPA). Factor #7 was elevated during periods of ruffling suppression and stress fiber reorganization. Cells showing high #7 values were examined by scanning electron microscopy (SEM) and found to exhibit strap-shaped and cupola-shaped projections. Because RhoA regulates stress fiber formation, we sought to perturb #7 features by introducing dominant-acting negative and positive constructs of RhoA, RhoA-N19, and RhoA-V14. Neither affected #7 values. Although overexpression of the kinase inhibitory domain of p21-activated kinase 1 (PAK) had no effect on #7 values, they were affected by overexpression of a domain binding PAK-interacting guanine nucleotide exchange factor (PIX). Because a PAK-PIX complex is implicated in the remodeling of focal complexes (FCs) and recycling of PAK to the cytoplasm, the results implicate a component of FCs in the formation of #7 features. The data suggested that feature formation is driven by activated Cdc42-binding kinase (ACK) and Rac. Moreover, they suggested that the #7 protrusions are neurite-like structures and that their development involves FC regulation.