Fogarty, Zachary C

Abstract Multiplexed Immunofluorescence (MxIF) enables detailed immune cell phenotyping, providing critical insights into cell behavior within the tumor immune microenvironment (TIME). However, signal integrity can be compromised due to the complex cyclic staining processes inherent to MxIF. Hematoxylin and Eosin (H&E) staining, on the other hand, offers complementary information through its depiction of cell morphology and texture patterns and is often visually cross-referenced with MxIF in clinical settings. In this study, we proposed a novel framework to align H&E and MxIF images for precise cross-modal cell feature validation. Using cell detection outputs from each modality as anchors, we formulated the multimodal image registration problem as point set alignment. Coherent Point Drift (CPD) is employed for initial alignment, followed by Graph Matching (GM) for refinement. Evaluations on ovarian cancer tissue microarrays (TMAs) demonstrate that our method achieves high alignment accuracy, enabling reliable validation of cell-level features across modalities for both restained and serial sections. Our results indicate that restained H&E enhances confidence in findings derived from MxIF. Additionally, we demonstrated the feasibility of generating high-quality virtual H&E images from MxIF data when restained H&E is unavailable, offering a viable alternative for integrated multimodal analysis.

Abstract Advances in genetic technology have led to the increasing use of genomic panels in precision oncology practice, with panels ranging from a couple to hundreds of genes. However, the clinical utilization and utility of oncology genomic panels, especially among vulnerable populations, is unclear. We examined the association of panel size with socioeconomic status and clinical trial matching. We retrospectively identified 9,886 eligible adult subjects in the Mayo Clinic Health System who underwent genomic testing between January 1, 2016 and June 30, 2020. Patient data were retrieved from structured and unstructured data sources of institutional collections, including cancer registries, clinical data warehouses, and clinical notes. Socioeconomic surrogates were approximated using the Area Deprivation Index (ADI) corresponding to primary residence addresses. Logistic regression was performed to analyze relationships between ADI or rural/urban status and (i) use of genomic test by panel size; (ii) clinical trial matching status. Compared with patients from the most affluent areas, patients had a lower odds of receiving a panel test (vs. a single-gene test) if from areas of higher socioeconomic deprivation [OR (95% confidence interval (CI): 0.71 (0.61–0.83), P < 0.01] or a rural area [OR (95% CI): 0.85 (0.76–0.96), P < 0.01]. Patients in areas of higher socioeconomic deprivation were less likely to be matched to clinical trials if receiving medium panel tests [(OR) (95% CI): 0.69 (0.49–0.97), P = 0.03]; however, there was no difference among patients receiving large panel tests (P > 0.05) and rural patients were almost 2x greater odds of being matched if receiving a large panel test [(OR) (95% CI): 1.76 (1.21–2.55), P < 0.01]. Significance: We identified socioeconomic and rurality disparities in the use of genomic tests and trial matching by panel size, which may have implications for equal access to targeted therapies. The lack of association between large panel tests and clinical trial matching by socioeconomic status, suggests a potential health equity impact, while removing barriers in access to large panels for rural patients may improve access to trials. However, further research is needed.

Hypoplastic left heart syndrome (HLHS) is a severe congenital heart disease (CHD) with a likely oligogenic etiology, but our understanding of the genetic complexities and pathogenic mechanisms leading to HLHS is limited. We performed whole genome sequencing (WGS) on 183 HLHS patient-parent trios to identify candidate genes, which were functionally tested in the Drosophila heart model. Bioinformatic analysis of WGS data from an index family of a HLHS proband born to consanguineous parents prioritized 9 candidate genes with rare, predicted damaging homozygous variants. Of them, cardiac-specific knockdown (KD) of mitochondrial MICOS complex subunit dCHCHD3/6 resulted in drastically compromised heart contractility, diminished levels of sarcomeric actin and myosin, reduced cardiac ATP levels, and mitochondrial fission-fusion defects. These defects were similar to those inflicted by cardiac KD of ATP synthase subunits of the electron transport chain (ETC), consistent with the MICOS complex’s role in maintaining cristae morphology and ETC assembly. Five additional HLHS probands harbored rare, predicted damaging variants in CHCHD3 or CHCHD6 . Hypothesizing an oligogenic basis for HLHS, we tested 60 additional prioritized candidate genes from these patients for genetic interactions with CHCHD3/6 in sensitized fly hearts. Moderate KD of CHCHD3/6 in combination with Cdk12 (activator of RNA polymerase II), RNF149 ( goliath , E3 ubiquitin ligase), or SPTBN1 ( β-Spectrin, scaffolding protein) caused synergistic heart defects, suggesting the likely involvement of diverse pathways in HLHS. Further elucidation of novel candidate genes and genetic interactions of potentially disease-contributing pathways is expected to lead to a better understanding of HLHS and other CHDs.

Background: Deciphering hypoplastic left heart syndrome (HLHS) pathogenesis is confounded by its genetic heterogeneity and oligogenic underpinnings. Methods: Whole genome sequences were analyzed by 3 independent strategies to identify HLHS gene candidates, ranked by variant, gene, and disease-level metrics. Results: First, a genome-wide association study of 174 cases and 853 controls revealed suggestive association with a MYO18B intron 33 variant (rs2269628-G; frequency=0.55 versus 0.39; OR, 1.97 [95% CI, 1.54–2.52]; P =6.70×10 − 8 ). Second, transmission disequilibrium testing of 161 HLHS proband-parent trios revealed overrepresentation of a MYO18B intron 42 variant (rs73154186-A; frequency=0.05; OR, 24 [95% CI, 3.2–177.4]; P =4.23×10 − 6 ). Third, rare, predicted-damaging variants were filtered in 2 multiplex families. In 141H, 2 fifth-degree relatives with HLHS shared a paternally-inherited MYO5A missense variant (p.Arg801Trp; frequency=0.00003; combined annotation-dependent depletion score=29), each with a maternally-inherited or de novo candidate modifier variant in a MYO5A-interacting conventional myosin. In 442H, a HLHS proband was compound heterozygous for MYO15A variants—a maternally-inherited pathogenic stop-gain variant co-segregating with tetralogy of Fallot and bicuspid aortic valve in maternal relatives (p.Tyr2819Ter; frequency=0.00003) and a paternally-inherited intronic variant altering a canonical transcription factor binding site (rs1277068603; frequency=0.00001; position weight matrix score=0.98). Conclusions: Collectively, these findings suggest that common and rare alleles within unconventional myosin genes are associated with HLHS susceptibility. The identified candidate MYO18B regulates cardiac sarcomerogenesis, supporting the hypothesis of intrinsic myogenic perturbation in arrested left heart development.

Background: Whole-genome sequencing in families enables deciphering of congenital heart disease causes. A shared genetic basis for familial bicuspid aortic valve (BAV) and hypoplastic left heart syndrome (HLHS) was postulated. Methods: Whole-genome sequencing was performed in affected members of 6 multiplex BAV families, an HLHS cohort of 197 probands and 546 relatives, and 813 controls. Data were filtered for rare, predicted-damaging variants that cosegregated with familial BAV and disrupted genes associated with congenital heart disease in humans and mice. Candidate genes were further prioritized by rare variant burden testing in HLHS cases versus controls. Modifier variants in HLHS proband-parent trios were sought to account for the severe developmental phenotype. Results: In 5 BAV families, missense variants in 6 ontologically diverse genes for structural ( SPTBN1 , PAXIP1 , and FBLN1 ) and signaling ( CELSR1 , PLXND1 , and NOS3 ) proteins fulfilled filtering metrics. CELSR1 , encoding cadherin epidermal growth factor laminin G seven-pass G-type receptor, was identified as a candidate gene in 2 families and was the only gene demonstrating rare variant enrichment in HLHS probands ( P =0.003575). HLHS-associated CELSR1 variants included 16 missense, one splice site, and 3 noncoding variants predicted to disrupt canonical transcription factor binding sites, most of which were inherited from a parent without congenital heart disease. Filtering whole-genome sequencing data for rare, predicted-damaging variants inherited from the other parent revealed 2 cases of CELSR1 compound heterozygosity, one case of CELSR1 - CELSR3 synergistic heterozygosity, and 4 cases of CELSR1 - MYO15A digenic heterozygosity. Conclusions: CELSR1 is a susceptibility gene for familial BAV and HLHS, further implicating planar cell polarity pathway perturbation in congenital heart disease.

Prostate‐specific membrane antigen (PSMA) is a validated target for molecular diagnostics and targeted radionuclide therapy. Our purpose was to evaluate PSMA expression in hepatocellular carcinoma (HCC), cholangiocarcinoma (CCA), and hepatic adenoma (HCA); investigate the genetic pathways in HCC associated with PSMA expression; and evaluate HCC detection rate with 68Ga‐PSMA‐11 positron emission tomography (PET). In phase 1, PSMA immunohistochemistry (IHC) on HCC (n = 148), CCA (n = 111), and HCA (n = 78) was scored. In a subset (n = 30), messenger RNA (mRNA) data from the Cancer Genome Atlas HCC RNA sequencing were correlated with PSMA expression. In phase 2, 68Ga‐PSMA‐11 PET was prospectively performed in patients with treatment‐naïve HCC on a digital PET scanner using cyclotron‐produced 68Ga. Uptake was graded qualitatively and semi‐quantitatively using standard metrics. On IHC, PSMA expression was significantly higher in HCC compared with CCA and HCA (P < 0.0001); 91% of HCCs (n = 134) expressed PSMA, which principally localized to tumor‐associated neovasculature. Higher tumor grade was associated with PSMA expression (P = 0.012) but there was no association with tumor size (P = 0.14), fibrosis (P = 0.35), cirrhosis (P = 0.74), hepatitis B virus (P = 0.31), or hepatitis C virus (P = 0.15). Overall survival tended to be longer in patients without versus with PSMA expression (median overall survival: 4.2 vs. 1.9 years; P = 0.273). FGF14 (fibroblast growth factor 14) mRNA expression correlated positively (rho = 0.70; P = 1.70 × 10‐5) and MAD1L1 (Mitotic spindle assembly checkpoint protein MAD1) correlated negatively with PSMA expression (rho = −0.753; P = 1.58 × 10‐6). Of the 190 patients who met the eligibility criteria, 31 patients with 39 HCC lesions completed PET; 64% (n = 25) lesions had pronounced 68Ga‐PSMA‐11 standardized uptake value: SUVmax (median [range] 9.2 [4.9‐28.4]), SUVmean 4.7 (2.4‐12.7), and tumor‐to‐liver background ratio 2 (1.1‐11). Conclusion: Ex vivo expression of PSMA in neovasculature of HCC translates to marked tumor avidity on 68Ga‐PSMA‐11 PET, which suggests that PSMA has the potential as a theranostic target in patients with HCC.

Congenital heart diseases (CHDs), including hypoplastic left heart syndrome (HLHS), are genetically complex and poorly understood. Here, a multidisciplinary platform was established to functionally evaluate novel CHD gene candidates, based on whole-genome and iPSC RNA sequencing of a HLHS family-trio. Filtering for rare variants and altered expression in proband iPSCs prioritized 10 candidates. siRNA/RNAi-mediated knockdown in healthy human iPSC-derived cardiomyocytes (hiPSC-CM) and in developing Drosophila and zebrafish hearts revealed that LDL receptor-related protein LRP2 is required for cardiomyocyte proliferation and differentiation. Consistent with hypoplastic heart defects, compared to parents the proband’s iPSC-CMs exhibited reduced proliferation. Interestingly, rare, predicted-damaging LRP2 variants were enriched in a HLHS cohort; however, understanding their contribution to HLHS requires further investigation. Collectively, we have established a multi-species high-throughput platform to rapidly evaluate candidate genes and their interactions during heart development, which are crucial first steps toward deciphering oligogenic underpinnings of CHDs, including hypoplastic left hearts.

AbstractX chromosome inactivation (XCI) is a key epigenetic gene expression regulatory process, which may play a role in women’s cancer. In particular tissues, some genes are known to escape XCI, yet patterns of XCI in ovarian cancer (OC) and their clinical associations are largely unknown. To examine XCI in OC, we integrated germline genotype with tumor copy number, gene expression and DNA methylation information from 99 OC patients. Approximately 10% of genes showed different XCI status (either escaping or being subject to XCI) compared with the studies of other tissues. Many of these genes are known oncogenes or tumor suppressors (e.g. DDX3X, TRAPPC2 and TCEANC). We also observed strong association between cis promoter DNA methylation and allele-specific expression imbalance (P = 2.0 × 10−10). Cluster analyses of the integrated data identified two molecular subgroups of OC patients representing those with regulated (N = 47) and dysregulated (N = 52) XCI. This XCI cluster membership was associated with expression of X inactive specific transcript (P = 0.002), a known driver of XCI, as well as age, grade, stage, tumor histology and extent of rl disease following surgical debulking. Patients with dysregulated XCI (N = 52) had shorter time to recurrence (HR = 2.34, P = 0.001) and overall survival time (HR = 1.87, P = 0.02) than those with regulated XCI, although results were attenuated after covariate adjustment. Similar findings were observed when restricted to high-grade serous tumors. We found evidence of a unique OC XCI profile, suggesting that XCI may play an important role in OC biology. Additional studies to examine somatic changes with paired tumor-normal tissue are needed.

ABSTRACTX‐chromosome inactivation (XCI) epigenetically silences transcription of an X chromosome in females; patterns of XCI are thought to be aberrant in women's cancers, but are understudied due to statistical challenges. We develop a two‐stage statistical framework to assess skewed XCI and evaluate gene‐level patterns of XCI for an individual sample by integration of RNA sequence, copy number alteration, and genotype data. Our method relies on allele‐specific expression (ASE) to directly measure XCI and does not rely on male samples or paired normal tissue for comparison. We model ASE using a two‐component mixture of beta distributions, allowing estimation for a given sample of the degree of skewness (based on a composite likelihood ratio test) and the posterior probability that a given gene escapes XCI (using a Bayesian beta‐binomial mixture model). To illustrate the utility of our approach, we applied these methods to data from tumors of ovarian cancer patients. Among 99 patients, 45 tumors were informative for analysis and showed evidence of XCI skewed toward a particular parental chromosome. For 397 X‐linked genes, we observed tumor XCI patterns largely consistent with previously identified consensus states based on multiple normal tissue types. However, 37 genes differed in XCI state between ovarian tumors and the consensus state; 17 genes aberrantly escaped XCI in ovarian tumors (including many oncogenes), whereas 20 genes were unexpectedly inactivated in ovarian tumors (including many tumor suppressor genes). These results provide evidence of the importance of XCI in ovarian cancer and demonstrate the utility of our two‐stage analysis.