Report of the third conference on next-generation sequencing for adventitious virus detection in biologics for humans and animals.

Next-generation sequencing (NGS) has been proven to address some of the limitations of the current testing methods for adventitious virus detection in biologics. The International Alliance for Biological Standardization (IABS), the U.S. Food and Drug Administration (FDA), and the European Directorate for the Quality of Medicines and Healthcare (EDQM) co-organized the "3rd Conference on Next-generation Sequencing for Adventitious Virus Detection in Biologics for Humans and Animals", which was held on September 27-28, 2022, in Rockville, Maryland, U.S.A. The meeting gathered international representatives from regulatory and public health authorities and other government agencies, industry, contract research organizations, and academia to present the current status of NGS applications and the progress on NGS standardization and validation for detection of viral adventitious agents in biologics, including human and animal vaccines, gene therapies, and biotherapeutics. Current regulatory expectations were discussed for developing a scientific consensus regarding using NGS for detection of adventitious viruses. Although there are ongoing improvements in the NGS workflow, the development of reference materials for facilitating method qualification and validation support the current use of NGS for adventitious virus detection.

Report of the 2019 NIST-FDA workshop on standards for next generation sequencing detection of viral adventitious agents in biologics and biomanufacturing.

Adventitious virus testing assures product safety by demonstrating the absence of viruses that could be unintentionally introduced during the manufacturing process. The capabilities of next-generation sequencing (NGS) for broad virus detection in biologics have been demonstrated by the detection of known and novel viruses that were previously missed using the recommended routine assays for adventitious agent testing. A meeting was co-organized by the National Institute of Standards and Technology and the U.S. Food and Drug Administration on September 18-19, 2019 in Gaithersburg, Maryland, USA, to facilitate standardization of NGS technologies for applications of adventitious virus testing in biologics. The goal was to assess the currently used standards for virus detection by NGS and their public availability, and to identify additional needs for different types of reference materials and standards (natural and synthetic). The meeting focused on the NGS processes from sample preparation through sequencing but did not thoroughly cover bioinformatics, since this was considered to be the topic of a separate meeting.

Elevated coding mutation rate during the reprogramming of human somatic cells into induced pluripotent stem cells.

Mutations in human induced pluripotent stem cells (iPSCs) pose a risk for their clinical use due to preferential reprogramming of mutated founder cell and selection of mutations during maintenance of iPSCs in cell culture. It is unknown, however, if mutations in iPSCs are due to stress associated with oncogene expression during reprogramming. We performed whole exome sequencing of human foreskin fibroblasts and their derived iPSCs at two different passages. We found that in vitro passaging contributed 7% to the iPSC coding point mutation load, and ultradeep amplicon sequencing revealed that 19% of the mutations preexist as rare mutations in the parental fibroblasts suggesting that the remaining 74% of the mutations were acquired during cellular reprogramming. Simulation suggests that the mutation intensity during reprogramming is ninefold higher than the background mutation rate in culture. Thus the factor induced reprogramming stress contributes to a significant proportion of the mutation load of iPSCs.

The MEF2A transcription factor interactome in cardiomyocytes.

Transcriptional regulators encoded by the Myocyte Enhancer Factor 2 (MEF2) gene family play a fundamental role in cardiac development, homeostasis and pathology. Previous studies indicate that MEF2A protein-protein interactions serve as a network hub in several cardiomyocyte cellular processes. Based on the idea that interactions with regulatory protein partners underly the diverse roles of MEF2A in cardiomyocyte gene expression, we undertook a systematic unbiased screen of the MEF2A protein interactome in primary cardiomyocytes using an affinity purification-based quantitative mass spectrometry approach. Bioinformatic processing of the MEF2A interactome revealed protein networks involved in the regulation of programmed cell death, inflammatory responses, actin dynamics and stress signaling in primary cardiomyocytes. Further biochemical and functional confirmation of specific protein-protein interactions documented a dynamic interaction between MEF2A and STAT3 proteins. Integration of transcriptome level data from MEF2A and STAT3-depleted cardiomyocytes reveals that the balance between MEF2A and STAT3 activity exerts a level of executive control over the inflammatory response and cardiomyocyte cell survival and experimentally ameliorates Phenylephrine induced cardiomyocyte hypertrophy. Lastly, we identified several MEF2A/STAT3 co-regulated genes, including the MMP9 gene. Herein, we document the cardiomyocyte MEF2A interactome, which furthers our understanding of protein networks involved in the hierarchical control of normal and pathophysiological cardiomyocyte gene expression in the mammalian heart.

Trans-regulation of mouse meiotic recombination hotspots by Rcr1.

Meiotic recombination is required for the orderly segregation of chromosomes during meiosis and for providing genetic diversity among offspring. Among mammals, as well as yeast and higher plants, recombination preferentially occurs at highly delimited chromosomal sites 1-2 kb long known as hotspots. Although considerable progress has been made in understanding the roles various proteins play in carrying out the molecular events of the recombination process, relatively little is understood about the factors controlling the location and relative activity of mammalian recombination hotspots. To search for trans-acting factors controlling the positioning of recombination events, we compared the locations of crossovers arising in an 8-Mb segment of a 100-Mb region of mouse Chromosome 1 (Chr 1) when the longer region was heterozygous C57BL/6J (B6) x CAST/EiJ (CAST) and the remainder of the genome was either similarly heterozygous or entirely homozygous B6. The lack of CAST alleles in the remainder of the genome resulted in profound changes in hotspot activity in both females and males. Recombination activity was lost at several hotspots; new, previously undetected hotspots appeared; and still other hotspots remained unaffected, indicating the presence of distant trans-acting gene(s) whose CAST allele(s) activate or suppress the activity of specific hotspots. Testing the activity of three activated hotspots in sperm samples from individual male progeny of two genetic crosses, we identified a single trans-acting regulator of hotspot activity, designated Rcr1, that is located in a 5.30-Mb interval (11.74-17.04 Mb) on Chr 17. Using an Escherichia coli cloning assay to characterize the molecular products of recombination at two of these hotspots, we found that Rcr1 controls the appearance of both crossover and noncrossover gene conversion events, indicating that it likely controls the sites of the double-strand DNA breaks that initiate the recombination process.

The recombinational anatomy of a mouse chromosome.

Among mammals, genetic recombination occurs at highly delimited sites known as recombination hotspots. They are typically 1-2 kb long and vary as much as a 1,000-fold or more in recombination activity. Although much is known about the molecular details of the recombination process itself, the factors determining the location and relative activity of hotspots are poorly understood. To further our understanding, we have collected and mapped the locations of 5,472 crossover events along mouse Chromosome 1 arising in 6,028 meioses of male and female reciprocal F1 hybrids of C57BL/6J and CAST/EiJ mice. Crossovers were mapped to a minimum resolution of 225 kb, and those in the telomere-proximal 24.7 Mb were further mapped to resolve individual hotspots. Recombination rates were evolutionarily conserved on a regional scale, but not at the local level. There was a clear negative-exponential relationship between the relative activity and abundance of hotspot activity classes, such that a small number of the most active hotspots account for the majority of recombination. Females had 1.2x higher overall recombination than males did, although the sex ratio showed considerable regional variation. Locally, entirely sex-specific hotspots were rare. The initiation of recombination at the most active hotspot was regulated independently on the two parental chromatids, and analysis of reciprocal crosses indicated that parental imprinting has subtle effects on recombination rates. It appears that the regulation of mammalian recombination is a complex, dynamic process involving multiple factors reflecting species, sex, individual variation within species, and the properties of individual hotspots.

Colocalization of somatic and meiotic double strand breaks near the Myc oncogene on mouse chromosome 15.

Both somatic and meiotic recombinations involve the repair of DNA double strand breaks (DSBs) that occur at preferred locations in the genome. Improper repair of DSBs during either mitosis or meiosis can lead to mutations, chromosomal aberration such as translocations, cancer, and/or cell death. Currently, no model exists that explains the locations of either spontaneous somatic DSBs or programmed meiotic DSBs or relates them to each other. One common class of tumorigenic translocations arising from DSBs is chromosomal rearrangements near the Myc oncogene. Myc translocations have been associated with Burkitt lymphoma in humans, plasmacytoma in mice, and immunocytoma in rats. Comparing the locations of somatic and meiotic DSBs near the mouse Myc oncogene, we demonstrated that the placement of these DSBs is not random and that both events clustered in the same short discrete region of the genome. Our work shows that both somatic and meiotic DSBs tend to occur in proximity to each other within the Myc region, suggesting that they share common originating features. It is likely that some regions of the genome are more susceptible to both somatic and meiotic DSBs, and the locations of meiotic hotspots may be an indicator of genomic regions more susceptible to DNA damage.

Sensitivity and breadth of detection of high-throughput sequencing for adventitious virus detection.

High-throughput sequencing (HTS) is capable of broad virus detection encompassing both known and unknown adventitious viruses in a variety of sample matrices. We describe the development of a general-purpose HTS-based method for the detection of adventitious viruses. Performance was evaluated using 16 viruses equivalent to well-characterized National Institutes of Health (NIH) virus stocks and another six viruses of interest. A viral vaccine crude harvest and a cell substrate matrix were spiked with 22 viruses. Specificity was demonstrated for all 22 viruses at the species level. Our method was capable of detecting and identifying adventitious viruses spiked at 104 genome copies per milliliter in a viral vaccine crude harvest and 0.01 viral genome copies spiked per cell in a cell substrate matrix. Moreover, 9 of the 11 NIH model viruses with published in vivo data were detected by HTS with an equivalent or better sensitivity (in a viral vaccine crude harvest). Our general-purpose HTS method is unbiased and highly sensitive for the detection of adventitious viruses, and has a large breadth of detection, which may obviate the need to perform in vivo testing.

Parental origin of chromosomes influences crossover activity within the Kcnq1 transcriptionally imprinted domain of Mus musculus.

BACKGROUND: Among the three functions of DNA, mammalian replication and transcription can be subject to epigenetic imprinting specified by the parental origin of chromosomes, and although there is suggestive indication that this is also true for meiotic recombination, no definitive evidence has yet been reported. RESULTS: We have now obtained such evidence on mouse chromosome 7 by assaying meiotic recombination as it occurs in reciprocal F1 mice. A 166 kb region near the Kcnq1 transcriptionally imprinted domain showed significantly higher recombination activity in the CAST x B6 parental direction (p < 0.03). Characterizing hotspots within this domain revealed a cluster of three hotspots lying within a 100 kb span, among these hotspots, Slc22a18 showed a definitive parent of origin effect on recombination frequency (p < 0.02). Comparing recombination activity in the mouse Kcnq1 and neighboring H19-Igf2 imprinted domains with their human counterparts, we found that elevated recombination activity in these domains is a consequence of their chromosomal position relative to the telomere and not an intrinsic characteristic of transcriptionally imprinted domains as has been previously suggested. CONCLUSION: Similar to replication and transcription, we demonstrate that meiotic recombination can be subjected to epigenetic imprinting and hotspot activity can be influenced by the parental origin of chromosomes. Furthermore, transcriptionally imprinted regions exhibiting elevated recombination activity are likely a consequence of their chromosomal location rather than their transcriptional characteristic.

A quantitative assay for crossover and noncrossover molecular events at individual recombination hotspots in both male and female gametes.

Meiotic recombination is a fundamental process in all eukaryotes. Among organisms in which recombination initiates prior to synapsis, recombination preferentially occurs in short 1-to 2-kb regions, known as recombination hotspots. Among mammals, genotyping sperm DNA has provided a means of monitoring recombination events at specific hotspots in male meiosis. To complement these current techniques, we developed an assay for amplifying all copies of a hotspot from the DNA of male and female germ cells, cloning the products into Escherichia coli, and SNP genotyping the resulting colonies using fluorescence technology. This approach examines the molecular details of crossover and noncrossover events of individual meioses directly at active hotspots while retaining the simplicity of using pooled DNA. Using this technique, we analyzed recombination events at the Hlx1 hotspot located on mouse chromosome 1, finding that the results agree well with a prior genetic characterization of 3026 male and 3002 female meioses.

Selection and evaluation of an efficient method for the recovery of viral nucleic acids from complex biologicals [corrected].

There is a need for a broad and efficient testing strategy for the detection of both known and novel viral adventitious agents in vaccines and biologicals. High-throughput sequencing (HTS) is an approach for such testing; however, an optimized testing method is one with a sample-processing pipeline that can help detect any viral adventitious agent that may be present. In this study, 11 commercial methods were assessed for efficient extraction of nucleic acids from a panel of viruses. An extraction strategy with two parallel arms, consisting of both the Invitrogen PureLink™ Virus RNA/DNA kit for total nucleic acid extraction and the Wako DNA Extractor® kit with an RNase A digestion for enrichment of double-stranded nucleic acid, was selected as the strategy for the extraction of all viral nucleic acid types (ssRNA, dsRNA, and dsDNA). Downstream processes, such as double-strand DNA synthesis and whole-genome amplification (WGA), were also assessed for the retrieval of viral sequences. Double-stranded DNA synthesis yielded larger numbers of viral reads, whereas WGA exhibited a strong bias toward amplification of double-stranded DNA, including host cellular DNA. The final sample-processing strategy consisted of the dual extraction approach followed by double-stranded DNA synthesis, which yielded a viral population with increased detection of some viruses by 8600-fold. Here we describe an efficient extraction procedure to support viral adventitious agent detection in cell substrates used for biological products using HTS.

Erratum: Author Correction: Selection and evaluation of an efficient method for the recovery of viral nucleic acids from complex biologicals.

[This corrects the article DOI: 10.1038/s41541-018-0067-3.].

Current Perspectives on High-Throughput Sequencing (HTS) for Adventitious Virus Detection: Upstream Sample Processing and Library Preparation.

A key step for broad viral detection using high-throughput sequencing (HTS) is optimizing the sample preparation strategy for extracting viral-specific nucleic acids since viral genomes are diverse: They can be single-stranded or double-stranded RNA or DNA, and can vary from a few thousand bases to over millions of bases, which might introduce biases during nucleic acid extraction. In addition, viral particles can be enveloped or non-enveloped with variable resistance to pre-treatment, which may influence their susceptibility to extraction procedures. Since the identity of the potential adventitious agents is unknown prior to their detection, efficient sample preparation should be unbiased toward all different viral types in order to maximize the probability of detecting any potential adventitious viruses using HTS. Furthermore, the quality assessment of each step for sample processing is also a critical but challenging aspect. This paper presents our current perspectives for optimizing upstream sample processing and library preparation as part of the discussion in the Advanced Virus Detection Technologies Interest group (AVDTIG) The topics include: use of nuclease treatment to enrich for encapsidated nucleic acids, techniques for amplifying low amounts of virus nucleic acids, selection of different extraction methods, relevant controls, the use of spike recovery experiments, and quality control measures during library preparation.

A Multicenter Study To Evaluate the Performance of High-Throughput Sequencing for Virus Detection.

The capability of high-throughput sequencing (HTS) for detection of known and unknown viruses makes it a powerful tool for broad microbial investigations, such as evaluation of novel cell substrates that may be used for the development of new biological products. However, like any new assay, regulatory applications of HTS need method standardization. Therefore, our three laboratories initiated a study to evaluate performance of HTS for potential detection of viral adventitious agents by spiking model viruses in different cellular matrices to mimic putative materials for manufacturing of biologics. Four model viruses were selected based upon different physical and biochemical properties and commercial availability: human respiratory syncytial virus (RSV), Epstein-Barr virus (EBV), feline leukemia virus (FeLV), and human reovirus (REO). Additionally, porcine circovirus (PCV) was tested by one laboratory. Independent samples were prepared for HTS by spiking intact viruses or extracted viral nucleic acids, singly or mixed, into different HeLa cell matrices (resuspended whole cells, cell lysate, or total cellular RNA). Data were obtained using different sequencing platforms (Roche 454, Illumina HiSeq1500 or HiSeq2500). Bioinformatic analyses were performed independently by each laboratory using available tools, pipelines, and databases. The results showed that comparable virus detection was obtained in the three laboratories regardless of sample processing, library preparation, sequencing platform, and bioinformatic analysis: between 0.1 and 3 viral genome copies per cell were detected for all of the model viruses used. This study highlights the potential for using HTS for sensitive detection of adventitious viruses in complex biological samples containing cellular background. IMPORTANCE Recent high-throughput sequencing (HTS) investigations have resulted in unexpected discoveries of known and novel viruses in a variety of sample types, including research materials, clinical materials, and biological products. Therefore, HTS can be a powerful tool for supplementing current methods for demonstrating the absence of adventitious or unwanted viruses in biological products, particularly when using a new cell line. However, HTS is a complex technology with different platforms, which needs standardization for evaluation of biologics. This collaborative study was undertaken to investigate detection of different virus types using two different HTS platforms. The results of the independently performed studies demonstrated a similar sensitivity of virus detection, regardless of the different sample preparation and processing procedures and bioinformatic analyses done in the three laboratories. Comparable HTS detection of different virus types supports future development of reference virus materials for standardization and validation of different HTS platforms.

Advanced Virus Detection Technologies Interest Group (AVDTIG): Efforts on High Throughput Sequencing (HTS) for Virus Detection.

Several nucleic-acid based technologies have recently emerged with capabilities for broad virus detection. One of these, high throughput sequencing, has the potential for novel virus detection because this method does not depend upon prior viral sequence knowledge. However, the use of high throughput sequencing for testing biologicals poses greater challenges as compared to other newly introduced tests due to its technical complexities and big data bioinformatics. Thus, the Advanced Virus Detection Technologies Users Group was formed as a joint effort by regulatory and industry scientists to facilitate discussions and provide a forum for sharing data and experiences using advanced new virus detection technologies, with a focus on high throughput sequencing technologies. The group was initiated as a task force that was coordinated by the Parenteral Drug Association and subsequently became the Advanced Virus Detection Technologies Interest Group to continue efforts for using new technologies for detection of adventitious viruses with broader participation, including international government agencies, academia, and technology service providers.

Cataloguing the taxonomic origins of sequences from a heterogeneous sample using phylogenomics: applications in adventitious agent detection.

We have designed and implemented a software system, named PhyloID™, that can be used to detect putative adventitious agents in biological samples characterized by next-generation sequencing. PhyloID is run in two steps, each being a self-contained automated process amenable to GMP validation. The first module, MiLY, is responsible for assembling individual sequence reads into contigs, and annotating all sequences with a unique sequence identifier, the number of reads in each contig, and the length of the sequence. The trimmed, assembled and annotated data are then processed by PhyloID's second module, NGmapper. NGmapper takes the FASTA-formatted output from MiLY and identifies the taxonomic origins of the contigs and singletons therein. It compares each sequence's BLASTN hit profile against the patterns of evolutionary relationships described within phylogenomic distance matrices for all of the various taxonomic groups, in order to find the best fit. NGmapper then produces lists of taxonomic assignments in both summarized and detailed form, and tree files for viewing results graphically. We optimized PhyloID's parameters and measured its performance using simulated metagenomic data and subsets of the reference phylogenies. PhyloID's precision and recall in identifying simulated sequences were measured by information retrieval analysis, focusing on read length, read number, sequence accuracy, background complexity, taxonomy and reference data coverage. We found PhyloID to be highly accurate and quantitative in its taxonomic mapping of sequences, with excellent precision, sensitivity and robustness. The degree of taxonomic representation available in publicly available databases remains an issue, as expected, for any sequence classifier, but community sequencing efforts are poised to overcome this problem. In order to illustrate real-world usage of the application, we also describe some simple spike-recovery experiments as well as a multi-site comparative characterization of a viral suspension. These data help to illustrate, to corroborate, and to extend results using simulated data. LAY ABSTRACT: In order to address gaps in the detection of contaminating viruses and microorganisms in vaccines and other biologicals, manufacturers are exploring the use of new technologies that promise greater sensitivity and breadth of coverage. One challenge in implementing such new methods is the complexity of analysis of the "big data" generated by these new instruments: hundreds of millions of sequence reads (segments of genetic material from viruses and cells) need to be compared against a vast and growing number of entries in genetic databases, in order to come up with a confident identification. These large-scale analyses must furthermore be carried out within the strict regulatory environment that governs the industry. We have developed an automated software pipeline named PhyloID™ that is capable of identifying viruses and microorganisms from large-scale sequence data. Using simulated data as well as real samples, we show that PhyloID is both sensitive and accurate in identifying any type of potential contaminant. Such a powerful new assay will be an important addition to the adventitious agent testing package, providing further assurance about product safety.

Preliminary Evaluation of Next-Generation Sequencing Performance Relative to qPCR and In Vitro Cell Culture Tests for Human Cytomegalovirus.

To compare the performances of conventional in vitro indicator cell culture, quantitative polymerase chain reaction (qPCR), and next-generation sequencing (NGS) as detection methods for adventitious agents, a preliminary assessment was performed using human cytomegalovirus (HCMV) as a model virus. HCMV was spiked into a crude viral harvest at various concentrations and inoculated onto MRC-5 cell monolayers. The cultures were observed for cytopathic effects (CPEs) as per the compendial method requirements, and samples were taken at various time points for analysis by qPCR or NGS. When using NGS, the detection of HCMV is 10 fold more sensitive than observed using the conventional CPE endpoint method. It may be possible for qPCR to achieve the detection level demonstrated by NGS, but further optimization of the technique would be required. In addition, NGS was able to detect HCMV in the initial inoculum when it was spiked in at 10 CCID50/mL, suggesting the potential to eliminate cell culture amplification with an NGS-based assay. This study highlights the advantage of NGS as a surrogate broad-spectrum technology for the detection of adventitious agents in biologics. LAY ABSTRACT: Human cytomeglovirus (HCMV) is highly prevalent in the general population and can lead to serious health issues in both immumocompromised individuals and/or newborns. The testing of HCMV in biological materials is stipulated by multiple regulatory agencies where HCMV is a potential risk. This test involves inoculating cell lines that are susceptible to HCMV infection, incubating the cultures for 28 days, and observing the cells for signs of viral infection. Next-generation sequencing and quantitative polymerase chain reaction (qPCR) are two technologies that can potentially shorten the extended 28 day cell culture incubation. In this study, we compared the sensitivity of the compendial cell culture assay with NGS and qPCR for the detection of HCMV. Our results show that NGS can potentially achieve a detection limit that is 10 times more sensitive than the cell culture assay. In addition, NGS was able to detect HCMV in the initial inoculum, potentially eliminating the need for cell culture amplification of the virus. Finally, sequence data generated by NGS directly demonstrate the presence of HCMV, and such information can serve as the foundation for any follow-up investigation.

A physical map of the genome of Atlantic salmon, Salmo salar.

A physical map of the Atlantic salmon (Salmo salar) genome was generated based on HindIII fingerprints of a publicly available BAC (bacterial artificial chromosome) library constructed from DNA isolated from a Norwegian male. Approximately 11.5 haploid genome equivalents (185,938 clones) were successfully fingerprinted. Contigs were first assembled via FPC using high-stringency (1e-16), and then end-to-end joins yielded 4354 contigs and 37,285 singletons. The accuracy of the contig assembly was verified by hybridization and PCR analysis using genetic markers. A subset of the BACs in the library contained few or no HindIII recognition sites in their insert DNA. BglI digestion fragment patterns of these BACs allowed us to identify three classes: (1) BACs containing histone genes, (2) BACs containing rDNA-repeating units, and (3) those that do not have BglI recognition sites. End-sequence analysis of selected BACs representing these three classes confirmed the identification of the first two classes and suggested that the third class contained highly repetitive DNA corresponding to tRNAs and related sequences.

Development and application of a salmonid EST database and cDNA microarray: data mining and interspecific hybridization characteristics.

We report 80,388 ESTs from 23 Atlantic salmon (Salmo salar) cDNA libraries (61,819 ESTs), 6 rainbow trout (Oncorhynchus mykiss) cDNA libraries (14,544 ESTs), 2 chinook salmon (Oncorhynchus tshawytscha) cDNA libraries (1317 ESTs), 2 sockeye salmon (Oncorhynchus nerka) cDNA libraries (1243 ESTs), and 2 lake whitefish (Coregonus clupeaformis) cDNA libraries (1465 ESTs). The majority of these are 3' sequences, allowing discrimination between paralogs arising from a recent genome duplication in the salmonid lineage. Sequence assembly reveals 28,710 different S. salar, 8981 O. mykiss, 1085 O. tshawytscha, 520 O. nerka, and 1176 C. clupeaformis putative transcripts. We annotate the submitted portion of our EST database by molecular function. Higher- and lower-molecular-weight fractions of libraries are shown to contain distinct gene sets, and higher rates of gene discovery are associated with higher-molecular weight libraries. Pyloric caecum library group annotations indicate this organ may function in redox control and as a barrier against systemic uptake of xenobiotics. A microarray is described, containing 7356 salmonid elements representing 3557 different cDNAs. Analyses of cross-species hybridizations to this cDNA microarray indicate that this resource may be used for studies involving all salmonids.