Multiomic analysis of uterine leiomyomas in self-described Black and White women: molecular insights into health disparities.

BACKGROUND: Black women are at an increased risk of developing uterine leiomyomas and experiencing worse disease prognosis than White women. Epidemiologic and molecular factors have been identified as underlying these disparities, but there remains a paucity of deep, multiomic analysis investigating molecular differences in uterine leiomyomas from Black and White patients. OBJECTIVE: To identify molecular alterations within uterine leiomyoma tissues correlating with patient race by multiomic analyses of uterine leiomyomas collected from cohorts of Black and White women. STUDY DESIGN: We performed multiomic analysis of uterine leiomyomas from Black (42) and White (47) women undergoing hysterectomy for symptomatic uterine leiomyomata. In addition, our analysis included the application of orthogonal methods to evaluate fibroid biomechanical properties, such as second harmonic generation microscopy, uniaxial compression testing, and shear-wave ultrasonography analyses. RESULTS: We found a greater proportion of MED12 mutant uterine leiomyomas from Black women (>35% increase; Mann-Whitney U, P<.001). MED12 mutant tumors exhibited an elevated abundance of extracellular matrix proteins, including several collagen isoforms, involved in the regulation of the core matrisome. Histologic analysis of tissue fibrosis using trichrome staining and secondary harmonic generation microscopy confirmed that MED12 mutant tumors are more fibrotic than MED12 wild-type tumors. Using shear-wave ultrasonography in a prospectively collected cohort, Black patients had fibroids that were firmer than White patients, even when similar in size. In addition, these analyses uncovered ancestry-linked expression quantitative trait loci with altered allele frequencies in African and European populations correlating with differential abundance of several proteins in uterine leiomyomas independently of MED12 mutation status, including tetracoidpeptide repeat protein 38. CONCLUSION: Our study shows that Black women have a higher prevalence of uterine leiomyomas harboring mutations in MED12 and that this mutational status correlates with increased tissue fibrosis compared with wild-type uterine leiomyomas. Our study provides insights into molecular alterations correlating with racial disparities in uterine leiomyomas and improves our understanding of the molecular etiology underlying uterine leiomyoma development within these populations.

Identification of germline population variants misclassified as cancer-associated somatic variants.

Introduction: Databases used for clinical interpretation in oncology rely on genetic data derived primarily from patients of European ancestry, leading to biases in cancer genetics research and clinical practice. One practical issue that arises in this context is the potential misclassification of multi-ancestral population variants as tumor-associated because they are not represented in reference genomes against which tumor sequencing data is aligned. Methods: To systematically find misclassified variants, we compared somatic variants in census genes from the Catalogue of Somatic Mutations in Cancer (COSMIC) V99 with multi-ancestral population variants from the Genome Aggregation Databases' Linkage Disequilibrium (GnomAD). By comparing genomic coordinates, reference, and alternate alleles, we could identify misclassified variants in genes associated with cancer. Results: We found 192 of 208 genes in COSMIC's cancer-associated census genes (92.31%) to be associated with variant misclassifications. Among the 1,906,732 variants in COSMIC, 6,957 variants (0.36%) aligned with normal population variants in GnomAD, concerning for misclassification. The African / African American ancestral population included the greatest number of misclassified variants and also had the greatest number of unique misclassified variants. Conclusion: The direct, systematic comparison of variants from COSMIC for co-occurrence in GnomAD supports a more accurate interpretation of tumor sequencing data and reduces bias related to genomic ancestry.

Immune response profiles from humans experimentally exposed to Phlebotomus duboscqi bites.

Introduction: Cutaneous leishmaniasis is a neglected vector-borne parasitic disease prevalent in 92 countries with approximately one million new infections annually. Interactions between vector saliva and the human host alter the response to infection and outcome of disease. Methods: To characterize the human immunological responses developed against saliva of Phlebotomus duboscqi, a Leishmania major (L. major) vector, we repeatedly exposed the arms of 14 healthy U.S volunteers to uninfected P. duboscqi bites. Blood was collected a week after each exposure and used to assess total IgG antibodies against the proteins of P. duboscqi salivary gland homogenate (SGH) and the levels of IFN-gamma and IL-10 from peripheral blood mononuclear cells (PBMCs) stimulated with SGH or recombinant sand fly proteins. We analyzed skin punch biopsies of the human volunteer arms from the insect bite site and control skin site after multiple P. duboscqi exposures (four volunteers) using immunohistochemical staining. Results: A variety of immediate insect bite skin reactions were observed. Late skin reactions to insect bites were characterized by macular hyperpigmentation and/or erythematous papules. Hematoxylin and eosin staining showed moderate mononuclear skin infiltrate with eosinophils in those challenged recently (within 2 months), eosinophils were not seen in biopsies with recall challenge (6 month post bites). An increase in plasma antigen-specific IgG responses to SGH was observed over time. Western Blot results showed strong plasma reactivity to five P. duboscqi salivary proteins. Importantly, volunteers developed a cellular immunity characterized by the secretion of IFN-gamma upon PBMC stimulation with P. duboscqi SGH and recombinant antigens. Discussion: Our results demonstrate that humans mounted a local and systemic immune response against P. duboscqi salivary proteins. Specifically, PduM02/SP15-like and PduM73/adenosine deaminase recombinant salivary proteins triggered a Th1 type immune response that might be considered in future development of a potential Leishmania vaccine.

Genome-Wide Analysis of Exertional Rhabdomyolysis in Sickle Cell Trait Positive African Americans.

Sickle cell trait (SCT), although generally a benign carrier state of hemoglobin S (HbAS), is a risk factor for exertional rhabdomyolysis (ERM), a rare but potentially fatal consequence of highly intense physical exercise, particularly among active-duty military personnel and high-performance athletes. The association between SCT and ERM is poorly understood. The objective of this study was to elucidate the genetic basis of ERM in an SCT-positive African American cohort. SCT-positive African Americans with a personal history of ERM (cases, n = 30) and without history of ERM (controls, n = 53) were enrolled in this study. Whole-genome sequencing was performed on DNA samples isolated from peripheral white blood cells. Participants' demographic, behavioral, and medical history information was obtained. An additional 131 controls were extracted from SCT-positive subjects of African descent from the 1000 Genomes Project. SCT carriers with ERM were characterized by myotoxicity features, significant muscle involvement dominated by muscle weakness, and severe pain and substantial increase in serum creatine kinase, with a mean value of 50,480 U/L. A distinctive feature of the SCT individuals with ERM was exertional collapse, which was reported in 53.3% of the cases in the study cohort. An important factor for the development of ERM was the duration and frequency of strenuous physical activity in the cases compared to the controls. Whole-genome sequencing identified 79,696 protein-coding variants. Genome-wide association analysis revealed that the p.C477R, rs115958260 variant in the SLC44A3 gene was significantly associated with ERM event in SCT-positive African Americans. The study results suggest that a combination of vigorous exercise and a genetic predisposing factor is involved in ERM.

ProteoMixture: A cell type deconvolution tool for bulk tissue proteomic data.

Numerous multi-omic investigations of cancer tissue have documented varying and poor pairwise transcript:protein quantitative correlations, and most deconvolution tools aiming to predict cell type proportions (cell admixture) have been developed and credentialed using transcript-level data alone. To estimate cell admixture using protein abundance data, we analyzed proteome and transcriptome data generated from contrived admixtures of tumor, stroma, and immune cell models or those selectively harvested from the tissue microenvironment by laser microdissection from high grade serous ovarian cancer (HGSOC) tumors. Co-quantified transcripts and proteins performed similarly to estimate stroma and immune cell admixture (r ≥ 0.63) in two commonly used deconvolution algorithms, ESTIMATE or ConsensusTME. We further developed and optimized protein-based signatures estimating cell admixture proportions and benchmarked these using bulk tumor proteomic data from over 150 patients with HGSOC. The optimized protein signatures supporting cell type proportion estimates from bulk tissue proteomic data are available at https://lmdomics.org/ProteoMixture/.

The influence of microbial colonization on inflammatory versus pro-healing trajectories in combat extremity wounds.

A combination of improved body armor, medical transportation, and treatment has led to the increased survival of warfighters from combat extremity injuries predominantly caused by blasts in modern conflicts. Despite advances, a high rate of complications such as wound infections, wound failure, amputations, and a decreased quality of life exist. To study the molecular underpinnings of wound failure, wound tissue biopsies from combat extremity injuries had RNA extracted and sequenced. Wounds were classified by colonization (colonized vs. non-colonized) and outcome (healed vs. failed) status. Differences in gene expression were investigated between timepoints at a gene level, and longitudinally by multi-gene networks, inferred proportions of immune cells, and expression of healing-related functions. Differences between wound outcomes in colonized wounds were more apparent than in non-colonized wounds. Colonized/healed wounds appeared able to mount an adaptive immune response to infection and progress beyond the inflammatory stage of healing, while colonized/failed wounds did not. Although, both colonized and non-colonized failed wounds showed increasing inferred immune and inflammatory programs, non-colonized/failed wounds progressed beyond the inflammatory stage, suggesting different mechanisms of failure dependent on colonization status. Overall, these data reveal gene expression profile differences in healing wounds that may be utilized to improve clinical treatment paradigms.

Proteogenomic analysis of enriched HGSOC tumor epithelium identifies prognostic signatures and therapeutic vulnerabilities.

We performed a deep proteogenomic analysis of bulk tumor and laser microdissection enriched tumor cell populations from high-grade serous ovarian cancer (HGSOC) tissue specimens spanning a broad spectrum of purity. We identified patients with longer progression-free survival had increased immune-related signatures and validated proteins correlating with tumor-infiltrating lymphocytes in 65 tumors from an independent cohort of HGSOC patients, as well as with overall survival in an additional 126 HGSOC patient cohort. We identified that homologous recombination deficient (HRD) tumors are enriched in pathways associated with metabolism and oxidative phosphorylation that we validated in independent patient cohorts. We further identified that polycomb complex protein BMI-1 is elevated in HR proficient (HRP) tumors, that elevated BMI-1 correlates with poor overall survival in HRP but not HRD HGSOC patients, and that HRP HGSOC cells are uniquely sensitive to BMI-1 inhibition.

Uncommon Protein-Coding Variants Associated With Suicide Attempt in a Diverse Sample of U.S. Army Soldiers.

BACKGROUND: Suicide is a societal and public health concern of global scale. Identifying genetic risk factors for suicide attempt can characterize underlying biology and enable early interventions to prevent deaths. Recent studies have described common genetic variants for suicide-related behaviors. Here, we advance this search for genetic risk by analyzing the association between suicide attempt and uncommon variation exome-wide in a large, ancestrally diverse sample. METHODS: We sequenced whole genomes of 13,584 soldiers from the Army STARRS (Army Study to Assess Risk and Resilience in Servicemembers), including 979 individuals with a history of suicide attempt. Uncommon, nonsilent protein-coding variants were analyzed exome-wide for association with suicide attempt using gene-collapsed and single-variant analyses. RESULTS: We identified 19 genes with variants enriched in individuals with history of suicide attempt, either through gene-collapsed or single-variant analysis (Bonferroni padjusted < .05). These genes were CIB2, MLF1, HERC1, YWHAE, RCN2, VWA5B1, ATAD3A, NACA, EP400, ZNF585A, LYST, RC3H2, PSD3, STARD9, SGMS1, ACTR6, RGS7BP, DIRAS2, and KRTAP10-1. Most genes had variants across multiple genomic ancestry groups. Seventeen of these genes were expressed in healthy brain tissue, with 9 genes expressed at the highest levels in the brain versus other tissues. Brains from individuals deceased from suicide aberrantly expressed RGS7BP (padjusted = .035) in addition to nominally significant genes including YWHAE and ACTR6, all of which have reported associations with other mental disorders. CONCLUSIONS: These results advance the molecular characterization of suicide attempt behavior and support the utility of whole-genome sequencing for complementing the findings of genome-wide association studies in suicide research.