Cryptic genetic variation (CGV) is the hidden genetic variation that can be unlocked by perturbing normal conditions. CGV can drive the emergence of novel complex phenotypes through changes in gene expression. Although our theoretical understanding of CGV has thoroughly increased over the past decade, insight into polymorphic gene expression regulation underlying CGV is scarce. Here we investigated the transcriptional architecture of CGV in response to rapid temperature changes in the nematode Caenorhabditis elegans. We analyzed regulatory variation in gene expression (and mapped eQTL) across the course of a heat stress and recovery response in a recombinant inbred population. We measured gene expression over three temperature treatments: i) control, ii) heat stress, and iii) recovery from heat stress. Compared to control, exposure to heat stress affected the transcription of 3305 genes, whereas 942 were affected in recovering animals. These affected genes were mainly involved in metabolism and reproduction. The gene expression pattern in recovering animals resembled both the control and the heat-stress treatment. We mapped eQTL using the genetic variation of the recombinant inbred population and detected 2626 genes with an eQTL in the heat-stress treatment, 1797 in the control, and 1880 in the recovery. The cis-eQTL were highly shared across treatments. A considerable fraction of the trans-eQTL (40–57%) mapped to 19 treatment specific trans-bands. In contrast to cis-eQTL, trans-eQTL were highly environment specific and thus cryptic. Approximately 67% of the trans-eQTL were only induced in a single treatment, with heat-stress showing the most unique trans-eQTL. These results illustrate the highly dynamic pattern of CGV across three different environmental conditions that can be evoked by a stress response over a relatively short time-span (2 h) and that CGV is mainly determined by response related trans regulatory eQTL.

West Nile virus (WNV) and dengue virus (DENV) are highly pathogenic, mosquito-borne flaviviruses (family Flaviviridae ) that cause severe disease and death in humans. WNV and DENV actively replicate in mosquitoes and human hosts and thus encounter different host immune responses. RNA interference (RNAi) is the predominant antiviral response against invading RNA viruses in insects and plants. As a countermeasure, plant and insect RNA viruses encode RNA silencing suppressor (RSS) proteins to block the generation/activity of small interfering RNA (siRNA). Enhanced flavivirus replication in mosquitoes depleted for RNAi factors suggests an important biological role for RNAi in restricting virus replication, but it has remained unclear whether or not flaviviruses counteract RNAi via expression of an RSS. First, we established that flaviviral RNA replication suppressed siRNA-induced gene silencing in WNV and DENV replicon-expressing cells. Next, we showed that none of the WNV encoded proteins displayed RSS activity in mammalian and insect cells and in plants by using robust RNAi suppressor assays. In contrast, we found that the 3′-untranslated region-derived RNA molecule known as subgenomic flavivirus RNA (sfRNA) efficiently suppressed siRNA- and miRNA-induced RNAi pathways in both mammalian and insect cells. We also showed that WNV sfRNA inhibits in vitro cleavage of double-stranded RNA by Dicer. The results of the present study suggest a novel role for sfRNA, i.e., as a nucleic acid-based regulator of RNAi pathways, a strategy that may be conserved among flaviviruses.

The Hawaiian strain (CB4856) of Caenorhabditis elegans is one of the most divergent from the canonical laboratory strain N2 and has been widely used in developmental, population, and evolutionary studies. To enhance the utility of the strain, we have generated a draft sequence of the CB4856 genome, exploiting a variety of resources and strategies. When compared against the N2 reference, the CB4856 genome has 327,050 single nucleotide variants (SNVs) and 79,529 insertion–deletion events that result in a total of 3.3 Mb of N2 sequence missing from CB4856 and 1.4 Mb of sequence present in CB4856 but not present in N2. As previously reported, the density of SNVs varies along the chromosomes, with the arms of chromosomes showing greater average variation than the centers. In addition, we find 61 regions totaling 2.8 Mb, distributed across all six chromosomes, which have a greatly elevated SNV density, ranging from 2 to 16% SNVs. A survey of other wild isolates show that the two alternative haplotypes for each region are widely distributed, suggesting they have been maintained by balancing selection over long evolutionary times. These divergent regions contain an abundance of genes from large rapidly evolving families encoding F-box, MATH, BATH, seven-transmembrane G-coupled receptors, and nuclear hormone receptors, suggesting that they provide selective advantages in natural environments. The draft sequence makes available a comprehensive catalog of sequence differences between the CB4856 and N2 strains that will facilitate the molecular dissection of their phenotypic differences. Our work also emphasizes the importance of going beyond simple alignment of reads to a reference genome when assessing differences between genomes.

Arboviruses are transmitted by distantly related arthropod vectors such as mosquitoes (class Insecta) and ticks (class Arachnida). RNA interference (RNAi) is the major antiviral mechanism in arthropods against arboviruses. Unlike in mosquitoes, tick antiviral RNAi is not understood, although this information is important to compare arbovirus/host interactions in different classes of arbovirus vectos. Using an Ixodes scapularis-derived cell line, key Argonaute proteins involved in RNAi and the response against tick-borne Langat virus (Flaviviridae) replication were identified and phylogenetic relationships characterized. Analysis of small RNAs in infected cells showed the production of virus-derived small interfering RNAs (viRNAs), which are key molecules of the antiviral RNAi response. Importantly, viRNAs were longer (22 nucleotides) than those from other arbovirus vectors and mapped at highest frequency to the termini of the viral genome, as opposed to mosquito-borne flaviviruses. Moreover, tick-borne flaviviruses expressed subgenomic flavivirus RNAs that interfere with tick RNAi. Our results characterize the antiviral RNAi response in tick cells including phylogenetic analysis of genes encoding antiviral proteins, and viral interference with this pathway. This shows important differences in antiviral RNAi between the two major classes of arbovirus vectors, and our data broadens our understanding of arthropod antiviral RNAi.