Rapid evolution of piRNA clusters in the Drosophila melanogaster ovary.

The piRNA pathway is a highly conserved mechanism to repress transposable element (TE) activity in the animal germline via a specialized class of small RNAs called piwi-interacting RNAs (piRNAs). piRNAs are produced from discrete genomic regions called piRNA clusters (piCs). Although the molecular processes by which piCs function are relatively well understood in Drosophila melanogaster, much less is known about the origin and evolution of piCs in this or any other species. To investigate piC origin and evolution, we use a population genomic approach to compare piC activity and sequence composition across eight geographically distant strains of D. melanogaster with high-quality long-read genome assemblies. We perform annotations of ovary piCs and genome-wide TE content in each strain. Our analysis uncovers extensive variation in piC activity across strains and signatures of rapid birth and death of piCs. Most TEs inferred to be recently active show an enrichment of insertions into old and large piCs, consistent with the previously proposed "trap" model of piC evolution. In contrast, a small subset of active LTR families is enriched for the formation of new piCs, suggesting that these TEs have higher proclivity to form piCs. Thus, our findings uncover processes leading to the origin of piCs. We propose that piC evolution begins with the emergence of piRNAs from individual insertions of a few select TE families prone to seed new piCs that subsequently expand by accretion of insertions from most other TE families during evolution to form larger "trap" clusters. Our study shows that TEs themselves are the major force driving the rapid evolution of piCs.

Theoretical microwave spectra of interstellar nitrogen-containing PAHs.

The recent discovery of naphthalene (C10H8) in cyano-substituted polycyclic aromatic hydrocarbon (CN-PAH) form in the Taurus molecular cloud (TMC-1) has sparked curiosity regarding the search for other nitrogen-containing naphthalenes in similar interstellar environments. In this light, naphthalenes having N atoms in the structure are promising candidates to be searched for in cold, dark molecular clouds such as TMC-1. Since obtaining data on such samples in the laboratory is complicated, the present work reports theoretical microwave spectra of naphthalene in all N-substituted forms. Density functional theory (DFT) calculations are employed to calculate the spectroscopic constants and simulate the rotational spectra with hyperfine splitting. For cold temperature regions such as TMC-1 (about 5 K), the considered N-naphthalene species show the strongest transition around centimetre wavelengths, a typical range for PAH-related species in dark molecular clouds. Accurate rotational data provided here may act as a guide for laboratory experiments and astronomical searches.

Rapid evolution of piRNA clusters in the Drosophila melanogaster ovary.

Animal genomes are parasitized by a horde of transposable elements (TEs) whose mutagenic activity can have catastrophic consequences. The piRNA pathway is a conserved mechanism to repress TE activity in the germline via a specialized class of small RNAs associated with effector Piwi proteins called piwi-associated RNAs (piRNAs). piRNAs are produced from discrete genomic regions called piRNA clusters (piCs). While piCs are generally enriched for TE sequences and the molecular processes by which they are transcribed and regulated are relatively well understood in Drosophila melanogaster, much less is known about the origin and evolution of piCs in this or any other species. To investigate piC evolution, we use a population genomics approach to compare piC activity and sequence composition across 8 geographically distant strains of D. melanogaster with high quality long-read genome assemblies. We perform extensive annotations of ovary piCs and TE content in each strain and test predictions of two proposed models of piC evolution. The 'de novo' model posits that individual TE insertions can spontaneously attain the status of a small piC to generate piRNAs silencing the entire TE family. The 'trap' model envisions large and evolutionary stable genomic clusters where TEs tend to accumulate and serves as a long-term "memory" of ancient TE invasions and produce a great variety of piRNAs protecting against related TEs entering the genome. It remains unclear which model best describes the evolution of piCs. Our analysis uncovers extensive variation in piC activity across strains and signatures of rapid birth and death of piCs in natural populations. Most TE families inferred to be recently or currently active show an enrichment of strain-specific insertions into large piCs, consistent with the trap model. By contrast, only a small subset of active LTR retrotransposon families is enriched for the formation of strain-specific piCs, suggesting that these families have an inherent proclivity to form de novo piCs. Thus, our findings support aspects of both 'de novo' and 'trap' models of piC evolution. We propose that these two models represent two extreme stages along an evolutionary continuum, which begins with the emergence of piCs de novo from a few specific LTR retrotransposon insertions that subsequently expand by accretion of other TE insertions during evolution to form larger 'trap' clusters. Our study shows that piCs are evolutionarily labile and that TEs themselves are the major force driving the formation and evolution of piCs.

A comprehensive rotational study of astronomical iso-pentane within 84 to 111 GHz.

The rotational line survey by ALMA (Atacama Large Millimeter/submillimeter Array) recently revealed the presence of i-C3H7CN (i-PrCN) and n-C3H7CN (n-PrCN) in 3-mm atmospheric window between 84 to 111 GHz towards the hot core region Sagittarius B2(N) (Sgr B2(N)). This was the first interstellar detection of a linear straight chain molecule. In this light, we report the rotational spectra of C5H12 isomeric group in the same frequency range. We performed quantum chemical calculations for spectroscopic parameters. The pure rotational spectrum of the species has been simulated using the PGOPHER program. The rotational spectrum of this molecule makes it a good candidate for future astronomical detections since the radio lines can be calculated to very high accuracy in mm/sub-mm wave region.

Genetic variation in P-element dysgenic sterility is associated with double-strand break repair and alternative splicing of TE transcripts.

The germline mobilization of transposable elements (TEs) by small RNA mediated silencing pathways is conserved across eukaryotes and critical for ensuring the integrity of gamete genomes. However, genomes are recurrently invaded by novel TEs through horizontal transfer. These invading TEs are not targeted by host small RNAs, and their unregulated activity can cause DNA damage in germline cells and ultimately lead to sterility. Here we use hybrid dysgenesis-a sterility syndrome of Drosophila caused by transposition of invading P-element DNA transposons-to uncover host genetic variants that modulate dysgenic sterility. Using a panel of highly recombinant inbred lines of Drosophila melanogaster, we identified two linked quantitative trait loci (QTL) that determine the severity of dysgenic sterility in young and old females, respectively. We show that ovaries of fertile genotypes exhibit increased expression of splicing factors that suppress the production of transposase encoding transcripts, which likely reduces the transposition rate and associated DNA damage. We also show that fertile alleles are associated with decreased sensitivity to double-stranded breaks and enhanced DNA repair, explaining their ability to withstand high germline transposition rates. Together, our work reveals a diversity of mechanisms whereby host genotype modulates the cost of an invading TE, and points to genetic variants that were likely beneficial during the P-element invasion.

Identification and Characterization of Small RNA Markers of Age in the Blow Fly Cochliomyia macellaria (Fabricius) (Diptera: Calliphoridae).

Blow fly development is important in decomposition ecology, agriculture, and forensics. Much of the impact of these species is from immature samples, thus knowledge of their development is important to enhance or ameliorate their effects. One application of this information is the estimation of immature insect age to provide temporal information for death investigations. While traditional markers of age such as stage and size are generally accurate, they lack precision in later developmental stages. We used miRNA sequencing to measure miRNA expression, throughout development, of the secondary screwworm, Cochliomyia macellaria (Fabricius) (Diptera: Calliphoridae) and identified 217 miRNAs present across the samples. Ten were identified to be significantly differentially expressed in larval samples and seventeen were found to be significantly differentially expressed in intrapuparial samples. Twenty-eight miRNAs were identified to be differentially expressed between sexes. Expression patterns of two miRNAs, miR-92b and bantam, were qPCR-validated in intrapuparial samples; these and likely food-derived miRNAs appear to be stable markers of age in C. macellaria. Our results support the use of miRNAs for developmental markers of age and suggest further investigations across species and under a range of abiotic and biotic conditions.

Transposable element landscapes in aging Drosophila.

Genetic mechanisms that repress transposable elements (TEs) in young animals decline during aging, as reflected by increased TE expression in aged animals. Does increased TE expression during aging lead to more genomic TE copies in older animals? To address this question, we quantified TE Landscapes (TLs) via whole genome sequencing of young and aged Drosophila strains of wild-type and mutant backgrounds. We quantified TLs in whole flies and dissected brains and validated the feasibility of our approach in detecting new TE insertions in aging Drosophila genomes when small RNA and RNA interference (RNAi) pathways are compromised. We also describe improved sequencing methods to quantify extra-chromosomal DNA circles (eccDNAs) in Drosophila as an additional source of TE copies that accumulate during aging. Lastly, to combat the natural progression of aging-associated TE expression, we show that knocking down PAF1, a conserved transcription elongation factor that antagonizes RNAi pathways, may bolster suppression of TEs during aging and extend lifespan. Our study suggests that in addition to a possible influence by different genetic backgrounds, small RNA and RNAi mechanisms may mitigate genomic TL expansion despite the increase in TE transcripts during aging.

A mosquito small RNA genomics resource reveals dynamic evolution and host responses to viruses and transposons.

Although mosquitoes are major transmission vectors for pathogenic arboviruses, viral infection has little impact on mosquito health. This immunity is caused in part by mosquito RNA interference (RNAi) pathways that generate antiviral small interfering RNAs (siRNAs) and Piwi-interacting RNAs (piRNAs). RNAi also maintains genome integrity by potently repressing mosquito transposon activity in the germline and soma. However, viral and transposon small RNA regulatory pathways have not been systematically examined together in mosquitoes. Therefore, we developed an integrated mosquito small RNA genomics (MSRG) resource that analyzes the transposon and virus small RNA profiles in mosquito cell cultures and somatic and gonadal tissues across four medically important mosquito species. Our resource captures both somatic and gonadal small RNA expression profiles within mosquito cell cultures, and we report the evolutionary dynamics of a novel Mosquito-Conserved piRNA Cluster Locus (MCpiRCL) made up of satellite DNA repeats. In the larger culicine mosquito genomes we detected highly regular periodicity in piRNA biogenesis patterns coinciding with the expansion of Piwi pathway genes. Finally, our resource enables detection of cross talk between piRNA and siRNA populations in mosquito cells during a response to virus infection. The MSRG resource will aid efforts to dissect and combat the capacity of mosquitoes to tolerate and spread arboviruses.

Diverse Defenses: A Perspective Comparing Dipteran Piwi-piRNA Pathways.

Animals face the dual threat of virus infections hijacking cellular function and transposons proliferating in germline genomes. For insects, the deeply conserved RNA interference (RNAi) pathways and other chromatin regulators provide an important line of defense against both viruses and transposons. For example, this innate immune system displays adaptiveness to new invasions by generating cognate small RNAs for targeting gene silencing measures against the viral and genomic intruders. However, within the Dipteran clade of insects, Drosophilid fruit flies and Culicids mosquitoes have evolved several unique mechanistic aspects of their RNAi defenses to combat invading transposons and viruses, with the Piwi-piRNA arm of the RNAi pathways showing the greatest degree of novel evolution. Whereas central features of Piwi-piRNA pathways are conserved between Drosophilids and Culicids, multiple lineage-specific innovations have arisen that may reflect distinct genome composition differences and specific ecological and physiological features dividing these two branches of Dipterans. This perspective review focuses on the most recent findings illuminating the Piwi/piRNA pathway distinctions between fruit flies and mosquitoes, and raises open questions that need to be addressed in order to ameliorate human diseases caused by pathogenic viruses that mosquitoes transmit as vectors.

Effect of Phenotype Selection on Genome Size Variation in Two Species of Diptera.

Genome size varies widely across organisms yet has not been found to be related to organismal complexity in eukaryotes. While there is no evidence for a relationship with complexity, there is evidence to suggest that other phenotypic characteristics, such as nucleus size and cell-cycle time, are associated with genome size, body size, and development rate. However, what is unknown is how the selection for divergent phenotypic traits may indirectly affect genome size. Drosophila melanogaster were selected for small and large body size for up to 220 generations, while Cochliomyia macellaria were selected for 32 generations for fast and slow development. Size in D. melanogaster significantly changed in terms of both cell-count and genome size in isolines, but only the cell-count changed in lines which were maintained at larger effective population sizes. Larger genome sizes only occurred in a subset of D. melanogaster isolines originated from flies selected for their large body size. Selection for development time did not change average genome size yet decreased the within-population variation in genome size with increasing generations of selection. This decrease in variation and convergence on a similar mean genome size was not in correspondence with phenotypic variation and suggests stabilizing selection on genome size in laboratory conditions.

Facultative Viviparity in a Flesh Fly (Diptera: Sarcophagidae): Forensic Implications of High Variability in Rates of Oviparity in Blaesoxipha plinthopyga (Diptera: Sarcophagidae).

Flesh flies are major primary consumers of carrion and are commonly found on human remains. Due to this latter feeding habit, their development rates can be used to provide temporal information in forensic investigations. This is usually done by referencing published flesh fly development datasets. Flesh flies are typically assumed to be strictly viviparous and datasets reporting their development rates therefore start at the first larval instar. However, an increasing number of studies has identified oviposition by flesh flies, including the forensically relevant species Blaesoxipha plinthopyga Wiedemann. To assess the impact of egg-laying behavior on casework, oviparity rates and time before larval hatching were assessed under controlled laboratory conditions that reflect common casework conditions in Harris County, Texas. We demonstrated systematic deposition of viable eggs but at a very variable rate between samples. Similarly, the duration between oviposition and larval hatching was highly variable, with some eggs taking more than a day to hatch after deposition. These results highlight the need to account for embryonic development in forensic investigations including B. plinthopyga and advocates for the re-evaluation of the assumed strict viviparity of the Sarcophagidae.

Har-P, a short P-element variant, weaponizes P-transposase to severely impair Drosophila development.

Without transposon-silencing Piwi-interacting RNAs (piRNAs), transposition causes an ovarian atrophy syndrome in Drosophila called gonadal dysgenesis (GD). Harwich (Har) strains with P-elements cause severe GD in F1 daughters when Har fathers mate with mothers lacking P-element-piRNAs (i.e. ISO1 strain). To address the mystery of why Har induces severe GD, we bred hybrid Drosophila with Har genomic fragments into the ISO1 background to create HISR-D or HISR-N lines that still cause Dysgenesis or are Non-dysgenic, respectively. In these lines, we discovered a highly truncated P-element variant we named 'Har-P' as the most frequent de novo insertion. Although HISR-D lines still contain full-length P-elements, HISR-N lines lost functional P-transposase but retained Har-P's that when crossed back to P-transposase restores GD induction. Finally, we uncovered P-element-piRNA-directed repression on Har-P's transmitted paternally to suppress somatic transposition. The Drosophila short Har-P's and full-length P-elements relationship parallels the MITEs/DNA-transposase in plants and SINEs/LINEs in mammals.

Paternal Induction of Hybrid Dysgenesis in Drosophila melanogaster Is Weakly Correlated with Both P-Element and hobo Element Dosage.

Transposable elements (TEs) are virtually ubiquitous components of genomes, yet they often impose significant fitness consequences on their hosts. In addition to producing specific deleterious mutations by insertional inactivation, TEs also impose general fitness costs by inducing DNA damage and participating in ectopic recombination. These latter fitness costs are often assumed to be dosage-dependent, with stronger effects occurring in the presence of higher TE copy numbers. We test this assumption in Drosophila melanogaster by considering the relationship between the copy number of two active DNA transposons, P-element and hobo element, and the incidence of hybrid dysgenesis, a sterility syndrome associated with transposon activity in the germline. By harnessing a subset of the Drosophila Genetic Reference Panel (DGRP), a group of fully-sequenced D. melanogaster strains, we describe quantitative and structural variation in P-elements and hobo elements among wild-derived genomes and associate these factors with hybrid dysgenesis. We find that the incidence of hybrid dysgenesis is associated with both P-element and hobo element copy number in a dosage-dependent manner. However, the relationship is weak for both TEs, suggesting that dosage alone explains only a small part of TE-associated fitness costs.