Search

The 2022 multi-country monkeypox outbreak concurrent with the ongoing COVID-19 pandemic has further highlighted the need for genomic surveillance and pathogen whole genome sequencing. While metagenomic sequencing approaches have been used to sequence many of the early human monkeypox virus infections, these methods are resource intensive and require samples with high viral DNA concentrations. Given the atypical clinical presentation of cases associated with the current outbreak and uncertainty regarding viral load across both the course of infection and anatomical body sites, there is an urgent need for a more sensitive and broadly applicable sequencing approach. Amplicon-based sequencing (PrimalSeq) was initially developed for sequencing of Zika virus, and later adapted as the main sequencing approach for SARS-CoV-2. Here, we used PrimalScheme to develop a primer scheme for human monkeypox virus that can be used with many sequencing and bioinformatics pipelines implemented during the COVID-19 pandemic. We sequenced clinical samples that tested presumptive positive for monkeypox virus with amplicon-based and metagenomic sequencing approaches. Upon comparison, we found notably higher genome coverage across the virus genome, particularly in higher PCR cycle threshold (lower DNA titer) samples, with minimal amplicon drop-outs, in using the amplicon-based sequencing approach. By sending out primer pool aliquots to laboratories across the United States and internationally, we validated the primer scheme in 12 public health laboratories with their established Illumina or Oxford Nanopore Technologies sequencing workflows and with different sample types across a range of Ct values. Our findings suggest that amplicon-based sequencing increases the success rate across a wider range of viral DNA concentrations, with the PCR Ct value threshold at which laboratories were able to achieve 80% genome coverage at 10X ranging between Ct 25-33. Therefore, it increases the number of samples where near-complete genomes can be generated and it provides a cost-effective and widely applicable alternative to metagenomics for continued human monkeypox virus genomic surveillance. Importantly, we show that the human monkeypox virus primer scheme can be used with currently implemented amplicon-based SARS-CoV-2 sequencing workflows, with minimal change to the protocol.

The 2022 multi-country monkeypox (mpox) outbreak concurrent with the ongoing COVID-19 pandemic has further highlighted the need for genomic surveillance and rapid pathogen whole genome sequencing. While metagenomic sequencing approaches have been used to sequence many of the early mpox infections, these methods are resource intensive and require samples with high viral DNA concentrations. Given the atypical clinical presentation of cases associated with the outbreak and uncertainty regarding viral load across both the course of infection and anatomical body sites, there was an urgent need for a more sensitive and broadly applicable sequencing approach. Highly multiplexed amplicon-based sequencing (PrimalSeq) was initially developed for sequencing of Zika virus, and later adapted as the main sequencing approach for SARS-CoV-2. Here, we used PrimalScheme to develop a primer scheme for human monkeypox virus that can be used with many sequencing and bioinformatics pipelines implemented in public health laboratories during the COVID-19 pandemic. We sequenced clinical samples that tested presumptive positive for human monkeypox virus with amplicon-based and metagenomic sequencing approaches. We found notably higher genome coverage across the virus genome, with minimal amplicon drop-outs, in using the amplicon-based sequencing approach, particularly in higher PCR cycle threshold (lower DNA titer) samples. Further testing demonstrated that Ct value correlated with the number of sequencing reads and influenced the percent genome coverage. To maximize genome coverage when resources are limited, we recommend selecting samples with a PCR cycle threshold below 31 Ct and generating 1 million sequencing reads per sample. To support national and international public health genomic surveillance efforts, we sent out primer pool aliquots to 10 laboratories across the United States, United Kingdom, Brazil, and Portugal. These public health laboratories successfully implemented the human monkeypox virus primer scheme in various amplicon sequencing workflows and with different sample types across a range of Ct values. Thus, we show that amplicon based sequencing can provide a rapidly deployable, cost-effective, and flexible approach to pathogen whole genome sequencing in response to newly emerging pathogens. Importantly, through the implementation of our primer scheme into existing SARS-CoV-2 workflows and across a range of sample types and sequencing platforms, we further demonstrate the potential of this approach for rapid outbreak response.

Transposable elements are powerful agents of evolution that can diversify transcriptional programs by distributing transcription factor DNA-binding sites throughout genomes. To investigate the extent that transposable elements provide transcription factor-binding motifs in Caenorhabditis elegans, we determined the genomic positions of DNA-binding motifs for 201 different transcription factors. Surprisingly, we found that almost all examined transcription factors have binding motifs that reside within transposable elements, and all types of transposable elements have at least 1 instance of a transcription factor motif, demonstrating that transposable elements provide previously unappreciated numbers of transcription factor-binding motifs to the C. elegans genome. After determining the occurrence of transcription factor motifs in transposable elements relative to the rest of the genome, we identified DNA-binding motifs for 45 different transcription factors that are greater than 20-fold enriched within transposable elements compared to what would be expected by chance. Consistent with potential functional roles for these transposable element-enriched transcription factor-binding sequences, we determined that all transcription factor motif types found in transposable elements have instances of residing within accessible chromatin sites associated with transcription factor binding. The overwhelming majority of transcription factor-binding motifs located within transposable elements associate with their cognate transcription factors, suggesting extensive binding of transcription factors to sequences within transposable elements. In addition, transposable elements with accessible or transcription factor-bound motifs reside in the putative promoter regions of approximately 12% of all protein-coding genes, providing widespread possibilities for influencing gene expression. This work represents the first comprehensive analysis of transposable element–transcription factor interactions in C. elegans and demonstrates that transposable element-provided transcription factor-binding sites are prevalent in this important model organism.

Elevated temperatures activate a heat shock response (HSR) to protect cells from the pathological effects of protein mis-folding, cellular mis-organization, organelle dysfunction and altered membrane fluidity. This response includes activation of the conserved transcription factor heat shock factor 1 (HSF-1), which binds heat shock elements (HSEs) in the promoters of genes induced by heat shock (HS). The upregulation of protein-coding genes (PCGs), such as heat shock proteins and cytoskeletal regulators, is critical for cellular survival during elevated temperatures. While the transcriptional response of PCGs to HS has been comprehensively analyzed in a variety of organisms, the effect of this stress on the expression of non-coding RNAs (ncRNAs) has not been systematically examined. Here we show that in Caenorhabditis elegans HS induces up- and downregulation of specific ncRNAs from multiple classes, including miRNA, piRNA, lincRNA, pseudogene and repeat elements. Moreover, some ncRNA genes appear to be direct targets of the HSR, as they contain HSF-1 bound HSEs in their promoters and their expression is regulated by this factor during HS. These results demonstrate that multiple ncRNA genes respond to HS, some as direct HSF-1 targets, providing new candidates that may contribute to organismal survival during this stress.

Heat Shock Factor 1 (HSF-1) is a key regulator of the heat shock response (HSR). Upon heat shock, HSF-1 binds well-conserved motifs, called Heat Shock Elements (HSEs), and drives expression of genes important for cellular protection during this stress. Remarkably, we found that substantial numbers of HSEs in multiple Caenorhabditis species reside within Helitrons, a type of DNA transposon. Consistent with Helitron-embedded HSEs being functional, upon heat shock they display increased HSF-1 and RNA polymerase II occupancy and up-regulation of nearby genes in C. elegans. Interestingly, we found that different genes appear to be incorporated into the HSR by species-specific Helitron insertions in C. elegans and C. briggsae and by strain-specific insertions among different wild isolates of C. elegans. Our studies uncover previously unidentified targets of HSF-1 and show that Helitron insertions are responsible for rewiring and diversifying the Caenorhabditis HSR.

Elevated temperatures activate a Heat Shock Response (HSR) to protect cells from the pathological effects of protein mis-folding, cellular mis-organization, organelle dysfunction and altered membrane fluidity. This response includes activation of the conserved transcription factor Heat Shock Factor 1 (HSF-1), which binds Heat Shock Elements (HSEs) in the promoters of genes induced by heat shock (HS). The up-regulation of protein-coding genes (PCGs), such as Heat Shock Proteins (HSPs) and cytoskeletal regulators, is critical for cellular survival during elevated temperatures. While the transcriptional response of PCGs to heat shock has been comprehensively analyzed in a variety of organisms, the effect of this stress on the expression of non-coding RNAs (ncRNAs) has not been systematically examined. Here we show that inCaenorhabditis elegansHS induces up- and down-regulation of specific ncRNAs from multiple classes, including miRNA, piRNA, lincRNA, pseudogene, and repeat elements. Moreover, some ncRNA genes appear to be direct targets of the HSR, as they contain HSF-1 bound HSEs in their promoters and their expression is regulated by this factor during HS. These results demonstrate that multiple ncRNA genes respond to HS, some as direct HSF-1 targets, providing new candidates that may contribute to organismal survival during this stress.

Formation of heterochromatin serves a critical role in organizing the genome and regulating gene expression. In most organisms, heterochromatin flanks centromeres and telomeres. To identify heterochromatic regions in the heavily studied model C. elegans, which possesses holocentric chromosomes with dispersed centromeres, we analyzed the genome-wide distribution of the heterochromatin protein 1 (HP1) ortholog HPL-2 and compared its distribution to other features commonly associated with heterochromatin. HPL-2 binding highly correlates with histone H3 mono- and dimethylated at lysine 9 (H3K9me1 and H3K9me2) and forms broad domains on autosomal arms. Although HPL-2, like other HP1 orthologs, binds H3K9me peptides in vitro, the distribution of HPL-2 in vivo appears relatively normal in mutant embryos that lack H3K9me, demonstrating that the chromosomal distribution of HPL-2 can be achieved in an H3K9me-independent manner. Consistent with HPL-2 serving roles independent of H3K9me, hpl-2 mutant worms display more severe defects than mutant worms lacking H3K9me. HPL-2 binding is enriched for repetitive sequences, and on chromosome arms is anticorrelated with centromeres. At the genic level, HPL-2 preferentially associates with well-expressed genes, and loss of HPL-2 results in up-regulation of some binding targets and down-regulation of others. Our work defines heterochromatin in an important model organism and uncovers both shared and distinctive properties of heterochromatin relative to other systems.

Genome function is dynamically regulated in part by chromatin, which consists of the histones, non-histone proteins and RNA molecules that package DNA. Studies in Caenorhabditis elegans and Drosophila melanogaster have contributed substantially to our understanding of molecular mechanisms of genome function in humans, and have revealed conservation of chromatin components and mechanisms. Nevertheless, the three organisms have markedly different genome sizes, chromosome architecture and gene organization. On human and fly chromosomes, for example, pericentric heterochromatin flanks single centromeres, whereas worm chromosomes have dispersed heterochromatin-like regions enriched in the distal chromosomal ‘arms’, and centromeres distributed along their lengths. To systematically investigate chromatin organization and associated gene regulation across species, we generated and analysed a large collection of genome-wide chromatin data sets from cell lines and developmental stages in worm, fly and human. Here we present over 800 new data sets from our ENCODE and modENCODE consortia, bringing the total to over 1,400. Comparison of combinatorial patterns of histone modifications, nuclear lamina-associated domains, organization of large-scale topological domains, chromatin environment at promoters and enhancers, nucleosome positioning, and DNA replication patterns reveals many conserved features of chromatin organization among the three organisms. We also find notable differences in the composition and locations of repressive chromatin. These data sets and analyses provide a rich resource for comparative and species-specific investigations of chromatin composition, organization and function., National Science Foundation (U.S.) (1122374)