Fogarty, Zachary C.

Abstract Background Multiplexed immunofluorescence (MxIF) enables high-dimensional immune cell phenotyping and detailed characterization of the tumor immune microenvironment (TIME), but complex cyclic staining can introduce signal degradation and uncertainty in cell-level measurements. In routine practice, pathologists rely on hematoxylin and eosin (H&E) staining as the primary reference for tissue morphology and frequently cross-reference H&E to interpret MxIF findings. Methods We present a framework for aligning H&E and MxIF images at the tissue microarray (TMA) core level to support cross-modal analysis. Cell nucleus detections from each modality are used as anchor points, formulating the task as a point-set alignment problem. Coherent Point Drift (CPD) is applied for global alignment, followed by graph-matching–based refinement to improve local consistency. Evaluation on ovarian TMAs demonstrates consistent alignment performance for both restained and serial sections. We further explore the use of restained H&E as a reference for MxIF interpretation and investigate virtual H&E generation from MxIF data as a complementary approach when restained H&E is unavailable. Results The aligned images enable quantitative assessment of cross-modal feature concordance between MxIF-derived measurements and H&E morphology. The generated virtual H&E demonstrated similar cell population compared to the real H&E. Conclusions Cell-centric alignment can facilitate integrative multimodal histopathology analysis for TMA cores in spatial co-localization, cell feature validation, and virtual H&E generation.

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.

Rare familial spontaneous coronary artery dissection (SCAD) kindreds implicate genetic disease predisposition and provide a unique opportunity for candidate gene discovery. Whole-genome sequencing was performed in fifteen probands with non-syndromic SCAD who had a relative with SCAD, eight of whom had a second relative with extra-coronary arteriopathy. Co-segregating variants and associated genes were prioritized by quantitative variant, gene, and disease-level metrics. Curated public databases were queried for functional relationships among encoded proteins. Fifty-four heterozygous coding variants in thirteen families co-segregated with disease and fulfilled primary filters of rarity, gene variation constraint, and predicted-deleterious protein effect. Secondary filters yielded 11 prioritized candidate genes in 12 families, with high arterial tissue expression (n = 7), high-confidence protein-level interactions with genes associated with SCAD previously (n = 10), and/or previous associations with connective tissue disorders and aortopathies (n = 3) or other vascular phenotypes in mice or humans (n = 11). High-confidence associations were identified among 10 familial SCAD candidate-gene-encoded proteins. A collagen-encoding gene was identified in five families, two with distinct variants in COL4A2. Familial SCAD is genetically heterogeneous, yet perturbations of extracellular matrix, cytoskeletal, and cell–cell adhesion proteins implicate common disease-susceptibility pathways. Incomplete penetrance and variable expression suggest genetic or environmental modifiers.

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.