Balancing selection and the functional effects of shared polymorphism in cryptic Daphnia species.

The patterns of genetic variation within and between related taxa represent the genetic history of a species. Shared polymorphisms, loci with identical alleles across species, are of unique interest as they may represent cases of ancient selection maintaining functional variation post-speciation. In this study, we investigate the abundance of shared polymorphism in the Daphnia pulex species complex. We test whether shared mutations are consistent with the action of balancing selection or alternative hypotheses such as hybridization, incomplete lineage sorting, or convergent evolution. We analyzed over 2,000 genomes from North American and European D. pulex and several outgroup species to examine the prevalence and distribution of shared alleles between the focal species pair, North American and European D. pulex. We show that while North American and European D. pulex diverged over ten million years ago, they retained tens of thousands of shared alleles. We found that the number of shared polymorphisms between North American and European D. pulex cannot be explained by hybridization or incomplete lineage sorting alone. Instead, we show that most shared polymorphisms could be the product of convergent evolution, that a limited number appear to be old trans-specific polymorphisms, and that balancing selection is affecting young and ancient mutations alike. Finally, we provide evidence that a blue wavelength opsin gene with trans-specific polymorphisms has functional effects on behavior and fitness in the wild. Ultimately, our findings provide insights into the genetic basis of adaptation and the maintenance of genetic diversity between species.

DNA fragility at topologically associated domain boundaries is promoted by alternative DNA secondary structure and topoisomerase II activity.

CCCTC-binding factor (CTCF) binding sites are hotspots of genome instability. Although many factors have been associated with CTCF binding site fragility, no study has integrated all fragility-related factors to understand the mechanism(s) of how they work together. Using an unbiased, genome-wide approach, we found that DNA double-strand breaks (DSBs) are enriched at strong, but not weak, CTCF binding sites in five human cell types. Energetically favorable alternative DNA secondary structures underlie strong CTCF binding sites. These structures coincided with the location of topoisomerase II (TOP2) cleavage complex, suggesting that DNA secondary structure acts as a recognition sequence for TOP2 binding and cleavage at CTCF binding sites. Furthermore, CTCF knockdown significantly increased DSBs at strong CTCF binding sites and at CTCF sites that are located at topologically associated domain (TAD) boundaries. TAD boundary-associated CTCF sites that lost CTCF upon knockdown displayed increased DSBs when compared to the gained sites, and those lost sites are overrepresented with G-quadruplexes, suggesting that the structures act as boundary insulators in the absence of CTCF, and contribute to increased DSBs. These results model how alternative DNA secondary structures facilitate recruitment of TOP2 to CTCF binding sites, providing mechanistic insight into DNA fragility at CTCF binding sites.

Role of ERα and aromatase in juvenile gigantomastia.

CONTEXT: Approximately 150 patients with juvenile gigantomastia have been reported in the literature but the underlying biologic mechanisms remain unknown. OBJECTIVE: Conduct extensive clinical, biochemical, immunochemical, and genetic studies in three patients with juvenile gigantomastia to determine causative biologic factors. DESIGN: We examined clinical effects of estrogen by blockading estrogen synthesis or its action. Breast tissue aromatase expression and activity were quantitated in one patient and five controls. Other biochemical markers including ERα, Cyclin D1 and E, p-RB, p-MAPK, p-AKT, BCL-2, EGF-R, IGF-IR ß, and p-EGFR. Immunohistochemical analyses for aromatase, ERα and ß, PGR, Ki67, sulfotransferase, estrone sulfatase, and 17ßHD were performed in all three patients. The entire genome of the mother, father, and patient in the three families was sequenced. RESULTS: Blockade of estrogen synthesis or action in patients resulted in demonstrable clinical effects. Biochemical studies on fresh frozen tissue revealed no differences between patients and controls, presumably due to tissue dilution from the large proportion of stroma. However, Immunohistochemically analysis of ductal breast cells in the three patients revealed a high percent of ERα (i.e. 64.1 ± 7.8% vs normal women 9.6% -range 2.3-15%); aromatase score of 4 (76-100% of cells positive) versus 30.4 ± 5.6%); PgR (69.5 ± 15.2% versus 6.0% -range 2.7-11.9%) and Ki67 (23.7 ± 0.54% versus 4.2%). Genetic studies were inconclusive though some intriguing variants were identified. CONCLUSION: The data implicate an important biologic role for ERα to increase tissue sensitivity to estrogen and aromatase to enhance local tissue production as biologic factors involved in juvenile gigantomastia.

The T Cell Immunoscore as a Reference for Biomarker Development Utilizing Real-World Data from Patients with Advanced Malignancies Treated with Immune Checkpoint Inhibitors.

BACKGROUND: We aimed to determine the prognostic value of an immunoscore reflecting CD3+ and CD8+ T cell density estimated from real-world transcriptomic data of a patient cohort with advanced malignancies treated with immune checkpoint inhibitors (ICIs) in an effort to validate a reference for future machine learning-based biomarker development. METHODS: Transcriptomic data was collected under the Total Cancer Care Protocol (NCT03977402) Avatar® project. The real-world immunoscore for each patient was calculated based on the estimated densities of tumor CD3+ and CD8+ T cells utilizing CIBERSORTx and the LM22 gene signature matrix. Then, the immunoscore association with overall survival (OS) was estimated using Cox regression and analyzed using Kaplan-Meier curves. The OS predictions were assessed using Harrell's concordance index (C-index). The Youden index was used to identify the optimal cut-off point. Statistical significance was assessed using the log-rank test. RESULTS: Our study encompassed 522 patients with four cancer types. The median duration to death was 10.5 months for the 275 participants who encountered an event. For the entire cohort, the results demonstrated that transcriptomics-based immunoscore could significantly predict patients at risk of death (p-value < 0.001). Notably, patients with an intermediate-high immunoscore achieved better OS than those with a low immunoscore. In subgroup analysis, the prediction of OS was significant for melanoma and head and neck cancer patients but did not reach significance in the non-small cell lung cancer or renal cell carcinoma cohorts. CONCLUSIONS: Calculating CD3+ and CD8+ T cell immunoscore using real-world transcriptomic data represents a promising signature for estimating OS with ICIs and can be used as a reference for future machine learning-based biomarker development.

Induction of PARP7 Creates a Vulnerability for Growth Inhibition by RBN2397 in Prostate Cancer Cells.

The ADP-ribosyltransferase PARP7 modulates protein function by conjugating ADP-ribose to the side chains of acceptor amino acids. PARP7 has been shown to affect gene expression in prostate cancer cells and certain other cell types by mechanisms that include transcription factor ADP-ribosylation. Here, we use a recently developed catalytic inhibitor to PARP7, RBN2397, to study the effects of PARP7 inhibition in androgen receptor (AR)-positive and AR-negative prostate cancer cells. We find that RBN2397 has nanomolar potency for inhibiting androgen-induced ADP-ribosylation of the AR. RBN2397 inhibits the growth of prostate cancer cells in culture when cells are treated with ligands that activate the AR, or the aryl hydrocarbon receptor, and induce PARP7 expression. We show that the growth-inhibitory effects of RBN2397 are distinct from its enhancement of IFN signaling recently shown to promote tumor immunogenicity. RBN2397 treatment also induces trapping of PARP7 in a detergent-resistant fraction within the nucleus, which is reminiscent of how inhibitors such as talazoparib affect PARP1 compartmentalization. Because PARP7 is expressed in AR-negative metastatic tumors and RBN2397 can affect cancer cells through multiple mechanisms, PARP7 may be an actionable target in advanced prostate cancer. Significance: RBN2397 is a potent and selective inhibitor of PARP7 that reduces the growth of prostate cancer cells, including a model for treatment-emergent neuroendocrine prostate cancer. RBN2397 induces PARP7 trapping on chromatin, suggesting its mechanism of action might be similar to clinically used PARP1 inhibitors.

Prehistoric human migration between Sundaland and South Asia was driven by sea-level rise.

Rapid sea-level rise between the Last Glacial Maximum (LGM) and the mid-Holocene transformed the Southeast Asian coastal landscape, but the impact on human demography remains unclear. Here, we create a paleogeographic map, focusing on sea-level changes during the period spanning the LGM to the present-day and infer the human population history in Southeast and South Asia using 763 high-coverage whole-genome sequencing datasets from 59 ethnic groups. We show that sea-level rise, in particular meltwater pulses 1 A (MWP1A, ~14,500-14,000 years ago) and 1B (MWP1B, ~11,500-11,000 years ago), reduced land area by over 50% since the LGM, resulting in segregation of local human populations. Following periods of rapid sea-level rises, population pressure drove the migration of Malaysian Negritos into South Asia. Integrated paleogeographic and population genomic analysis demonstrates the earliest documented instance of forced human migration driven by sea-level rise.

Whole genome sequencing identifies structural variants contributing to hematologic traits in the NHLBI TOPMed program.

Genome-wide association studies have identified thousands of single nucleotide variants and small indels that contribute to variation in hematologic traits. While structural variants are known to cause rare blood or hematopoietic disorders, the genome-wide contribution of structural variants to quantitative blood cell trait variation is unknown. Here we utilized whole genome sequencing data in ancestrally diverse participants of the NHLBI Trans Omics for Precision Medicine program (N = 50,675) to detect structural variants associated with hematologic traits. Using single variant tests, we assessed the association of common and rare structural variants with red cell-, white cell-, and platelet-related quantitative traits and observed 21 independent signals (12 common and 9 rare) reaching genome-wide significance. The majority of these associations (N = 18) replicated in independent datasets. In genome-editing experiments, we provide evidence that a deletion associated with lower monocyte counts leads to disruption of an S1PR3 monocyte enhancer and decreased S1PR3 expression.

Genetic Variation in Reproductive Investment Across an Ephemerality Gradient in Daphnia pulex.

Species across the tree of life can switch between asexual and sexual reproduction. In facultatively sexual species, the ability to switch between reproductive modes is often environmentally dependent and subject to local adaptation. However, the ecological and evolutionary factors that influence the maintenance and turnover of polymorphism associated with facultative sex remain unclear. We studied the ecological and evolutionary dynamics of reproductive investment in the facultatively sexual model species, Daphnia pulex. We found that patterns of clonal diversity, but not genetic diversity varied among ponds consistent with the predicted relationship between ephemerality and clonal structure. Reconstruction of a multi-year pedigree demonstrated the coexistence of clones that differ in their investment into male production. Mapping of quantitative variation in male production using lab-generated and field-collected individuals identified multiple putative quantitative trait loci (QTL) underlying this trait, and we identified a plausible candidate gene. The evolutionary history of these QTL suggests that they are relatively young, and male limitation in this system is a rapidly evolving trait. Our work highlights the dynamic nature of the genetic structure and composition of facultative sex across space and time and suggests that quantitative genetic variation in reproductive strategy can undergo rapid evolutionary turnover.

CCL22 mutations drive natural killer cell lymphoproliferative disease by deregulating microenvironmental crosstalk.

Chronic lymphoproliferative disorder of natural killer cells (CLPD-NK) is characterized by clonal expansion of natural killer (NK) cells where the underlying genetic mechanisms are incompletely understood. In the present study, we report somatic mutations in the chemokine gene CCL22 as the hallmark of a distinct subset of CLPD-NK. CCL22 mutations were enriched at highly conserved residues, mutually exclusive of STAT3 mutations and associated with gene expression programs that resembled normal CD16dim/CD56bright NK cells. Mechanistically, the mutations resulted in ligand-biased chemokine receptor signaling, with decreased internalization of the G-protein-coupled receptor (GPCR) for CCL22, CCR4, via impaired ß-arrestin recruitment. This resulted in increased cell chemotaxis in vitro, bidirectional crosstalk with the hematopoietic microenvironment and enhanced NK cell proliferation in vivo in transgenic human IL-15 mice. Somatic CCL22 mutations illustrate a unique mechanism of tumor formation in which gain-of-function chemokine mutations promote tumorigenesis by biased GPCR signaling and dysregulation of microenvironmental crosstalk.

Integrated bioinformatic pipeline using whole-exome and RNAseq data to identify germline variants correlated with cancer.

Germline Variants (GVs) are effective in predicting cancer risk and may be relevant in predicting patient outcomes. Here we provide a bioinformatic pipeline to identify GVs from the TCGA lower grade glioma cohort in Genomics Data Commons. We integrate paired whole exome sequences from normal and tumor samples and RNA sequences from tumor samples to determine a patient's GV status. We then identify the subset of GVs that are predictive of patient outcomes by Cox regression. For complete details on the use and execution of this protocol, please refer to Chatrath et al. (2019) and Chatrath et al. (2020).

Telomere length is not a main factor for the development of islet autoimmunity and type 1 diabetes in the TEDDY study.

The Environmental Determinants of Diabetes in the Young (TEDDY) study enrolled 8676 children, 3-4 months of age, born with HLA-susceptibility genotypes for islet autoimmunity (IA) and type 1 diabetes (T1D). Whole-genome sequencing (WGS) was performed in 1119 children in a nested case-control study design. Telomere length was estimated from WGS data using five tools: Computel, Telseq, Telomerecat, qMotif and Motif_counter. The estimated median telomere length was 5.10 kb (IQR 4.52-5.68 kb) using Computel. The age when the blood sample was drawn had a significant negative correlation with telomere length (P = 0.003). European children, particularly those from Finland (P = 0.041) and from Sweden (P = 0.001), had shorter telomeres than children from the U.S.A. Paternal age (P = 0.019) was positively associated with telomere length. First-degree relative status, presence of gestational diabetes in the mother, and maternal age did not have a significant impact on estimated telomere length. HLA-DR4/4 or HLA-DR4/X children had significantly longer telomeres compared to children with HLA-DR3/3 or HLA-DR3/9 haplogenotypes (P = 0.008). Estimated telomere length was not significantly different with respect to any IA (P = 0.377), IAA-first (P = 0.248), GADA-first (P = 0.248) or T1D (P = 0.861). These results suggest that telomere length has no major impact on the risk for IA, the first step to develop T1D. Nevertheless, telomere length was shorter in the T1D high prevalence populations, Finland and Sweden.

Detecting molecular subtypes from multi-omics datasets using SUMO.

We present a data integration framework that uses non-negative matrix factorization of patient-similarity networks to integrate continuous multi-omics datasets for molecular subtyping. It is demonstrated to have the capability to handle missing data without using imputation and to be consistently among the best in detecting subtypes with differential prognosis and enrichment of clinical associations in a large number of cancers. When applying the approach to data from individuals with lower-grade gliomas, we identify a subtype with a significantly worse prognosis. Tumors assigned to this subtype are hypomethylated genome wide with a gain of AP-1 occupancy in demethylated distal enhancers. The tumors are also enriched for somatic chromosome 7 (chr7) gain, chr10 loss, and other molecular events that have been suggested as diagnostic markers for "IDH wild type, with molecular features of glioblastoma" by the cIMPACT-NOW consortium but have yet to be included in the World Health Organization (WHO) guidelines.

Protocol for integrative subtyping of lower-grade gliomas using the SUMO pipeline.

Grouping patients into subtypes with homogeneous molecular features can guide diagnosis and therapeutic interventions. SUMO is a computational pipeline that uses nonnegative matrix factorization of patient-similarity networks to integrate continuous multi-omic datasets for molecular subtyping of a disease. Here, we present a detailed protocol to demonstrate its use in determining subtypes of lower-grade gliomas by integrating gene expression, DNA methylation, and miRNA expression data from the TCGA-LGG cohort. For complete details on the use and execution of this profile, please refer to Sienkiewicz et al. (2022).

Genomic landscape of TCRαß and TCRγδ T-large granular lymphocyte leukemia.

Large granular lymphocyte (LGL) leukemia comprises a group of rare lymphoproliferative disorders whose molecular landscape is incompletely defined. We leveraged paired whole-exome and transcriptome sequencing in the largest LGL leukemia cohort to date, which included 105 patients (93 T-cell receptor αß [TCRαß] T-LGL and 12 TCRγδ T-LGL). Seventy-six mutations were observed in 3 or more patients in the cohort, and out of those, STAT3, KMT2D, PIK3R1, TTN, EYS, and SULF1 mutations were shared between both subtypes. We identified ARHGAP25, ABCC9, PCDHA11, SULF1, SLC6A15, DDX59, DNMT3A, FAS, KDM6A, KMT2D, PIK3R1, STAT3, STAT5B, TET2, and TNFAIP3 as recurrently mutated putative drivers using an unbiased driver analysis approach leveraging our whole-exome cohort. Hotspot mutations in STAT3, PIK3R1, and FAS were detected, whereas truncating mutations in epigenetic modifying enzymes such as KMT2D and TET2 were observed. Moreover, STAT3 mutations co-occurred with mutations in chromatin and epigenetic modifying genes, especially KMT2D and SETD1B (P < .01 and P < .05, respectively). STAT3 was mutated in 50.5% of the patients. Most common Y640F STAT3 mutation was associated with lower absolute neutrophil count values, and N647I mutation was associated with lower hemoglobin values. Somatic activating mutations (Q160P, D170Y, L287F) in the STAT3 coiled-coil domain were characterized. STAT3-mutant patients exhibited increased mutational burden and enrichment of a mutational signature associated with increased spontaneous deamination of 5-methylcytosine. Finally, gene expression analysis revealed enrichment of interferon-γ signaling and decreased phosphatidylinositol 3-kinase-Akt signaling for STAT3-mutant patients. These findings highlight the clinical and molecular heterogeneity of this rare disorder.

Genomic analyses of the metastasis-derived prostate cancer cell lines LNCaP, VCaP, and PC3-AR.

BACKGROUND: The lymph node metastasis-derived LNCaP, the bone metastasis-derived PC3 (skull), and VCaP (vertebral) cell lines are widely used as preclinical models of human prostate cancer (CaP) and have been described in more than 19,000 publications. Here, we report on short-read whole-genome sequencing and genomic analyses of LNCaP, VCaP, and PC3 cells stably transduced with WT AR (PC3-AR). METHODS: LNCaP, VCaP, and PC3-AR cell lines were sequenced to an average depth of more than 30-fold using Illumina short-read sequencing. Using various computational methods, we identified and compared the single-nucleotide variants, copy-number profiles, and the structural variants observed in the three cell lines. RESULTS: LNCaP cells are composed of multiple subpopulations, which results in nonintegral copy number states and a high mutational load when the data is analyzed in bulk. All three cell lines contain pathogenic mutations and homozygous deletions in genes involved in DNA mismatch repair, along with deleterious mutations in cell-cycle, Wnt signaling, and other critical cellular processes. PC3-AR cells have a truncating mutation in TP53 and do not express the p53 protein. The VCaP cells contain a homozygous gain-of-function mutation in TP53 (p.R248W) that promotes cancer invasion, metastasis, and progression and has also been observed in prostate adenocarcinomas. In addition, we detect the signatures of chromothripsis of the q arms of chromosome 5 in both PC3-AR and VCaP cells, strengthening the association of TP53 inactivation with chromothripsis reported in other systems. CONCLUSIONS: Our work provides a resource for genetic, genomic, and biological studies employing these commonly-used prostate cancer cell lines.

Pangenomics enables genotyping of known structural variants in 5202 diverse genomes.

We introduce Giraffe, a pangenome short-read mapper that can efficiently map to a collection of haplotypes threaded through a sequence graph. Giraffe maps sequencing reads to thousands of human genomes at a speed comparable to that of standard methods mapping to a single reference genome. The increased mapping accuracy enables downstream improvements in genome-wide genotyping pipelines for both small variants and larger structural variants. We used Giraffe to genotype 167,000 structural variants, discovered in long-read studies, in 5202 diverse human genomes that were sequenced using short reads. We conclude that pangenomics facilitates a more comprehensive characterization of variation and, as a result, has the potential to improve many genomic analyses.

Human ACE2 receptor polymorphisms and altered susceptibility to SARS-CoV-2.

COVID-19 is a respiratory illness caused by a novel coronavirus called SARS-CoV-2. The viral spike (S) protein engages the human angiotensin-converting enzyme 2 (ACE2) receptor to invade host cells with ~10-15-fold higher affinity compared to SARS-CoV S-protein, making it highly infectious. Here, we assessed if ACE2 polymorphisms can alter host susceptibility to SARS-CoV-2 by affecting this interaction. We analyzed over 290,000 samples representing >400 population groups from public genomic datasets and identified multiple ACE2 protein-altering variants. Using reported structural data, we identified natural ACE2 variants that could potentially affect virus-host interaction and thereby alter host susceptibility. These include variants S19P, I21V, E23K, K26R, T27A, N64K, T92I, Q102P and H378R that were predicted to increase susceptibility, while variants K31R, N33I, H34R, E35K, E37K, D38V, Y50F, N51S, M62V, K68E, F72V, Y83H, G326E, G352V, D355N, Q388L and D509Y were predicted to be protective variants that show decreased binding to S-protein. Using biochemical assays, we confirmed that K31R and E37K had decreased affinity, and K26R and T92I variants showed increased affinity for S-protein when compared to wildtype ACE2. Consistent with this, soluble ACE2 K26R and T92I were more effective in blocking entry of S-protein pseudotyped virus suggesting that ACE2 variants can modulate susceptibility to SARS-CoV-2.

Germline variants predictive of tumor mutational burden and immune checkpoint inhibitor efficacy.

High tumor mutational burden (TMB) is associated with response to checkpoint blockade in several cancers. We identify pathogenic germline variants associated with increased TMB (GVITMB). GVITMB were found in 7 genes using a pan-cancer approach (APC, FANCL, SLC25A13, ERCC3, MSH6, PMS2, and TP53) and 38 gene sets (e.g., those involved in DNA repair and programmed cell death). GVITMB were also associated with mutational signatures related to the dysfunction of the gene carrying the variant, somatic mutations that further affect the gene or pathway with the variant, or transcriptomic changes concordant with the expected effect of the variant. In a validation cohort of 140 patients with cutaneous melanoma, we found that patients with GVITMB had prolonged progression-free survival (p = 0.0349, hazard ratio = 0.688) and responded favorably (p = 0.0341, odds = 1.842) when treated with immune checkpoint inhibitors. Our results suggest that germline variants can influence the molecular phenotypes of tumors and predict the response to immune checkpoint inhibitors.

Frequent somatic TET2 mutations in chronic NK-LGL leukemia with distinct patterns of cytopenias.

Chronic natural killer large granular lymphocyte (NK-LGL) leukemia, also referred to as chronic lymphoproliferative disorder of NK cells, is a rare disorder defined by prolonged expansion of clonal NK cells. Similar prevalence of STAT3 mutations in chronic T-LGL and NK-LGL leukemia is suggestive of common pathogenesis. We undertook whole-genome sequencing to identify mutations unique to NK-LGL leukemia. The results were analyzed to develop a resequencing panel that was applied to 58 patients. Phosphatidylinositol 3-kinase pathway gene mutations (PIK3CD/PIK3AP1) and TNFAIP3 mutations were seen in 5% and 10% of patients, respectively. TET2 was exceptional in that mutations were present in 16 (28%) of 58 patient samples, with evidence that TET2 mutations can be dominant and exclusive to the NK compartment. Reduced-representation bisulfite sequencing revealed that methylation patterns were significantly altered in TET2 mutant samples. The promoter of TET2 and that of PTPRD, a negative regulator of STAT3, were found to be methylated in additional cohort samples, largely confined to the TET2 mutant group. Mutations in STAT3 were observed in 19 (33%) of 58 patient samples, 7 of which had concurrent TET2 mutations. Thrombocytopenia and resistance to immunosuppressive agents were uniquely observed in those patients with only TET2 mutation (Games-Howell post hoc test, P = .0074; Fisher's exact test, P = .00466). Patients with STAT3 mutation, inclusive of those with TET2 comutation, had lower hematocrit, hemoglobin, and absolute neutrophil count compared with STAT3 wild-type patients (Welch's t test, P ≤ .015). We present the discovery of TET2 mutations in chronic NK-LGL leukemia and evidence that it identifies a unique molecular subtype.

Upregulation of human endogenous retrovirus-K (HML-2) mRNAs in hepatoblastoma: Identification of potential new immunotherapeutic targets and biomarkers.

PURPOSE: Hepatoblastoma is the most common liver malignancy in children. In order to advance therapy against hepatoblastoma, novel immunologic targets and biomarkers are needed. Our purpose in this investigation is to examine hepatoblastoma transcriptomes for the expression of a class of genomic elements known as Human Endogenous Retrovirus (HERVs). HERVs are abundant in the human genome and are biologically active elements that have been associated with multiple malignancies and proposed as immunologic targets in a subset of tumors. A sub-family of HERVs, HERV-K(HML-2) (HERV-K), have been shown to be tightly regulated in fetal development, making investigation of these elements in pediatric tumors paramount. METHODS: We first created a HERVK-FASTA file utilizing 91 previously described HML-2 proviruses. We then concatenated the file onto the GRCh38.95 cDNA library from Ensembl. We used this reference database to evaluate existing RNA-seq data from 10 hepatoblastoma tumors and 3 normal liver controls (GEO accession ID: GSE8977575). Quantification and differential proviral expression analysis between hepatoblastoma and normal liver controls was performed using the pseudo-alignment program Salmon and DESeq2, respectively. RESULTS: HERV-K mRNA was expressed in hepatoblastoma from multiple proviral loci. All expressed HERV-K proviral loci were upregulated in hepatoblastoma compared to normal liver controls. Five HERV-K proviruses (1q21.3, 3q27.2, 7q22.2, 12q24.33 and 17p13.1) were significantly differentially expressed (p-adjusted value <0.05, |log2 fold change| > 1.5) across conditions. The provirus at 17p13.1 had an approximately 300-fold increased expression in hepatoblastoma as compared to normal liver. This was in part due to the near absence of HERV-K mRNA at the 17p13.1 locus in fully differentiated liver samples. CONCLUSIONS: Our investigation demonstrates that HERV-K is expressed from multiple loci in hepatoblastoma and that the expression is increased for several proviruses compared to normal liver controls. Our results suggest that HERV-K mRNA expression may be useful as a biomarker in hepatoblastoma, given the large differential expression profiles in hepatoblastoma, with very low mRNA levels in liver control samples.