Detection of Mpox Virus Using Microbial Cell-Free DNA: The Potential of Pathogen-Agnostic Sequencing for Rapid Identification of Emerging Pathogens

BACKGROUND

The 2022 global outbreak of Monkeypox virus (MPXV) highlighted challenges with polymerase chain reaction detection as divergent strains emerged and atypical presentations limited the applicability of swab sampling. Recommended testing in the United States requires a swab of lesions, which arise late in infection and may be unrecognized. We present MPXV detections using plasma microbial cell-free DNA (mcfDNA) sequencing.

METHODS

Fifteen plasma samples from 12 case-patients were characterized through mcfDNA sequencing. Assay performance was confirmed through in silico inclusivity and exclusivity assessments. MPXV isolates were genotyped using mcfDNA, and phylodynamic information was imputed using publicly available sequences.

RESULTS

MPXV mcfDNA was detected in 12 case-patients. Mpox was not suspected in 5, with 1 having documented resolution of mpox >6 months previously. Six had moderate to severe mpox, supported by high MPXV mcfDNA concentrations; 4 died. In 7 case-patients, mcfDNA sequencing detected coinfections. Genotyping by mcfDNA sequencing identified 22 MPXV mutations at 10 genomic loci in 9 case-patients. Consistent with variation observed in the 2022 outbreak, 21 of 22 variants were G > A/C > T. Phylogenetic analyses imputed isolates to sublineages arising at different time points and from different geographic locations.

CONCLUSIONS

We demonstrate the potential of plasma mcfDNA sequencing to detect, quantify, and, for acute infections with high sequencing coverage, subtype MPXV using a single noninvasive test. Sequencing plasma mcfDNA may augment existing mpox testing in vulnerable patient populations or in patients with atypical symptoms or unrecognized mpox. Strain type information may supplement disease surveillance and facilitate tracking emerging pathogens.

Cataloging Human PRDM9 Allelic Variation Using Long-Read Sequencing Reveals PRDM9 Population Specificity and Two Distinct Groupings of Related Alleles

The PRDM9 protein determines sites of meiotic recombination in humans by directing meiotic DNA double-strand breaks to specific loci. Targeting specificity is encoded by a long array of C₂H₂ zinc fingers that bind to DNA. This zinc finger array is hypervariable, and the resulting alleles each have a potentially different DNA binding preference. The assessment of PRDM9 diversity is important for understanding the complexity of human population genetics, inheritance linkage patterns, and predisposition to genetic disease. Due to the repetitive nature of the PRDM9 zinc finger array, the large-scale sequencing of human PRDM9 is challenging. We, therefore, developed a long-read sequencing strategy to infer the diploid PRDM9 zinc finger array genotype in a high-throughput manner. From an unbiased study of PRDM9 allelic diversity in 720 individuals from seven human populations, we detected 69 PRDM9 alleles. Several alleles differ in frequency among human populations, and 32 alleles had not been identified by previous studies, which were heavily biased to European populations. PRDM9 alleles are distinguished by their DNA binding site preferences and fall into two major categories related to the most common PRDM9-A and PRDM9-C alleles. We also found that it is likely that inter-conversion between allele types is rare. By mapping meiotic double-strand breaks (DSBs) in the testis, we found that small variations in PRDM9 can substantially alter the meiotic recombination landscape, demonstrating that minor PRDM9 variants may play an under-appreciated role in shaping patterns of human recombination. In summary, our data greatly expands knowledge of PRDM9 diversity in humans.

Rat PRDM9 shapes recombination landscapes, duration of meiosis, gametogenesis, and age of fertility

BACKGROUND

Vertebrate meiotic recombination events are concentrated in regions (hotspots) that display open chromatin marks, such as trimethylation of lysines 4 and 36 of histone 3 (H3K4me3 and H3K36me3). Mouse and human PRDM9 proteins catalyze H3K4me3 and H3K36me3 and determine hotspot positions, whereas other vertebrates lacking PRDM9 recombine in regions with chromatin already opened for another function, such as gene promoters. While these other vertebrate species lacking PRDM9 remain fertile, inactivation of the mouse Prdm9 gene, which shifts the hotspots to the functional regions (including promoters), typically causes gross fertility reduction; and the reasons for these species differences are not clear.

RESULTS

We introduced Prdm9 deletions into the Rattus norvegicus genome and generated the first rat genome-wide maps of recombination-initiating double-strand break hotspots. Rat strains carrying the same wild-type Prdm9 allele shared 88% hotspots but strains with different Prdm9 alleles only 3%. After Prdm9 deletion, rat hotspots relocated to functional regions, about 40% to positions corresponding to Prdm9-independent mouse hotspots, including promoters. Despite the hotspot relocation and decreased fertility, Prdm9-deficient rats of the SHR/OlaIpcv strain produced healthy offspring. The percentage of normal pachytene spermatocytes in SHR-Prdm9 mutants was almost double than in the PWD male mouse oligospermic sterile mutants. We previously found a correlation between the crossover rate and sperm presence in mouse Prdm9 mutants. The crossover rate of SHR is more similar to sperm-carrying mutant mice, but it did not fully explain the fertility of the SHR mutants. Besides mild meiotic arrests at rat tubular stages IV (mid-pachytene) and XIV (metaphase), we also detected postmeiotic apoptosis of round spermatids. We found delayed meiosis and age-dependent fertility in both sexes of the SHR mutants.

CONCLUSIONS

We hypothesize that the relative increased fertility of rat versus mouse Prdm9 mutants could be ascribed to extended duration of meiotic prophase I. While rat PRDM9 shapes meiotic recombination landscapes, it is unnecessary for recombination. We suggest that PRDM9 has additional roles in spermatogenesis and speciation-spermatid development and reproductive age-that may help to explain male-specific hybrid sterility.

After the break: DSB end processing in mouse meiosis

This Outlook discusses findings by Yamada et al. that provide some of the first genome-wide insights into double-strand break resection and repair in a mammalian genome.

The exchange of genetic information between parental chromosomes in meiosis is an integral process for the creation of gametes. To generate a crossover, hundreds of DNA double-strand breaks (DSBs) are introduced in the genome of each meiotic cell by the SPO11 protein. The nucleolytic resection of DSB-adjacent DNA is a key step in meiotic DSB repair, but this process has remained understudied. In this issue of Genes & Development, Yamada and colleagues (pp. 806–818) capture some of the first details of resection and DSB repair intermediates in mouse meiosis using a method that maps blunt-ended DNA after ssDNA digestion. This yields some of the first genome-wide insights into DSB resection and repair in a mammalian genome and offers a tantalizing glimpse of how to quantitatively dissect this difficult to study, yet integral, nuclear process.

Cell-type-specific genomics reveals histone modification dynamics in mammalian meiosis

Meiosis is the specialized cell division during which parental genomes recombine to create genotypically unique gametes. Despite its importance, mammalian meiosis cannot be studied in vitro, greatly limiting mechanistic studies. In vivo, meiocytes progress asynchronously through meiosis and therefore the study of specific stages of meiosis is a challenge. Here, we describe a method for isolating pure sub-populations of nuclei that allows for detailed study of meiotic substages. Interrogating the H3K4me3 landscape revealed dynamic chromatin transitions between substages of meiotic prophase I, both at sites of genetic recombination and at gene promoters. We also leveraged this method to perform the first comprehensive, genome-wide survey of histone marks in meiotic prophase, revealing a heretofore unappreciated complexity of the epigenetic landscape at meiotic recombination hotspots. Ultimately, this study presents a straightforward, scalable framework for interrogating the complexities of mammalian meiosis.

Meiotic DSB formation, repair and recombination occur in a continuum of substages termed leptonema, zygonema, pachynema, and diplonema. Here, authors develop a method for isolating pure sub-populations of nuclei that allows for detailed study of meiotic substages.

Histone methyltransferase PRDM9 is not essential for meiosis in male mice

A hallmark of meiosis is the rearrangement of parental alleles to ensure genetic diversity in the gametes. These chromosome rearrangements are mediated by the repair of programmed DNA double-strand breaks (DSBs) as genetic crossovers between parental homologs. In mice, humans, and many other mammals, meiotic DSBs occur primarily at hotspots, determined by sequence-specific binding of the PRDM9 protein. Without PRDM9, meiotic DSBs occur near gene promoters and other functional sites. Studies in a limited number of mouse strains showed that functional PRDM9 is required to complete meiosis, but despite its apparent importance, Prdm9 has been repeatedly lost across many animal lineages. Both the reason for mouse sterility in the absence of PRDM9 and the mechanism by which Prdm9 can be lost remain unclear. Here, we explore whether mice can tolerate the loss of Prdm9. By generating Prdm9 functional knockouts in an array of genetic backgrounds, we observe a wide range of fertility phenotypes and ultimately demonstrate that PRDM9 is not required for completion of male meiosis. Although DSBs still form at a common subset of functional sites in all mice lacking PRDM9, meiotic outcomes differ substantially. We speculate that DSBs at functional sites are difficult to repair as a crossover and that by increasing the efficiency of crossover formation at these sites, genetic modifiers of recombination rates can allow for meiotic progression. This model implies that species with a sufficiently high recombination rate may lose Prdm9 yet remain fertile.

REC114 Partner ANKRD31 Controls Number, Timing, and Location of Meiotic DNA Breaks

Double-strand breaks (DSBs) initiate the homologous recombination that is crucial for meiotic chromosome pairing and segregation. Here, we unveil mouse ANKRD31 as a lynchpin governing multiple aspects of DSB formation. Spermatocytes lacking ANKRD31 have altered DSB locations and fail to target DSBs to the pseudoautosomal regions (PARs) of sex chromosomes. They also have delayed and/or fewer recombination sites but, paradoxically, more DSBs, suggesting DSB dysregulation. Unrepaired DSBs and pairing failures-stochastic on autosomes, nearly absolute on X and Y-cause meiotic arrest and sterility in males. Ankrd31-deficient females have reduced oocyte reserves. A crystal structure defines a pleckstrin homology (PH) domain in REC114 and its direct intermolecular contacts with ANKRD31. In vivo, ANKRD31 stabilizes REC114 association with the PAR and elsewhere. Our findings inform a model in which ANKRD31 is a scaffold anchoring REC114 and other factors to specific genomic locations, thereby regulating DSB formation.

Analysis of Meiotic Double-Strand Break Initiation in Mammals

The repair of programmed DNA double-strand breaks (DSBs) physically tethers homologous chromosomes in meiosis to allow for accurate segregation through meiotic cell divisions. This process, known as recombination, also results in the exchange of alleles between parental chromosomes and contributes to genetic diversity. In mammals, meiotic DSBs occur predominantly in a small fraction of the genome, at sites known as hotspots. Studies of the formation and repair of meiotic DSBs in mammals are challenging, because few cells undergo meiotic DSB formation at a given time. To better understand the initiation and control of meiotic recombination in mammals, we have devised a highly sensitive method to map the sites of meiotic DSBs genome wide. Our method first isolates DNA bound to DSB repair proteins and then specifically sequences the associated single-stranded DNA. This protocol has generated the first meiotic DSB maps in several mammals and the only map of meiotic DSBs in humans.

Erratum: The evolutionary turnover of recombination hot spots contributes to speciation in mice

Due to a technical error in processing the figures in the above-mentioned article, Figures 3, A and B; 4B; 5B; and 6, A and C contained errors or missing elements. The errors have been corrected in both the PDF and full-text HTML files online.

The evolutionary turnover of recombination hot spots contributes to speciation in mice

In this study, Smagulova et al investigate how the speciation gene in mammals, Prdm9, regulates hot spots of genetic recombination and how this contributes to speciation. The authors generated a comprehensive panel of high-quality maps of recombination hot spots across major mouse subspecies and their hybrids and provide novel insights into how the action of PRDM9 and the evolutionary turnover of recombination hot spots alter the genome and lead to incompatibilities in hybrids.

Meiotic recombination is required for the segregation of homologous chromosomes and is essential for fertility. In most mammals, the DNA double-strand breaks (DSBs) that initiate meiotic recombination are directed to a subset of genomic loci (hot spots) by sequence-specific binding of the PRDM9 protein. Rapid evolution of the DNA-binding specificity of PRDM9 and gradual erosion of PRDM9-binding sites by gene conversion will alter the recombination landscape over time. To better understand the evolutionary turnover of recombination hot spots and its consequences, we mapped DSB hot spots in four major subspecies of Mus musculus with different Prdm9 alleles and in their F1 hybrids. We found that hot spot erosion governs the preferential usage of some Prdm9 alleles over others in hybrid mice and increases sequence diversity specifically at hot spots that become active in the hybrids. As crossovers are disfavored at such hot spots, we propose that sequence divergence generated by hot spot turnover may create an impediment for recombination in hybrids, potentially leading to reduced fertility and, eventually, speciation.

DNA recombination. Recombination initiation maps of individual human genomes

DNA double-strand breaks (DSBs) are introduced in meiosis to initiate recombination and generate crossovers, the reciprocal exchanges of genetic material between parental chromosomes. Here, we present high-resolution maps of meiotic DSBs in individual human genomes. Comparing DSB maps between individuals shows that along with DNA binding by PRDM9, additional factors may dictate the efficiency of DSB formation. We find evidence for both GC-biased gene conversion and mutagenesis around meiotic DSB hotspots, while frequent colocalization of DSB hotspots with chromosome rearrangement breakpoints implicates the aberrant repair of meiotic DSBs in genomic disorders. Furthermore, our data indicate that DSB frequency is a major determinant of crossover rate. These maps provide new insights into the regulation of meiotic recombination and the impact of meiotic recombination on genome function.

Suppression of genetic recombination in the pseudoautosomal region and at subtelomeres in mice with a hypomorphic Spo11 allele

BACKGROUND

Homologous recombination is the key process that generates genetic diversity and drives evolution. SPO11 protein triggers recombination by introducing DNA double stranded breaks at discreet areas of the genome called recombination hotspots. The hotspot locations are largely determined by the DNA binding specificity of the PRDM9 protein in human, mice and most other mammals. In budding yeast Saccharomyces cerevisae, which lacks a Prdm9 gene, meiotic breaks are formed opportunistically in the regions of accessible chromatin, primarily at gene promoters. The genome-wide distribution of hotspots in this organism can be altered by tethering Spo11 protein to Gal4 recognition sequences in the strain expressing Spo11 attached to the DNA binding domain of the Gal4 transcription factor. To establish whether similar re-targeting of meiotic breaks can be achieved in PRDM9-containing organisms we have generated a Gal4BD-Spo11 mouse that expresses SPO11 protein joined to the DNA binding domain of yeast Gal4.

RESULTS

We have mapped the genome-wide distribution of the recombination initiation sites in the Gal4BD-Spo11 mice. More than two hundred of the hotspots in these mice were novel and were likely defined by Gal4BD, as the Gal4 consensus motif was clustered around the centers in these hotspots. Surprisingly, meiotic DNA breaks in the Gal4BD-Spo11 mice were significantly depleted near the ends of chromosomes. The effect is particularly striking at the pseudoautosomal region of the X and Y chromosomes - normally the hottest region in the genome.

CONCLUSIONS

Our data suggest that specific, yet-unidentified factors influence the initiation of meiotic recombination at subtelomeric chromosomal regions.

Genome-wide analysis reveals novel molecular features of mouse recombination hotspots

Meiotic recombination predominantly occurs at discrete genomic loci called recombination hotspots, but the features defining these areas are still largely unknown (reviewed in refs 1-5). To allow a comprehensive analysis of hotspot-associated DNA and chromatin characteristics, we developed a direct molecular approach for mapping meiotic DNA double-strand breaks that initiate recombination. Here we present the genome-wide distribution of recombination initiation sites in the mouse genome. Hotspot centres are mapped with approximately 200-nucleotide precision, which allows analysis of the fine structural details of the preferred recombination sites. We determine that hotspots share a centrally distributed consensus motif, possess a nucleotide skew that changes polarity at the centres of hotspots and have an intrinsic preference to be occupied by a nucleosome. Furthermore, we find that the vast majority of recombination initiation sites in mouse males are associated with testis-specific trimethylation of lysine 4 on histone H3 that is distinct from histone H3 lysine 4 trimethylation marks associated with transcription. The recombination map presented here has been derived from a homogeneous mouse population with a defined genetic background and therefore lends itself to extensive future experimental exploration. We note that the mapping technique developed here does not depend on the availability of genetic markers and hence can be easily adapted to other species with complex genomes. Our findings uncover several fundamental features of mammalian recombination hotspots and underline the power of the new recombination map for future studies of genetic recombination, genome stability and evolution.

Mediation of CTCF transcriptional insulation by DEAD-box RNA-binding protein p68 and steroid receptor RNA activator SRA

CCCTC-binding factor (CTCF) is a DNA-binding protein that plays important roles in chromatin organization, although the mechanism by which CTCF carries out these functions is not fully understood. Recent studies show that CTCF recruits the cohesin complex to insulator sites and that cohesin is required for insulator activity. Here we showed that the DEAD-box RNA helicase p68 (DDX5) and its associated noncoding RNA, steroid receptor RNA activator (SRA), form a complex with CTCF that is essential for insulator function. p68 was detected at CTCF sites in the IGF2/H19 imprinted control region (ICR) as well as other genomic CTCF sites. In vivo depletion of SRA or p68 reduced CTCF-mediated insulator activity at the IGF2/H19 ICR, increased levels of IGF2 expression, and increased interactions between the endodermal enhancer and IGF2 promoter. p68/SRA also interacts with members of the cohesin complex. Depletion of either p68 or SRA does not affect CTCF binding to its genomic sites, but does reduce cohesin binding. The results suggest that p68/SRA stabilizes the interaction of cohesin with CTCF by binding to both, and is required for proper insulator function.

Revisiting the Plasmodium falciparum RIFIN family: from comparative genomics to 3D-model prediction

Background

Subtelomeric RIFIN genes constitute the most abundant multigene family in Plasmodium falciparum. RIFIN products are targets for the human immune response and contribute to the antigenic variability of the parasite. They are transmembrane proteins grouped into two sub-families (RIF_A and RIF_B). Although recent data show that RIF_A and RIF_B have different sub-cellular localisations and possibly different functions, the same structural organisation has been proposed for members of the two sub-families. Despite recent advances, our knowledge of the regulation of RIFIN gene expression is still poor and the biological role of the protein products remain obscure.

Results

Comparative studies on RIFINs in three clones of P. falciparum (3D7, HB3 and Dd2) by Multidimensional scaling (MDS) showed that gene sequences evolve differently in the 5'upstream, coding, and 3'downstream regions, and suggested a possible role of highly conserved 3' downstream sequences. Despite the expected polymorphism, we found that the overall structure of RIFIN repertoires is conserved among clones suggesting a balance between genetic drift and homogenisation mechanisms which guarantees emergence of novel variants but preserves the functionality of genes. Protein sequences from a bona fide set of 3D7 RIFINs were submitted to predictors of secondary structure elements. In contrast with the previously proposed structural organisation, no signal peptide and only one transmembrane helix were predicted for the majority of RIF_As. Finally, we developed a strategy to obtain a reliable 3D-model for RIF_As. We generated 265 possible structures from 53 non-redundant sequences, from which clustering and quality assessments selected two models as the most representative for putative RIFIN protein structures.

Conclusion

First, comparative analyses of RIFIN repertoires in different clones of P. falciparum provide insights on evolutionary mechanisms shaping the multigene family. Secondly, we found that members of the two sub-families RIF_As and RIF_Bs have different structural organization in accordance with recent experimental results. Finally, representative models for RIF_As have an "Armadillo-like" fold which is known to promote protein-protein interactions in diverse contexts.

Core promoters are predicted by their distinct physicochemical properties in the genome of Plasmodium falciparum

A method is presented to computationally identify core promoters in the Plasmodium falciparum genome using only DNA physicochemical properties.

Little is known about the structure and distinguishing features of core promoters in Plasmodium falciparum. In this work, we describe the first method to computationally identify core promoters in this AT-rich genome. This prediction algorithm uses solely DNA physicochemical properties as descriptors. Our results add to a growing body of evidence that a physicochemical code for eukaryotic genomes plays a crucial role in core promoter recognition.

A novel series of compositionally biased substitution matrices for comparing Plasmodium proteins

Background

The most common substitution matrices currently used (BLOSUM and PAM) are based on protein sequences with average amino acid distributions, thus they do not represent a fully accurate substitution model for proteins characterized by a biased amino acid composition. This problem has been addressed recently by adjusting existing matrices, however, to date, no empirical approach has been taken to build matrices which offer a substitution model for comparing proteins sharing an amino acid compositional bias. Here, we present a novel procedure to construct series of symmetrical substitution matrices to align proteins from similarly biased Plasmodium proteomes.

Results

We generated substitution matrices by selecting from the BLOCKS database those multiple alignments with a compositional bias similar to that of P. falciparum and P. yoelii proteins. A novel 'fuzzy' clustering method was adopted to group sequences within these alignments, showing that this method retains more complete information on the amino acid substitutions when compared to hierarchical clustering. We assessed the performance against the BLOSUM62 series and showed that the usage of our matrices results in an improvement in the performance of BLAST database searches, greatly reducing the number of false positive hits. We then demonstrated applications of the use of novel matrices to improve the annotation of homologs between the two Plasmodium species and to classify members of the P. falciparum RIFIN/STEVOR family.

Conclusion

We confirmed that in the case of compositionally biased proteins, standard BLOSUM matrices are not suited for optimal alignments, and specific substitution matrices are required. In addition, we showed that the usage of these matrices leads to a reduction of false positive hits, facilitating the automatic annotation process.