Early Detection of Severe Acute Respiratory Syndrome Coronavirus 2 Variants Using Traveler-based Genomic Surveillance at 4 US Airports, September 2021-January 2022.

We enrolled arriving international air travelers in a severe acute respiratory syndrome coronavirus 2 genomic surveillance program. We used molecular testing of pooled nasal swabs and sequenced positive samples for sublineage. Traveler-based surveillance provided early-warning variant detection, reporting the first US Omicron BA.2 and BA.3 in North America.

Lineage abundance estimation for SARS-CoV-2 in wastewater using transcriptome quantification techniques.

Effectively monitoring the spread of SARS-CoV-2 mutants is essential to efforts to counter the ongoing pandemic. Predicting lineage abundance from wastewater, however, is technically challenging. We show that by sequencing SARS-CoV-2 RNA in wastewater and applying algorithms initially used for transcriptome quantification, we can estimate lineage abundance in wastewater samples. We find high variability in signal among individual samples, but the overall trends match those observed from sequencing clinical samples. Thus, while clinical sequencing remains a more sensitive technique for population surveillance, wastewater sequencing can be used to monitor trends in mutant prevalence in situations where clinical sequencing is unavailable.

Variant abundance estimation for SARS-CoV-2 in wastewater using RNA-Seq quantification.

Effectively monitoring the spread of SARS-CoV-2 variants is essential to efforts to counter the ongoing pandemic. Wastewater monitoring of SARS-CoV-2 RNA has proven an effective and efficient technique to approximate COVID-19 case rates in the population. Predicting variant abundances from wastewater, however, is technically challenging. Here we show that by sequencing SARS-CoV-2 RNA in wastewater and applying computational techniques initially used for RNA-Seq quantification, we can estimate the abundance of variants in wastewater samples. We show by sequencing samples from wastewater and clinical isolates in Connecticut U.S.A. between January and April 2021 that the temporal dynamics of variant strains broadly correspond. We further show that this technique can be used with other wastewater sequencing techniques by expanding to samples taken across the United States in a similar timeframe. We find high variability in signal among individual samples, and limited ability to detect the presence of variants with clinical frequencies <10%; nevertheless, the overall trends match what we observed from sequencing clinical samples. Thus, while clinical sequencing remains a more sensitive technique for population surveillance, wastewater sequencing can be used to monitor trends in variant prevalence in situations where clinical sequencing is unavailable or impractical.

Designing and Implementing NGS Tests for Inherited Disorders: A Practical Framework with Step-by-Step Guidance for Clinical Laboratories.

Comprehensive next-generation sequencing (NGS) tests are increasingly used as first-line tests in the evaluation of patients with suspected heritable disease. Despite major technical simplifications, these assays still pose significant challenges for molecular testing laboratories. Existing professional guidelines and recommendations provide a framework for laboratories implementing such tests, but in-depth, concrete guidance is generally not provided. Consequently, there is variability in how laboratories interpret and subsequently implement these regulatory frameworks. To address the need for more detailed guidance, the College of American Pathologists with representation from the Association for Molecular Pathologists assembled a working group to create a practical resource for clinical laboratories. This initial work is focused on variant detection in the setting of inherited disease and provides structured worksheets that guide the user through the entire life cycle of an NGS test, including design, optimization, validation, and quality management with additional guidance for clinical bioinformatics. This resource is designed to be a living document that is publicly available and will be updated with user and expert feedback as the wet bench and bioinformatic landscapes continue to evolve. It is intended to facilitate the standardization of NGS testing across laboratories and therefore to improve patient care.

A Follow-Up of the Multicenter Collaborative Study on HIV-1 Drug Resistance and Tropism Testing Using 454 Ultra Deep Pyrosequencing.

BACKGROUND: Ultra deep sequencing is of increasing use not only in research but also in diagnostics. For implementation of ultra deep sequencing assays in clinical laboratories for routine diagnostics, intra- and inter-laboratory testing are of the utmost importance. METHODS: A multicenter study was conducted to validate an updated assay design for 454 Life Sciences' GS FLX Titanium system targeting protease/reverse transcriptase (RTP) and env (V3) regions to identify HIV-1 drug-resistance mutations and determine co-receptor use with high sensitivity. The study included 30 HIV-1 subtype B and 6 subtype non-B samples with viral titers (VT) of 3,940-447,400 copies/mL, two dilution series (52,129-1,340 and 25,130-734 copies/mL), and triplicate samples. Amplicons spanning PR codons 10-99, RT codons 1-251 and the entire V3 region were generated using barcoded primers. Analysis was performed using the GS Amplicon Variant Analyzer and geno2pheno for tropism. For comparison, population sequencing was performed using the ViroSeq HIV-1 genotyping system. RESULTS: The median sequencing depth across the 11 sites was 1,829 reads per position for RTP (IQR 592-3,488) and 2,410 for V3 (IQR 786-3,695). 10 preselected drug resistant variants were measured across sites and showed high inter-laboratory correlation across all sites with data (P<0.001). The triplicate samples of a plasmid mixture confirmed the high inter-laboratory consistency (mean% ± stdev: 4.6 ±0.5, 4.8 ±0.4, 4.9 ±0.3) and revealed good intra-laboratory consistency (mean% range ± stdev range: 4.2-5.2 ± 0.04-0.65). In the two dilutions series, no variants >20% were missed, variants 2-10% were detected at most sites (even at low VT), and variants 1-2% were detected by some sites. All mutations detected by population sequencing were also detected by UDS. CONCLUSIONS: This assay design results in an accurate and reproducible approach to analyze HIV-1 mutant spectra, even at variant frequencies well below those routinely detectable by population sequencing.

Validation of a next-generation-sequencing cancer panel for use in the clinical laboratory.

CONTEXT: Next-generation sequencing allows for high-throughput processing and sensitive variant detection in multiple genes from small samples. For many diseases, including cancer, a comprehensive mutational profile of a targeted list of genes can be used to simultaneously inform patient care, establish eligibility for ongoing clinical trials, and further research. OBJECTIVE: To validate a pan-cancer, next-generation-sequencing assay for use in the clinical laboratory. DESIGN: DNA was extracted from 68 clinical specimens (formalin-fixed, paraffin-embedded; fine-needle aspirates; peripheral blood; or bone marrow) and 5 normal controls. Sixty-four DNA samples (94%; 64 of 68) were successfully processed with the TruSeq Amplicon Cancer Panel (Illumina Inc, San Diego, California) and sequenced in 4 sequencing runs. The data were analyzed at 4 different filter settings for sequencing coverage and variant frequency cutoff. RESULTS: Libraries created from 40 specimens could be successfully sequenced in a single run and still yield sufficient coverage for robust data analysis of individual samples. Sensitivity for mutation detection down to 5% was demonstrated using dilutions of clinical specimens and control samples. The test was highly repeatable and reproducible and showed 100% concordance with clinically validated Sanger sequencing results. Comparison to an alternate next-generation sequencing technology was performed by also processing 9 of the specimens with the AmpliSeq Cancer Hotspot Panel (version 2; Life Technologies, Grand Island, New York). Thirty of the 31 (97%) TruSeq-detected variants covered by the designs of both panels were confirmed. CONCLUSIONS: A sensitive, high-throughput, pan-cancer mutation panel for sequencing of cancer hot-spot mutations in 42 genes was validated for routine use in clinical testing.

Deep sequencing of HIV-1 variants from paired plasma and cerebrospinal fluid during primary HIV infection.

BACKGROUND: Limited data exist comparing viral quasispecies between cerebrospinal fluid (CSF) and plasma compartments during primary HIV infection. Deep sequencing is a new method to examine the HIV plasma and CSF quasispecies. METHODS: In this pilot study, deep sequencing of protease (PR) and reverse transcriptase (RT) was performed in plasma and CSF from participants during primary HIV infection. Estimated mutational load was calculated by mutant variant frequency multiplied by HIV-RNA level. RESULTS: Paired plasma and CSF samples were studied from five antiretroviral therapy-naïve male participants with median 109 days post estimated transmission, age 32 years, CD4 cell count 580 cells/µL, HIV-RNA 5.18 log10 copies/mL in plasma and 3.67 log10 copies/mL in CSF. Plasma samples averaged 7,124 reads of PR and 2,448 reads of RT, whereas CSF samples averaged 7,082 and 2,792 reads, respectively. A distinct drug-resistance pattern with linked mutations present at significant levels (5-10%) was detected in one participant in CSF. Other low abundance variants (>0.2%) were detected in plasma and CSF of four out of five participants. CONCLUSIONS: Deep sequencing of CSF HIV is technically possible with sufficient HIV-RNA levels. Differences between the quasispecies in the two compartments detected in one participant, which were present with a high mutational load in CSF at an estimated 3.6 months after HIV infection, suggest that early CNS compartmentalisation may be revealed by sensitive deep-sequencing methods. The presence of distinct low abundance (<1%) resistance variants in plasma and CSF of three other subjects may be significant, but further investigation is needed.

Low-frequency NNRTI-resistant HIV-1 variants and relationship to mutational load in antiretroviral-naïve subjects.

Low-frequency HIV variants possessing resistance mutations against non­nucleoside reverse transcriptase inhibitors (NNRTI), especially at HIV reverse transcriptase (RT) amino acid (aa) positions K103 and Y181, have been shown to adversely affect treatment response. Therapeutic failure correlates with both the mutant viral variant frequency and the mutational load. We determined the prevalence of NNRTI resistance mutations at several RT aa positions in viruses from 204 antiretroviral (ARV)-naïve HIV-infected individuals using deep sequencing, and examined the relationship between mutant variant frequency and mutational load for those variants. Deep sequencing to ≥0.4% levels found variants with major NNRTI-resistance mutations having a Stanford-HIVdb algorithm value ≥30 for efavirenz and/or nevirapine in 52/204 (25.5%) ARV-naïve HIV-infected persons. Eighteen different major NNRTI mutations were identified at 11 different positions, with the majority of variants being at frequency >1%. The frequency of these variants correlated strongly with the mutational load, but this correlation weakened at low frequencies. Deep sequencing detected additional major NNRTI-resistant viral variants in treatment-naïve HIV-infected individuals. Our study suggests the significance of screening for mutations at all RT aa positions (in addition to K103 and Y181) to estimate the true burden of pre-treatment NNRTI-resistance. An important finding was that variants at low frequency had a wide range of mutational loads (>100-fold) suggesting that frequency alone may underestimate the impact of specific NNRTI-resistant variants. We recommend further evaluation of all low-frequency NNRTI-drug resistant variants with special attention given to the impact of mutational loads of these variants on treatment outcomes.

Human responses to influenza vaccination show seroconversion signatures and convergent antibody rearrangements.

B cells produce a diverse antibody repertoire by undergoing gene rearrangements. Pathogen exposure induces the clonal expansion of B cells expressing antibodies that can bind the infectious agent. To assess human B cell responses to trivalent seasonal influenza and monovalent pandemic H1N1 vaccination, we sequenced gene rearrangements encoding the immunoglobulin heavy chain, a major determinant of epitope recognition. The magnitude of B cell clonal expansions correlates with an individual's secreted antibody response to the vaccine, and the expanded clones are enriched with those expressing influenza-specific monoclonal antibodies. Additionally, B cell responses to pandemic influenza H1N1 vaccination and infection in different people show a prominent family of convergent antibody heavy chain gene rearrangements specific to influenza antigens. These results indicate that microbes can induce specific signatures of immunoglobulin gene rearrangements and that pathogen exposure can potentially be assessed from B cell repertoires.

An international multicenter study on HIV-1 drug resistance testing by 454 ultra-deep pyrosequencing.

The detection of mutant spectra within the viral quasispecies is critical for therapeutic management of HIV-1 infections. Routine clinical application of ultrasensitive genotyping requires reproducibility and concordance within and between laboratories. The goal of the study was to evaluate a new protocol on HIV-1 drug resistance testing by 454 ultra-deep pyrosequencing (454-UDS) in an international multicenter study. Sixteen blinded HIV-1 subtype B samples were provided for 454-UDS as both RNA and cDNA with viral titers of 88,600-573,000 HIV-1 RNA copies/ml. Eight overlapping amplicons spanning protease (PR) codons 10-99 and reverse transcriptase (RT) codons 1-251 were generated using molecular barcoded primers. 454-UDS was performed using the 454 Life Sciences/Roche GS FLX platform. PR and RT sequences were analyzed using 454 Life Sciences Amplicon Variant Analyzer (AVA) software. Quantified variation data were analyzed for intra-laboratory reproducibility and inter-laboratory concordance. Routine population sequencing was performed using the ViroSeq HIV-1 genotyping system. Eleven laboratories and the reference laboratory 454 Life Sciences sequenced the HIV-1 sample set. Data presented are derived from seven laboratories and the reference laboratory since severe study protocol execution errors occurred in four laboratories leading to exclusion. The median sequencing depth across all sites was 1364 reads per position (IQR=809-2065). 100% of the ViroSeq-reported mutations were also detected by 454-UDS. Minority HIV-1 drug resistance mutations, defined as HIV-1 drug resistance mutations identified at frequencies of 1-25%, were only detected by 454-UDS. Analysis of 10 preselected majority and minority mutations were consistently found across sites. The analysis of drug-resistance mutations detected between 1 and 10% demonstrated high intra- and inter-laboratory consistency in frequency estimates for both RNA and prepared cDNA samples, indicating robustness of the method. HIV-1 drug resistance testing using 454 ultra-deep pyrosequencing results in an accurate and highly reproducible, albeit complex, approach to the analysis of HIV-1 mutant spectra, even at frequencies well below those detected by routine population sequencing.

High-resolution antibody dynamics of vaccine-induced immune responses.

The adaptive immune system confers protection by generating a diverse repertoire of antibody receptors that are rapidly expanded and contracted in response to specific targets. Next-generation DNA sequencing now provides the opportunity to survey this complex and vast repertoire. In the present work, we describe a set of tools for the analysis of antibody repertoires and their application to elucidating the dynamics of the response to viral vaccination in human volunteers. By analyzing data from 38 separate blood samples across 2 y, we found that the use of the germ-line library of V and J segments is conserved between individuals over time. Surprisingly, there appeared to be no correlation between the use level of a particular VJ combination and degree of expansion. We found the antibody RNA repertoire in each volunteer to be highly dynamic, with each individual displaying qualitatively different response dynamics. By using combinatorial phage display, we screened selected VH genes paired with their corresponding VL library for affinity against the vaccine antigens. Altogether, this work presents an additional set of tools for profiling the human antibody repertoire and demonstrates characterization of the fast repertoire dynamics through time in multiple individuals responding to an immune challenge.

The effects of somatic hypermutation on neutralization and binding in the PGT121 family of broadly neutralizing HIV antibodies.

Broadly neutralizing HIV antibodies (bnAbs) are typically highly somatically mutated, raising doubts as to whether they can be elicited by vaccination. We used 454 sequencing and designed a novel phylogenetic method to model lineage evolution of the bnAbs PGT121-134 and found a positive correlation between the level of somatic hypermutation (SHM) and the development of neutralization breadth and potency. Strikingly, putative intermediates were characterized that show approximately half the mutation level of PGT121-134 but were still capable of neutralizing roughly 40-80% of PGT121-134 sensitive viruses in a 74-virus panel at median titers between 15- and 3-fold higher than PGT121-134. Such antibodies with lower levels of SHM may be more amenable to elicitation through vaccination while still providing noteworthy coverage. Binding characterization indicated a preference of inferred intermediates for native Env binding over monomeric gp120, suggesting that the PGT121-134 lineage may have been selected for binding to native Env at some point during maturation. Analysis of glycan-dependent neutralization for inferred intermediates identified additional adjacent glycans that comprise the epitope and suggests changes in glycan dependency or recognition over the course of affinity maturation for this lineage. Finally, patterns of neutralization of inferred bnAb intermediates suggest hypotheses as to how SHM may lead to potent and broad HIV neutralization and provide important clues for immunogen design.

Convergent antibody signatures in human dengue.

Dengue is the most prevalent mosquito-borne viral disease in humans, and the lack of early prognostics, vaccines, and therapeutics contributes to immense disease burden. To identify patterns that could be used for sequence-based monitoring of the antibody response to dengue, we examined antibody heavy-chain gene rearrangements in longitudinal peripheral blood samples from 60 dengue patients. Comparing signatures between acute dengue, postrecovery, and healthy samples, we found increased expansion of B cell clones in acute dengue patients, with higher overall clonality in secondary infection. Additionally, we observed consistent antibody sequence features in acute dengue in the highly variable major antigen-binding determinant, complementarity-determining region 3 (CDR3), with specific CDR3 sequences highly enriched in acute samples compared to postrecovery, healthy, or non-dengue samples. Dengue thus provides a striking example of a human viral infection where convergent immune signatures can be identified in multiple individuals. Such signatures could facilitate surveillance of immunological memory in communities.

Next generation sequencing reveals the association of DRB3*02:02 with type 1 diabetes.

The primary associations of the HLA class II genes, HLA-DRB1 and HLA-DQB1, and the class I genes, HLA-A and HLA-B, with type 1 diabetes (T1D) are well established. However, the role of polymorphism at the HLA-DRB3, HLA-DRB4, and HLA-DRB5 loci remains unclear. In two separate studies, one of 500 subjects and 500 control subjects and one of 366 DRB1*03:01-positive samples from selected multiplex T1D families, we used Roche 454 sequencing with Conexio Genomics ASSIGN ATF 454 HLA genotyping software analysis to analyze sequence variation at these three HLA-DRB loci. Association analyses were performed on the two HLA-DRB loci haplotypes (DRB1-DRB3, -DRB4, or -DRB5). Three common HLA-DRB3 alleles (*01:01, *02:02, *03:01) were observed. DRB1*03:01 haplotypes carrying DRB3*02:02 conferred a higher T1D risk than did DRB1*03:01 haplotypes carrying DRB3*01:01 in DRB1*03:01/*03:01 homozygotes with two DRB3*01:01 alleles (odds ratio [OR] 3.4 [95% CI 1.46-8.09]), compared with those carrying one or two DRB3*02:02 alleles (OR 25.5 [3.43-189.2]) (P = 0.033). For DRB1*03:01/*04:01 heterozygotes, however, the HLA-DRB3 allele did not significantly modify the T1D risk of the DRB1*03:01 haplotype (OR 7.7 for *02:02; 6.8 for *01:01). These observations were confirmed by sequence analysis of HLA-DRB3 exon 2 in a targeted replication study of 281 informative T1D family members and 86 affected family-based association control (AFBAC) haplotypes. The frequency of DRB3*02:02 was 42.9% in the DRB1*03:01/*03:01 patients and 27.6% in the DRB1*03:01/*04 (P = 0.005) compared with 22.6% in AFBAC DRB1*03:01 chromosomes (P = 0.001). Analysis of T1D-associated alleles at other HLA loci (HLA-A, HLA-B, and HLA-DPB1) on DRB1*03:01 haplotypes suggests that DRB3*02:02 on the DRB1*03:01 haplotype can contribute to T1D risk.

Therapy failure resulting from superinfection by a drug-resistant HIV variant.

BACKGROUND: HIV-1-infected patients can be superinfected with additional HIV-1 variants. Therapy failure can be the consequence of an infection with a resistant strain. METHODS: A patient was diagnosed with a recent HIV-1 infection in April 2005 and subsequently clinically monitored. HIV-1 evolution was studied by population sequencing of the first 984 bases of the pol gene as well as 454 ultra-deep pyrosequencing (UDPS) of parts of the pol and env genes. RESULTS: The patient was diagnosed with a wild-type HIV-1 strain, but experienced rapid virological failure after initiating a non-nucleoside reverse transcriptase inhibitor (NNRTI)-based treatment regimen 3 years later. Population sequencing and UDPS revealed the presence of a second HIV-1 strain with a Y188L NNRTI resistance mutation in a sample obtained shortly prior to initiation of therapy. Phylogenetic analyses showed that the two HIV-1 strains were genetically distinct, providing evidence for superinfection. CONCLUSIONS: The virological treatment failure in this patient was probably due to the superinfection with an NNRTI-resistant HIV-1 variant. Superinfection with drug-resistant strains can undermine HIV-1 treatment regimens selected on the basis of resistance testing at diagnosis. Patients, especially in high-risk groups, as well as their clinicians, should be aware of the risks and dangers of superinfections.

Initial antibodies binding to HIV-1 gp41 in acutely infected subjects are polyreactive and highly mutated.

The initial antibody response to HIV-1 is targeted to envelope (Env) gp41, and is nonneutralizing and ineffective in controlling viremia. To understand the origins and characteristics of gp41-binding antibodies produced shortly after HIV-1 transmission, we isolated and studied gp41-reactive plasma cells from subjects acutely infected with HIV-1. The frequencies of somatic mutations were relatively high in these gp41-reactive antibodies. Reverted unmutated ancestors of gp41-reactive antibodies derived from subjects acutely infected with HIV-1 frequently did not react with autologous HIV-1 Env; however, these antibodies were polyreactive and frequently bound to host or bacterial antigens. In one large clonal lineage of gp41-reactive antibodies, reactivity to HIV-1 Env was acquired only after somatic mutations. Polyreactive gp41-binding antibodies were also isolated from uninfected individuals. These data suggest that the majority of gp41-binding antibodies produced after acute HIV-1 infection are cross-reactive responses generated by stimulating memory B cells that have previously been activated by non-HIV-1 antigens.

High-throughput sequencing of PCR products tagged with universal primers using 454 life sciences systems.

This unit describes a method to convert PCR products (amplicons) flanked by universal M13 primers into a library for use on all 454 Sequencing Systems (454 Life Sciences, a Roche Company). This is especially useful for simultaneous sequencing and analysis of large numbers of amplicons or for the detection of minor variations within the amplified products. The method described here involves preparing a library of DNA with specific primers containing adaptor sequences recognized by the GS Junior System sequencing process. The data from the sequencing run are processed and analyzed by 454 Life Sciences software. This approach allows for multiplexing of a large number of amplicons to streamline processing and analysis. Any pre-existing universally tagged amplicon, primer set, or plasmid with M13 sequences flanking the cloning site can be used in this protocol.

Prevalence of low-level HIV-1 variants with reverse transcriptase mutation K65R and the effect of antiretroviral drug exposure on variant levels.

BACKGROUND: It has been reported that treatment-naive individuals infected with HIV-1 subtype C may be more likely to harbour viral variants possessing a K65R reverse transcriptase gene mutation. The objectives of this study were to determine the prevalence of low-level K65R variants within different HIV-1 subtypes and to assess the effects of antiretroviral exposure on K65R variant levels. METHODS: Treatment-naive individuals infected with different HIV-1 subtypes were genotyped by ultra-deep sequencing. Samples were evaluated for low-level variants to 0.4% or 1% levels depending upon viral load. Estimated mutational load was calculated by multiplying the percentage of the variant by the plasma viral load. RESULTS: A total of 411 treatment-naive individuals were evaluated by ultra-deep sequencing to 1% levels; 4 subjects (0.97%) had K65R variants at ≥1% or had a very high mutation load. All four subjects had variants with linked drug resistance mutations suggesting transmitted resistant variants. 147 ARV-naive subjects were sequenced to 0.4% levels; 8.8% (13/147) had K65R low-level variants identified: 2.2% (2/92) in subtype B, 35.7% (10/28) in subtype C (P<0.001 for B versus C) and 3.7% (1/27) in non-B/C subtypes. The 13 ARV-naive subjects with K65R variants at <1% received tenofovir plus emtricitabine plus a ritonavir-boosted protease inhibitor (TDF+FTC+PI/r) and 5 subsequently experienced virological failure. There was no enhancement in K65R levels by percentage or mutational load compared to pre-therapy levels. CONCLUSIONS: Low-level K65R variants were more frequently identified in subtype C. K65R variants at >1% levels likely represent transmitted resistant variants. The clinical implication of low-level K65R variants below 1% in treatment-naive subjects who receive TDF+FTC+PI/r remains to be determined as the majority are very low-level and did not increase after antiretroviral exposure.