Search

Introduction: Circulating tumor DNA (ctDNA) offers the ability to repeatedly interrogate tumor genomic information, providing an opportunity for real-time monitoring of tumor genomic dynamics. In this study, we deeply analyzed multiple ctDNA samples collected over narrow time frames (days-to-weeks) from seven patients with metastatic triple-negative breast cancer (mTNBC), a cancer type known to have high ctDNA content. Methods: Patients with mTNBC were enrolled in a clinical trial of multi-kinase inhibitor cabozantinib, providing uniform and targeted treatment, and samples were collected day 1, day 8, then every 21 to 42 days on study. ctDNA was extracted from each plasma sample and underwent ultra-low pass whole genome sequencing (ULP-WGS; average depth 0.1x; n=42 samples), deep targeted panel sequencing (TPS) of 402 cancer-related genes with unique molecular identifier indexing (depth 10,000x; n=42 samples), and samples with tumor fraction (TFx) >10% underwent whole exome sequencing (WES; depth 200x; n=31 samples), with whole blood germline sequencing of both TPS and WES for subsequent analyses. Somatic copy number alterations (SCNAs) were identified from ULP-WGS and WES. PyClone with TPS was employed for clonal dynamic analyses. Predicted neoantigens were determined from WES using HLAthena. Results: A total of 42 total plasma samples from 7 patients (range 4-8 samples per patient) were collected at narrow time intervals, median 21 days (range 6 to 42 days) between samples. The median TFx across all samples was 18.1% (range 2.5% to 44.3%). TFx estimates were concordant when comparing orthogonal sequencing approaches (ULP-WGS, WES) and tumor fraction estimation algorithms (ichorCNA, FACETS). Despite all seven patients having ‘stable disease’ as best objective response, TFx dynamics were widely variable with TFx declining to lower limit of detection in three of seven patients. Of all samples, 31/42 (73.8%) had tumor fraction >10% and underwent WES; each patient had at least 3 samples that underwent both WES and TPS. There was strong agreement between TPS and WES: across all 31 shared samples, mutation recall in TPS versus WES (gold standard) was 95.5%. Variant allele frequency across all mutations detected in both TPS and WES was highly concordant (Pearson’s r=0.949). Clonal mutations were consistently detected across multiple samples within patients. When comparing genome-wide copy number from last to first available sample within each patient, copy number log ratios were largely stable within patients (union Pearson’s r=0.924) and there were not recurrent shifts in SCNAs across patients. Through statistical modeling of TPS data, we tracked distinct clonal populations for each patient over their sampling windows. Modeled clonal architecture in most patients revealed stable, polyclonal profiles, with important breast cancer driver alterations (e.g. TP53 and PIK3CA) recurrently presenting at high prevalence. Infrequently, we also detected emergence and expansion of clones over narrow time frames (weeks) containing acquired alterations poorly annotated in the breast cancer literature. We successfully predicted neoantigens from ctDNA WES at multiple time points in each patient, with evidence that patients acquired new mutations predicted to be ‘strong binder’ neoantigens over time on therapy. Conclusions: Analysis of serial ctDNA samples collected at narrow time intervals (days-to-weeks) provides unique insight into the dynamics of ctDNA. We demonstrate strong concordance across ctDNA sequencing appraoches. Evolving genomic features of tumor populations can be identified via ctDNA while on treatment, potentially providing real time insight for clinical decision-making. Citation Format: Daniel G Stover, Katharine A Collier, David Tallman, Juliet Forman, Sachet Shukla, Sarah Asad, Justin Rhoades, Samuel Freeman, Mathew Cherian, Sagar Sardesai, Romualdo Barroso-Sousa, Carrie Cibulskis, Niall Lennon, Gavin Ha, Sara M Tolaney, Viktor A Adalsteinsson, Zachary Weber. Modeling clonal structure over narrow time frames via circulating tumor DNA in metastatic breast cancer [abstract]. In: Proceedings of the 2020 San Antonio Breast Cancer Virtual Symposium; 2020 Dec 8-11; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2021;81(4 Suppl):Abstract nr PD9-08.

BRAF mutations occur in approximately 10% of colorectal cancers. Although RAF inhibitor monotherapy is highly effective in BRAF-mutant melanoma, response rates in BRAF-mutant colorectal cancer are poor. Recent clinical trials of combined RAF/EGFR or RAF/MEK inhibition have produced improved efficacy, but patients ultimately develop resistance. To identify molecular alterations driving clinical acquired resistance, we performed whole-exome sequencing on paired pretreatment and postprogression tumor biopsies from patients with BRAF-mutant colorectal cancer treated with RAF inhibitor combinations. We identified alterations in MAPK pathway genes in resistant tumors not present in matched pretreatment tumors, including KRAS amplification, BRAF amplification, and a MEK1 mutation. These alterations conferred resistance to RAF/EGFR or RAF/MEK combinations through sustained MAPK pathway activity, but an ERK inhibitor could suppress MAPK activity and overcome resistance. Identification of MAPK pathway reactivating alterations upon clinical acquired resistance underscores the MAPK pathway as a critical target in BRAF-mutant colorectal cancer and suggests therapeutic options to overcome resistance. Significance: RAF inhibitor combinations represent promising approaches in clinical development for BRAF-mutant colorectal cancer. Initial characterization of clinical acquired resistance mechanisms to these regimens identified several MAPK pathway alterations driving resistance by reactivating MAPK signaling, highlighting the critical dependence of BRAF-mutant colorectal cancers on MAPK signaling and offering potential strategies to overcome resistance. Cancer Discov; 5(4); 358–67. ©2015 AACR. See related commentary by Meador and Pao, p. 348 This article is highlighted in the In This Issue feature, p. 333

Analogous to observations in RNA viruses such as human immunodeficiency virus, genetic variation associated with intrahost dengue virus (DENV) populations has been postulated to influence viral fitness and disease pathogenesis. Previous attempts to investigate intrahost genetic variation in DENV characterized only a few viral genes or a limited number of full-length genomes. We developed a whole-genome amplification approach coupled with deep sequencing to capture intrahost diversity across the entire coding region of DENV-2. Using this approach, we sequenced DENV-2 genomes from the serum of 22 Nicaraguan individuals with secondary DENV infection and captured ∼75% of the DENV genome in each sample (range, 40 to 98%). We identified and quantified variants using a highly sensitive and specific method and determined that the extent of diversity was considerably lower than previous estimates. Significant differences in intrahost diversity were detected between genes and also between antigenically distinct domains of the Envelope gene. Interestingly, a strong association was discerned between the extent of intrahost diversity in a few genes and viral clade identity. Additionally, the abundance of viral variants within a host, as well as the impact of viral mutations on amino acid encoding and predicted protein function, determined whether intrahost variants were observed at the interhost level in circulating Nicaraguan DENV-2 populations, strongly suggestive of purifying selection across transmission events. Our data illustrate the value of high-coverage genome-wide analysis of intrahost diversity for high-resolution mapping of the relationship between intrahost diversity and clinical, epidemiological, and virological parameters of viral infection.

A variety of microbial communities and their genes (the microbiome) exist throughout the human body, with fundamental roles in human health and disease. The National Institutes of Health (NIH)-funded Human Microbiome Project Consortium has established a population-scale framework to develop metagenomic protocols, resulting in a broad range of quality-controlled resources and data including standardized methods for creating, processing and interpreting distinct types of high-throughput metagenomic data available to the scientific community. Here we present resources from a population of 242 healthy adults sampled at 15 or 18 body sites up to three times, which have generated 5,177 microbial taxonomic profiles from 16S ribosomal RNA genes and over 3.5 terabases of metagenomic sequence so far. In parallel, approximately 800 reference strains isolated from the human body have been sequenced. Collectively, these data represent the largest resource describing the abundance and variety of the human microbiome, while providing a framework for current and future studies.

BACKGROUND: The RTS,S/AS01 vaccine targets the circumsporozoite protein of Plasmodium falciparum and has partial protective efficacy against clinical and severe malaria disease in infants and children. We investigated whether the vaccine efficacy was specific to certain parasite genotypes at the circumsporozoite protein locus. METHODS: We used polymerase chain reaction-based next-generation sequencing of DNA extracted from samples from 4985 participants to survey circumsporozoite protein polymorphisms. We evaluated the effect that polymorphic positions and haplotypic regions within the circumsporozoite protein had on vaccine efficacy against first episodes of clinical malaria within 1 year after vaccination. RESULTS: In the per-protocol group of 4577 RTS,S/AS01-vaccinated participants and 2335 control-vaccinated participants who were 5 to 17 months of age, the 1-year cumulative vaccine efficacy was 50.3% (95% confidence interval [CI], 34.6 to 62.3) against clinical malaria in which parasites matched the vaccine in the entire circumsporozoite protein C-terminal (139 infections), as compared with 33.4% (95% CI, 29.3 to 37.2) against mismatched malaria (1951 infections) (P=0.04 for differential vaccine efficacy). The vaccine efficacy based on the hazard ratio was 62.7% (95% CI, 51.6 to 71.3) against matched infections versus 54.2% (95% CI, 49.9 to 58.1) against mismatched infections (P=0.06). In the group of infants 6 to 12 weeks of age, there was no evidence of differential allele-specific vaccine efficacy. CONCLUSIONS: These results suggest that among children 5 to 17 months of age, the RTS,S vaccine has greater activity against malaria parasites with the matched circumsporozoite protein allele than against mismatched malaria. The overall vaccine efficacy in this age category will depend on the proportion of matched alleles in the local parasite population; in this trial, less than 10% of parasites had matched alleles. (Funded by the National Institutes of Health and others.).

West Nile virus (WNV) has become firmly established in northeastern US, reemerging every summer since its introduction into North America in 1999. To determine whether WNV overwinters locally or is reseeded annually, we examined the patterns of viral lineage persistence and replacement in Connecticut over 10 consecutive transmission seasons by phylogenetic analysis. In addition, we compared the full protein coding sequence among WNV isolates to search for evidence of convergent and adaptive evolution. Viruses sampled from Connecticut segregated into a number of well-supported subclades by year of isolation with few clades persisting ≥ 2 years. Similar viral strains were dispersed in different locations across the state and divergent strains appeared within a single location during a single transmission season, implying widespread movement and rapid colonization of virus. Numerous amino acid substitutions arose in the population but only one change, V → A at position 159 of the envelope protein, became permanently fixed. Several instances of parallel evolution were identified in independent lineages, including one amino acid change in the NS4A protein that appears to be positively selected. Our results suggest that annual reemergence of WNV is driven by both reintroduction and local-overwintering of virus. Despite ongoing evolution of WNV, most amino acid variants occurred at low frequencies and were transient in the virus population.

Resistance mutations to hepatitis C virus (HCV) nonstructural protein 3 (NS3) protease inhibitors in1% of the viral quasispecies may still allow1000-fold viral load reductions upon treatment, consistent with their reported reduced replicative fitness in vitro. Recently, however, an R155K protease mutation was reported as the dominant quasispecies in a treatment-naïve individual, raising concerns about possible full drug resistance. To investigate the prevalence of dominant resistance mutations against specifically targeted antiviral therapy for HCV (STAT-C) in the population, we analyzed HCV genome sequences from 507 treatment-naïve patients infected with HCV genotype 1 from the United States, Germany, and Switzerland. Phylogenetic sequence analysis and viral load data were used to identify the possible spread of replication-competent, drug-resistant viral strains in the population and to infer the consequences of these mutations upon viral replication in vivo. Mutations described to confer resistance to the protease inhibitors Telaprevir, BILN2061, ITMN-191, SCH6 and Boceprevir; the NS5B polymerase inhibitor AG-021541; and to the NS4A antagonist ACH-806 were observed mostly as sporadic, unrelated cases, at frequencies between 0.3% and 2.8% in the population, including two patients with possible multidrug resistance. Collectively, however, 8.6% of the patients infected with genotype 1a and 1.4% of those infected with genotype 1b carried at least one dominant resistance mutation. Viral loads were high in the majority of these patients, suggesting that drug-resistant viral strains might achieve replication levels comparable to nonresistant viruses in vivo.Naturally occurring dominant STAT-C resistance mutations are common in treatment-naïve patients infected with HCV genotype 1. Their influence on treatment outcome should further be characterized to evaluate possible benefits of drug resistance testing for individual tailoring of drug combinations when treatment options are limited due to previous nonresponse to peginterferon and ribavirin.

High throughput sequencing methods are widely used in analyses of microbial diversity, but are generally applied to small numbers of samples, which precludes characterization of patterns of microbial diversity across space and time. We have designed a primer-tagging approach that allows pooling and subsequent sorting of numerous samples, which is directed to amplification of a region spanning the nuclear ribosomal internal transcribed spacers and partial large subunit from fungi in environmental samples. To test the method for phylogenetic biases, we constructed a controlled mixture of four taxa representing the Chytridiomycota, Zygomycota, Ascomycota and Basidiomycota. Following cloning and colony restriction fragment length polymorphism analysis, we found no significant difference in representation in 19 of the 23 tested primers. We also generated a clone library from two soil DNA extracts using two primers for each extract and compared 456 clone sequences. Community diversity statistics and contingency table tests applied to counts of operational taxonomic units revealed that the two DNA extracts differed significantly, while the pairs of tagged primers from each extract were indistinguishable. Similar results were obtained using UniFrac phylogenetic comparisons. Together, these results suggest that the pig-tagged primers can be used to increase ecological inference in high throughput sequencing projects on fungi.

The rapid spread of dengue is a worldwide public health problem. In two clinical studies of dengue in Managua, Nicaragua, we observed an abrupt increase in disease severity across several epidemic seasons of dengue virus serotype 2 (DENV-2) transmission. Waning DENV-1 immunity appeared to increase the risk of severe disease in subsequent DENV-2 infections after a period of cross-protection. The increase in severity coincided with replacement of the Asian/American DENV-2 NI-1 clade with a new virus clade, NI-2B. In vitro analyses of viral isolates from the two clades and analysis of viremia in patient blood samples support the emergence of a fitter virus in later, relative to earlier, epidemic seasons. In addition, the NI-1 clade of viruses was more virulent specifically in children who were immune to DENV-1, whereas DENV-3 immunity was associated with more severe disease among NI-2B infections. Our data demonstrate that the complex interaction between viral genetics and population dynamics of serotype-specific immunity contributes to the risk of severe dengue disease. Furthermore, this work provides insights into viral evolution and the interaction between viral and immunological determinants of viral fitness and virulence.

Infection with genotype 4 of the Hepatitis C virus is common in Africa and the Mediterranean area, but has also been found at increasing frequencies in injection drug users in Europe and North America. Full length viral sequences to characterize viral diversity and structure have recently become available mostly for subtype 4a, and studies in Egypt and Saudi Arabia, where high proportions of subtype 4a infected patients exist, have begun to establish optimized treatment regimens. However knowledge about other subtype variants of genotype 4 present in less developed African states is lacking. In this study the full coding region from so far poorly characterized variants of HCV genotype 4 was amplified and sequenced using a long range PCR technique. Sequences were analyzed with respect to phylogenetic relationship, possible recombination and prominent sequence characteristics compared to other known HCV strains. We present for the first time two full-length sequences from the HCV genotype 4k, in addition to five strains from HCV genotypes 4d and 4f. Reference sequences for accurate HCV genotyping are required for optimized treatment, and a better knowledge of the global viral sequence diversity is needed to guide vaccines or new drugs effective in the world wide epidemic.

Mutations in the dysferlin gene cause limb girdle muscular dystrophy type 2B and Miyoshi myopathy. We report here the results of expression profile analyses and in vitro investigations that point to an interaction between dysferlin and the Ca2+ and lipid-binding proteins, annexins A1 and A2, and define a role for dysferlin in Ca2+-dependent repair of sarcolemmal injury through a process of vesicle fusion. Expression profiling identified a network of genes that are co-regulated in dysferlinopathic mice. Co-immunofluorescence, co-immunoprecipitation, and fluorescence lifetime imaging microscopy revealed that dysferlin normally associates with both annexins A1 and A2 in a Ca2+ and membrane injury-dependent manner. The distribution of the annexins and the efficiency of sarcolemmal wound-healing are significantly disrupted in dysferlin-deficient muscle. We propose a model of muscle membrane healing mediated by dysferlin that is relevant to both normal and dystrophic muscle and defines the annexins as potential muscular dystrophy genes.

Familial amyotrophic lateral sclerosis (ALS)-linked mutations in the copper-zinc superoxide dismutase (SOD1) gene cause motor neuron death in about 3% of ALS cases. While the wild-type (wt) protein is anti-apoptotic, mutant SOD1 promotes apoptosis. We now demonstrate that both wt and mutant SOD1 bind the anti-apoptotic protein Bcl-2, providing evidence of a direct link between SOD1 and an apoptotic pathway. This interaction is evident in vitro and in vivo in mouse and human spinal cord. We also demonstrate that in mice and humans, Bcl-2 binds to high molecular weight SDS-resistant mutant SOD1 containing aggregates that are present in mitochondria from spinal cord but not liver. These findings provide new insights into the anti-apoptotic function of SOD1 and suggest that entrapment of Bcl-2 by large SOD1 aggregates may deplete motor neurons of this anti-apoptotic protein.