McDonough, Samantha J.

ABSTRACT Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has had an unprecedented impact on human health and highlights the need for genomic epidemiology studies to increase our understanding of virus evolution and spread and to inform policy decisions. We sequenced viral genomes from over 22,000 patient samples tested at Mayo Clinic Laboratories between 2020 and 2022 and used Bayesian phylodynamics to describe county and regional spread in Minnesota. The earliest calculated introduction into Minnesota was to Hennepin County from a domestic source around 22 January 2020; 6 weeks before the first confirmed case in the state. This led to the virus spreading to Northern Minnesota and, eventually, the rest of the state. International introductions were most abundant in Hennepin (home to the Minneapolis/St. Paul International Airport) totaling 45 (out of 107) over the 2-year period. Southern Minnesota counties were most common for domestic introductions, with 19 (out of 64), potentially driven by bordering states such as Iowa and Wisconsin as well as Illinois, which is nearby. Hennepin also was, by far, the most dominant source of in-state transmissions to other Minnesota locations ( n = 772) over the 2-year period. We also analyzed the diversity of the location source of SARS-CoV-2 viruses in each county and noted the timing of state-wide policies as well as trends in clinical cases. Neither the number of clinical cases nor the major policy decisions, such as the end of the lockdown period in 2020 or the end of all restrictions in 2021, appeared to have an impact on virus diversity across each individual county. IMPORTANCE We analyzed over 22,000 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes of patient samples tested at Mayo Clinic Laboratories during a 2-year period in the COVID-19 pandemic, which included Alpha, Delta, and Omicron variants of concern to examine the roles and relationships of Minnesota virus transmission. We found that Hennepin County, the most populous county, drove the transmission of SARS-CoV-2 viruses in the state after including the formation of earlier clades including 20A, 20C, and 20G, as well as variants of concern Alpha and Delta. We also found that Hennepin County was the source for most of the county-to-county introductions after an initial predicted introduction with the virus in early 2020 from an international source, while other counties acted as transmission “sinks.” In addition, major policies, such as the end of the lockdown period in 2020 or the end of all restrictions in 2021, did not appear to have an impact on virus diversity across individual counties.

Abstract BackgroundFormalin-fixed, paraffin-embedded (FFPE) tissues have many advantages for identification of risk biomarkers, including wide availability and potential for extended follow-up endpoints. However, RNA derived from archival FFPE samples has limited quality. Here we identified parameters that determine which FFPE samples have the potential for successful RNA extraction, library preparation, and generation of usable RNAseq data.MethodsWe optimized library preparation protocols designed for use with FFPE samples using seven FFPE and Fresh Frozen replicate pairs, and tested optimized protocols using a study set of 130 FFPE biopsies from women with benign breast disease. Metrics from RNA extraction and preparation procedures were collected and compared with bioinformatics sequencing summary statistics. Finally, a decision tree model was built to learn the relationship between pre-sequencing lab metrics and qc pass/fail status as determined by bioinformatics metrics.ResultsSamples that failed bioinformatics qc tended to have low median sample-wise correlation within the cohort (Spearman correlation < 0.75), low number of reads mapped to gene regions (< 25 million), or low number of detectable genes (11,400 # of detected genes with CPM > 2). The median RNA concentration and pre-capture library Qubit values for qc failed samples were 18.7 ng/ul and 1.68 ng/ul respectively, which were significantly lower than those of qc pass samples (40.8 ng/ul and 5.58 ng/ul). We built a decision tree model based on input RNA concentration, input library qubit values, and achieved an F score of 0.803 in predicting QC status (pass/fail) of FFPE samples. ConclusionsWe provide a bioinformatics quality control recommendation for FFPE samples from breast tissue by evaluating bioinformatic and sample metrics. Our results suggest a minimum concentration of 25 ng/ul FFPE-extracted RNA for library preparation and 1.8 ng/ul pre-capture library output to achieve adequate RNA-seq data for downstream bioinformatics analysis.

Abstract Benign breast disease (BBD) is an established risk factor for development of breast cancer (BC). However, little is known about the prevalence of somatic genetic alterations at the time of BBD diagnosis and relationship to subsequent BC. To characterize DNA variants in BBD and investigate association with BC risk, we isolated DNA from formalin-fixed paraffin-embedded (FFPE) BBD biopsies and assessed mutations using a panel of 93 breast cancer predisposition genes (Qiagen GeneRead). The study set consisted of 120 BBD patients from three age-matched groups defined based on 16 years of BC follow-up: BBD-controls (cancer-free after 16 years), BBD-ER+, or BBD-ER- breast cancer (within 16 years). Bioinformatics quality checks included comparisons of paired fresh-frozen and FFPE replicates to reduce false variant calls induced by FFPE artifacts. Gene-level variant burden differences among study groups were assessed using 12 statistical approaches, adjusted for patient age, year of biopsy, and histology. To investigate possible immune implications of increased variant burden, extent of leukocytes in the BBD biopsies was assessed using CD45 immunohistochemistry. To assess variant profile changes during the progression from BBD to BC, we also performed targeted sequencing of 12 BBD-matched subsequent tumor samples (six ER+ and six ER-) developed within 5 years of BBD diagnosis. Gene-level association results showed consistently skewed distributions suggesting a higher variant burden in BBD controls: among the top 10 genes associated with case-status (p&lt;0.05), all had BC-risk odds ratio estimates less than 1.0. Immune infiltrates, assessed by CD45 H-scores, were positively associated with overall variant burden (r=0.42, p=0.0031), and noticeably higher in cancer-free controls (r=0.50, p=0.005). Comparison of mutation variants in 12 BBD cases that developed cancer within 5 years of diagnosis and their subsequent tumors revealed that more than half of the non-SNP variants were shared between the BBD and the subsequent tumor (54.6%), whereas less than one quarter were unique to the tumor (20.6%) or unique to the BBD (24.9%). Through BBD-BC paired variant analysis, three subjects were found with newly emergent clones uniquely appearing in the tumors and four subjects had loss of clones in the tumor that were present in the BBD. These results show that somatic variant burden is lower in patients with BBD that progress to BC, mutation profiles of BC progression are heterogeneous, and BBD tissues with reduced variant burden also show reduced immune infiltrate. This implicates a potential protective role for immune activation in BBD. Future directions will define the characteristics of the potentially protective immune response and its role in variant heterogeneity in a larger sample set, integrated with paired RNA sequencing to better understand the pathways involved and the relationship to immune infiltration. Citation Format: Stacey J. Winham, Julie M. Cunningham, Yuanhang Liu, Aditya V. Bhagwate, Ethan P. Heinzen, Samantha J. McDonough, Melody L. Stallings-Mann, Robert A. Vierkant, Jodi M. Carter, Mark E. Sherman, Derek C. Radisky, Amy C. Degnim, Chen Wang. Somatic variant burdenin benign breast disease and association with risk of subsequent breast cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 2274.

Abstract Background: Benign breast disease (BBD) is a risk factor for breast cancer (BC); however, little is known about the genetic alterations present at the time of BBD diagnosis, and how these relate to risk of incident BC. Methods: A subset of a long-term BBD cohort was selected to examine DNA variation across three BBD groups (42 future estrogen receptor-positive (ER+) BC, 36 future estrogen receptor-negative (ER-) BC and 42 controls cancer-free for at least 16 years post-BBD). DNA extracted from archival formalin fixed, paraffin-embedded (FFPE) tissue blocks was tested and analyzed for presence of DNA alterations using a targeted panel of 93 BC-associated genes. To address artifacts frequently observed in FFPE tissues (e.g., C>T changes), we applied three filtering strategies based on alternative allele frequencies and nucleotide substitution context. Gene-level associations were performed using two burden tests, SKAT-O and logistic regression, and adjusted for clinical and technical covariates including epithelial percentage, histologic impression, patient age and year of BBD biopsy. Results: After filtering, the variant frequency of SNPs in our samples was highly consistent with population allele frequencies reported in 1KG/ExAC (0.986, p <1e-16). Ten genes were found to be associated with later cancer status by four of 12 association methods, with nominal p-values <0.05: MED12, MSH2, BRIP1, PMS1, GATA3, MUC16, FAM175A, EXT2, MLH1 and TGFB1. Additionally, analysis revealed 11 and 19 gene-level associations with p-values <0.05 and OR<1 for all cancer cases and ER- cases, respectively. Comparing between the three case groups, BBD ER+ cases were closer to controls in mutation profile, while BBD ER- cases were distinct. Notably, the variant burden was significantly higher in controls than in either ER+ or ER- cases. CD45 expression was associated with mutational burden (p=0.003).Conclusions: Somatic mutations were more frequent in benign breast tissue from women who did not develop cancer, opening questions of clonal diversity or immune-mediated restraint on future cancer development. CD45 expression was positively associated with mutational burden, most strongly in controls. Further studies in both normal and premalignant tissues are needed to better understand the role of somatic gene mutations and their contribution to future cancer development.

Abstract Background Next Generation sequencing (NGS) is a powerful tool to identify somatic mutations associated with tumor onset and drug response. While it is well suited for high quality fresh/frozen samples, NGS is not proven for FFPE tissue which is the most common type of clinical specimen. Since the nucleic acids can be readily extracted from FFPE samples for a variety of genomic analyses, a comparative mutational analysis of paired frozen and FFPE tissues is urgently needed. Our long term goal is to establish a lab protocol to detect mutations in FFPE tumors using a targeted capture and sequencing approach for genes of interest. This pilot study focuses on the comparison of FFPE and frozen samples to test the validity of using FFPE tissues in such application. Methods Gene Selection: 128 genes associated with known pathogenic mutations in lymphoma Sample Selection: 9 diffuse large B-cell lymphoma (DLBCL) cases with FFPE, frozen and germline samples, as well as 10 frozen normal lymphatic tissues as references for CNV detections Capture Probe Design: We targeted coding exons and UTR, as well as the evolutionarily conserved intronic regions. The capture probes were designed using the Agilent eArray tool. The titling density of the probes was set to 3 probes overlapping with every base in the target region to improve the capture efficiency in FFPE samples. The least stringent masking of the repeat regions was allowed to include regions with small repeats that are shorter than the length of the sequencing reads (100-bp). In addition, boosting parameters were picked to set various levels of probe replication in different regions in order to minimize the local coverage differences (e.g. between regions of different GC contents) Sequencing and Bioinformatics: The target capture and sequencing were performed by the Mayo Clinic Medical Genome Facility. The reads were mapped to Human Reference Genome Build 37 using Novalign, and SNVs were called using GATK. The CNVs were identified using an in-house developed algorithm, patternCNV. Results The designed probes covered 99.65937% of the target regions. We generated 2.2-6.7 Gbp of reads per sample, 57.4-71.5% of which were on target. This equalled an average coverage of 2100-6700 folds which is 10-30 times higher than the minimal coverage recommended by Agilent. Due to this high coverage, we observed duplicate reads that accounted for 7.7-73.5% of the total reads. When we analysed the data with and without the duplicated reads, the concordance of the called SNVs was between 84-93% out of 207-249 mutated positions per trio-sample. There were 7.8-8.9% and 1.1-2.2% unique SNVs per sample by excluding or including duplicate reads, respectively. The dis-concordances were mostly missed calls, where a SNV was observed in only 1 or 2 of the trio samples. The missed calls from frozen samples ranged from 0-10.4% compared to 1.4-10.4% from the FFPE tissues, with 0.88-2.4% more SNVs missed in FFPE. Further analyses showed that all of the missing calls came from the lack of or low coverage of the corresponding positions. There were also differences of the called SNVs between the trio samples. However, this was extremely rare. Only 2 out of the 9 trio samples at a total of 3 positions had disagreements in called SNVs between FFPE and frozen tissues, all due to the allelic imbalance where the percentage of reads supporting the alternative alleles were below 20%. Therefore, this dis-concordance can be removed by back-filling of the read-level information for each position. Unfortunately only 11.9-47.4% of the CNVs called in frozen tissues were identified in FFPE samples, due to the widely various coverage in FFPE samples. The consequent large noises of the log ratio values between the FFPEs and normal references significantly reduced the sensitivity for CNV calling. Conclusions This pilot study compared the performance of SNV and CNV detection in FFPE and paired frozen tissues using a target capture and sequencing approach. With a capture probe design strategized to benefit FFPE samples, we observed SNV detection rates in FFPE that were only slightly lower (0.88-2.4%) than those of frozen tissues due to poor coverage of some positions in FFPE samples. With a proper back-filling step, there was no dis-concordance of the called SNVs between FFPE and frozen samples. However, CNV detections in FFPE were more problematic due to the un-predictable regional coverage in FFPE samples. Disclosures: No relevant conflicts of interest to declare.

Abstract Background: Archived formalin fixed paraffin embedded (FFPE) samples are valuable clinical resources to examine clinically relevant morphology features and also to study genetic changes. However, DNA quality and quantity of FFPE samples are sub-optimal, and resulting NGS-based genetics variant detections are prone to false positives. Evaluations of wet-lab and bioinformatics approaches are needed to optimize variant detection from FFPE samples. Results: As a pilot study, we designed within-subject triplicate samples of DNA derived from paired FFPE and fresh frozen breast tissues to highlight FFPE-specific artifacts. For FFPE samples, we tested two FFPE DNA extraction methods to determine impact of wet-lab procedures on variant calling: QIAGEN QIAamp DNA Mini Kit ("QA"), and QIAGEN GeneRead DNA FFPE Kit ("QGR"). We also used negative-control (NA12891) and positive control samples (Horizon Discovery Reference Standard FFPE). All DNA sample libraries were prepared for NGS according to the QIAseq Human Breast Cancer Targeted DNA Panel protocol and sequenced on the HiSeq 4000. Variant calling and filtering were performed using QIAGEN Gene Globe Data Portal. Detailed variant concordance comparisons and mutational signature analysis were performed to investigate effects of FFPE samples compared to paired fresh frozen samples, along with different DNA extraction methods. In this study, we found that five times or more variants were called with FFPE samples, compared to their paired fresh-frozen tissue samples even after applying molecular barcoding error-correction and default bioinformatics filtering recommended by the vendor. We also found that QGR as an optimized FFPE-DNA extraction approach leads to much fewer discordant variants between paired fresh frozen and FFPE samples. Approximately 92 % of the uniquely called FFPE variants were of low allelic frequency range (<5%), and collectively shared a “C>T|G>A” mutational signature known to be representative of FFPE artifacts resulting from cytosine deamination. Based on control samples and FFPE-frozen replicates, we derived an effective filtering strategy with associated empirical false-discovery estimates. Conclusions: Through this study, we demonstrated feasibility of calling and filtering genetic variants from FFPE tissue samples using a combined strategy with molecular barcodes, optimized DNA extraction, and bioinformatics methods incorporating genomics context such as mutational signature and variant allelic frequency.

Abstract Background Archived formalin fixed paraffin embedded (FFPE) samples are valuable clinical resources to examine clinically relevant morphology features and also to study genetic changes. However, DNA quality and quantity of FFPE samples are often sub-optimal, and resulting NGS-based genetics variant detections are prone to false positives. Evaluations of wet-lab and bioinformatics approaches are needed to optimize variant detection from FFPE samples. Results As a pilot study, we designed within-subject triplicate samples of DNA derived from paired FFPE and fresh frozen breast tissues to highlight FFPE-specific artifacts. For FFPE samples, we tested two FFPE DNA extraction methods to determine impact of wet-lab procedures on variant calling: QIAGEN QIAamp DNA Mini Kit ("QA"), and QIAGEN GeneRead DNA FFPE Kit ("QGR"). We also used negative-control (NA12891) and positive control samples (Horizon Discovery Reference Standard FFPE). All DNA sample libraries were prepared for NGS according to the QIAseq Human Breast Cancer Targeted DNA Panel protocol and sequenced on the HiSeq 4000. Variant calling and filtering were performed using QIAGEN Gene Globe Data Portal. Detailed variant concordance comparisons and mutational signature analysis were performed to investigate effects of FFPE samples compared to paired fresh frozen samples, along with different library preparations. In this study, we found that five times or more variants were called with FFPE samples, compared to their paired fresh-frozen tissue samples even after applying molecular barcoding error-correction and default bioinformatics filtering recommended by the vendor. We also found that QGR as an optimized FFPE-DNA extraction approach leads to much fewer discordant variants between paired fresh frozen and FFPE samples. Approximately 92% of the uniquely called FFPE variants were of low allelic frequency range (<5%), and collectively shared a “C>T|G>A” mutational signature known to be representative of FFPE artifacts resulting fromcytosine deamination. Based on control samples and FFPE-frozen replicates, we derived an effective filtering strategy with associated empirical false-discovery estimates. Conclusions Through this study, we demonstrated feasibility of calling and filtering genetic variants from FFPE tissue samples using a combined strategy with molecular barcodes, optimized DNA extraction, and bioinformatics methods incorporating genomics context such as mutational signature and variant allelic frequency.

Abstract Background Archived formalin fixed paraffin embedded (FFPE) samples are valuable clinical resources to examine clinically relevant morphology features and also to study genetic changes. However, DNA quality and quantity of FFPE samples are often sub-optimal, and resulting NGS-based genetics variant detections are prone to false positives. Evaluations of wet-lab and bioinformatics approaches are needed to optimize variant detection from FFPE samples. Results As a pilot study, we designed within-subject triplicate samples of DNA derived from paired FFPE and fresh frozen breast tissues to highlight FFPE-specific artifacts. For FFPE samples, we tested two FFPE DNA extraction methods to determine impact of wet-lab procedures on variant calling: QIAGEN QIAamp DNA Mini Kit ("QA"), and QIAGEN GeneRead DNA FFPE Kit ("QGR"). We also used negative-control (NA12891) and positive control samples (Horizon Discovery Reference Standard FFPE). All DNA sample libraries were prepared for NGS according to the QIAseq Human Breast Cancer Targeted DNA Panel protocol and sequenced on the HiSeq 4000. Variant calling and filtering were performed using QIAGEN Gene Globe Data Portal. Detailed variant concordance comparisons and mutational signature analysis were performed to investigate effects of FFPE samples compared to paired fresh frozen samples, along with different library preparations. In this study, we found that five times or more variants were called with FFPE samples, compared to their paired fresh-frozen tissue samples even after applying molecular barcoding error-correction and default bioinformatics filtering recommended by the vendor. We also found that QGR as an optimized FFPE-DNA extraction approach leads to much fewer discordant variants between paired fresh frozen and FFPE samples. Approximately 92% of the uniquely called FFPE variants were of low allelic frequency range (<5%), and collectively shared a “C>T|G>A” mutational signature known to be representative of FFPE artifacts resulting fromcytosine deamination. Based on control samples and FFPE-frozen replicates, we derived an effective filtering strategy with associated empirical false-discovery estimates. Conclusions Through this study, we demonstrated feasibility of calling and filtering genetic variants from FFPE tissue samples using a combined strategy with molecular barcodes, optimized DNA extraction, and bioinformatics methods incorporating genomics context such as mutational signature and variant allelic frequency.