Assessing nanobody interaction with SARS-CoV-2 Nsp9.

The interaction between SARS-CoV-2 non-structural protein Nsp9 and the nanobody 2NSP90 was investigated by NMR spectroscopy using the paramagnetic perturbation methodology PENELOP (Paramagnetic Equilibrium vs Nonequilibrium magnetization Enhancement or LOss Perturbation). The Nsp9 monomer is an essential component of the replication and transcription complex (RTC) that reproduces the viral gRNA for subsequent propagation. Therefore preventing Nsp9 recruitment in RTC would represent an efficient antiviral strategy that could be applied to different coronaviruses, given the Nsp9 relative invariance. The NMR results were consistent with a previous characterization suggesting a 4:4 Nsp9-to-nanobody stoichiometry with the occurrence of two epitope pairs on each of the Nsp9 units that establish the inter-dimer contacts of Nsp9 tetramer. The oligomerization state of Nsp9 was also analyzed by molecular dynamics simulations and both dimers and tetramers resulted plausible. A different distribution of the mapped epitopes on the tetramer surface with respect to the former 4:4 complex could also be possible, as well as different stoichiometries of the Nsp9-nanobody assemblies such as the 2:2 stoichiometry suggested by the recent crystal structure of the Nsp9 complex with 2NSP23 (PDB ID: 8dqu), a nanobody exhibiting essentially the same affinity as 2NSP90. The experimental NMR evidence, however, ruled out the occurrence in liquid state of the relevant Nsp9 conformational change observed in the same crystal structure.

Evolution of a biological thermocouple by adaptation of cytochrome c oxidase in a subterrestrial metazoan.

In this study we report a naturally evolved temperature-sensing electrical regulator in the cytochrome c oxidase of the Devil Worm, Halicephalobus mephisto. This extremophile metazoan was isolated 1.3 km underground in a South African goldmine, where it adapted to heat and potentially to hypoxia, making its mitochondrial sequence a likely target of adaptational change. We obtained the full mitochondrial genome sequence of this organism, and show through dN/dS analysis statistically robust evidence of positive selection in H. mephisto cytochrome c oxidase subunits. Seventeen of these positively-selected amino acid substitutions were localized in proximity to the H- and K-pathway proton channels of the complex. Surprisingly, the H. mephisto cytochrome c oxidase proton pump completely shuts down at low temperatures (20°C) leading to approximately a 4.8-fold reduction in the transmembrane proton gradient voltage (ΔΨm) compared to optimal temperature (37°C). Direct measurement of oxygen consumption found a corresponding 4.7-fold drop at 20°C compared to 37°C. Correspondingly, the lifecycle of H. mephisto takes four-fold longer at the low temperature compared to higher. This elegant evolutionary adaptation creates a finely-tuned mitochondrial temperature sensor, allowing this ectothermic organism to maximize its reproductive success in varying environmental temperatures. Our study shows that evolutionary innovation may remodel core metabolism to make it more accurately map onto environmental variation.

ß-actin mediated H3K27ac changes demonstrate the link between compartment switching and enhancer-dependent transcriptional regulation.

BACKGROUND: Recent work has demonstrated that three-dimensional genome organization is directly affected by changes in the levels of nuclear cytoskeletal proteins such as ß-actin. The mechanisms which translate changes in 3D genome structure into changes in transcription, however, are not fully understood. Here, we use a comprehensive genomic analysis of cells lacking nuclear ß-actin to investigate the mechanistic links between compartment organization, enhancer activity, and gene expression. RESULTS: Using HiC-Seq, ATAC-Seq, and RNA-Seq, we first demonstrate that transcriptional and chromatin accessibility changes observed upon ß-actin loss are highly enriched in compartment-switching regions. Accessibility changes within compartment switching genes, however, are mainly observed in non-promoter regions which potentially represent distal regulatory elements. Our results also show that ß-actin loss induces widespread accumulation of the enhancer-specific epigenetic mark H3K27ac. Using the ABC model of enhancer annotation, we then establish that these epigenetic changes have a direct impact on enhancer activity and underlie transcriptional changes observed upon compartment switching. A complementary analysis of fibroblasts undergoing reprogramming into pluripotent stem cells further confirms that this relationship between compartment switching and enhancer-dependent transcriptional change is not specific to ß-actin knockout cells but represents a general mechanism linking compartment-level genome organization to gene expression. CONCLUSIONS: We demonstrate that enhancer-dependent transcriptional regulation plays a crucial role in driving gene expression changes observed upon compartment-switching. Our results also reveal a novel function of nuclear ß-actin in regulating enhancer function by influencing H3K27 acetylation levels.

A statistical framework for high-content phenotypic profiling using cellular feature distributions.

High-content screening (HCS) uses microscopy images to generate phenotypic profiles of cell morphological data in high-dimensional feature space. While HCS provides detailed cytological information at single-cell resolution, these complex datasets are usually aggregated into summary statistics that do not leverage patterns of biological variability within cell populations. Here we present a broad-spectrum HCS analysis system that measures image-based cell features from 10 cellular compartments across multiple assay panels. We introduce quality control measures and statistical strategies to streamline and harmonize the data analysis workflow, including positional and plate effect detection, biological replicates analysis and feature reduction. We also demonstrate that the Wasserstein distance metric is superior over other measures to detect differences between cell feature distributions. With this workflow, we define per-dose phenotypic fingerprints for 65 mechanistically diverse compounds, provide phenotypic path visualizations for each compound and classify compounds into different activity groups.

The Caenorhabditis elegans TDRD5/7-like protein, LOTR-1, interacts with the helicase ZNFX-1 to balance epigenetic signals in the germline.

LOTUS and Tudor domain containing proteins have critical roles in the germline. Proteins that contain these domains, such as Tejas/Tapas in Drosophila, help localize the Vasa helicase to the germ granules and facilitate piRNA-mediated transposon silencing. The homologous proteins in mammals, TDRD5 and TDRD7, are required during spermiogenesis. Until now, proteins containing both LOTUS and Tudor domains in Caenorhabditis elegans have remained elusive. Here we describe LOTR-1 (D1081.7), which derives its name from its LOTUS and Tudor domains. Interestingly, LOTR-1 docks next to P granules to colocalize with the broadly conserved Z-granule helicase, ZNFX-1. The Tudor domain of LOTR-1 is required for its Z-granule retention. Like znfx-1 mutants, lotr-1 mutants lose small RNAs from the 3' ends of WAGO and mutator targets, reminiscent of the loss of piRNAs from the 3' ends of piRNA precursor transcripts in mouse Tdrd5 mutants. Our work shows that LOTR-1 acts with ZNFX-1 to bring small RNA amplifying mechanisms towards the 3' ends of its RNA templates.

NMR-Based Analysis of Nanobodies to SARS-CoV-2 Nsp9 Reveals a Possible Antiviral Strategy Against COVID-19.

Following the entry into the host cell, SARS-CoV-2 replication is mediated by the replication transcription complex (RTC) assembled through a number of nonstructural proteins (Nsps). A monomeric form of Nsp9 is particularly important for RTC assembly and function. In the present study, 136 unique nanobodies targeting Nsp9 are generated. Several nanobodies belonging to different B-cell lineages are expressed, purified, and characterized. Results from immunoassays applied to purified Nsp9 and neat saliva from coronavirus disease (COVID-19) patients show that these nanobodies effectively and specifically recognize both recombinant and endogenous Nsp9. Nuclear magnetic resonance analyses supported by molecular dynamics reveal a composite Nsp9 oligomerization pattern and demonstrate that both nanobodies stabilize the tetrameric form of wild-type Nsp9 also identifying the epitopes on the tetrameric assembly. These results can have important implications in the potential use of these nanobodies to combat viral replication.

ß-actin dependent chromatin remodeling mediates compartment level changes in 3D genome architecture.

ß-actin is a crucial component of several chromatin remodeling complexes that control chromatin structure and accessibility. The mammalian Brahma-associated factor (BAF) is one such complex that plays essential roles in development and differentiation by regulating the chromatin state of critical genes and opposing the repressive activity of polycomb repressive complexes (PRCs). While previous work has shown that ß-actin loss can lead to extensive changes in gene expression and heterochromatin organization, it is not known if changes in ß-actin levels can directly influence chromatin remodeling activities of BAF and polycomb proteins. Here we conduct a comprehensive genomic analysis of ß-actin knockout mouse embryonic fibroblasts (MEFs) using ATAC-Seq, HiC-seq, RNA-Seq and ChIP-Seq of various epigenetic marks. We demonstrate that ß-actin levels can induce changes in chromatin structure by affecting the complex interplay between chromatin remodelers such as BAF/BRG1 and EZH2. Our results show that changes in ß-actin levels and associated chromatin remodeling activities can not only impact local chromatin accessibility but also induce reversible changes in 3D genome architecture. Our findings reveal that ß-actin-dependent chromatin remodeling plays a role in shaping the chromatin landscape and influences the regulation of genes involved in development and differentiation.

Novel LOTUS-domain proteins are organizational hubs that recruit C. elegans Vasa to germ granules.

We describe MIP-1 and MIP-2, novel paralogous C. elegans germ granule components that interact with the intrinsically disordered MEG-3 protein. These proteins promote P granule condensation, form granules independently of MEG-3 in the postembryonic germ line, and balance each other in regulating P granule growth and localization. MIP-1 and MIP-2 each contain two LOTUS domains and intrinsically disordered regions and form homo- and heterodimers. They bind and anchor the Vasa homolog GLH-1 within P granules and are jointly required for coalescence of MEG-3, GLH-1, and PGL proteins. Animals lacking MIP-1 and MIP-2 show temperature-sensitive embryonic lethality, sterility, and mortal germ lines. Germline phenotypes include defects in stem cell self-renewal, meiotic progression, and gamete differentiation. We propose that these proteins serve as scaffolds and organizing centers for ribonucleoprotein networks within P granules that help recruit and balance essential RNA processing machinery to regulate key developmental transitions in the germ line.

Pheniqs 2.0: accurate, high-performance Bayesian decoding and confidence estimation for combinatorial barcode indexing.

BACKGROUND: Systems biology increasingly relies on deep sequencing with combinatorial index tags to associate biological sequences with their sample, cell, or molecule of origin. Accurate data interpretation depends on the ability to classify sequences based on correct decoding of these combinatorial barcodes. The probability of correct decoding is influenced by both sequence quality and the number and arrangement of barcodes. The rising complexity of experimental designs calls for a probability model that accounts for both sequencing errors and random noise, generalizes to multiple combinatorial tags, and can handle any barcoding scheme. The needs for reproducibility and community benchmark standards demand a peer-reviewed tool that preserves decoding quality scores and provides tunable control over classification confidence that balances precision and recall. Moreover, continuous improvements in sequencing throughput require a fast, parallelized and scalable implementation. RESULTS AND DISCUSSION: We developed a flexible, robustly engineered software that performs probabilistic decoding and supports arbitrarily complex barcoding designs. Pheniqs computes the full posterior decoding error probability of observed barcodes by consulting basecalling quality scores and prior distributions, and reports sequences and confidence scores in Sequence Alignment/Map (SAM) fields. The product of posteriors for multiple independent barcodes provides an overall confidence score for each read. Pheniqs achieves greater accuracy than minimum edit distance or simple maximum likelihood estimation, and it scales linearly with core count to enable the classification of > 11 billion reads in 1 h 15 m using < 50 megabytes of memory. Pheniqs has been in production use for seven years in our genomics core facility. CONCLUSION: We introduce a computationally efficient software that implements both probabilistic and minimum distance decoders and show that decoding barcodes using posterior probabilities is more accurate than available methods. Pheniqs allows fine-tuning of decoding sensitivity using intuitive confidence thresholds and is extensible with alternative decoders and new error models. Any arbitrary arrangement of barcodes is easily configured, enabling computation of combinatorial confidence scores for any barcoding strategy. An optimized multithreaded implementation assures that Pheniqs is faster and scales better with complex barcode sets than existing tools. Support for POSIX streams and multiple sequencing formats enables easy integration with automated analysis pipelines.

FRET-Based Probe for High-Throughput DNA Intercalator Drug Discovery and In Vivo Imaging.

Molecules that bind DNA by intercalating its bases remain among the most potent cancer therapies and antimicrobials due to their interference with DNA-processing proteins. To accelerate the discovery of novel intercalating drugs, we designed a fluorescence resonance energy transfer (FRET)-based probe that reports on DNA intercalation, allowing rapid and sensitive screening of chemical libraries in a high-throughput format. We demonstrate that the method correctly identifies known DNA intercalators in approved drug libraries and discover previously unreported intercalating compounds. When introduced in cells, the oligonucleotide-based probe rapidly distributes in the nucleus, allowing direct imaging of the dynamics of drug entry and its interaction with DNA in its native environment. This enabled us to directly correlate the potency of intercalators in killing cultured cancer cells with the ability of the drug to penetrate the cell membrane. The combined capability of the single probe to identify intercalators in vitro and follow their function in vivo can play a valuable role in accelerating the discovery of novel DNA-intercalating drugs or repurposing approved ones.

Structural Modeling of the SARS-CoV-2 Spike/Human ACE2 Complex Interface can Identify High-Affinity Variants Associated with Increased Transmissibility.

The COVID-19 pandemic has triggered concerns about the emergence of more infectious and pathogenic viral strains. As a public health measure, efficient screening methods are needed to determine the functional effects of new sequence variants. Here we show that structural modeling of SARS-CoV-2 Spike protein binding to the human ACE2 receptor, the first step in host-cell entry, predicts many novel variant combinations with enhanced binding affinities. By focusing on natural variants at the Spike-hACE2 interface and assessing over 700 mutant complexes, our analysis reveals that high-affinity Spike mutations (including N440K, S443A, G476S, E484R, G502P) tend to cluster near known human ACE2 recognition sites (K31 and K353). These Spike regions are structurally flexible, allowing certain mutations to optimize interface interaction energies. Although most human ACE2 variants tend to weaken binding affinity, they can interact with Spike mutations to generate high-affinity double mutant complexes, suggesting variation in individual susceptibility to infection. Applying structural analysis to highly transmissible variants, we find that circulating point mutations S477N, E484K and N501Y form high-affinity complexes (~40% more than wild-type). By combining predicted affinities and available antibody escape data, we show that fast-spreading viral variants exploit combinatorial mutations possessing both enhanced affinity and antibody resistance, including S477N/E484K, E484K/N501Y and K417T/E484K/N501Y. Thus, three-dimensional modeling of the Spike/hACE2 complex predicts changes in structure and binding affinity that correlate with transmissibility and therefore can help inform future intervention strategies.

The Microbiome of the Lebanese Wild Apple, Malus trilobata, is a Rich Source of Potential Biocontrol Agents for Fungal Post-harvest Pathogens of Apples.

The widespread use of harmful fungicides in the agricultural sector has led to a demand for safer alternatives to protect against crop pathogens. The domestic apple is the second most highly consumed fruit in the world and encounters several pre- and post-harvest fungal and bacterial phytopathogens. The goal of this study was to explore the uncharacterized microbiome of a wild apple, Malus trilobata, as a potential source of novel biocontrol agents for two post-harvest fungi that affect commercial apples: Botrytis cinerea and Penicillium expansum. We sampled microflora associated with the leaves, bulk soil, and roots of Malus trilobata in two regions of Lebanon: Ehden reserve in the north and Dhour EL Choueir near Beirut. The two regions have different soil types Dhour EL Choueir and samples from the two regions showed very different microbial compositions, with greater microbial diversity among those from Ehden reserve. Molecular characterization revealed a wide variety of genera displaying activity against the two fungal pathogens, including several with previously unknown antifungal activity: Bosea, Microlunatus, Microbacterium, Mycetecola, Rhizobium and Paraphoma. In total, 92 strains inhibited Penicillium expansum (39%) and 87 strains inhibited Botrytis cinerea (38%) out of 237 screened. Further chemical and genetic characterization of one or more selected strains could pave the way for future development of new biocontrol agents for post-harvest applications.

Loss of ß-Actin Leads to Accelerated Mineralization and Dysregulation of Osteoblast-Differentiation Genes during Osteogenic Reprogramming.

Actin plays fundamental roles in both the cytoplasm and the cell nucleus. In the nucleus, ß-actin regulates neuronal reprogramming by consolidating a heterochromatin landscape required for transcription of neuronal gene programs, yet it remains unknown whether it has a role in other differentiation models. To explore the potential roles of ß-actin in osteogenesis, ß-actin wild-type (WT) and ß-actin knockout (KO) mouse embryonic fibroblasts (MEFs) are reprogrammed to osteoblast-like cells using small molecules in vitro. It is discovered that loss of ß-actin leads to an accelerated mineralization phenotype (hypermineralization), accompanied with enhanced formation of extracellular hydroxyapatite microcrystals, which originate in the mitochondria in the form of microgranules. This phenotype is a consequence of rapid upregulation of mitochondrial genes including those involved in oxidative phosphorylation (OXPHOS) in reprogrammed KO cells. It is further found that osteogenic gene programs are differentially regulated between WT and KO cells, with clusters of genes exhibiting different temporal expression patterns. A novel function for ß-actin in osteogenic reprogramming through a mitochondria-based mechanism that controls cell-mediated mineralization is proposed.

Diatom modulation of select bacteria through use of two unique secondary metabolites.

Unicellular eukaryotic phytoplankton, such as diatoms, rely on microbial communities for survival despite lacking specialized compartments to house microbiomes (e.g., animal gut). Microbial communities have been widely shown to benefit from diatom excretions that accumulate within the microenvironment surrounding phytoplankton cells, known as the phycosphere. However, mechanisms that enable diatoms and other unicellular eukaryotes to nurture specific microbiomes by fostering beneficial bacteria and repelling harmful ones are mostly unknown. We hypothesized that diatom exudates may tune microbial communities and employed an integrated multiomics approach using the ubiquitous diatom Asterionellopsis glacialis to reveal how it modulates its naturally associated bacteria. We show that A. glacialis reprograms its transcriptional and metabolic profiles in response to bacteria to secrete a suite of central metabolites and two unusual secondary metabolites, rosmarinic acid and azelaic acid. While central metabolites are utilized by potential bacterial symbionts and opportunists alike, rosmarinic acid promotes attachment of beneficial bacteria to the diatom and simultaneously suppresses the attachment of opportunists. Similarly, azelaic acid enhances growth of beneficial bacteria while simultaneously inhibiting growth of opportunistic ones. We further show that the bacterial response to azelaic acid is numerically rare but globally distributed in the world's oceans and taxonomically restricted to a handful of bacterial genera. Our results demonstrate the innate ability of an important unicellular eukaryotic group to modulate select bacteria in their microbial consortia, similar to higher eukaryotes, using unique secondary metabolites that regulate bacterial growth and behavior inversely across different bacterial populations.

NASQAR: a web-based platform for high-throughput sequencing data analysis and visualization.

BACKGROUND: As high-throughput sequencing applications continue to evolve, the rapid growth in quantity and variety of sequence-based data calls for the development of new software libraries and tools for data analysis and visualization. Often, effective use of these tools requires computational skills beyond those of many researchers. To ease this computational barrier, we have created a dynamic web-based platform, NASQAR (Nucleic Acid SeQuence Analysis Resource). RESULTS: NASQAR offers a collection of custom and publicly available open-source web applications that make extensive use of a variety of R packages to provide interactive data analysis and visualization. The platform is publicly accessible at http://nasqar.abudhabi.nyu.edu/ . Open-source code is on GitHub at https://github.com/nasqar/NASQAR , and the system is also available as a Docker image at https://hub.docker.com/r/aymanm/nasqarall . NASQAR is a collaboration between the core bioinformatics teams of the NYU Abu Dhabi and NYU New York Centers for Genomics and Systems Biology. CONCLUSIONS: NASQAR empowers non-programming experts with a versatile and intuitive toolbox to easily and efficiently explore, analyze, and visualize their Transcriptomics data interactively. Popular tools for a variety of applications are currently available, including Transcriptome Data Preprocessing, RNA-seq Analysis (including Single-cell RNA-seq), Metagenomics, and Gene Enrichment.

Phenotyping of the thrashing forces exerted by partially immobilized C. elegans using elastomeric micropillar arrays.

As a simple model organism, C. elegans plays an important role in gaining insight into the relationship between bodily thrashing forces and biological effects, such as disease and aging, or physical stimuli, like touch and light. Due to their similar length scale, microfluidic chips have been extensively explored for use in various biological studies involving C. elegans. However, a formidable challenge still exists due to the complexity of integrating external stimuli (chemical, mechanical or optical) with free-moving worms and subsequent imaging on the chip. In this report, we use a microfluidic device to partially immobilize a worm, which allows for measurements of the relative changes in the thrashing force under different assay conditions. Using a device adapted to the natural escape-like coiling response of a worm to immobilization, we have quantified the relative changes in the thrashing force during different developmental stages (L1, L3, L4, and young adult) and in response to various glucose concentrations and drug treatment. Our findings showed a loss of thrashing force following the introduction of glucose into a wild type worm culture that could be reversed upon treatment with the type 2 diabetes drug metformin. A morphological study of the actin filament structures in the body wall muscles provided supporting evidence for the force measurement data. Finally, we demonstrated the multiplexing capabilities of our device through recording the thrashing activities of eight worms simultaneously. The multiplexing capabilities and facile imaging available using our device open the door for high-throughput neuromuscular studies using C. elegans.

The Role of Tertiary Structure in MicroRNA Target Recognition.

Translational repression and degradation of transcripts by microRNAs (miRNAs) is mediated by a ribonucleoprotein complex called the miRNA-induced silencing complex (miRISC, or RISC). Advances in experimental determination of RISC structures have enabled detailed analysis and modeling of known miRNA targets, yet a full appreciation of the structural factors influencing target recognition remains a challenge, primarily because target recognition involves a combination of RNA-RNA and RNA-protein interactions that can vary greatly among different miRNA-target pairs. In this chapter, we review progress toward understanding the role of tertiary structure in miRNA target recognition using computational approaches to assemble RISC complexes at known targets and physics-based methods for computing target interactions. Using this framework to examine RISC structures and dynamics, we describe how the conformational flexibility of Argonautes plays an important role in accommodating the diversity of miRNA-target duplexes formed at canonical and noncanonical target sites. We then discuss applications of tertiary structure-based approaches to emerging topics, including the structural effects of SNPs in miRNA targets and cooperative interactions involving Argonaute-Argonaute complexes. We conclude by assessing the prospects for genome-scale modeling of RISC structures and modeling of higher-order Argonaute complexes associated with miRNA biogenesis, mRNA regulation, and other functions.

miR-122 and Ago interactions with the HCV genome alter the structure of the viral 5' terminus.

Hepatitis C virus (HCV) is a positive-sense RNA virus that interacts with the liver-specific microRNA, miR-122. miR-122 binds to two sites in the 5' untranslated region (UTR) and this interaction promotes HCV RNA accumulation, although the precise role of miR-122 in the HCV life cycle remains unclear. Using biophysical analyses and Selective 2' Hydroxyl Acylation analyzed by Primer Extension (SHAPE) we investigated miR-122 interactions with the 5' UTR. Our data suggests that miR-122 binding results in alteration of nucleotides 1-117 to suppress an alternative secondary structure and promote functional internal ribosomal entry site (IRES) formation. Furthermore, we demonstrate that two hAgo2:miR-122 complexes are able to bind to the HCV 5' terminus simultaneously and SHAPE analyses revealed further alterations to the structure of the 5' UTR to accommodate these complexes. Finally, we present a computational model of the hAgo2:miR-122:HCV RNA complex at the 5' terminus of the viral genome as well as hAgo2:miR-122 interactions with the IRES-40S complex that suggest hAgo2 is likely to form additional interactions with SLII which may further stabilize the HCV IRES. Taken together, our results support a model whereby hAgo2:miR-122 complexes alter the structure of the viral 5' terminus and promote formation of the HCV IRES.

Tissue- and sex-specific small RNAomes reveal sex differences in response to the environment.

RNA interference (RNAi) related pathways are essential for germline development and fertility in metazoa and can contribute to inter- and trans-generational inheritance. In the nematode Caenorhabditis elegans, environmental double-stranded RNA provided by feeding can lead to heritable changes in phenotype and gene expression. Notably, transmission efficiency differs between the male and female germline, yet the underlying mechanisms remain elusive. Here we use high-throughput sequencing of dissected gonads to quantify sex-specific endogenous piRNAs, miRNAs and siRNAs in the C. elegans germline and the somatic gonad. We identify genes with exceptionally high levels of secondary 22G RNAs that are associated with low mRNA expression, a signature compatible with silencing. We further demonstrate that contrary to the hermaphrodite germline, the male germline, but not male soma, is resistant to environmental RNAi triggers provided by feeding, in line with previous work. This sex-difference in silencing efficacy is associated with lower levels of gonadal RNAi amplification products. Moreover, this tissue- and sex-specific RNAi resistance is regulated by the germline, since mutant males with a feminized germline are RNAi sensitive. This study provides important sex- and tissue-specific expression data of miRNA, piRNA and siRNA as well as mechanistic insights into sex-differences of gene regulation in response to environmental cues.