Global genomics of Aedes aegypti unveils widespread and novel infectious viruses capable of triggering a small RNA response.

The mosquito Aedes aegypti is a prominent vector for arboviruses, but the breadth of mosquito viruses that infects this specie is not fully understood. In the broadest global survey to date of over 200 Ae. aegypti small RNA samples, we detected viral small interfering RNAs (siRNAs) and Piwi interacting RNAs (piRNAs) arising from mosquito viruses. We confirmed that most academic laboratory colonies of Ae. aegypti lack persisting viruses, yet two commercial strains were infected by a novel tombus-like virus. Ae. aegypti from North to South American locations were also teeming with multiple insect viruses, with Anphevirus and a bunyavirus displaying geographical boundaries from the viral small RNA patterns. Asian Ae. aegypti small RNA patterns indicate infections by similar mosquito viruses from the Americas and reveal the first wild example of dengue virus infection generating viral small RNAs. African Ae. aegypti also contained various viral small RNAs including novel viruses only found in these African substrains. Intriguingly, viral long RNA patterns can differ from small RNA patterns, indicative of viral transcripts evading the mosquitoes' RNA interference (RNAi) machinery. To determine whether the viruses we discovered via small RNA sequencing were replicating and transmissible, we infected C6/36 and Aag2 cells with Ae. aegypti homogenates. Through blind passaging, we generated cell lines stably infected by these mosquito viruses which then generated abundant viral siRNAs and piRNAs that resemble the native mosquito viral small RNA patterns. This mosquito small RNA genomics approach augments surveillance approaches for emerging infectious diseases.

Myc Dysregulation in Activated Macrophages Initiates Iron-Mediated Lipid Peroxidation that Fuels Type I Interferon and Compromises TB Resistance.

A quarter of human population is infected with Mycobacterium tuberculosis, but less than 10% of those infected develop clinical, mostly pulmonary, TB. To dissect mechanisms of susceptibility in immunocompetent individuals, we developed a genetically defined sst1-susceptible mouse model that uniquely reproduces a defining feature of human TB: development of necrotic lung lesions after infection with virulent Mtb. In this study, we explored the connectivity of the sst1-regulated pathways during prolonged macrophage activation with TNF. We determined that the aberrant response of the sst1-susceptible macrophages to TNF was primarily driven by conflicting Myc and antioxidant response pathways that resulted in a coordinated failure to properly sequester intracellular iron and activate ferroptosis inhibitor enzymes. Consequently, iron-mediated lipid peroxidation fueled IFNß superinduction and sustained the Type I Interferon (IFN-I) pathway hyperactivity that locked the sst1-susceptible macrophages in a state of unresolving stress and compromised their resistance to Mtb. The accumulation of the aberrantly activated, stressed, macrophages within granuloma microenvironment led to the local failure of anti-tuberculosis immunity and tissue necrosis. Our findings suggest a novel link between metabolic dysregulation in macrophages and susceptibility to TB, offering insights into potential therapeutic targets aimed at modulating macrophage function and improving TB control.

Dengue Virus-2 Infection Affects Fecundity and Elicits Specific Transcriptional Changes in the Ovaries of Aedes aegypti Mosquitoes.

Dengue fever (DF), caused by the dengue virus (DENV), is the most burdensome arboviral disease in the world, with an estimated 400 million infections each year. The Aedes aegypti mosquito is the main vector of DENV and transmits several other human pathogens, including Zika, yellow fever, and chikungunya viruses. Previous studies have shown that the pathogen infection of mosquitoes can alter reproductive fitness, revealing specific vector-pathogen interactions that are key determinants of vector competence. However, only a handful of studies have examined the effect of DENV infection in A. aegypti, showing a reduction in lifespan and fecundity over multiple blood meals. To provide a more comprehensive analysis of the impact of DENV infection on egg laying and fecundity, we assessed egg laying timing in DENV-2 blood-fed mosquitoes (infected group) compared to mock blood-fed mosquitoes (control group). We confirmed a significant decrease in fecundity during the first gonadotrophic cycle. To further investigate this phenotype and the underlying DENV-2 infection-dependent changes in gene expression, we conducted a transcriptomic analysis for differentially expressed genes in the ovaries of A. aegypti infected with DENV-2 vs. mock-infected mosquitoes. This analysis reveals several DENV-2-regulated genes; among them, we identified a group of 12 metabolic genes that we validated using reverse transcription-quantitative PCR (RT-qPCR). Interestingly, two genes found to be upregulated in DENV-infected mosquito ovaries exhibited an antiviral role for DENV-2 in an Aedes cell line. Altogether, this study offers useful insights into the virus-vector interface, highlighting the importance of gene expression changes in the mosquito's ovary during DENV-2 infection in the first gonadotrophic  cycle,  triggering  antiviral  responses  that  may  possibly  interfere  with mosquito reproduction. This information is extremely relevant for further investigation of A. aegypti's ability to tolerate viruses since virally infected mosquitoes in nature constitute a powerful source of supporting viruses during intra-epidemic periods, causing a huge burden on the public health system.

Extending and Running the Mosquito Small RNA Genomics Resource Pipeline.

The Mosquito Small RNA Genomics (MSRG) resource is a repository of analyses on the small RNA transcriptomes of mosquito cell cultures and somatic and gonadal tissues. This resource allows for comparing the regulation dynamics of small RNAs generated from transposons and viruses across mosquito species. This chapter covers the procedures to set up the MSRG resource pipeline as a new installation by detailing the necessary collection of genome reference and annotation files and lists of microRNAs (miRNAs) hairpin sequences, transposon repeats consensus sequences, and virus genome sequences. Proper execution of the MSRG resource pipeline yields outputs amenable to biologists to further analyze with desktop and spreadsheet software to gain insights into the balance between arthropod endogenous small RNA populations and the proportions of virus-derived small RNAs that include Piwi-interacting RNAs (piRNAs) and endogenous small interfering RNAs (siRNAs).

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.

DNA templates with blocked long 3' end single-stranded overhangs (BL3SSO) promote bona fide Cas9-stimulated homology-directed repair of long transgenes into endogenous gene loci.

Knock-in of large transgenes by Cas9-mediated homology-directed repair (HDR) is an extremely inefficient process. Although the use of single-stranded oligonucleotides (ssODN) as an HDR donor has improved the integration of smaller transgenes, they do not support efficient insertion of large DNA sequences. In an effort to gain insights into the mechanism(s) governing the HDR-mediated integration of larger transgenes and to improve the technology, we conducted knock-in experiments targeting the human EMX1 locus and applied rigorous genomic PCR analyses in the human HEK293 cell line. This exercise revealed an unexpected molecular complication arising from the transgene HDR being initiated at the single homology arm and the subsequent genomic integration of plasmid backbone sequences. To pivot around this problem, we devised a novel PCR-constructed template containing blocked long 3' single-stranded overhangs (BL3SSO) that greatly improved the efficiency of bona fide Cas9-stimulated HDR at the EMX1 locus. We further refined BL3SSO technology and successfully used it to insert GFP transgenes into two important interferon-stimulated genes (ISGs) loci, Viperin/RSAD2, and ISG15. This study demonstrates the utility of the BL3SSO platform for inserting long DNA sequences into both constitutive and inducible endogenous loci to generate novel human cell lines for the study of important biological processes.

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.

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.

Klotho Is Neuroprotective in the Superoxide Dismutase (SOD1G93A) Mouse Model of ALS.

Amyotrophic lateral sclerosis (ALS) is a progressive neurodegenerative disorder characterized by the loss of motor neurons in the brain and spinal cord. ALS neuropathology is associated with increased oxidative stress, excitotoxicity, and inflammation. We and others reported that the anti-aging and cognition-enhancing protein Klotho is a neuroprotective, antioxidative, anti-inflammatory, and promyelinating protein. In mice, its absence leads to an extremely shortened life span and to multiple phenotypes resembling human aging, including motor and hippocampal neurodegeneration and cognitive impairment. In contrast, its overexpression extends life span, enhances cognition, and confers resistance against oxidative stress; it also reduces premature mortality and cognitive and behavioral abnormalities in an animal model for Alzheimer's disease (AD). These pleiotropic beneficial properties of Klotho suggest that Klotho could be a potent therapeutic target for preventing neurodegeneration in ALS. Klotho overexpression in the SOD1 mouse model of ALS resulted in delayed onset and progression of the disease and extended survival that was more prominent in females than in males. Klotho reduced the expression of neuroinflammatory markers and prevented neuronal loss with the more profound effect in the spinal cord than in the motor cortex. The effect of Klotho was accompanied by reduced expression of proinflammatory cytokines and enhanced the expression of antioxidative and promyelinating factors in the motor cortex and spinal cord of Klotho × SOD1 compared to SOD1 mice. Our study provides evidence that increased levels of Klotho alleviate ALS-associated pathology in the SOD1 mouse model and may serve as a basis for developing Klotho-based therapeutic strategies for ALS.

Interplay between small RNA pathways shapes chromatin landscapes in C. elegans.

The nematode Caenorhabditis elegans contains several types of endogenous small interfering RNAs (endo-siRNAs) produced by RNA-dependent RNA polymerase (RdRP) complexes. Both 'silencing' siRNAs bound by Worm-specific Argonautes (WAGO) and 'activating' siRNAs bound by the CSR-1 Argonaute require the DRH-3 helicase, an RdRP component. Here, we show that, in the drh-3(ne4253) mutant deficient in RdRP-produced secondary endo-siRNAs, the silencing histone mark H3K9me3 is largely depleted, whereas in the csr-1 partially rescued null mutant strain (WM193), this mark is ectopically deposited on CSR-1 target genes. Moreover, we observe ectopic H3K9me3 at enhancer elements and an increased number of small RNAs that match enhancers in both drh-3 and csr-1 mutants. Finally, we detect accumulation of H3K27me3 at highly expressed genes in the drh-3(ne4253) mutant, which correlates with their reduced transcription. Our study shows that when abundant RdRP-produced siRNAs are depleted, there is ectopic elevation of noncoding RNAs linked to sites with increased silencing chromatin marks. Moreover, our results suggest that enhancer small RNAs may guide local H3K9 methylation.

Drosophila PAF1 Modulates PIWI/piRNA Silencing Capacity.

To test the directness of factors in initiating PIWI-directed gene silencing, we employed a Piwi-interacting RNA (piRNA)-targeted reporter assay in Drosophila ovary somatic sheet (OSS) cells [1]. This assay confirmed direct silencing roles for piRNA biogenesis factors and PIWI-associated factors [2-12] but suggested that chromatin-modifying proteins may act downstream of the initial silencing event. Our data also revealed that RNA-polymerase-II-associated proteins like PAF1 and RTF1 antagonize PIWI-directed silencing. PAF1 knockdown enhances PIWI silencing of reporters when piRNAs target the transcript region proximal to the promoter. Loss of PAF1 suppresses endogenous transposable element (TE) transcript maturation, whereas a subset of gene transcripts and long-non-coding RNAs adjacent to TE insertions are affected by PAF1 knockdown in a similar fashion to piRNA-targeted reporters. Additionally, transcription activation at specific TEs and TE-adjacent loci during PIWI knockdown is suppressed when PIWI and PAF1 levels are both reduced. Our study suggests a mechanistic conservation between fission yeast PAF1 repressing AGO1/small interfering RNA (siRNA)-directed silencing [13, 14] and Drosophila PAF1 opposing PIWI/piRNA-directed silencing.

Xenopus Piwi proteins interact with a broad proportion of the oocyte transcriptome.

Piwi proteins utilize small RNAs (piRNAs) to recognize target transcripts such as transposable elements (TE). However, extensive piRNA sequence diversity also suggests that Piwi/piRNA complexes interact with many transcripts beyond TEs. To determine Piwi target RNAs, we used ribonucleoprotein-immunoprecipitation (RIP) and cross-linking and immunoprecipitation (CLIP) to identify thousands of transcripts associated with the Piwi proteins XIWI and XILI (Piwi-protein-associated transcripts, PATs) from early stage oocytes of X. laevis and X. tropicalis Most PATs associate with both XIWI and XILI and include transcripts of developmentally important proteins in oogenesis and embryogenesis. Only a minor fraction of PATs in both frog species displayed near perfect matches to piRNAs. Since predicting imperfect pairing between all piRNAs and target RNAs remains intractable, we instead determined that PAT read counts correlate well with the lengths and expression levels of transcripts, features that have also been observed for oocyte mRNAs associated with Drosophila Piwi proteins. We used an in vitro assay with exogenous RNA to confirm that XIWI associates with RNAs in a length- and concentration-dependent manner. In this assay, noncoding transcripts with many perfectly matched antisense piRNAs were unstable, whereas coding transcripts with matching piRNAs were stable, consistent with emerging evidence that Piwi proteins both promote the turnover of TEs and other RNAs, and may also regulate mRNA localization and translation. Our study suggests that Piwi proteins play multiple roles in germ cells and establishes a tractable vertebrate system to study the role of Piwi proteins in transcript regulation.

Small RNAs from a Big Genome: The piRNA Pathway and Transposable Elements in the Salamander Species Desmognathus fuscus.

Most of the largest vertebrate genomes are found in salamanders, a clade of amphibians that includes 686 species. Salamander genomes range in size from 14 to 120 Gb, reflecting the accumulation of large numbers of transposable element (TE) sequences from all three TE classes. Although DNA loss rates are slow in salamanders relative to other vertebrates, high levels of TE insertion are also likely required to explain such high TE loads. Across the Tree of Life, novel TE insertions are suppressed by several pathways involving small RNA molecules. In most known animals, TE activity in the germline is primarily regulated by the Piwi-interacting RNA (piRNA) pathway. In this study, we test the hypothesis that salamanders' unusually high TE loads reflect the loss of the ancestral piRNA-mediated TE-silencing machinery. We characterized the small RNA pool in the female and male adult gonads, testing for the presence of small RNA molecules that bear the characteristics of TE-targeting piRNAs. We also analyzed the amino acid sequences of piRNA pathway proteins from salamanders and other vertebrates, testing whether the overall patterns of sequence divergence are consistent with conserved pathway function across the vertebrate clade. Our results do not support the hypothesis of piRNA pathway loss; instead, they suggest that the piRNA pathway is expressed in salamanders. Given these results, we propose hypotheses to explain how the extraordinary TE loads in salamander genomes could have accumulated, despite the expression of TE-silencing machinery.

Unique transposon landscapes are pervasive across Drosophila melanogaster genomes.

To understand how transposon landscapes (TLs) vary across animal genomes, we describe a new method called the Transposon Insertion and Depletion AnaLyzer (TIDAL) and a database of >300 TLs in Drosophila melanogaster (TIDAL-Fly). Our analysis reveals pervasive TL diversity across cell lines and fly strains, even for identically named sub-strains from different laboratories such as the ISO1 strain used for the reference genome sequence. On average, >500 novel insertions exist in every lab strain, inbred strains of the Drosophila Genetic Reference Panel (DGRP), and fly isolates in the Drosophila Genome Nexus (DGN). A minority (<25%) of transposon families comprise the majority (>70%) of TL diversity across fly strains. A sharp contrast between insertion and depletion patterns indicates that many transposons are unique to the ISO1 reference genome sequence. Although TL diversity from fly strains reaches asymptotic limits with increasing sequencing depth, rampant TL diversity causes unsaturated detection of TLs in pools of flies. Finally, we show novel transposon insertions negatively correlate with Piwi-interacting RNA (piRNA) levels for most transposon families, except for the highly-abundant roo retrotransposon. Our study provides a useful resource for Drosophila geneticists to understand how transposons create extensive genomic diversity in fly cell lines and strains.

Conserved piRNA Expression from a Distinct Set of piRNA Cluster Loci in Eutherian Mammals.

The Piwi pathway is deeply conserved amongst animals because one of its essential functions is to repress transposons. However, many Piwi-interacting RNAs (piRNAs) do not base-pair to transposons and remain mysterious in their targeting function. The sheer number of piRNA cluster (piC) loci in animal genomes and infrequent piRNA sequence conservation also present challenges in determining which piC loci are most important for development. To address this question, we determined the piRNA expression patterns of piC loci across a wide phylogenetic spectrum of animals, and reveal that most genic and intergenic piC loci evolve rapidly in their capacity to generate piRNAs, regardless of known transposon silencing function. Surprisingly, we also uncovered a distinct set of piC loci with piRNA expression conserved deeply in Eutherian mammals. We name these loci Eutherian-Conserved piRNA cluster (ECpiC) loci. Supporting the hypothesis that conservation of piRNA expression across ~100 million years of Eutherian evolution implies function, we determined that one ECpiC locus generates abundant piRNAs antisense to the STOX1 transcript, a gene clinically associated with preeclampsia. Furthermore, we confirmed reduced piRNAs in existing mouse mutations at ECpiC-Asb1 and -Cbl, which also display spermatogenic defects. The Asb1 mutant testes with strongly reduced Asb1 piRNAs also exhibit up-regulated gene expression profiles. These data indicate ECpiC loci may be specially adapted to support Eutherian reproduction.

The capacity of target silencing by Drosophila PIWI and piRNAs.

Although Piwi proteins and Piwi-interacting RNAs (piRNAs) genetically repress transposable elements (TEs), it is unclear how the highly diverse piRNA populations direct Piwi proteins to silence TE targets without silencing the entire transcriptome. To determine the capacity of piRNA-mediated silencing, we introduced reporter genes into Drosophila OSS cells, which express microRNAs (miRNAs) and piRNAs, and compared the Piwi pathway to the Argonaute pathway in gene regulation. Reporter constructs containing several target sites that were robustly silenced by miRNAs were not silenced to the same degrees by piRNAs. However, another set of reporters we designed to enable a large number of both TE-directed and genic piRNAs to bind were robustly silenced by the PIWI/piRNA complex in OSS cells. These reporters show that a bulk of piRNAs are required to pair to the reporter's transcripts and not the reporter's DNA sequence to engage PIWI-mediated silencing. Following our genome-wide study of PIWI-regulated targets in OSS cells, we assessed candidate gene elements with our reporter platform. These results suggest TE sequences are the most direct of PIWI regulatory targets while coding genes are less directly affected by PIWI targeting. Finally, our study suggests that the PIWI transcriptional silencing mechanism triggers robust chromatin changes on targets with sufficient piRNA binding, and preferentially regulates TE transcripts because protein-coding transcripts lack a threshold of targeting by piRNA populations. This reporter platform will facilitate future dissections of the PIWI-targeting mechanism.

Transposable element dynamics and PIWI regulation impacts lncRNA and gene expression diversity in Drosophila ovarian cell cultures.

Piwi proteins and Piwi-interacting RNAs (piRNAs) repress transposable elements (TEs) from mobilizing in gonadal cells. To determine the spectrum of piRNA-regulated targets that may extend beyond TEs, we conducted a genome-wide survey for transcripts associated with PIWI and for transcripts affected by PIWI knockdown in Drosophila ovarian somatic sheet (OSS) cells, a follicle cell line expressing the Piwi pathway. Despite the immense sequence diversity among OSS cell piRNAs, our analysis indicates that TE transcripts are the major transcripts associated with and directly regulated by PIWI. However, several coding genes were indirectly regulated by PIWI via an adjacent de novo TE insertion that generated a nascent TE transcript. Interestingly, we noticed that PIWI-regulated genes in OSS cells greatly differed from genes affected in a related follicle cell culture, ovarian somatic cells (OSCs). Therefore, we characterized the distinct genomic TE insertions across four OSS and OSC lines and discovered dynamic TE landscapes in gonadal cultures that were defined by a subset of active TEs. Particular de novo TEs appeared to stimulate the expression of novel candidate long noncoding RNAs (lncRNAs) in a cell lineage-specific manner, and some of these TE-associated lncRNAs were associated with PIWI and overlapped PIWI-regulated genes. Our analyses of OSCs and OSS cells demonstrate that despite having a Piwi pathway to suppress endogenous mobile elements, gonadal cell TE landscapes can still dramatically change and create transcriptome diversity.