Genetic variation across and within individuals.

Germline variation and somatic mutation are intricately connected and together shape human traits and disease risks. Germline variants are present from conception, but they vary between individuals and accumulate over generations. By contrast, somatic mutations accumulate throughout life in a mosaic manner within an individual due to intrinsic and extrinsic sources of mutations and selection pressures acting on cells. Recent advancements, such as improved detection methods and increased resources for association studies, have drastically expanded our ability to investigate germline and somatic genetic variation and compare underlying mutational processes. A better understanding of the similarities and differences in the types, rates and patterns of germline and somatic variants, as well as their interplay, will help elucidate the mechanisms underlying their distinct yet interlinked roles in human health and biology.

Clinical Sequencing Uncovers Origins and Evolution of Lassa Virus.

The 2013-2015 West African epidemic of Ebola virus disease (EVD) reminds us of how little is known about biosafety level 4 viruses. Like Ebola virus, Lassa virus (LASV) can cause hemorrhagic fever with high case fatality rates. We generated a genomic catalog of almost 200 LASV sequences from clinical and rodent reservoir samples. We show that whereas the 2013-2015 EVD epidemic is fueled by human-to-human transmissions, LASV infections mainly result from reservoir-to-human infections. We elucidated the spread of LASV across West Africa and show that this migration was accompanied by changes in LASV genome abundance, fatality rates, codon adaptation, and translational efficiency. By investigating intrahost evolution, we found that mutations accumulate in epitopes of viral surface proteins, suggesting selection for immune escape. This catalog will serve as a foundation for the development of vaccines and diagnostics. VIDEO ABSTRACT.

Managing differential performance of polygenic risk scores across groups: Real-world experience of the eMERGE Network.

The differential performance of polygenic risk scores (PRSs) by group is one of the major ethical barriers to their clinical use. It is also one of the main practical challenges for any implementation effort. The social repercussions of how people are grouped in PRS research must be considered in communications with research participants, including return of results. Here, we outline the decisions faced and choices made by a large multi-site clinical implementation study returning PRSs to diverse participants in handling this issue of differential performance. Our approach to managing the complexities associated with the differential performance of PRSs serves as a case study that can help future implementers of PRSs to plot an anticipatory course in response to this issue.

Neoantigen vaccine generates intratumoral T cell responses in phase Ib glioblastoma trial.

Neoantigens, which are derived from tumour-specific protein-coding mutations, are exempt from central tolerance, can generate robust immune responses1,2 and can function as bona fide antigens that facilitate tumour rejection3. Here we demonstrate that a strategy that uses multi-epitope, personalized neoantigen vaccination, which has previously been tested in patients with high-risk melanoma4-6, is feasible for tumours such as glioblastoma, which typically have a relatively low mutation load1,7 and an immunologically 'cold' tumour microenvironment8. We used personalized neoantigen-targeting vaccines to immunize patients newly diagnosed with glioblastoma following surgical resection and conventional radiotherapy in a phase I/Ib study. Patients who did not receive dexamethasone-a highly potent corticosteroid that is frequently prescribed to treat cerebral oedema in patients with glioblastoma-generated circulating polyfunctional neoantigen-specific CD4+ and CD8+ T cell responses that were enriched in a memory phenotype and showed an increase in the number of tumour-infiltrating T cells. Using single-cell T cell receptor analysis, we provide evidence that neoantigen-specific T cells from the peripheral blood can migrate into an intracranial glioblastoma tumour. Neoantigen-targeting vaccines thus have the potential to favourably alter the immune milieu of glioblastoma.

Corrigendum: An immunogenic personal neoantigen vaccine for patients with melanoma.

This corrects the article DOI: 10.1038/nature22991.

Returning integrated genomic risk and clinical recommendations: The eMERGE study.

PURPOSE: Assessing the risk of common, complex diseases requires consideration of clinical risk factors as well as monogenic and polygenic risks, which in turn may be reflected in family history. Returning risks to individuals and providers may influence preventive care or use of prophylactic therapies for those individuals at high genetic risk. METHODS: To enable integrated genetic risk assessment, the eMERGE (electronic MEdical Records and GEnomics) network is enrolling 25,000 diverse individuals in a prospective cohort study across 10 sites. The network developed methods to return cross-ancestry polygenic risk scores, monogenic risks, family history, and clinical risk assessments via a genome-informed risk assessment (GIRA) report and will assess uptake of care recommendations after return of results. RESULTS: GIRAs include summary care recommendations for 11 conditions, education pages, and clinical laboratory reports. The return of high-risk GIRA to individuals and providers includes guidelines for care and lifestyle recommendations. Assembling the GIRA required infrastructure and workflows for ingesting and presenting content from multiple sources. Recruitment began in February 2022. CONCLUSION: Return of a novel report for communicating monogenic, polygenic, and family history-based risk factors will inform the benefits of integrated genetic risk assessment for routine health care.

An immunogenic personal neoantigen vaccine for patients with melanoma.

Effective anti-tumour immunity in humans has been associated with the presence of T cells directed at cancer neoantigens, a class of HLA-bound peptides that arise from tumour-specific mutations. They are highly immunogenic because they are not present in normal tissues and hence bypass central thymic tolerance. Although neoantigens were long-envisioned as optimal targets for an anti-tumour immune response, their systematic discovery and evaluation only became feasible with the recent availability of massively parallel sequencing for detection of all coding mutations within tumours, and of machine learning approaches to reliably predict those mutated peptides with high-affinity binding of autologous human leukocyte antigen (HLA) molecules. We hypothesized that vaccination with neoantigens can both expand pre-existing neoantigen-specific T-cell populations and induce a broader repertoire of new T-cell specificities in cancer patients, tipping the intra-tumoural balance in favour of enhanced tumour control. Here we demonstrate the feasibility, safety, and immunogenicity of a vaccine that targets up to 20 predicted personal tumour neoantigens. Vaccine-induced polyfunctional CD4+ and CD8+ T cells targeted 58 (60%) and 15 (16%) of the 97 unique neoantigens used across patients, respectively. These T cells discriminated mutated from wild-type antigens, and in some cases directly recognized autologous tumour. Of six vaccinated patients, four had no recurrence at 25 months after vaccination, while two with recurrent disease were subsequently treated with anti-PD-1 (anti-programmed cell death-1) therapy and experienced complete tumour regression, with expansion of the repertoire of neoantigen-specific T cells. These data provide a strong rationale for further development of this approach, alone and in combination with checkpoint blockade or other immunotherapies.

Durvalumab plus tremelimumab alone or in combination with low-dose or hypofractionated radiotherapy in metastatic non-small-cell lung cancer refractory to previous PD(L)-1 therapy: an open-label, multicentre, randomised, phase 2 trial.

BACKGROUND: Patients with non-small-cell lung cancer (NSCLC) that is resistant to PD-1 and PD-L1 (PD[L]-1)-targeted therapy have poor outcomes. Studies suggest that radiotherapy could enhance antitumour immunity. Therefore, we investigated the potential benefit of PD-L1 (durvalumab) and CTLA-4 (tremelimumab) inhibition alone or combined with radiotherapy. METHODS: This open-label, multicentre, randomised, phase 2 trial was done by the National Cancer Institute Experimental Therapeutics Clinical Trials Network at 18 US sites. Patients aged 18 years or older with metastatic NSCLC, an Eastern Cooperative Oncology Group performance status of 0 or 1, and progression during previous PD(L)-1 therapy were eligible. They were randomly assigned (1:1:1) in a web-based system by the study statistician using a permuted block scheme (block sizes of three or six) without stratification to receive either durvalumab (1500 mg intravenously every 4 weeks for a maximum of 13 cycles) plus tremelimumab (75 mg intravenously every 4 weeks for a maximum of four cycles) alone or with low-dose (0·5 Gy delivered twice per day, repeated for 2 days during each of the first four cycles of therapy) or hypofractionated radiotherapy (24 Gy total delivered over three 8-Gy fractions during the first cycle only), 1 week after initial durvalumab-tremelimumab administration. Study treatment was continued until 1 year or until progression. The primary endpoint was overall response rate (best locally assessed confirmed response of a partial or complete response) and, along with safety, was analysed in patients who received at least one dose of study therapy. The trial is registered with ClinicalTrials.gov, NCT02888743, and is now complete. FINDINGS: Between Aug 24, 2017, and March 29, 2019, 90 patients were enrolled and randomly assigned, of whom 78 (26 per group) were treated. This trial was stopped due to futility assessed in an interim analysis. At a median follow-up of 12·4 months (IQR 7·8-15·1), there were no differences in overall response rates between the durvalumab-tremelimumab alone group (three [11·5%, 90% CI 1·2-21·8] of 26 patients) and the low-dose radiotherapy group (two [7·7%, 0·0-16·3] of 26 patients; p=0·64) or the hypofractionated radiotherapy group (three [11·5%, 1·2-21·8] of 26 patients; p=0·99). The most common grade 3-4 adverse events were dyspnoea (two [8%] in the durvalumab-tremelimumab alone group; three [12%] in the low-dose radiotherapy group; and three [12%] in the hypofractionated radiotherapy group) and hyponatraemia (one [4%] in the durvalumab-tremelimumab alone group vs two [8%] in the low-dose radiotherapy group vs three [12%] in the hypofractionated radiotherapy group). Treatment-related serious adverse events occurred in one (4%) patient in the durvalumab-tremelimumab alone group (maculopapular rash), five (19%) patients in the low-dose radiotherapy group (abdominal pain, diarrhoea, dyspnoea, hypokalemia, and respiratory failure), and four (15%) patients in the hypofractionated group (adrenal insufficiency, colitis, diarrhoea, and hyponatremia). In the low-dose radiotherapy group, there was one death from respiratory failure potentially related to study therapy. INTERPRETATION: Radiotherapy did not increase responses to combined PD-L1 plus CTLA-4 inhibition in patients with NSCLC resistant to PD(L)-1 therapy. However, PD-L1 plus CTLA-4 therapy could be a treatment option for some patients. Future studies should refine predictive biomarkers in this setting. FUNDING: The US National Institutes of Health and the Dana-Farber Cancer Institute.

Cross-Sectional Assessment of SARS-CoV-2 Viral Load by Symptom Status in Massachusetts Congregate Living Facilities.

Transmission of coronavirus disease 2019 (COVID-19) from people without symptoms confounds societal mitigation strategies. From April to June 2020, we tested nasopharyngeal swabs by reverse transcriptase quantitative polymerase chain reaction (RT-qPCR) from 15 514 staff and 16 966 residents of nursing homes and assisted living facilities in Massachusetts. Cycle threshold (Ct) distributions were very similar between populations with (n = 739) and without (n = 2179) symptoms at the time of sampling (mean Ct, 25.7 vs 26.4; ranges 12-38). However, as local cases waned, those without symptoms shifted towards higher Ct. With such similar viral load distributions, existing testing modalities should perform comparably regardless of symptoms, contingent upon time since infection.

Implementation of SARS-CoV2 Screening in K-12 Schools Using In-School Pooled Molecular Testing and Deconvolution by Rapid Antigen Test.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV2) testing is one component of a multilayered mitigation strategy to enable safe in-person school attendance for the K-12 school population. However, costs, logistics, and uncertainty about effectiveness are potential barriers to implementation. We assessed early data from the Massachusetts K-12 public school pooled SARS-CoV2 testing program, which incorporates two novel design elements: in-school "pod pooling" for assembling pools of dry anterior nasal swabs from 5 to 10 individuals and positive pool deconvolution using the BinaxNOW antigen rapid diagnostic test (Ag RDT), to assess the operational and analytical feasibility of this approach. Over 3 months, 187,597 individual swabs were tested across 39,297 pools from 738 schools. The pool positivity rate was 0.8%; 98.2% of pools tested negative and 0.2% inconclusive, and 0.8% of pools submitted could not be tested. Of 310 positive pools, 70.6% had an N1 or N2 probe cycle threshold (CT) value of ≤30. In reflex testing (performed on specimens newly collected from members of the positive pool), 92.5% of fully deconvoluted pools with an N1 or N2 target CT of ≤30 identified a positive individual using the BinaxNOW test performed 1 to 3 days later. However, of 124 positive pools with full reflex testing data available for analysis, 32 (25.8%) of BinaxNOW pool deconvolution testing attempts did not identify a positive individual, requiring additional reflex testing. With sufficient staffing support and low pool positivity rates, pooled sample collection and reflex testing were feasible for schools. These early program findings confirm that screening for K-12 students and staff is achievable at scale with a scheme that incorporates in-school pooling, primary testing by reverse transcription-PCR (RT-PCR), and Ag RDT reflex/deconvolution testing.

Performance and Implementation Evaluation of the Abbott BinaxNOW Rapid Antigen Test in a High-Throughput Drive-Through Community Testing Site in Massachusetts.

Rapid diagnostic tests (RDTs) for SARS-CoV-2 antigens (Ag) that can be performed at point of care (POC) can supplement molecular testing and help mitigate the COVID-19 pandemic. Deployment of an Ag RDT requires an understanding of its operational and performance characteristics under real-world conditions and in relevant subpopulations. We evaluated the Abbott BinaxNOW COVID-19 Ag card in a high-throughput, drive-through, free community testing site in Massachusetts using anterior nasal (AN) swab reverse transcriptase PCR (RT-PCR) for clinical testing. Individuals presenting for molecular testing in two of seven lanes were offered the opportunity to also receive BinaxNOW testing. Dual AN swabs were collected from symptomatic and asymptomatic children (≤18 years of age) and adults. BinaxNOW testing was performed in a testing pod with temperature/humidity monitoring. One individual performed testing and official result reporting for each test, but most tests had a second independent reading to assess interoperator agreement. Positive BinaxNOW results were scored as faint, medium, or strong. Positive BinaxNOW results were reported to patients by phone, and they were instructed to isolate pending RT-PCR results. The paired RT-PCR result was the reference for sensitivity and specificity calculations. Of 2,482 participants, 1,380 adults and 928 children had paired RT-PCR/BinaxNOW results and complete symptom data. In this study, 974/1,380 (71%) adults and 829/928 (89%) children were asymptomatic. BinaxNOW had 96.5% (95% confidence interval [CI], 90.0 to 99.3) sensitivity and 100% (95% CI, 98.6 to 100.0) specificity in adults within 7 days of symptoms and 84.6% (95% CI, 65.1 to 95.6) sensitivity and 100% (95% CI, 94.5 to 100.0) specificity in children within 7 days of symptoms. Sensitivity and specificity in asymptomatic adults were 70.2% (95% CI, 56.6 to 81.6) and 99.6% (95% CI, 98.9 to 99.9), respectively, and in asymptomatic children, they were 65.4% (95% CI, 55.6 to 74.4) and 99.0% (95% CI, 98.0 to 99.6), respectively. By cycle threshold (CT ) value cutoff, sensitivity in all subgroups combined (n = 292 RT-PCR-positive individuals) was 99.3% with CT values of ≤25, 95.8% with CT values of ≤30, and 81.2% with CT values of ≤35. Twelve false-positive BinaxNOW results (out of 2,308 tests) were observed; in all 12, the test bands were faint but otherwise normal and were noted by both readers. One invalid BinaxNOW result was identified. Interoperator agreement (positive versus negative BinaxNOW result) was 100% (n = 2,230/2,230 double reads). Each operator was able to process 20 RDTs per hour. In a separate set of 30 specimens (from individuals with symptoms ≤7 days) run at temperatures below the manufacturer's recommended range (46 to 58.5°F), sensitivity was 66.7% and specificity 95.2%. BinaxNOW had very high specificity in both adults and children and very high sensitivity in newly symptomatic adults. Overall, 95.8% sensitivity was observed with CT values of ≤30. These data support public health recommendations for use of the BinaxNOW test in adults with symptoms for ≤7 days without RT-PCR confirmation. Excellent interoperator agreement indicates that an individual can perform and read the BinaxNOW test alone. A skilled laboratorian can perform and read 20 tests per hour. Careful attention to temperature is critical.

DeTiN: overcoming tumor-in-normal contamination.

Comparison of sequencing data from a tumor sample with data from a matched germline control is a key step for accurate detection of somatic mutations. Detection sensitivity for somatic variants is greatly reduced when the matched normal sample is contaminated with tumor cells. To overcome this limitation, we developed deTiN, a method that estimates the tumor-in-normal (TiN) contamination level and, in cases affected by contamination, improves sensitivity by reclassifying initially discarded variants as somatic.

Plasma Mucin-1 (CA15-3) Levels in Autosomal Dominant Tubulointerstitial Kidney Disease due to MUC1 Mutations.

INTRODUCTION: Patients with ADTKD-MUC1 have one allele producing normal mucin-1 (MUC1) and one allele producing mutant MUC1, which remains intracellular. We hypothesized that ADTKD-MUC1 patients, who have only 1 secretory-competent wild-type MUC1 allele, should exhibit decreased plasma mucin-1 (MUC1) levels. To test this hypothesis, we repurposed the serum CA15-3 assay used to measure MUC1 in breast cancer to measure plasma MUC1 levels in ADTKD-MUC1. METHODS: This cross-sectional study analyzed CA15-3 levels in a reference population of 6,850 individuals, in 85 individuals with ADTKD-MUC1, and in a control population including 135 individuals with ADTKD-UMOD and 114 healthy individuals. RESULTS: Plasma CA15-3 levels (mean ± standard deviation) were 8.6 ± 4.3 U/mL in individuals with ADTKD-MUC1 and 14.6 ± 5.6 U/mL in controls (p < 0.001). While there was a significant difference in mean CA15-3 levels, there was substantial overlap between the 2 groups. Plasma CA15-3 levels were <5 U/mL in 22% of ADTKD-MUC1 patients, in 0/249 controls, and in 1% of the reference population. Plasma CA15-3 levels were >20 U/mL in 1/85 ADTKD-MUC1 patients, in 18% of control individuals, and in 25% of the reference population. Segregation of plasma CA15-3 levels by the rs4072037 genotype did not significantly improve differentiation between affected and unaffected individuals. CA15-3 levels were minimally affected by gender and estimated glomerular filtration rate. DISCUSSION/CONCLUSIONS: Plasma CA15-3 levels in ADTKD-MUC1 patients are approximately 40% lower than levels in healthy individuals, though there is significant overlap between groups. Further investigations need to be performed to see if plasma CA15-3 levels would be useful in diagnosis, prognosis, or assessing response to new therapies in this disorder.

Characterization and remediation of sample index swaps by non-redundant dual indexing on massively parallel sequencing platforms.

BACKGROUND: Here we present an in-depth characterization of the mechanism of sequencer-induced sample contamination due to the phenomenon of index swapping that impacts Illumina sequencers employing patterned flow cells with Exclusion Amplification (ExAmp) chemistry (HiSeqX, HiSeq4000, and NovaSeq). We also present a remediation method that minimizes the impact of such swaps. RESULTS: Leveraging data collected over a two-year period, we demonstrate the widespread prevalence of index swapping in patterned flow cell data. We calculate mean swap rates across multiple sample preparation methods and sequencer models, demonstrating that different library methods can have vastly different swapping rates and that even non-ExAmp chemistry instruments display trace levels of index swapping. We provide methods for eliminating sample data cross contamination by utilizing non-redundant dual indexing for complete filtering of index swapped reads, and share the sequences for 96 non-combinatorial dual indexes we have validated across various library preparation methods and sequencer models. Finally, using computational methods we provide a greater insight into the mechanism of index swapping. CONCLUSIONS: Index swapping in pooled libraries is a prevalent phenomenon that we observe at a rate of 0.2 to 6% in all sequencing runs on HiSeqX, HiSeq 4000/3000, and NovaSeq. Utilizing non-redundant dual indexing allows for the removal (flagging/filtering) of these swapped reads and eliminates swapping induced sample contamination, which is critical for sensitive applications such as RNA-seq, single cell, blood biopsy using circulating tumor DNA, or clinical sequencing.

Differences in the Selection Bottleneck between Modes of Sexual Transmission Influence the Genetic Composition of the HIV-1 Founder Virus.

Due to the stringent population bottleneck that occurs during sexual HIV-1 transmission, systemic infection is typically established by a limited number of founder viruses. Elucidation of the precise forces influencing the selection of founder viruses may reveal key vulnerabilities that could aid in the development of a vaccine or other clinical interventions. Here, we utilize deep sequencing data and apply a genetic distance-based method to investigate whether the mode of sexual transmission shapes the nascent founder viral genome. Analysis of 74 acute and early HIV-1 infected subjects revealed that 83% of men who have sex with men (MSM) exhibit a single founder virus, levels similar to those previously observed in heterosexual (HSX) transmission. In a metadata analysis of a total of 354 subjects, including HSX, MSM and injecting drug users (IDU), we also observed no significant differences in the frequency of single founder virus infections between HSX and MSM transmissions. However, comparison of HIV-1 envelope sequences revealed that HSX founder viruses exhibited a greater number of codon sites under positive selection, as well as stronger transmission indices possibly reflective of higher fitness variants. Moreover, specific genetic "signatures" within MSM and HSX founder viruses were identified, with single polymorphisms within gp41 enriched among HSX viruses while more complex patterns, including clustered polymorphisms surrounding the CD4 binding site, were enriched in MSM viruses. While our findings do not support an influence of the mode of sexual transmission on the number of founder viruses, they do demonstrate that there are marked differences in the selection bottleneck that can significantly shape their genetic composition. This study illustrates the complex dynamics of the transmission bottleneck and reveals that distinct genetic bottleneck processes exist dependent upon the mode of HIV-1 transmission.

Medulloblastoma exome sequencing uncovers subtype-specific somatic mutations.

Medulloblastomas are the most common malignant brain tumours in children. Identifying and understanding the genetic events that drive these tumours is critical for the development of more effective diagnostic, prognostic and therapeutic strategies. Recently, our group and others described distinct molecular subtypes of medulloblastoma on the basis of transcriptional and copy number profiles. Here we use whole-exome hybrid capture and deep sequencing to identify somatic mutations across the coding regions of 92 primary medulloblastoma/normal pairs. Overall, medulloblastomas have low mutation rates consistent with other paediatric tumours, with a median of 0.35 non-silent mutations per megabase. We identified twelve genes mutated at statistically significant frequencies, including previously known mutated genes in medulloblastoma such as CTNNB1, PTCH1, MLL2, SMARCA4 and TP53. Recurrent somatic mutations were newly identified in an RNA helicase gene, DDX3X, often concurrent with CTNNB1 mutations, and in the nuclear co-repressor (N-CoR) complex genes GPS2, BCOR and LDB1. We show that mutant DDX3X potentiates transactivation of a TCF promoter and enhances cell viability in combination with mutant, but not wild-type, ß-catenin. Together, our study reveals the alteration of WNT, hedgehog, histone methyltransferase and now N-CoR pathways across medulloblastomas and within specific subtypes of this disease, and nominates the RNA helicase DDX3X as a component of pathogenic ß-catenin signalling in medulloblastoma.

Accurate Estimation of Fungal Diversity and Abundance through Improved Lineage-Specific Primers Optimized for Illumina Amplicon Sequencing.

While high-throughput sequencing methods are revolutionizing fungal ecology, recovering accurate estimates of species richness and abundance has proven elusive. We sought to design internal transcribed spacer (ITS) primers and an Illumina protocol that would maximize coverage of the kingdom Fungi while minimizing nontarget eukaryotes. We inspected alignments of the 5.8S and large subunit (LSU) ribosomal genes and evaluated potential primers using PrimerProspector. We tested the resulting primers using tiered-abundance mock communities and five previously characterized soil samples. We recovered operational taxonomic units (OTUs) belonging to all 8 members in both mock communities, despite DNA abundances spanning 3 orders of magnitude. The expected and observed read counts were strongly correlated (r = 0.94 to 0.97). However, several taxa were consistently over- or underrepresented, likely due to variation in rRNA gene copy numbers. The Illumina data resulted in clustering of soil samples identical to that obtained with Sanger sequence clone library data using different primers. Furthermore, the two methods produced distance matrices with a Mantel correlation of 0.92. Nonfungal sequences comprised less than 0.5% of the soil data set, with most attributable to vascular plants. Our results suggest that high-throughput methods can produce fairly accurate estimates of fungal abundances in complex communities. Further improvements might be achieved through corrections for rRNA copy number and utilization of standardized mock communities. IMPORTANCE: Fungi play numerous important roles in the environment. Improvements in sequencing methods are providing revolutionary insights into fungal biodiversity, yet accurate estimates of the number of fungal species (i.e., richness) and their relative abundances in an environmental sample (e.g., soil, roots, water, etc.) remain difficult to obtain. We present improved methods for high-throughput Illumina sequencing of the species-diagnostic fungal ribosomal marker gene that improve the accuracy of richness and abundance estimates. The improvements include new PCR primers and library preparation, validation using a known mock community, and bioinformatic parameter tuning.

Fine-mapping classical HLA variation associated with durable host control of HIV-1 infection in African Americans.

A small proportion of human immunodeficiency virus-1 (HIV-1) infected individuals, termed HIV-1 controllers, suppress viral replication to very low levels in the absence of therapy. Genetic investigations of this phenotype have strongly implicated variation in the class I major histocompatibility complex (MHC) region as key to HIV-1 control. We collected sequence-based classical class I HLA genotypes at 4-digit resolution in HIV-1-infected African American controllers and progressors (n = 1107), and tested them for association with host control using genome-wide single nucleotide polymorphism data to account for population structure. Several classical alleles at HLA-B were associated with host control, including B*57:03 [odds ratio (OR) = 5.1; P= 3.4 × 10(-18)] and B*81:01 (OR = 4.8; P= 1.3 × 10(-9)). Analysis of variable amino acid positions demonstrates that HLA-B position 97 is the most significant association with host control in African Americans (omnibus P = 1.2 × 10(-21)) and explains the signal of several HLA-B alleles, including B*57:03. Within HLA-B, we also identified independent effects at position 116 (omnibus P= 2.8 × 10(-15)) in the canonical F pocket, position 63 in the B pocket (P= 1.5 × 10(-3)) and the non-pocket position 245 (P= 8.8 × 10(-10)), which is thought to influence CD8-binding kinetics. Adjusting for these HLA-B effects, there is evidence for residual association in the MHC region. These results underscore the key role of HLA-B in affecting HIV-1 replication, likely through the molecular interaction between HLA-B and viral peptides presented by infected cells, and suggest that sites outside the peptide-binding pocket also influence HIV-1 control.

Whole genome deep sequencing of HIV-1 reveals the impact of early minor variants upon immune recognition during acute infection.

Deep sequencing technologies have the potential to transform the study of highly variable viral pathogens by providing a rapid and cost-effective approach to sensitively characterize rapidly evolving viral quasispecies. Here, we report on a high-throughput whole HIV-1 genome deep sequencing platform that combines 454 pyrosequencing with novel assembly and variant detection algorithms. In one subject we combined these genetic data with detailed immunological analyses to comprehensively evaluate viral evolution and immune escape during the acute phase of HIV-1 infection. The majority of early, low frequency mutations represented viral adaptation to host CD8+ T cell responses, evidence of strong immune selection pressure occurring during the early decline from peak viremia. CD8+ T cell responses capable of recognizing these low frequency escape variants coincided with the selection and evolution of more effective secondary HLA-anchor escape mutations. Frequent, and in some cases rapid, reversion of transmitted mutations was also observed across the viral genome. When located within restricted CD8 epitopes these low frequency reverting mutations were sufficient to prime de novo responses to these epitopes, again illustrating the capacity of the immune response to recognize and respond to low frequency variants. More importantly, rapid viral escape from the most immunodominant CD8+ T cell responses coincided with plateauing of the initial viral load decline in this subject, suggestive of a potential link between maintenance of effective, dominant CD8 responses and the degree of early viremia reduction. We conclude that the early control of HIV-1 replication by immunodominant CD8+ T cell responses may be substantially influenced by rapid, low frequency viral adaptations not detected by conventional sequencing approaches, which warrants further investigation. These data support the critical need for vaccine-induced CD8+ T cell responses to target more highly constrained regions of the virus in order to ensure the maintenance of immunodominant CD8 responses and the sustained decline of early viremia.

Whole genome pyrosequencing of rare hepatitis C virus genotypes enhances subtype classification and identification of naturally occurring drug resistance variants.

BACKGROUND: Infection with hepatitis C virus (HCV) is a burgeoning worldwide public health problem, with 170 million infected individuals and an estimated 20 million deaths in the coming decades. While 6 main genotypes generally distinguish the global geographic diversity of HCV, a multitude of closely related subtypes within these genotypes are poorly defined and may influence clinical outcome and treatment options. Unfortunately, the paucity of genetic data from many of these subtypes makes time-consuming primer walking the limiting step for sequencing understudied subtypes. METHODS: Here we combined long-range polymerase chain reaction amplification with pyrosequencing for a rapid approach to generate the complete viral coding region of 31 samples representing poorly defined HCV subtypes. RESULTS: Phylogenetic classification based on full genome sequences validated previously identified HCV subtypes, identified a recombinant sequence, and identified a new distinct subtype of genotype 4. Unlike conventional sequencing methods, use of deep sequencing also facilitated characterization of minor drug resistance variants within these uncommon or, in some cases, previously uncharacterized HCV subtypes. CONCLUSIONS: These data aid in the classification of uncommon HCV subtypes while also providing a high-resolution view of viral diversity within infected patients, which may be relevant to the development of therapeutic regimens to minimize drug resistance.