Exome Sequencing Identifies Carriers of the Autosomal Dominant Cancer Predisposition Disorders Beyond Current Practice Guideline Recommendations.

PURPOSE: The autosomal dominant cancer predisposition disorders hereditary breast and ovarian cancer (HBOC) and Lynch syndrome (LS) are genetic conditions for which early identification and intervention have a positive effect on the individual and public health. The goals of this study were to determine whether germline genetic screening using exome sequencing could be used to efficiently identify carriers of HBOC and LS. METHODS: Participants were recruited from three geographically and racially diverse sites in the United States (Rochester, MN; Phoenix, AZ; Jacksonville, FL). Participants underwent Exome+ sequencing (Helix Inc, San Mateo, CA) and return of results for specific genetic findings: HBOC (BRCA1 and BRCA1) and LS (MLH1, MSH2, MSH6, PMS2, and EPCAM). Chart review was performed to collect demographics and personal and family cancer history. RESULTS: To date, 44,306 participants have enrolled in Tapestry. Annotation and interpretation of all variants in genes for HBOC and LS resulted in the identification of 550 carriers (prevalence, 1.24%), which included 387 with HBOC (27.2% BRCA1, 42.8% BRCA2) and 163 with LS (12.3% MSH6, 8.8% PMS2, 4.5% MLH1, 3.8% MSH2, and 0.2% EPCAM). More than half of these participants (52.1%) were newly diagnosed carriers with HBOC and LS. In all, 39.2% of HBOC/LS carriers did not satisfy National Comprehensive Cancer Network (NCCN) criteria for genetic evaluation. NCCN criteria were less commonly met in underrepresented minority populations versus self-reported White race (51.5% v 37.5%, P = .028). CONCLUSION: Our results emphasize the need for wider utilization of germline genetic sequencing for enhanced screening and detection of individuals who have LS and HBOC cancer predisposition syndromes.

De novo and inherited monoallelic variants in TUBA4A cause ataxia and spasticity.

Alpha-tubulin 4A encoding gene (TUBA4A) has been associated with familial amyotrophic lateral sclerosis (fALS) and fronto-temporal dementia (FTD), based on identification of likely pathogenic variants in patients from distinct ALS and FTD cohorts. By screening a multicentric French cohort of 448 unrelated probands presenting with cerebellar ataxia, we identified ultra-rare TUBA4A missense variants, all being absent from public databases and predicted pathogenic by multiple in-silico tools. In addition, gene burden analyses in the 100,000 genomes project (100KGP) showed enrichment of TUBA4A rare variants in the inherited ataxia group compared to controls (OR: 57.0847 [10.2- 576.7]; p = 4.02 x10-07). Altogether, we report 12 patients presenting with spasticity and/or cerebellar ataxia and harboring a predicted pathogenic TUBA4A missense mutation, including 5 confirmed de novo cases and a mutation previously reported in a large family presenting with spastic ataxia. Cultured fibroblasts from 3 patients harboring distinct TUBA4A missense showed significant alterations in microtubule organisation and dynamics, providing insight of TUBA4A variants pathogenicity. Our data confirm the identification of a hereditary spastic ataxia disease gene with variable age of onset, expanding the clinical spectrum of TUBA4A associated phenotypes.

The expanding clinical and genetic spectrum of DYNC1H1-related disorders.

Intracellular trafficking involves an intricate machinery of motor complexes including the dynein complex to shuttle cargo for autophagolysosomal degradation. Deficiency in dynein axonemal chains as well as cytoplasmic light and intermediate chains have been linked with ciliary dyskinesia and skeletal dysplasia. The cytoplasmic dynein 1 heavy chain protein (DYNC1H1) serves as a core complex for retrograde trafficking in neuronal axons. Dominant pathogenic variants in DYNC1H1 have been previously implicated in peripheral neuromuscular disorders (NMD) and neurodevelopmental disorders (NDD). As heavy-chain dynein is ubiquitously expressed, the apparent selectivity of heavy-chain dyneinopathy for motor neuronal phenotypes remains currently unaccounted for. Here, we aimed to evaluate the full DYNC1H1-related clinical, molecular and imaging spectrum, including multisystem features and novel phenotypes presenting throughout life. We identified 47 cases from 43 families with pathogenic heterozygous variants in DYNC1H1 (aged 0-59 years) and collected phenotypic data via a comprehensive standardized survey and clinical follow-up appointments. Most patients presented with divergent and previously unrecognized neurological and multisystem features, leading to significant delays in genetic testing and establishing the correct diagnosis. Neurological phenotypes include novel autonomic features, previously rarely described behavioral disorders, movement disorders, and periventricular lesions. Sensory neuropathy was identified in nine patients (median age of onset 10.6 years), of which five were only diagnosed after the second decade of life, and three had a progressive age-dependent sensory neuropathy. Novel multisystem features included primary immunodeficiency, bilateral sensorineural hearing loss, organ anomalies, and skeletal manifestations, resembling the phenotypic spectrum of other dyneinopathies. We also identified an age-dependent biphasic disease course with developmental regression in the first decade and, following a period of stability, neurodegenerative progression after the second decade of life. Of note, we observed several cases in whom neurodegeneration appeared to be prompted by intercurrent systemic infections with double-stranded DNA viruses (Herpesviridae) or single-stranded RNA viruses (Ross-River fever, SARS-CoV-2). Moreover, the disease course appeared to be exacerbated by viral infections regardless of age and/or severity of NDD manifestations, indicating a role of dynein in anti-viral immunity and neuronal health. In summary, our findings expand the clinical, imaging, and molecular spectrum of pathogenic DYNC1H1 variants beyond motor neuropathy disorders and suggest a life-long continuum and age-related progression due to deficient intracellular trafficking. This study will facilitate early diagnosis and improve counselling and health surveillance of affected patients.

Novel protein-truncating variants of a chromatin-modifying gene MSL2 in syndromic neurodevelopmental disorders.

Numerous large scale genomic studies have uncovered rare but recurrent pathogenetic variants in a significant number of genes encoding epigenetic machinery in cases with neurodevelopmental disorders (NDD) especially autism spectrum disorder (ASD). These findings provide strong support for the functional importance of epigenetic regulators in neurodevelopment. After the clinical genomics evaluation of the patients using exome sequencing, we have identified, three novel protein-truncating variants (PTVs) in the MSL2 gene (OMIM: 614802) which encodes a chromatin modifying enzyme. MSL2 modifies chromatin through both mono-ubiquitination of histone 2B on lysine 34 (K34) and acetylation of histone H4 on lysine 16 (K16). We reported first time the detailed clinical features associated with 3 MSL2 PTVs. There are 15 PTVs (13 de novo) reported from the large genomics studies (12 cases) or ClinVar (3 cases) of NDD, ASD, and developmental disorders (DD) but the specific clinical features for these cases are not described. Taken together, our descriptions of dysmorphic face and other features support the causal role of MSL2 in a likely syndromic neurodevelopmental disorder and add MSL2 to a growing list of epigenetic genes implicated in ASD.

Diagnostic yield of exome and genome sequencing after non-diagnostic multi-gene panels in patients with single-system diseases.

BACKGROUND: Though next-generation sequencing (NGS) tests like exome sequencing (ES), genome sequencing (GS), and panels derived from exome and genome data (EGBP) are effective for rare diseases, the ideal diagnostic approach is debated. Limited research has explored reanalyzing raw ES and GS data post-negative EGBP results for diagnostics. RESULTS: We analyzed complete ES/GS raw sequencing data from Mayo Clinic's Program for Rare and Undiagnosed Diseases (PRaUD) patients to assess whether supplementary findings could augment diagnostic yield. ES data from 80 patients (59 adults) and GS data from 20 patients (10 adults), averaging 43 years in age, were analyzed. Most patients had renal (n=44) and auto-inflammatory (n=29) phenotypes. Ninety-six cases had negative findings and in four cases additional genetic variants were found, including a variant related to a recently described disease (RRAGD-related hypomagnesemia), a variant missed due to discordant inheritance pattern (COL4A3), a variant with high allelic frequency (NPHS2) in the general population, and a variant associated with an initially untargeted phenotype (HNF1A). CONCLUSION: ES and GS show diagnostic yields comparable to EGBP for single-system diseases. However, EGBP's limitations in detecting new disease-associated genes underscore the necessity for periodic updates.

The Natural History of Splenic Artery Aneurysms: Factors That Predict Aneurysm Growth.

PURPOSE: To examine the natural history of splenic artery aneurysms (SAAs) at a single institution and assess the effect of patient factors and aneurysm characteristics on aneurysm growth. MATERIALS AND METHODS: This single-center retrospective study included patients with SAAs who underwent serial imaging over 30 years (1990-2020). Data regarding patient demographics and aneurysm characteristics were collected. The variables contributing to aneurysm growth were assessed using nonparametric tests for continuous variables and chi-square test for categorical variables. Multivariable linear regression was performed using aneurysm growth rate as a continuous dependent variable. RESULTS: A total of 132 patients were included in this study. The median maximum diameter of the SAAs was 15.8 mm (range, 4.0-50.0 mm). Growth over time was observed in 39% of the aneurysms, whereas the remaining 61% were stable in size. Of aneurysms that increased in size, the median aneurysm growth rate was 0.60 mm/y (range, 0.03-5.00 mm/y). Maximum aneurysm diameter of >2 cm and the presence of >50% mural thrombus were significant positive predictors for aneurysm growth (P = .020 and P = .022, respectively). Greater than 50% rim calcification was a significant negative predictor for aneurysm growth (P = .009) in multivariate analysis. CONCLUSIONS: A larger baseline SAA size, presence of mural thrombus, and lack of rim calcification are associated with increased aneurysm growth rate.

OmicsFootPrint: a framework to integrate and interpret multi-omics data using circular images and deep neural networks.

The OmicsFootPrint framework addresses the need for advanced multi-omics data analysis methodologies by transforming data into intuitive two-dimensional circular images and facilitating the interpretation of complex diseases. Utilizing Deep Neural Networks and incorporating the SHapley Additive exPlanations (SHAP) algorithm, the framework enhances model interpretability. Tested with The Cancer Genome Atlas (TCGA) data, OmicsFootPrint effectively classified lung and breast cancer subtypes, achieving high Area Under Curve (AUC) scores - 0.98±0.02 for lung cancer subtype differentiation, 0.83±0.07 for breast cancer PAM50 subtypes, and successfully distinguishe between invasive lobular and ductal carcinomas in breast cancer, showcasing its robustness. It also demonstrated notable performance in predicting drug responses in cancer cell lines, with a median AUC of 0.74, surpassing existing algorithms. Furthermore, its effectiveness persists even with reduced training sample sizes. OmicsFootPrint marks an enhancement in multi-omics research, offering a novel, efficient, and interpretable approach that contributes to a deeper understanding of disease mechanisms.

Identification of skewed X chromosome inactivation using exome and transcriptome sequencing in patients with suspected rare genetic disease.

BACKGROUND: X-chromosome inactivation (XCI) is an epigenetic process that occurs during early development in mammalian females by randomly silencing one of two copies of the X chromosome in each cell. The preferential inactivation of either the maternal or paternal copy of the X chromosome in a majority of cells results in a skewed or non-random pattern of X inactivation and is observed in over 25% of adult females. Identifying skewed X inactivation is of clinical significance in patients with suspected rare genetic diseases due to the possibility of biased expression of disease-causing genes present on the active X chromosome. The current clinical test for the detection of skewed XCI relies on the methylation status of the methylation-sensitive restriction enzyme (Hpall) binding site present in proximity of short tandem polymorphic repeats on the androgen receptor (AR) gene. This approach using one locus results in uninformative or inconclusive data for 10-20% of tests. Further, recent studies have shown inconsistency between methylation of the AR locus and the state of inactivation of the X chromosome. Herein, we develop a method for estimating X inactivation status, using exome and transcriptome sequencing data derived from blood in 227 female samples. We built a reference model for evaluation of XCI in 135 females from the GTEx consortium. We tested and validated the model on 11 female individuals with different types of undiagnosed rare genetic disorders who were clinically tested for X-skew using the AR gene assay and compared results to our outlier-based analysis technique. RESULTS: In comparison to the AR clinical test for identification of X inactivation, our method was concordant with the AR method in 9 samples, discordant in 1, and provided a measure of X inactivation in 1 sample with uninformative clinical results. We applied this method on an additional 81 females presenting to the clinic with phenotypes consistent with different hereditary disorders without a known genetic diagnosis. CONCLUSIONS: This study presents the use of transcriptome and exome sequencing data to provide an accurate and complete estimation of X-inactivation and skew status in a cohort of female patients with different types of suspected rare genetic disease.

Bi-allelic variants in CELSR3 are implicated in central nervous system and urinary tract anomalies.

CELSR3 codes for a planar cell polarity protein. We describe twelve affected individuals from eleven independent families with bi-allelic variants in CELSR3. Affected individuals presented with an overlapping phenotypic spectrum comprising central nervous system (CNS) anomalies (7/12), combined CNS anomalies and congenital anomalies of the kidneys and urinary tract (CAKUT) (3/12) and CAKUT only (2/12). Computational simulation of the 3D protein structure suggests the position of the identified variants to be implicated in penetrance and phenotype expression. CELSR3 immunolocalization in human embryonic urinary tract and transient suppression and rescue experiments of Celsr3 in fluorescent zebrafish reporter lines further support an embryonic role of CELSR3 in CNS and urinary tract formation.

Semiautomated approach focused on new genomic information results in time and effort-efficient reannotation of negative exome data.

Most rare disease patients (75-50%) undergoing genomic sequencing remain unsolved, often due to lack of information about variants identified. Data review over time can leverage novel information regarding disease-causing variants and genes, increasing this diagnostic yield. However, time and resource constraints have limited reanalysis of genetic data in clinical laboratories setting. We developed RENEW, (REannotation of NEgative WES/WGS) an automated reannotation procedure that uses relevant new information in on-line genomic databases to enable rapid review of genomic findings. We tested RENEW in an unselected cohort of 1066 undiagnosed cases with a broad spectrum of phenotypes from the Mayo Clinic Center for Individualized Medicine using new information in ClinVar, HGMD and OMIM between the date of previous analysis/testing and April of 2022. 5741 variants prioritized by RENEW were rapidly reviewed by variant interpretation specialists. Mean analysis time was approximately 20 s per variant (32 h total time). Reviewed cases were classified as: 879 (93.0%) undiagnosed, 63 (6.6%) putatively diagnosed, and 4 (0.4%) definitively diagnosed. New strategies are needed to enable efficient review of genomic findings in unsolved cases. We report on a fast and practical approach to address this need and improve overall diagnostic success in patient testing through a recurrent reannotation process.

A supervised learning method for classifying methylation disorders.

BACKGROUND: DNA methylation is one of the most stable and well-characterized epigenetic alterations in humans. Accordingly, it has already found clinical utility as a molecular biomarker in a variety of disease contexts. Existing methods for clinical diagnosis of methylation-related disorders focus on outlier detection in a small number of CpG sites using standardized cutoffs which differentiate healthy from abnormal methylation levels. The standardized cutoff values used in these methods do not take into account methylation patterns which are known to differ between the sexes and with age. RESULTS: Here we profile genome-wide DNA methylation from blood samples drawn from within a cohort composed of healthy controls of different age and sex alongside patients with Prader-Willi syndrome (PWS), Beckwith-Wiedemann syndrome, Fragile-X syndrome, Angelman syndrome, and Silver-Russell syndrome. We propose a Generalized Additive Model to perform age and sex adjusted outlier analysis of around 700,000 CpG sites throughout the human genome. Utilizing z-scores among the cohort for each site, we deployed an ensemble based machine learning pipeline and achieved a combined prediction accuracy of 0.96 (Binomial 95% Confidence Interval 0.868[Formula: see text]0.995). CONCLUSION: We demonstrate a method for age and sex adjusted outlier detection of differentially methylated loci based on a large cohort of healthy individuals. We present a custom machine learning pipeline utilizing this outlier analysis to classify samples for potential methylation associated congenital disorders. These methods are able to achieve high accuracy when used with machine learning methods to classify abnormal methylation patterns.

Synergizing Thermal Ablation Modalities with Immunotherapy: Enough to Induce Systemic Antitumoral Immunity?

Thermal ablation modalities (cryoablation, radiofrequency ablation, and microwave ablation) have long been noted to occasionally induce a systemic antitumoral response. With the widespread use of checkpoint inhibitors, there is a significant interest in whether thermal ablation can promote immune system tumor recognition and increase checkpoint inhibitor response rates. In this review, we examine the current state of preclinical and clinical evidence examining the combination of checkpoint inhibitor therapies and thermal ablation modalities as well as discuss remaining the unanswered questions and directions for future research.

Nuclear Abnormalities in LMNA p.(Glu2Lys) Variant Segregating with LMNA-Associated Cardiocutaneous Progeria Syndrome.

The LMNA gene encodes lamin A and lamin C, which play important roles in nuclear organization. Pathogenic variants in LMNA cause laminopathies, a group of disorders with diverse phenotypes. There are two main groups of disease-causing variants: missense variants affecting dimerization and intermolecular interactions, and heterozygous substitutions activating cryptic splice sites. These variants lead to different disorders, such as dilated cardiomyopathy and Hutchinson-Gilford progeria (HGP). Among these, the phenotypic terms for LMNA-associated cardiocutaneous progeria syndrome (LCPS), which does not alter lamin A processing and has an older age of onset, have been described. Here, we present the workup of an LMNA variant of uncertain significance, NM_170707.2 c. 4G>A, p.(Glu2Lys), in a 36-year-old female with severe calcific aortic stenosis, a calcified mitral valve, premature aging, and a family history of similar symptoms. Due to the uncertainty of in silico predictions for this variant, an assessment of nuclear morphology was performed using the immunocytochemistry of stable cell lines to indicate whether the p.(Glu2Lys) had a similar pathogenic mechanism as a previously described pathogenic variant associated with LCPS, p.Asp300Gly. Indirect immunofluorescence analysis of nuclei from stable cell lines showed abnormal morphology, including lobulation and occasional ringed nuclei. Relative to the controls, p.Glu2Lys and p.Asp300Gly nuclei had significantly (p < 0.001) smaller average nuclear areas than controls (mean = 0.10 units, SD = 0.06 for p.Glu2Lys; and mean = 0.09 units, SD = 0.05 for p.Asp300Gly versus mean = 0.12, SD = 0.05 for WT). After functional studies and segregation studies, this variant was upgraded to likely pathogenic. In summary, our findings suggest that p.Glu2Lys impacts nuclear morphology in a manner comparable to what was observed in p.Asp300Gly cells, indicating that the variant is the likely cause of the LCPS segregating within this family.

Complex balanced intrachromosomal rearrangement involving PITX2 identified as a cause of Axenfeld-Rieger Syndrome.

Axenfeld-Rieger Syndrome (ARS) type 1 is a rare autosomal dominant condition characterized by anterior chamber anomalies, umbilical defects, dental hypoplasia, and craniofacial anomalies, with Meckel's diverticulum in some individuals. Here, we describe a clinically ascertained female of childbearing age with ARS for whom clinical targeted sequencing and deletion/duplication analysis followed by clinical exome and genome sequencing resulted in no pathogenic variants or variants of unknown significance in PITX2 or FOXC1. Advanced bioinformatic analysis of the genome data identified a complex, balanced rearrangement disrupting PITX2. This case is the first reported intrachromosomal rearrangement leading to ARS, illustrating that for patients with compelling clinical phenotypes but negative genomic testing, additional bioinformatic analysis are essential to identify subtle genomic abnormalities in target genes.

Further delineation of the rare GDACCF (global developmental delay, absent or hypoplastic corpus callosum, dysmorphic facies syndrome): genotype and phenotype of 22 patients with ZNF148 mutations.

BACKGROUND: Pathogenic variants in the zinc finger protein coding genes are rare causes of intellectual disability and congenital malformations. Mutations in the ZNF148 gene causing GDACCF syndrome (global developmental delay, absent or hypoplastic corpus callosum, dysmorphic facies; MIM #617260) have been reported in five individuals so far. METHODS: As a result of an international collaboration using GeneMatcher Phenome Central Repository and personal communications, here we describe the clinical and molecular genetic characteristics of 22 previously unreported individuals. RESULTS: The core clinical phenotype is characterised by developmental delay particularly in the domain of speech development, postnatal growth retardation, microcephaly and facial dysmorphism. Corpus callosum abnormalities appear less frequently than suggested by previous observations. The identified mutations concerned nonsense or frameshift variants that were mainly located in the last exon of the ZNF148 gene. Heterozygous deletion including the entire ZNF148 gene was found in only one case. Most mutations occurred de novo, but were inherited from an affected parent in two families. CONCLUSION: The GDACCF syndrome is clinically diverse, and a genotype-first approach, that is, exome sequencing is recommended for establishing a genetic diagnosis rather than a phenotype-first approach. However, the syndrome may be suspected based on some recurrent, recognisable features. Corpus callosum anomalies were not as constant as previously suggested, we therefore recommend to replace the term 'GDACCF syndrome' with 'ZNF148-related neurodevelopmental disorder'.

Percutaneous Intratumoral Immunoadjuvant Gel Increases the Abscopal Effect of Cryoablation for Checkpoint Inhibitor Resistant Cancer.

Percutaneous cryoablation is a common clinical therapy for metastatic and primary cancer. There are rare clinical reports of cryoablation inducing regression of distant metastases, known as the "abscopal" effect. Intratumoral immunoadjuvants may be able to augment the abscopal rate of cryoablation, but existing intratumoral therapies suffer from the need for frequent injections and inability to confirm target delivery, leading to poor clinical trial outcomes. To address these shortcomings, an injectable thermoresponsive gel-based controlled release formulation is developed for the FDA-approved Toll-like-receptor 7 (TLR7) agonist imiquimod ("Imigel") that forms a tumor-resident depot upon injection and contains a contrast agent for visualization under computed tomography (CT). The poly-lactic-co-glycolic acid-polyethylene glycol-poly-lactic-co-glycolic acid (PLGA-PEG-PLGA)-based amphiphilic copolymer gel's underlying micellar nature enables high drug concentration and a logarithmic release profile that is additive with the neo-antigen release from cryoablation, requiring only a single injection. Rheological testing demonstrated the thermoresponsive increase in viscosity at body temperature and radio-opacity via microCT. Its ability to significantly augment the abscopal rate of cryoablation is demonstrated in otherwise immunotherapy resistant metastatic tumors in two aggressive colorectal and breast cancer dual tumor models with an all or nothing response, responders generally demonstrating complete regression of bilateral tumors in 90-day survival studies.

MIDOS: a novel stochastic model towards a treatment planning system for microsphere dosimetry in liver tumors.

PURPOSE: Transarterial radioembolization (TARE) procedures treat liver tumors by injecting radioactive microspheres into the hepatic artery. Currently, there is a critical need to optimize TARE towards a personalized dosimetry approach. To this aim, we present a novel microsphere dosimetry (MIDOS) stochastic model to estimate the activity delivered to the tumor(s), normal liver, and lung. METHODS: MIDOS incorporates adult male/female liver computational phantoms with the hepatic arterial, hepatic portal venous, and hepatic venous vascular trees. Tumors can be placed in both models at user discretion. The perfusion of microspheres follows cluster patterns, and a Markov chain approach was applied to microsphere navigation, with the terminal location of microspheres determined to be in either normal hepatic parenchyma, hepatic tumor, or lung. A tumor uptake model was implemented to determine if microspheres get lodged in the tumor, and a probability was included in determining the shunt of microspheres to the lung. A sensitivity analysis of the model parameters was performed, and radiation segmentectomy/lobectomy procedures were simulated over a wide range of activity perfused. Then, the impact of using different microspheres, i.e., SIR-Sphere®, TheraSphere®, and QuiremSphere®, on the tumor-to-normal ratio (TNR), lung shunt fraction (LSF), and mean absorbed dose was analyzed. RESULTS: Highly vascularized tumors translated into increased TNR. Treatment results (TNR and LSF) were significantly more variable for microspheres with high particle load. In our scenarios with 1.5 GBq perfusion, TNR was maximum for TheraSphere® at calibration time in segmentectomy/lobar technique, for SIR-Sphere® at 1-3 days post-calibration, and regarding QuiremSphere® at 3 days post-calibration. CONCLUSION: This novel approach is a decisive step towards developing a personalized dosimetry framework for TARE. MIDOS assists in making clinical decisions in TARE treatment planning by assessing various delivery parameters and simulating different tumor uptakes. MIDOS offers evaluation of treatment outcomes, such as TNR and LSF, and quantitative scenario-specific decisions.

Genomic epidemiology reveals the dominance of Hennepin County in the transmission of SARS-CoV-2 in Minnesota from 2020 to 2022.

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.