Pan-cancer profiling of tumor-infiltrating natural killer cells through transcriptional reference mapping.

The functional diversity of natural killer (NK) cell repertoires stems from differentiation, homeostatic, receptor-ligand interactions and adaptive-like responses to viral infections. In the present study, we generated a single-cell transcriptional reference map of healthy human blood- and tissue-derived NK cells, with temporal resolution and fate-specific expression of gene-regulatory networks defining NK cell differentiation. Transfer learning facilitated incorporation of tumor-infiltrating NK cell transcriptomes (39 datasets, 7 solid tumors, 427 patients) into the reference map to analyze tumor microenvironment (TME)-induced perturbations. Of the six functionally distinct NK cell states identified, a dysfunctional stressed CD56bright state susceptible to TME-induced immunosuppression and a cytotoxic TME-resistant effector CD56dim state were commonly enriched across tumor types, the ratio of which was predictive of patient outcome in malignant melanoma and osteosarcoma. This resource may inform the design of new NK cell therapies and can be extended through transfer learning to interrogate new datasets from experimental perturbations or disease conditions.

High-dimensional single-cell analysis of human natural killer cell heterogeneity.

Natural killer (NK) cells are innate lymphoid cells (ILCs) contributing to immune responses to microbes and tumors. Historically, their classification hinged on a limited array of surface protein markers. Here, we used single-cell RNA sequencing (scRNA-seq) and cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) to dissect the heterogeneity of NK cells. We identified three prominent NK cell subsets in healthy human blood: NK1, NK2 and NK3, further differentiated into six distinct subgroups. Our findings delineate the molecular characteristics, key transcription factors, biological functions, metabolic traits and cytokine responses of each subgroup. These data also suggest two separate ontogenetic origins for NK cells, leading to divergent transcriptional trajectories. Furthermore, we analyzed the distribution of NK cell subsets in the lung, tonsils and intraepithelial lymphocytes isolated from healthy individuals and in 22 tumor types. This standardized terminology aims at fostering clarity and consistency in future research, thereby improving cross-study comparisons.

HLA class I signal peptide polymorphism determines the level of CD94/NKG2-HLA-E-mediated regulation of effector cell responses.

Human leukocyte antigen (HLA)-E binds epitopes derived from HLA-A, HLA-B, HLA-C and HLA-G signal peptides (SPs) and serves as a ligand for CD94/NKG2A and CD94/NKG2C receptors expressed on natural killer and T cell subsets. We show that among 16 common classical HLA class I SP variants, only 6 can be efficiently processed to generate epitopes that enable CD94/NKG2 engagement, which we term 'functional SPs'. The single functional HLA-B SP, known as HLA-B/-21M, induced high HLA-E expression, but conferred the lowest receptor recognition. Consequently, HLA-B/-21M SP competes with other SPs for providing epitope to HLA-E and reduces overall recognition of target cells by CD94/NKG2A, calling for reassessment of previous disease models involving HLA-B/-21M. Genetic population data indicate a positive correlation between frequencies of functional SPs in humans and corresponding cytomegalovirus mimics, suggesting a means for viral escape from host responses. The systematic, quantitative approach described herein will facilitate development of prediction algorithms for accurately measuring the impact of CD94/NKG2-HLA-E interactions in disease resistance/susceptibility.

TREM2 macrophages drive NK cell paucity and dysfunction in lung cancer.

Natural killer (NK) cells are commonly reduced in human tumors, enabling many to evade surveillance. Here, we sought to identify cues that alter NK cell activity in tumors. We found that, in human lung cancer, the presence of NK cells inversely correlated with that of monocyte-derived macrophages (mo-macs). In a murine model of lung adenocarcinoma, we show that engulfment of tumor debris by mo-macs triggers a pro-tumorigenic program governed by triggering receptor expressed on myeloid cells 2 (TREM2). Genetic deletion of Trem2 rescued NK cell accumulation and enabled an NK cell-mediated regression of lung tumors. TREM2+ mo-macs reduced NK cell activity by modulating interleukin (IL)-18/IL-18BP decoy interactions and IL-15 production. Notably, TREM2 blockade synergized with an NK cell-activating agent to further inhibit tumor growth. Altogether, our findings identify a new axis, in which TREM2+ mo-macs suppress NK cell accumulation and cytolytic activity. Dual targeting of macrophages and NK cells represents a new strategy to boost antitumor immunity.

The CD94/NKG2A inhibitory receptor educates uterine NK cells to optimize pregnancy outcomes in humans and mice.

The conserved CD94/NKG2A inhibitory receptor is expressed by nearly all human and ∼50% of mouse uterine natural killer (uNK) cells. Binding human HLA-E and mouse Qa-1, NKG2A drives NK cell education, a process of unknown physiological importance influenced by HLA-B alleles. Here, we show that NKG2A genetic ablation in dams mated with wild-type males caused suboptimal maternal vascular responses in pregnancy, accompanied by perturbed placental gene expression, reduced fetal weight, greater rates of smaller fetuses with asymmetric growth, and abnormal brain development. These are features of the human syndrome pre-eclampsia. In a genome-wide association study of 7,219 pre-eclampsia cases, we found a 7% greater relative risk associated with the maternal HLA-B allele that does not favor NKG2A education. These results show that the maternal HLA-B→HLA-E→NKG2A pathway contributes to healthy pregnancy and may have repercussions on offspring health, thus establishing the physiological relevance for NK cell education. VIDEO ABSTRACT.

Immunology of COVID-19: Current State of the Science.

The coronavirus disease 2019 (COVID-19) pandemic, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has affected millions of people worldwide, igniting an unprecedented effort from the scientific community to understand the biological underpinning of COVID19 pathophysiology. In this Review, we summarize the current state of knowledge of innate and adaptive immune responses elicited by SARS-CoV-2 infection and the immunological pathways that likely contribute to disease severity and death. We also discuss the rationale and clinical outcome of current therapeutic strategies as well as prospective clinical trials to prevent or treat SARS-CoV-2 infection.

IL-15 synergizes with CD40 agonist antibodies to induce durable immunity against bladder cancer.

CD40 is a central costimulatory receptor implicated in productive antitumor immune responses across multiple cancers, including bladder cancer. Despite strong preclinical rationale, systemic administration of therapeutic agonistic antibodies targeting the CD40 pathway has demonstrated dose-limiting toxicities with minimal clinical activity, emphasizing an important need for optimized CD40-targeted approaches, including rational combination therapy strategies. Here, we describe a role for the endogenous IL-15 pathway in contributing to the therapeutic activity of CD40 agonism in orthotopic bladder tumors, with upregulation of transpresented IL-15/IL-15Rα surface complexes, particularly by cross-presenting conventional type 1 DCs (Dendritic Cells), and associated enrichment of activated CD8 T cells. In bladder cancer patient samples, we identify DCs as the primary source of IL-15, although they lack high levels of IL-15Rα at baseline. Using humanized immunocompetent orthotopic bladder tumor models, we demonstrate the ability to therapeutically augment this interaction through combined treatment with anti-CD40 agonist antibodies and exogenous IL-15, including the fully-human Fc-optimized antibody 2141-V11 currently in clinical development for the treatment of bladder cancer. Collectively, these data reveal an important role for IL-15 in mediating antitumor CD40 agonist responses in bladder cancer and provide key proof-of-concept for combined use of Fc-optimized anti-CD40 agonist antibodies and agents targeting the IL-15 pathway. These data support expansion of ongoing clinical studies evaluating anti-CD40 agonist antibodies and IL-15-based approaches to develop combinations of these promising therapeutics for the treatment of patients with bladder cancer.

Urine supernatant reveals a signature that predicts survival in clear-cell renal cell carcinoma.

OBJECTIVE: To profile the cell-free urine supernatant and plasma of a small cohort of clear-cell renal cell carcinoma (ccRCC) patients by measuring the relative concentrations of 92 proteins related to inflammation. Using The Cancer Genome Atlas (TCGA), we then performed a targeted mRNA analysis of genes encoding the above proteins and defined their effects on overall survival (OS). SUBJECTS/PATIENTS AND METHODS: Samples were collected prospectively from ccRCC patients. A multiplex proximity extension assay was used to measure the concentrations of 92 inflammation-related proteins in cell-free urine supernatants and plasma. Transcriptomic and clinical information from ccRCC patients was obtained from TCGA. Unsupervised clustering and differential protein expression analyses were performed on protein concentration data. Targeted mRNA analysis on genes encoding significant differentially expressed proteins was performed using TCGA. Backward stepwise regression analyses were used to build a nomogram. The performance of the nomogram and clinical benefit was assessed by discrimination and calibration, and a decision curve analysis, respectively. RESULTS: Unsupervised clustering analysis revealed inflammatory signatures in the cell-free urine supernatant of ccRCC patients. Backward stepwise regressions using TCGA data identified transcriptomic risk factors and risk groups associated with OS. A nomogram to predict 2-year and 5-year OS was developed using these risk factors. The decision curve analysis showed that our model was associated with a net benefit improvement compared to the treat-all/none strategies. CONCLUSION: We defined four novel biomarkers using proteomic and transcriptomic data that distinguish severity of prognosis in ccRCC. We showed that these biomarkers can be used in a model to predict 2-year and 5-year OS in ccRCC across different tumour stages. This type of analysis, if validated in the future, provides non-invasive prognostic information that could inform either management or surveillance strategies for patients.

MiDAS-Meaningful Immunogenetic Data at Scale.

Human immunogenetic variation in the form of HLA and KIR types has been shown to be strongly associated with a multitude of immune-related phenotypes. However, association studies involving immunogenetic loci most commonly involve simple analyses of classical HLA allelic diversity, resulting in limitations regarding the interpretability and reproducibility of results. We here present MiDAS, a comprehensive R package for immunogenetic data transformation and statistical analysis. MiDAS recodes input data in the form of HLA alleles and KIR types into biologically meaningful variables, allowing HLA amino acid fine mapping, analyses of HLA evolutionary divergence as well as experimentally validated HLA-KIR interactions. Further, MiDAS enables comprehensive statistical association analysis workflows with phenotypes of diverse measurement scales. MiDAS thus closes the gap between the inference of immunogenetic variation and its efficient utilization to make relevant discoveries related to immune and disease biology. It is freely available under a MIT license.

Once-Weekly Oral Dosing of MK-8591 Protects Male Rhesus Macaques From Intrarectal Challenge With SHIV109CP3.

BACKGROUND: MK-8591 (4'-ethynyl-2-fluoro-2'-deoxyadenosine [EFdA]) is a novel reverse transcriptase-translocation inhibitor. METHODS: We assessed MK-8591 as preexposure prophylaxis in the rhesus macaque model of intrarectal challenge with simian/human immunodeficiency virus (SHIV). In study 1, 8 rhesus macaques received 3.9 mg/kg of MK-8591 orally on day 0 and once weekly for the next 14 weeks. Eight controls were treated with vehicle. All rhesus macaques were challenged with SHIV109CP3 on day 6 and weekly for up to 12 challenges or until infection was confirmed. The dose of MK-8591 was reduced to 1.3 and 0.43 mg/kg/week in study 2 and further to 0.1 and 0.025 mg/kg/week in study 3. In studies 2 and 3, each dose was given up to 6 times once weekly, and animals were challenged 4 times once weekly with SHIV109CP3. RESULTS: Control macaques were infected after a median of 1 challenge (range, 1-4 challenges). All treated animals in studies 1 and 2 were protected, consistent with a 41.5-fold lower risk of infection (P < .0001, by the log-rank test). In study 3, at a 0.1-mg/kg dose, 2 rhesus macaques became infected, consistent with a 7.2-fold lower risk of infection (P = .0003, by the log-rank test). The 0.025-mg/kg dose offered no protection. CONCLUSIONS: These data support MK-8591's potential as a preexposure prophylaxis agent.

A Distinctive Cytoplasmic Tail Contributes to Low Surface Expression and Intracellular Retention of the Patr-AL MHC Class I Molecule.

Chimpanzees have orthologs of the six fixed, functional human MHC class I genes. But, in addition, the chimpanzee has a seventh functional gene, Patr-AL, which is not polymorphic but contributes substantially to population diversity by its presence on only 50% of MHC haplotypes. The ancestral AL gene emerged long before the separation of human and chimpanzee ancestors and then subsequently and specifically lost function during human evolution, but was maintained in chimpanzees. Patr-AL is an alloantigen that participates in negative and positive selection of the T cell repertoire. The three-dimensional structure and the peptide-binding repertoire of Patr-AL and HLA-A*02 are surprisingly similar. In contrast, the expression of these two molecules is very different, as shown using specific mAbs and polyclonal Abs made against Patr-AL. Peripheral blood cells and B cell lines express low levels of Patr-AL at the cell surface. Higher levels are seen for 221-cell transfectants expressing Patr-AL, but in these cells a large majority of Patr-AL molecules are retained in the early compartments of the secretory pathway: mainly the endoplasmic reticulum, but also cis-Golgi. Replacing the cytoplasmic tail of Patr-AL with that of HLA-A*02 increased the cell-surface expression of Patr-AL substantially. Four substitutions distinguish the Patr-AL and HLA-A*02 cytoplasmic tails. Systematic mutagenesis showed that each substitution contributes changes in cell-surface expression. The combination of residues present in Patr-AL appears unique, but each individual residue is present in other primate MHC class I molecules, notably MHC-E, the most ancient of the functional human MHC class I molecules.

Regulation of Adaptive NK Cells and CD8 T Cells by HLA-C Correlates with Allogeneic Hematopoietic Cell Transplantation and with Cytomegalovirus Reactivation.

Mass cytometry was used to investigate the effect of CMV reactivation on lymphocyte reconstitution in hematopoietic cell transplant patients. For eight transplant recipients (four CMV negative and four CMV positive), we studied PBMCs obtained 6 mo after unrelated donor hematopoietic cell transplantation (HCT). Forty cell-surface markers, distinguishing all major leukocyte populations in PBMC, were analyzed with mass cytometry. This group included 34 NK cell markers. Compared with healthy controls, transplant recipients had higher HLA-C expression on CD56(-)CD16(+) NK cells, B cells, CD33(bright) myeloid cells, and CD4CD8 T cells. The increase in HLA-C expression was greater for CMV-positive HCT recipients than for CMV negative recipients. Present in CMV-positive HCT recipients, but not in CMV-negative HCT recipients or controls, is a population of killer cell Ig-like receptor (KIR)-expressing CD8 T cells not previously described. These CD8 T cells coexpress CD56, CD57, and NKG2C. The HCT recipients also have a population of CD57(+)NKG2A(+) NK cells that preferentially express KIR2DL1. An inverse correlation was observed between the frequencies of CD57(+)NKG2C(+) NK cells and CD57(+)NKG2A(+) NK cells. Although CD57(+)NKG2A(+) NK cells are less abundant in CMV-positive recipients, their phenotype is of a more activated cell than the CD57(+)NKG2A(+) NK cells of controls and CMV-negative HCT recipients. These data demonstrate that HCT and CMV reactivation are associated with an increased expression of HLA-C. This could influence NK cell education during lymphocyte reconstitution. The increased inhibitory KIR expression by proliferating CMV-specific CD8 T cells suggests regulatory interactions between HLA-C and KIR might promote Graft-versus-Leukemia effects following transplantation.

Coordinated regulation of NK receptor expression in the maturing human immune system.

NK cells are responsible for recognizing and killing transformed, stressed, and infected cells. They recognize a set of non-Ag-specific features termed "altered self" through combinatorial signals from activating and inhibitory receptors. These NKRs are also expressed on CD4(+) and CD8(+) T cells, B cells, and monocytes, although a comprehensive inventory of NKR expression patterns across leukocyte lineages has never been performed. Using mass cytometry, we found that NKR expression patterns distinguish cell lineages in human peripheral blood. In individuals with high levels of CD57, indicative of a mature immune repertoire, NKRs are more likely to be expressed on non-NK cells, especially CD8(+) T cells. Mature NK and CD8(+) T cell populations show increased diversity of NKR surface expression patterns, but with distinct determinants: mature NK cells acquire primarily inhibitory receptors, whereas CD8(+) T cells attain a specific subset of both activating and inhibitory receptors, potentially imbuing them with a distinct functional role. Concurrently, monocytes show decreased expression of the generalized inhibitory receptor leukocyte Ig-like receptor subfamily b member 1, consistent with an increased activation threshold. Therefore, NKR expression is coordinately regulated as the immune system matures, resulting in the transfer of "altered self" recognition potential among leukocyte lineages. This likely reduces Ag specificity in the mature human immune system, and implies that vaccines and therapeutics that engage both its innate and adaptive branches may be more effective in the settings of aging and chronic infection.

Adaptive Natural Killer Cell and Killer Cell Immunoglobulin-Like Receptor-Expressing T Cell Responses are Induced by Cytomegalovirus and Are Associated with Protection against Cytomegalovirus Reactivation after Allogeneic Donor Hematopoietic Cell Transplantation.

Cytomegalovirus (CMV) reactivates in >30% of CMV-seropositive patients after allogeneic hematopoietic cell transplantation (HCT). Previously, we reported an increase of natural killer (NK) cells expressing NKG2C, CD57, and inhibitory killer cell immunoglobulin-like receptors (KIRs) in response to CMV reactivation after HCT. These NK cells persist after the resolution of infection and display "adaptive" or memory properties. Despite these findings, the differential impact of persistent/inactive versus reactivated CMV on NK versus T cell maturation after HCT from different graft sources has not been defined. We compared the phenotype of NK and T cells from 292 recipients of allogeneic sibling (n = 118) or umbilical cord blood (UCB; n = 174) grafts based on recipient pretransplantation CMV serostatus and post-HCT CMV reactivation. This cohort was utilized to evaluate CMV-dependent increases in KIR-expressing NK cells exhibiting an adaptive phenotype (NKG2C(+)CD57(+)). Compared with CMV-seronegative recipients, those who reactivated CMV had the highest adaptive cell frequencies, whereas intermediate frequencies were observed in CMV-seropositive recipients harboring persistent/nonreplicating CMV. The same effect was observed in T cells and CD56(+) T cells. These adaptive lymphocyte subsets were increased in CMV-seropositive recipients of sibling but not UCB grafts and were correlated with lower rates of CMV reactivation (sibling 33% versus UCB 51%; P < .01). These data suggest that persistent/nonreplicating recipient CMV induces rapid production of adaptive NK and T cells from mature cells from sibling but not UCB grafts. These adaptive lymphocytes are associated with protection from CMV reactivation.

Antigen-specific IL-2 secretion correlates with NK cell responses after immunization of Tanzanian children with the RTS,S/AS01 malaria vaccine.

RTS,S/AS01, a vaccine targeting pre-erythrocytic stages of Plasmodium falciparum, is undergoing clinical trials. We report an analysis of cellular immune response to component Ags of RTS,S-hepatitis B surface Ag (HBs) and P. falciparum circumsporozoite (CS) protein-among Tanzanian children in a phase IIb RTS,S/AS01(E) trial. RTS,S/AS01 (E) vaccinees make stronger T cell IFN-γ, CD69, and CD25 responses to HBs peptides than do controls, indicating that RTS,S boosts pre-existing HBs responses. T cell CD69 and CD25 responses to CS and CS-specific secreted IL-2 were augmented by RTS,S vaccination. Importantly, more than 50% of peptide-induced IFN-γ(+) lymphocytes were NK cells, and the magnitude of the NK cell CD69 response to HBs peptides correlated with secreted IL-2 concentration. CD69 and CD25 expression and IL-2 secretion may represent sensitive markers of RTS,S-induced, CS-specific T cells. The potential for T cell-derived IL-2 to augment NK cell activation in RTS,S-vaccinated individuals, and the relevance of this for protection, needs to be explored further.

C1QA is an invariant biomarker for tissue macrophages.

Macrophages play a pivotal role in immune responses, particularly in the context of combating microbial threats within tissues. The identification of reliable biomarkers associated with macrophage function is essential for understanding their diverse roles in host defense. This study investigates the potential of C1QA as an invariant biomarker for tissue macrophages, focusing on its correlation with the anti-microbial pathway. C1QA, a component of the complement system, has been previously implicated in various immune functions. Our research delves into the specific association of C1QA with tissue-resident macrophages and its implications in the context of anti-microbial responses. Through comprehensive systems biology and Boolean analysis of gene expression, we aim to establish C1QA as a consistent and reliable marker for identifying tissue macrophages. Furthermore, we explore the functional significance of C1QA in the anti-microbial pathway. This research seeks to provide valuable insights into the molecular mechanisms underlying the anti-microbial functions of tissue macrophages, with C1QA emerging as a potential key player in this intricate regulatory network. Understanding the relationship between C1QA, tissue macrophages, and the anti-microbial pathway could pave the way for the development of targeted therapeutic strategies aimed at enhancing the host's ability to combat infections. Ultimately, our findings contribute to the expanding knowledge of macrophage biology and may have implications for the diagnosis and treatment of infectious diseases.

A monoclonal antibody that recognizes a unique 13-residue epitope in the cytoplasmic tail of HLA-E.

The Class I MHC molecule (MHC-I) HLA-E presents peptides that are derived from the signal sequences, either those of other MHC-I products, or of viral type I membrane glycoproteins. Monoclonal antibodies with proven specificity for HLA-E, and with no cross-reactions with other MHC-I products, have yet to be described. To obtain anti-HLA-E-specific antibodies suitable for a range of applications, we generated monoclonal antibodies against a unique feature of HLA-E: its cytoplasmic tail. We created an immunogen by performing an enzymatically catalyzed transpeptidation reaction to obtain a fusion of the cytoplasmic tail of HLA-E with a nanobody that recognizes murine Class II MHC (MHC-II) products. We obtained a mouse monoclonal antibody that recognizes a 13-residue stretch in the HLA-E cytoplasmic tail. We cloned the genes that encode this antibody in expression vectors to place an LPETG sortase recognition motif at the C-terminus of the heavy and light chains. This arrangement allows the site-specific installation of fluorophores or biotin at these C-termini. The resulting immunoglobulin preparations, labeled with 4 equivalents of a fluorescent or biotinylated payload of choice, can then be used for direct immunofluorescence or detection of the tag by fluorescence or by streptavidin-based methods. We also show that the 13-residue sequence can serve as an epitope tag, independent of the site of its placement within a protein's sequence. The antibody can be used diagnostically to stain for HLA-E on patient tumor samples, it can be used as an antibody-epitope tag for extracellular proteins, and it enables research into the unique role of the cytoplasmic tail of HLA-E.

Urine scRNAseq reveals new insights into the bladder tumor immune microenvironment.

Due to bladder tumors' contact with urine, urine-derived cells (UDCs) may serve as a surrogate for monitoring the tumor microenvironment (TME) in bladder cancer (BC). However, the composition of UDCs and the extent to which they mirror the tumor remain poorly characterized. We generated the first single-cell RNA-sequencing of BC patient UDCs with matched tumor and peripheral blood mononuclear cells (PBMC). BC urine was more cellular than healthy donor (HD) urine, containing multiple immune populations including myeloid cells, CD4+ and CD8+ T cells, natural killer (NK) cells, B cells, and dendritic cells (DCs) in addition to tumor and stromal cells. Immune UDCs were transcriptionally more similar to tumor than blood. UDCs encompassed cytotoxic and activated CD4+ T cells, exhausted and tissue-resident memory CD8+ T cells, macrophages, germinal-center-like B cells, tissue-resident and adaptive NK cells, and regulatory DCs found in tumor but lacking or absent in blood. Our findings suggest BC UDCs may be surrogates for the TME and serve as therapeutic biomarkers.

Single-cell transcriptomic-informed deconvolution of bulk data identifies immune checkpoint blockade resistance in urothelial cancer.

Interactions within the tumor microenvironment (TME) significantly influence tumor progression and treatment responses. While single-cell RNA sequencing (scRNA-seq) and spatial genomics facilitate TME exploration, many clinical cohorts are assessed at the bulk tissue level. Integrating scRNA-seq and bulk tissue RNA-seq data through computational deconvolution is essential for obtaining clinically relevant insights. Our method, ProM, enables the examination of major and minor cell types. Through evaluation against existing methods using paired single-cell and bulk RNA sequencing of human urothelial cancer (UC) samples, ProM demonstrates superiority. Application to UC cohorts treated with immune checkpoint inhibitors reveals pre-treatment cellular features associated with poor outcomes, such as elevated SPP1 expression in macrophage/monocytes (MM). Our deconvolution method and paired single-cell and bulk tissue RNA-seq dataset contribute novel insights into TME heterogeneity and resistance to immune checkpoint blockade.

Characterization of renal masses with MRI-based radiomics: assessment of inter-package and inter-observer reproducibility in a prospective pilot study.

OBJECTIVES: To evaluate radiomics features' reproducibility using inter-package/inter-observer measurement analysis in renal masses (RMs) based on MRI and to employ machine learning (ML) models for RM characterization. METHODS: 32 Patients (23M/9F; age 61.8 ± 10.6 years) with RMs (25 renal cell carcinomas (RCC)/7 benign masses; mean size, 3.43 ± 1.73 cm) undergoing resection were prospectively recruited. All patients underwent 1.5 T MRI with T2-weighted (T2-WI), diffusion-weighted (DWI)/apparent diffusion coefficient (ADC), and pre-/post-contrast-enhanced T1-weighted imaging (T1-WI). RMs were manually segmented using volume of interest (VOI) on T2-WI, DWI/ADC, and T1-WI pre-/post-contrast imaging (1-min, 3-min post-injection) by two independent observers using two radiomics software packages for inter-package and inter-observer assessments of shape/histogram/texture features common to both packages (104 features; n = 26 patients). Intra-class correlation coefficients (ICCs) were calculated to assess inter-observer and inter-package reproducibility of radiomics measurements [good (ICC ≥ 0.8)/moderate (ICC = 0.5-0.8)/poor (ICC < 0.5)]. ML models were employed using reproducible features (between observers and packages, ICC > 0.8) to distinguish RCC from benign RM. RESULTS: Inter-package comparisons demonstrated that radiomics features from T1-WI-post-contrast had the highest proportion of good/moderate ICCs (54.8-58.6% for T1-WI-1 min), while most features extracted from T2-WI, T1-WI-pre-contrast, and ADC exhibited poor ICCs. Inter-observer comparisons found that radiomics measurements from T1-WI pre/post-contrast and T2-WI had the greatest proportion of features with good/moderate ICCs (95.3-99.1% T1-WI-post-contrast 1-min), while ADC measurements yielded mostly poor ICCs. ML models generated an AUC of 0.71 [95% confidence interval = 0.67-0.75] for diagnosis of RCC vs. benign RM. CONCLUSION: Radiomics features extracted from T1-WI-post-contrast demonstrated greater inter-package and inter-observer reproducibility compared to ADC, with fair accuracy for distinguishing RCC from benign RM. CLINICAL RELEVANCE: Knowledge of reproducibility of MRI radiomics features obtained on renal masses will aid in future study design and may enhance the diagnostic utility of radiomics models for renal mass characterization.