Secondary Logo

Journal Logo

Original Articles

HPV51-associated Leiomyosarcoma

A Novel Class of TP53/RB1-Wildtype Tumor With Predilection for the Female Lower Reproductive Tract

Williams, Erik A. MD*,†; Montesion, Meagan PhD; Lincoln, Vadim MD*; Tse, Julie Y. MD; Hiemenz, Matthew C. MD, MS; Mata, Douglas A. MD, MPH; Shah, Bhamini B. BS; Shoroye, Adebowale BS; Alexander, Brian M. MD, MPH; Werth, Adrienne J. MD; Foley-Peres, Kathleen PhD§; Milante, Riza R. MD, DPDS; Ross, Jeffrey S. MD†,¶; Ramkissoon, Shakti H. MD, PhD†,#; Williams, Kevin Jon MD**; Adhikari, Laura J. MD††; Zuna, Rosemary E. MD††; LeBoit, Philip E. MD*; Lin, Douglas I. MD; Elvin, Julia A. MD, PhD

Author Information
The American Journal of Surgical Pathology: June 2022 - Volume 46 - Issue 6 - p 729-741
doi: 10.1097/PAS.0000000000001862

Abstract

Leiomyosarcoma (LMS), a malignant neoplasm defined by smooth muscle differentiation, is the most common gynecologic (GYN) sarcoma and can also arise in soft tissues elsewhere.1,2 GYN LMS arises most frequently in the uterine wall and only rarely as a primary vaginal neoplasm.3–5 GYN LMS is aggressive and resists standard therapy, with high rates of recurrence and progression; five-year survival for women with metastatic disease at presentation is only 10% to 15%.6,7 Stage of disease, as defined by the International Federation of Gynecology and Obstetrics (FIGO)8 or the American Joint Committee on Cancer,9 is an important prognostic factor at the time of diagnosis.1 Surgery is the standard of care for localized tumors, with hormonal and cytotoxic chemotherapy reserved for advanced stages.10 Inactivating mutations in tumor suppressor genes TP53 and RB1 are considered central drivers in the majority of LMS.11

In high-risk (hr) human papillomavirus (HPV)-related tumors, a similar functional outcome is achieved through HPV oncoproteins E6 and E7, which inactivate the protein products p53 and RB1, respectively.12 Genotypes of hrHPV, such as HPV16 and HPV18, are associated with the development of a subset of carcinomas, including cervical, anogenital, and a proportion of oropharyngeal carcinomas.12–14 A potential role for hrHPV in the development of non-epithelial tumors, such as LMS, has not been defined.

Here, we hypothesized that HPV infection could provide an alternative mechanism for tumorigenesis in a distinct subset of TP53/RB1-wildtype LMS. To evaluate our hypothesis, we surveyed a large database of tumors previously diagnosed as LMS for recurrent HPV genome reads that occur in the absence of detectable mutations in TP53 or RB1. Strikingly, we discovered a novel recurrent HPV51-positive genomic signature that was restricted to LMSs from lower female genital tract that were both TP53- and RB1-wildtype.

METHODS

Clinical Cohort and Comprehensive Genomic Profiling

Comprehensive genomic profiling was performed in a Clinical Laboratory Improvement Amendments (CLIA)-certified, College of American Pathologists (CAP)-accredited laboratory (Foundation Medicine Inc., Cambridge, MA). Approval for this study, including a waiver of informed consent and a Health Insurance Portability and Accountability Act (HIPAA) waiver of individual authorization, was obtained from the Western Institutional Review Board (Protocol No. 20152817). Our cohort initially included all GYN and soft tissue LMS referred to Foundation Medicine for comprehensive genomic profiling as part of clinical care (n=2585 unique cases, including 12 of confirmed vaginal primary site). Board-certified pathologists on staff at Foundation Medicine reviewed the submitted pathologic diagnosis and examined routine hematoxylin and eosin (H&E)-stained slides of each case before sequencing. Molecular results were also reviewed by on-staff pathologists, and 16 cases were excluded from the final LMS cohort owing to the presence of molecular alterations that are considered defining for sarcomas other than LMS, meaning BCOR-ZC3H7B (n=6), PDGFB-COL1A1 (n=3), STAT6-NAB2 (n=2), GREB1-SS18 (n=1), SS18-SSX1 (n=1), SS18-SSX2 (n=1), and EWSR1-PATZ1 (n=1) rearrangements, and BCOR internal tandem duplication (n=1), for a final cohort of 2,569 LMS.

Before molecular analyses, sections were macrodissected to achieve >20% estimated percent tumor nuclei in each case, where percent tumor nuclei is defined as 100 times the number of tumor cells divided by total number of nucleated cells.

For genomic assessments, ≥60 ng DNAs was extracted from 40 μm sections cut from formalin-fixed, paraffin-embedded tissue blocks. The samples were assayed by adaptor ligation hybrid capture, performed for all coding exons of 236 (v1), 315 (v2), or 405 (v3) cancer‐related genes plus select introns from 19 (v1), 28 (v2), or 31 (v3) genes frequently rearranged in cancer (Supplemental Table 1, Supplemental Digital Content 1, https://links.lww.com/PAS/B306).15,16 For samples with available RNA (n=2,038), targeted RNA‐seq was performed for analysis of rearrangements in 265 genes.16 Sequencing of captured libraries was performed using the Illumina sequencing system to a mean exon coverage depth of targeted regions of >×500, and sequences were analyzed for genomic alterations (GAs), including short variant alterations (base substitutions, insertions, and deletions), copy number alterations (focal amplifications and homozygous deletions), and select gene fusions or rearrangements.15,17,18 To maximize mutation detection accuracy in impure clinical specimens, the test was previously optimized and validated to detect base substitutions at a ≥5% mutant allele frequency, indels with a ≥10% mutant allele frequency with ≥99% accuracy, and fusions occurring within baited introns/exons with >99% sensitivity.15 Tumor mutational burden (TMB, mutations/Mb) was determined on 0.8 to 1.1 Mbp of sequenced DNA.18 Microsatellite instability (MSI) was determined on up to 114 loci.19

HPV genome sequences were detected by de novo assembly of nonhuman sequencing reads followed by comparison against all viral nucleotide sequences in the National Center for Biotechnology Information (NCBI) RefSeq database using the Basic Local Alignment Search Tool for nucleotides (BLASTn). The RefSeq database is comprehensive and distinguishes HPV types including HPV6, 11, 16, 18, 26, 30, 31, 33, 35, 39, 40, 42, 43, 44, 45, 51, 52, 53, 54, 55, 56, 58, 59, 61, 62, 64, 66, 67, 68, 69, 70, 71, 72, 73, 81, 82, 83, 84, 87, CP6108, and IS39. HPV types identified in this study were stratified according to the HPV classification of Muñoz et al, with HPV16, 18, 31, 33, 51, and 59 labeled high-risk for causing cancer and HPV6 and 11 labeled low-risk.20 HPV30 and 87 were classified as undetermined risk.21 HPV positivity was designated for cases with assembled contiguous sequences (contigs) of viral genomic DNA ≥80 nucleotides in length with ≥97% sequence identity to the Refseq sequence identified by BLASTn. Within the HPV51 type, specific sublineages (A1-A4, B1-B2) were identified using short regions that contain lineage-specific and sublineage-specific single nucleotide polymorphisms.22

Patient ancestry was determined by classifying specific single nucleotide variations by genomic profiling based on their known variation among populations in the 1000 Genomes Project.23

Confirmatory Immunohistochemistry and In Situ Hybridization In-house

We performed confirmatory analyses in-house on all HPV51-positive LMS tumors with available tissue (n=3). These analyses included immunohistochemistry for desmin (Agilent Technologies, Lexington, MA) and caldesmon (Agilent Technologies), both key markers of smooth muscle lineage. In addition, we completed in situ hybridization for high-risk HPV E6/E7 mRNA using the RNAscope HPV kit (Advanced Cell Diagnostics Inc., Hayward, CA) following the manufacturer’s instructions. Formalin-fixed, paraffin-embedded tissue sections were hybridized with a cocktail to detect E6 or E7 mRNA from 18 high-risk HPV types: HPV16, 18, 26, 31, 33, 35, 39, 45, 51, 52, 53, 56, 58, 59, 66, 68, 73, and 82. Approval for immunohistochemistry and in situ hybridization for these targets, including a waiver of additional informed consent and a HIPAA waiver of authorization, was obtained from the University of California San Francisco Institutional Review Board (protocol No. 11-05569).

Clinicopathologic Assessments of HPV51-positive Leiomyosarcoma Cases

Within our overall cohort of LMS cases, we performed in-depth assessments of the subset harboring HPV51 genomic sequences (HPV51-positive). Clinicopathologic data on the HPV51-positive LMS cases included patient age, sex, tumor site, and FIGO stage or American Joint Committee on Cancer stage (8th edition), which were extracted from accompanying pathology reports.9,24 Primary site data were available for nearly all cases (91%), with the remaining case designated “indeterminant primary.” Histopathology of all HPV51-positive cases was re-assessed on routine H&E-stained slides of tissue sections by 2 board-certified pathologists (E.A.W., D.I.L.).

For categorical data, statistical comparisons between HPV51-positive versus HPV51-negative LMS cases were performed using the Fisher exact test owing to the size of the cohort. For continuous variables (age and TMB), comparisons between the 2 groups were performed using the nonparametric Mann-Whitney U test. A 2‐tailed P‐value of <0.05 was considered statistically significant.

RESULTS

Comprehensive Genomic Profiling

To address our hypothesis, we began with a pan-LMS analysis to find cases that were TP53/RB1-wildtype and harbored HPV genomic sequences. From the entire 2569 LMS cohort, 486 LMS were both TP53-wildtype and RB1-wildtype (486/2569; 18.9%). From the entire cohort, 48 unique LMS cases carried detectable HPV genomic reads (n=48/2569; 1.9%; Fig. 1). Overlap between these groups consisted of 18 LMS cases that were TP53-wildtype, RB1-wildtype, and HPV-positive (Figs. 1A–C). Surprisingly, most of these cases (n=11) were HPV51-positive (Fig. 1A), and these 11 represented all HPV51-positive tumors in our entire LMS database. HPV51-positive LMS comprised 2.3% of all TP53/RB1-wildtype LMS (n=11/486). A single HPV51-positive case had concurrent reads of additional HPV types 59 and 87 (Fig. 1A, rightmost case). All HPV51-positive LMS contained regions that matched to HPV51 sublineage A1, except the single case with concurrent non-HPV51 reads, which matched to HPV51 sublineage B1.

F1
FIGURE 1:
Summary tile plot of oncogenic variants in key genes in all leiomyosarcoma cases in which HPV genomic sequences were detected (n=48/2,569; 1.9%). Here, cases were split into 4 groups from left to right (A–D) based on genomic signature, with HPV the likely tumor driver in group A and potentially B, and HPV genomic reads presumed to be an incidental finding in groups C and D. All 11 cases of HPV51-positive leiomyosarcoma (A) were TP53/RB1-wildtype with no detected alterations in CDKN2C or FAF1, as were 3 additional cases harboring 3 specific non-HPV51 genotypes (B). Next, 4 HPV51-negative, FAF1-null cases are shown, 3 of which were also CDKN2C-null (C). The remaining 30 HPV-positive LMS cases were all HPV51-negative and harbored TP53 and/or RB1 mutations (D). Each column represents data for a single unique patient. Age, sex, HPV status, and tumor primary site are also provided for each case. The histogram on top shows tumor mutational burden (mutations/megabase).

The distribution of GAs in all HPV-positive LMS tumors in our database is displayed in Figure 1, and GAs for all LMS cases are shown in Table 1. HPV51-positive LMS cases (n=11) were genomically distinct from HPV-positive but HPV51-negative LMS cases (n=37; Fig. 1A vs. Figs. 1B–D, especially Fig. 1A vs. the 34 cases in Figs. 1C, D) and from all HPV51-negative cases (n=2558; Table 1), with most of this last group negative for any HPV sequences, as just noted. In particular, TP53 and RB1 GAs were strikingly absent from all of the HPV51-positive LMS tumors (Fig. 1A), in contrast to the high prevalence of TP53 and RB1 GAs in HPV-positive but HPV51-negative LMS cases (n=37; Figs. 1B–D; P<0.00001 and 0.001, respectively) and in all of the HPV51-negative LMS (n=2558; Table 1; P<0.00001 and 0.0002, respectively). The only recurrent pathogenic GAs in HPV51-positive LMS were in ATRX and TSC1, and these GAs were significantly enriched in HPV51-positive LMS compared with HPV51-negative LMS (Table 1; P=0.027 and 0.0047, respectively). No HPV-positive LMS contained any detectable pathogenic alterations in PIK3CA.

TABLE 1 - Demographics and Frequencies of Genomic Alterations in Leiomyosarcomas Stratified by HPV51 Status
HPV51-positive LMS HPV-negative LMS P
Number of cases 11 2558
% female 100% (11/11) 80% (2044/2558) 0.14
Median age (range; years) 54 (23-74) 58 (<1-89+) 0.16
TMB (Q1-Q3; mut/Mb) 1.6 (0.8-2.4) 2.4 (1.6-4.0) 0.50
% MSI high 0% (0/11) 0.2% (5/2144) >0.999
TP53 0% (0/11) 72% (1839/2558) <0.00001
RB1 0% (0/11) 53% (1359/2558) 0.0002
ATRX 55% (6/11) 24% (602/2538) 0.027
PTEN 0% (0/11) 16% (406/2558) 0.23
TSC1 18% (2/11) 0.6% (23/2558) 0.0047
For percent values, number of positive cases over number of evaluated cases are included in parentheses. Parameters significantly different between HPV51-positive versus HPV51-negative cases are indicated in bold (P<0.05).

Among the HPV51-negative LMS cases, ones with versus without other HPV sequences showed considerable genomic similarities, e.g., GAs in TP53 were 78% (29/37, Figs. 1B–D) and 72% (1810/2521), respectively; GAs in RB1 were 57% (21/37) and 53% (1338/2521), respectively; and cases that were both TP53-wildtype and RB1-wildtype were 18.9% (7/37, Figs. 1B–D) and 18.6% (468/2521), respectively. These similarities suggest that HPV sequences other than HPV51 in our LMS cohort may be merely incidental findings. Importantly, this pattern highlights the contrast between the distinctive genomic signatures of HPV51-positive (Fig. 1A, Table 1) versus HPV51-negative (Figs. 1B–D, Table 1) LMS cases.

Median TMB for HPV51-positive LMS was 1.6 mut/Mb (range <0.8 to 9.6; Q1-Q3=0.8 to 2.4), similar to the remainder of the LMS cohort (median=2.4 mut/Mb; range <0.8 to 203; Q1-Q3=1.6 to 4.0) (Table 1; P=0.50, Mann-Whitney U test). Median computational tumor purity for HPV51-positive LMS was 71.5% (range: 46.6% to 99.6%). No cases of MSI high were present in the HPV51-positive cohort, which is not statistically distinguishable from the rate of 0.2% in the HPV51-negative LMS group (Table 1).

Of the 37 LMS cases with non-HPV51 HPV reads, 26 contained HPV16, 6 contained HPV6, and 1 each contained HPV11, 18, 30, 31, and 33 (Figs. 1B–D). Of these 37 cases, only 7 (18.9%) were both TP53-wildtype and RB1-wildtype (Figs. 1B, C), as noted above. Four of these 7 cases harbored HPV16 reads, of which 3 were CDKN2C-null and the remaining case was null for an adjacent gene, FAF1, consistent with our prior report implicating a central oncogenic role for these gene deletions in a subset of TP53/RB1-wildtype LMS25 (Fig. 1C). The remaining three TP53/RB1-wildtype, HPV-positive, HPV51-negative cases were the ones harboring HPV30, 31, and 33 reads, respectively (Fig. 1B). The TP53/RB1-wildtype HPV30 and 33 cases were both ATRX-mutant (Fig. 1B). Of note, the TP53/RB1-wildtype HPV30 case was of vaginal primary site. These genomic signatures suggest that the HPV genotypes in Figure 1B could represent mimics of HPV51, but additional cases would be needed to support that inference.

All 6 vaginal primary LMS without HPV reads contained mutations in TP53; 4 of these 6 were also RB1-mutant, and 2 of the 6 were ATRX-mutant.

Clinicopathologic Features of HPV51-positive Leiomyosarcoma

From our cohort of 2569 LMS cases, HPV51-positive LMS accounted for 0.4% of the total (n=11/2569) and 42% of vaginal primary LMS (n=5/12), representing a significant enrichment for HPV51-positive cases at this site (42% [5/12] vs. 0.2% of HPV51-negative LMS [6/2557], P<0.0001). Clinical characteristics of HPV51-positive LMS are summarized in Tables 1 and 2. Patients were similar in age and in female-preponderance between HPV51-positive (n=11) and HPV51-negative LMS (n=2558; Table 1). Most patients with HPV51-positive LMS in our cohort had clinically advanced disease (Table 2). Anatomic locations and other characteristics of the HPV51-positive tumor specimens that we sequenced are listed in Table 3. Detailed past medical history was available for patient cases 1 and 2 and was unremarkable, other than a history of a cervical cone with cervical intraepithelial neoplasia 3 in patient case 2 that was anatomically discrete from the vaginal mass.

TABLE 2 - Clinical Characteristics of Patients With HPV51-positive Leiomyosarcoma
Characteristic N (%)
No of patients 11
Median age at diagnosis, years 54, range: 23-74
Sex
 Female 11 (100)
 Male 0
Ancestry*
 European 6 (55)
 African 3 (27)
 Ad-Mixed American 2 (18)
Primary site
 Vagina 5 (46)
 Uterine cervix 2 (18)
 Uterine, not otherwise specified 2 (18)
 Perineum/perirectal 1 (9)
 Indeterminant 1 (9)
FIGO staging (uterine primary)
 IVA 1 (25)
 IVB 2 (50)
 Unknown 1 (25)
AJCC staging (nonuterine or indeterminant primary)
 II 1 (14)
 IIIA 1 (14)
 IV 3 (43)
 Unknown 2 (29)
*Ancestry superpopulations as defined by the 1000 Genomes Project.23
The “Ad-Mixed American” superpopulation was defined by specific SNPs carried by individuals of Mexican-American, Puerto Rican, Colombian, or Peruvian ancestry.23
AJCC indicates American Joint Committee on Cancer.

TABLE 3 - Detailed Clinicopathologic and Molecular Features of All 11 HPV51-positive Leiomyosarcomas
Patient Case Age (y)/sex Primary Tumor Site Anatomic Location of Sequenced Tumor Primary Tumor Size (cm) Predominant Morphology of Sequenced Tumor Marked Nuclear Atypia Mitoses/10 HPF Necrosis Border of Infiltration Vascular Invasion Immunohistochemistry HPV E6/E7 mRNA In Situ Hybridization Genomic Alterations Disease and Survival Status
1 56/F Uterus Uterus 23.1×15.6×9.3 Spindle Present 12 Present Pushing Not identified Caldesmon: strongly positive Desmin: focally positive Smooth muscle myosin: strongly positive; CD34/CD117/DOG1: negative ER: negative Positive ATRX p.E1862fs*17 Dead of disease at 15 mo
2 54/F Vagina Vagina 3.0×2.4×1.7 Spindle, with focal myxoid areas Present 20 Absent Pushing Not identified Caldesmon: strongly positive Desmin: strongly positive Per report: Calponin, smooth muscle actin, muscle specific actin, ER, and PR: diffuse strong positive; Negative for DOG1, CD117, and pan-melanoma; P53: wild-type pattern Positive None detected Disease-free at 84 mo
3 30/F Vagina Vagina 11.0×9.0×5.0 Spindle, with focal myxoid areas Present 7 Absent CBD Not identified Caldesmon: focally positive Desmin: focally positive Per report: negative for CD34 Positive None detected Not available
4 69/F CBD Axillary soft tissue 4.0×3.0×3.0 (axillary mass) Epithelioid, with clear cell features Present 11 Absent Pushing Not identified Per report: Smooth muscle myosin, smooth muscle actin, and vimentin positive; Negative for S100, MART-1, HMB45, CAM5.2, and AE1/AE3 Not performed None detected Not available
5 40/F Vagina Vagina CBD Epithelioid, with focal spindle areas Present 20 Present CBD Not identified Per report: smooth muscle actin, desmin, vimentin, p16, ER, PR, and p63 positive; Negative for PIN, PAX8, S100, pancytokeratin, CK5/6, CK7, CK20, HMB-45, S100, myogenin, and MyoD-1 Not performed ATRX p.Q2242*, TSC2 homozygous loss Not available
6 47/F Vagina Lymph node, station 7 2.5×2.0×1.0 Spindle, with focal epithelioid areas Present 22 Absent CBD Not identified Per report: smooth muscle actin positive Not performed ATRX p.S541*, CCND3 amplification Not available
7 74/F Perineum/perirectal Perirectal CBD Epithelioid Present 5 Absent CBD Not identified Per report: smooth muscle actin, calponin, and caldesmon weakly positive; Negative for SOX-10 Not performed ATRX p.E1065* Not available
8 23/F Vagina Vagina 6.2×5.0×5.0 Spindle, with focal myxoid area Present 14 Absent CBD Not identified Per report: ER positive (50% moderate nuclear staining) and PR positive (80% strong nuclear staining) Not performed None detected Not available
9 57/F Uterine cervix Lung CBD Epithelioid, with focal spindle areas Present 15 Present CBD Not identified Per report: smooth muscle actin weakly positive, desmin focally positive, and vimentin positive; Negative for pancytokeratin and CK7 Not performed ATRX p.R654*, TSC2 homozygous loss Not available
10 66/F Uterus Pancreas CBD Epithelioid Present 10 Absent CBD Not identified Per report: negative for chromogranin and cytokeratin cocktail Not performed ATRX p.R444* Not available
11 41/F Uterine cervix Uterine cervix CBD Not evaluable owing to crush artifact Present Not evaluable owing to crush artifact Present CBD Not identified Not specified on available report Not performed None detected Not available
Key smooth muscle immunohistochemical markers are indicated in bold.
CBD indicates cannot be determined; F, female.

Histopathology

Histopathologic re-evaluation of all 11 HPV51-positive LMS tumors was performed on high-resolution digital pathology of H&E-stained slides. Histology was heterogeneous, as shown in Figures 2–5 and summarized in Table 3, which includes tumor size, predominant differentiation (spindle, epithelioid, or myxoid), nuclear atypia, number of mitoses, presence of tumor necrosis, border of infiltration and vascular invasion, if known. Overall, HPV51-positive LMS were mitotically active tumors (range: 5 to 22 mitoses per 10 high-power fields, median 14 per 10 high-power fields) with marked nuclear atypia, necrosis, and pushing borders. Five cases had predominant epithelioid morphology and 5 cases were predominantly spindle cell LMS. Predominant morphology of 1 case could not be determined owing to crush artifact (Table 3). Three cases, all vaginal primary cases, showed focal myxoid areas.

F2
FIGURE 2:
Radiology and pathology from patient case 1. A, Abdominopelvic computed tomography with contrast at initial presentation demonstrated a large heterogenously enhancing pelvic mass (white arrows) and a 2.6-cm hypodense mass in the right lobe of the liver (grey arrow). Multiple liver and lung masses were present on different abdominal and thoracic planes on initial imaging (not shown). Histopathologic examination of the pelvic tumor showed fascicles of atypical spindle cells with a pushing border of infiltration (blue arrow) (B; hematoxylin and eosin stain). The tumor also showed areas of geographic necrosis (C) and was mitotically active (C, inset). By immunohistochemistry, tumor cells showed diffuse caldesmon (D), diffuse smooth muscle myosin (E), and scattered desmin positivity (F). By in situ hybridization for high-risk HPV E6/E7 mRNA (detecting 18 HPV genotypes including HPV51), tumors cells showed prominent punctate brown staining, while adjacent nontumoral stroma was negative (blue arrow) (G, H).
F3
FIGURE 3:
Pathology from patient case 2. A, Histopathologic examination of the tumor showed fascicles of spindle cells with focal myxoid areas (inset) (hematoxylin and eosin stain). By immunohistochemistry, tumor cells showed diffuse caldesmon (B) and diffuse desmin positivity (C). By in situ hybridization for high-risk HPV E6/E7 mRNA (detecting 18 HPV genotypes including HPV51), tumor cell showed prominent punctate brown staining, while ad-mixed vasculature was negative (D).
F4
FIGURE 4:
Pathology from patient case 3. A, Histopathologic examination of the tumor showed fascicles of spindle cells with focal myxoid areas (inset) (hematoxylin and eosin stain). By immunohistochemistry, tumor cells showed focal caldesmon (B) and focal desmin positivity (C). By in situ hybridization for high-risk HPV E6/E7 mRNA (detecting 18 HPV genotypes including HPC 51), tumor cell showed prominent punctate brown staining (D).
F5
FIGURE 5:
Histopathology of the remaining HPV51-associated leiomyosarcomas showed fascicles and sheets of cells that ranged from epithelioid with clear cell features (patient case 4) (A), to epithelioid with focal spindled areas (case 5) (B; inset for spindled area), to predominantly spindled (case 6) (C). Cases also showed marked nuclear pleomorphism and hyperchromasia (case 7) (D) and a confirmed vaginal primary case showed focal myxoid histology (case 8) (E; inset for myxoid material). A subset of cases also showed areas of necrosis (case 9) (F). Higher power examination revealed readily identifiable mitotic figures (blue arrows) and significant nuclear pleomorphism (case 10) (G). A single case showed significant crush artifact, with large cells with markedly elevated nuclear-to-cytoplasmic ratio (case 11) (H) (A–F: H&E stains; G–H: H&E stains).

Three of 3 HPV51-positive LMS tumors for which in-house confirmatory immunohistochemistry was feasible showed variable positivity for both desmin and caldesmon (Patient cases 1 to 3, Table 3), indicating smooth muscle origin. For 10 of 11 HPV51-positive LMS, findings from immunohistochemistry reports are included in Table 3.

In-house in situ hybridization for high-risk HPV E6/E7 mRNA (detecting 18 HPV genotypes including HPV51) was performed on 3 cases (Table 3), all of which revealed diffuse punctate brown staining of the tumor in each case (Figs. 2G, H, 3D, 4D).

DISCUSSION

Within our cohort of 2569 cases of LMS, HPV51-positive tumors represent a genomically and clinically distinct class of LMS, characterized by the complete absence of detectable in TP53 or RB1 mutations, recurrent ATRX (55%) and TSC1 (18%) mutations, the A1 sublineage of HPV51, and a predilection for the female lower reproductive tract as the site of origin. The most common primary site was the vagina, representing a marked enrichment of HPV51-positive cases at this otherwise-rare primary site for LMS. The genomic signature and clinical characteristics of HPV51-positive LMS differed significantly from features of HPV51-negative LMS (Table 1) and even from HPV-positive but HPV51-negative LMS (Figs. 1B–D). Thus, although hrHPV types are generally considered pathogenic in carcinomas, our results implicate HPV51, sublineage A1, in the development of this distinct class of sarcoma.

Along similar lines, HPV-positive but HPV51-negative LMS and HPV-negative LMS resembled each other, particularly in their high—and nearly identical—prevalence of TP53 and RB1 mutations. Strikingly, no LMS tumor with either TP53 or RB1 mutations harbored any detectable HPV51 genomic sequences. Taken together, these data support our hypothesis that HPV can be an alternative tumorigenic mechanism to TP53 and RB1 mutations in LMS, specifically implicating HPV51 sublineage A1. Furthermore, we demonstrated direct evidence of high-risk HPV E6/E7 mRNA within the tumor cells of three of three unique HPV51-positive LMS by in situ hybridization (Figs. 2G, H, 3D, 4D), further supporting an oncogenic role of the virus in this new tumor type.

HPV infections have a striking tropism for epithelial cells, so the possibility that these unique lesions represent metaplasia from an epithelial intermediate cannot be entirely excluded. However, for several reasons, we do not favor that these tumors represent carcinomas with mesenchymal differentiation, that is, carcinosarcomas or sarcomatoid/metaplastic carcinomas. First, no areas of epithelial differentiation were identified in any HPV51-positive case on available tissue. Second, while >30% of HPV-related vulvar squamous cell carcinomas exhibit pathogenic alterations in PIK3CA,14 none of our HPV51-positive LMS contained detectable pathogenic PIK3CA alterations. Third, the frequent ATRX alterations that we identified in HPV51-positive LMS were not present in a large case series of HPV-related vulvar carcinomas,14 but have been reported as common recurrent genomic events in LMS.11

The current study adds a third class of oncogenic mechanisms that may participate in LMS. The first and most common mechanism in LMS involves inactivating mutations in the tumor suppressor genes TP53 and RB1 (Fig. 1D, Table 1).11 Second, we recently reported a role for deletions of CDKN2C and the adjacent FAF1 gene amongst TP53-wildtype and especially RB1-wildtype LMS, accounting for 3.0% of all LMS cases (Fig. 1C).25 Here, as a third apparent mechanism, we find a novel role for the HPV51 sublineage A1 in a subset of LMS that are wildtype for both TP53 and RB1 (Fig. 1A). This presumptive mechanism is through HPV oncoproteins E6 and E7, which are known to inactivate the p53 and RB1 proteins, respectively.12 Each of these three distinct classes of molecular participants in different LMS tumors may have therapeutic and prognostic implications that will require further study. A strength of the current study is that we were able to exclude cases with molecular alterations that are considered defining for sarcomas other than LMS.

HPV51 was first described in 198826 and is now classified as a hrHPV type.12,13,27 High-risk HPV genotypes, such as the more frequent HPV16 and HPV18, have been associated mainly with the development of carcinomas, particularly in the cervical, anogenital, and oropharyngeal regions.12–14 By convention, distinct HPV types are defined by the sequence of the open reading frame that encodes the major capsid protein L1, with a new HPV type defined as differing from all others by at least 10% within this sequence.22 While different hrHPV types can show different biological activities,28,29 their oncogenicity results primarily from oncoproteins E6 and E7, as noted above.12 High-risk oncogenic HPV genotypes are all clustered in one clade composed of species α5, α6, α7, α9, and α11; HPV51 is a member of α5, which consists of types HPV26, HPV51, HPV69, and HPV82.22 Specific HPV genotype lineages (defined as 1% to <10% genomic differences in the L1 open reading frame) and sublineages (defined as 0.5% to <1% differences), such as HPV16 sublineage D2,30 have been associated with increased risk of malignancy. HPV51 exhibits considerable genomic variability31–33 and contains sublineages A1-A4 and B1-B2.31 No previous study has linked any specific HPV51 lineages to increased cancer risk or to a specific cancer type. For 10 of 11 of the HPV51-positive LMS in our cohort, HPV51 sequences matched to sublineage A1. In a large study of HPV-positive cervical swabs, sublineage A1 was likewise the most common HPV51 sublineage, accounting for 73% of HPV51 infections,34 but was not linked to any subsequent cancers.

The discovery of HPV as a viral carcinogen led to the development of effective vaccines. Nevertheless, HPV51 is not directly targeted by current vaccines. Data remain inconclusive on cross-protection, that is, the development of immunity to HPV51 after immunization with vaccines that do not specifically target HPV51. One study analyzed several trial cohorts of the HPV 16/18 vaccine, and the data suggested cross-protection against HPV51-associated cervical intraepithelial neoplasia.35 Additional studies, however, readily identified new HPV51 infections in vaccinated women, with 1 series indicating that HPV51 was the most frequent HPV type identified in new infections (1.38/100 person-year).36,37

Importantly, the prevalence of HPV51 differs widely amongst different populations, with particularly high HPV51 infection rates in developing countries.38–40 Percentages of HPV51 infections range from 3.0% (Jilin province, China) to 31.4% (rural Brazil) and 30.7% (Adama Town, Ethiopia) of hrHPV-positive cervical smears, depending on the population.38–40 Unfortunately, high rates of nonvaccine hrHPV infections may translate into nonvaccine-related advanced tumors.41 Of note, 46% of HPV51-positive LMS in our cohort were in individuals of African or ad-mixed American ancestry (Table 2).23

HPV51 in routine cervical smears and in neoplasia has frequently been detected in co-infection with other HPV types.42–44 One recent study identified that biopsies of cervical intraepithelial neoplasia grade 2 to 3 that are infected with HPV51 often show co-infection with other HPV genotypes, with only 15/79 (19%) of HPV51-positive cases showing infection with no other HPV type.42 The authors of that study concluded that it was unknown whether the inclusion of HPV51 in a new vaccine would improve protection for cervical neoplasia.43 In contrast, in our cohort of HPV51-positive LMS, HPV51 was the sole HPV genotype identified in 10 of 11 cases (90.9%).

To our knowledge, although incidental HPV has been previously identified in some sarcoma biopsy specimens,45 no HPV type has been linked to sarcoma by in situ studies. In models of sarcoma in vitro and in animals, hrHPV types 16 and 18 have been shown to transform smooth muscle cells and fibroblasts.46–50 HPV E6 and E7 proteins have also been identified in smooth muscle cells of aortic tissue and coronary arteries with atherosclerosis51,52 although their pathogenic effects, if any, in those locations remain unknown. In a few reports, HPV genomic sequences were occasionally identified in rare carcinomas with mesenchymal differentiation, but these were not true sarcomas, and HPV51 was not identified in any of them.53–58 Uncommon hrHPV genotypes, particularly HPV33, have also been previously identified in a rare cancer type with a characteristic location and histology—namely, HPV-related multiphenotypic sinonasal carcinoma. This distinct sinonasal neoplasm shows characteristic myoepithelial differentiation, frequent surface epithelium involvement, and occasional descriptions of sarcomatoid differentiation,59 but not a true sarcoma.

Primary vulvovaginal LMS is rare, with an aggressive course and a high recurrence rate.3–5 A possible link to myxoid histology has been proposed,60–62 and in our study, HPV51-positive LMS with focal myxoid areas were from only vaginal primary sites and accounted for 3 of the 5 vaginal primary HPV51-positive LMS (Table 3). A substantial percentage of confirmed vulvovaginal primary GYN LMS in our cohort were HPV51-positive (5/12; 42%). An additional single separate vaginal primary case from our cohort was HPV30-positive, a rare HPV type of undetermined risk.21 Our study may provide insights into the tumor biology of LMS arising from this rare primary tumor site. A strength of our study is that we provide details of clinical, microscopic, and available immunohistochemical parameters, as well as photomicrographs of each of the 11 key cases, plus detailed genomic data, for comparison to guide future studies. If clinically indicated, HPV51 testing, such as by clinically available RNA in situ hybridization (Figs. 2G, H, 3D, 4D), may be warranted for malignant smooth muscle tumors of the lower female reproductive tract or in metastatic LMS of unknown origin, to aid the determination of site of origin.

We identified a few additional HPV genotypes in other rare LMS: in our cohort, nearly all cases with HPV reads that were not HPV51 harbored TP53/RB1 mutations (Fig. 1D) or were CDKN2C-null (Fig. 1C),25 suggesting that the viral reads identified in these cases were an incidental finding and not a driver of tumorigenesis. Supporting this interpretation, we found no statistically significant differences in oncogenic mutational signatures between HPV-positive but HPV51-negative LMS cases versus HPV-negative LMS cases. In addition, no predilection for the female lower reproductive tract was identified in HPV-positive but HPV51-negative LMS cases. In contrast, the single HPV30, 31, and 33 cases in our LMS cohort were TP53/RB1-wildtype (Fig. 1B) and may therefore represent mimics of HPV51, but additional LMS cases with those 3 HPV genotypes would be needed before inferring an oncogenic role.

Limitations of this study include its retrospective nature and the enrichment for aggressive tumors, mostly metastatic to distant sites, which presumably reflects collection bias from submission of specimens later in the course of disease. In addition, we reviewed H&E scanned slides and accompanying pathology reports from all cases, and we directly reviewed immunohistochemical slides on 3 of 3 evaluable cases shown in Figures 2–4. For other TP53/RB1-wildtype, HPV51-positive cases, we relied on submitted diagnoses without in-house immunohistochemical confirmation, which is a study limitation. Lastly, confirmatory detection of HPV mRNA by in situ hybridization was available for three tumors from three different patients of the 11 cases in which we identified HPV51 genomic sequences; 3 of 3 of those tumors contained detectable HPV E6/E7 mRNA.

Additional studies will be needed to determine if the presence of HPV51 in LMS affects prognosis and approaches or responses to treatment. Nevertheless, it is tantalizing to speculate that HPV51-associated GYN LMS may respond to immunotherapy despite low TMB and absence of MSI, similarly to HPV-positive cervical and head and neck carcinomas. If clinical relevance is confirmed, testing of all vaginal primary LMS by clinically available in situ hybridization to find HPV51 genomic sequences may be indicated. Future studies of interest also include characterizing the gene expression profiles of HPV51-associated LMS tumors,63,64 in comparison with LMS tumors that show mutations in TP53 and/or RB111 or deletion of CDKN2C.24 Taking our findings into consideration, it will also be important to resolve whether HPV51 should be included in new vaccines to prevent HPV51-associated LMS and possibly other tumors.

REFERENCES

1. Roberts ME, Aynardi JT, Chu CS. Uterine leiomyosarcoma: a review of the literature and update on management options. Gynecol Oncol. 2018;151:562–572.
2. Rahaman J, Cohen CJ Kufe DW, Pollock RE, Weichselbaum RR, Bast RC, Gansler TS, Holland JF, Frei E. Chapter 120: gynecologic sarcomas. Holland-Frei Cancer Medicine, 6th ed. Hamilton, ON: BC Decker; 2003.
3. Aartsen EJ, Albus-Lutter CE. Vulvar sarcoma: clinical implications. Eur J Obstet Gynecol Reprod Biol. 1994;56:181–189.
4. Keller NA, Godoy H. Leiomyosarcoma of the vagina: an exceedingly rare diagnosis. Case Rep Obstet Gynecol. 2015;2015:1–4.
5. González-Bugatto F, Añón-Requena MJ, López-Guerrero MA, et al. Vulvar leiomyosarcoma in Bartholin’s gland area: a case report and literature review. Arch Gynecol Obstet. 2009;279:171–174.
6. Seagle BLL, Sobecki-Rausch J, Strohl AE, et al. Prognosis and treatment of uterine leiomyosarcoma: A National Cancer Database study. Gynecol Oncol. 2017;145:61–70.
7. Hensley ML, Barrette BA, Baumann K, et al. Gynecologic Cancer InterGroup (GCIG) consensus review: uterine and ovarian leiomyosarcomas. Int J Gynecol Cancer. 2014;24:S61–S66.
8. Horn LC, Schmidt D, Fathke C, et al. New FIGO staging for uterine sarcomas. Pathologe. 2009;30:302–303.
9. Amin MB, Edge SB, Greene FL, et al. AJCC Cancer Staging System, 8th ed. New York, NY: Spinger; 2017.
10. Wang Z, Shi N, Naing A, et al. Survival of patients with metastatic leiomyosarcoma: the MD Anderson Clinical Center for targeted therapy experience. Cancer Med. 2016;5:3427–3444.
11. Chudasama P, Mughal SS, Sanders MA, et al. Integrative genomic and transcriptomic analysis of leiomyosarcoma. Nat Commun. 2018;9:144.
12. Schiffman M, Doorbar J, Wentzensen N, et al. Carcinogenic human papillomavirus infection. Nat Rev Dis Prim. 2016;2:16086.
13. Armstrong B, Baverstock K, Brenner D, et al. IARC Working Group on the Evaluation of Carcinogenic Risks to Humans. Lyon: International Agency for Research on Cancer; 2012.
14. Williams EA, Werth AJ, Sharaf R, et al. Vulvar squamous cell carcinoma: comprehensive genomic profiling of HPV+ versus HPV– forms reveals distinct sets of potentially actionable molecular targets. JCO Precis Oncol. 2020;4:647–661.
15. Frampton GM, Fichtenholtz A, Otto GA, et al. Development and validation of a clinical cancer genomic profiling test based on massively parallel DNA sequencing. Nat Biotechnol. 2013;31:1023–1031.
16. He J, Abdel-Wahab O, Nahas MK, et al. Integrated genomic DNA/RNA profiling of hematologic malignancies in the clinical setting. Blood. 2016;127:3004–3014.
17. Sun JX, He Y, Sanford E, et al. A computational approach to distinguish somatic vs. germline origin of genomic alterations from deep sequencing of cancer specimens without a matched normal. PLoS Pathog. 2018;14:e1005965.
18. Chalmers ZR, Connelly CF, Fabrizio D, et al. Analysis of 100,000 human cancer genomes reveals the landscape of tumor mutational burden. Genome Med. 2017;9:34.
19. Trabucco SE, Gowen K, Maund SL, et al. A novel next-generation sequencing approach to detecting microsatellite instability and pan-tumor characterization of 1000 microsatellite instability-high cases in 67,000 patient samples. J Mol Diagn. 2019;21:1053–1066.
20. Muñoz N, Bosch FX, De Sanjosé S, et al. Epidemiologic classification of human papillomavirus types associated with cervical cancer. N Engl J Med. 2003;348:518–527.
21. Silva J, Ribeiro J, Sousa H, et al. Oncogenic HPV types infection in adolescents and university women from North Portugal: from self-sampling to cancer prevention. J Oncol. 2011;2011:953469.
22. Chen Z, Schiffman M, Herrero R, et al. Classification and evolution of human papillomavirus genome variants: Alpha-5 (HPV26, 51, 69, 82), Alpha-6 (HPV30, 53, 56, 66), Alpha-11 (HPV34, 73), Alpha-13 (HPV54) and Alpha-3 (HPV61). Virology. 2018;516:86–101.
23. Carrot-Zhang J, Chambwe N, Damrauer JS, et al. Comprehensive analysis of genetic ancestry and its molecular correlates in cancer. Cancer Cell. 2020;37:639.e7–654.e6.
24. Prat J. FIGO staging for uterine sarcomas. Int J Gyn Obst. 2009;104:177–178.
25. Williams EA, Sharaf R, Decker B, et al. CDKN2C-null leiomyosarcoma: a novel, genomically distinct class of TP53/RB1–wild-type tumor with frequent CIC genomic alterations and 1p/19q-codeletion. JCO Precis Oncol. 2020;4:955–971.
26. Nuovo GJ, Crum CP, De Villiers EM, et al. Isolation of a novel human papillomavirus (type 51) from a cervical condyloma. J Virol. 1988;62:1452–1455.
27. Lungu O, Crum CP, Silverstein S. Biologic properties and nucleotide sequence analysis of human papillomavirus type 51. J Virol. 1991;65:4216–4225.
28. Boon SS, Tomaić V, Thomas M, et al. Cancer-causing human papillomavirus E6 proteins display major differences in the phospho-regulation of their PDZ interactions. J Virol. 2015;89:1579–1586.
29. Mischo A, Ohlenschläger O, Ramachandran R, et al. NMR assignment of a PDZ domain in complex with a HPV51 E6 derived N-terminally pyroglutamic acid modified peptide. Biomol NMR Assign. 2013;7:47–49.
30. Mirabello L, Yeager M, Cullen M, et al. HPV16 sublineage associations with histology-specific cancer risk using HPV whole-genome sequences in 3200 women. J Natl Cancer Inst. 2016;108:djw100.
31. Xu J, Tan L, Wang T, et al. Genetic variability of human papillomavirus type 51 E6, E7, L1 and L2 genes in Southwest China. Gene. 2019;690:99–122.
32. Godi A, Epifano I, Bissett SL, et al. Amino acid motifs in both the major and minor capsid proteins of HPV51 impact antigenicity and infectivity. J Gen Virol. 2015;96:1842–1849.
33. Kondo K, Ishii Y, Mori S, et al. Nuclear location of minor capsid protein L2 is required for expression of a reporter plasmid packaged in HPV51 pseudovirions. Virology. 2009;394:259–265.
34. Xi LF, Schiffman M, Koutsky LA, et al. Lineages of oncogenic human papillomavirus types other than type 16 and 18 and risk for cervical intraepithelial neoplasia. J Natl Cancer Inst. 2014;106:dju270.
35. Wheeler CM, Castellsagué X, Garland SM, et al. Cross-protective efficacy of HPV-16/18 AS04-adjuvanted vaccine against cervical infection and precancer caused by non-vaccine oncogenic HPV types: 4-year end-of-study analysis of the randomised, double-blind PATRICIA trial. Lancet Oncol. 2012;13:100–110.
36. McClymont E, Coutlée F, Lee M, et al. Brief report: persistence of non-vaccine oncogenic HPV genotypes in quadrivalent HPV-vaccinated women living with HIV. J Acquir Immune Defic Syndr. 2020;83:230–234.
37. Halfon P, Ravet S, Khiri H, et al. Incident HPV 51 infection after prophylactic quadrivalent human papillomavirus (types 6, 11, 16, and 18) L1 virus-like particle vaccine Gardasil/silgard®. Clin Med Insights Case Rep. 2010;3:69–71.
38. Hao S, Wang C, Liu S, et al. HPV genotypic spectrum in Jilin province, China, where non-vaccine-covered HPV53 and 51 are prevalent, exhibits a bimodal agespecific pattern. PLoS One. 2020;15:e0230640.
39. Haile EL, Cindy S, Ina B, et al. HPV testing on vaginal/cervical nurse-assisted self-samples versus clinician-taken specimens and the HPV prevalence, in Adama Town, Ethiopia. Med (Baltimore). 2019;98:e16970.
40. Lorenzi AT, Fregnani JH, Villa LL, et al. Diversity of human papillomavirus typing among women population living in rural and remote areas of Brazilian territory. Papillomavirus Res. 2019;98:e16970.
41. Williams EA, Newberg J, Williams KJ, et al. Prevalence of high-risk nonvaccine human papillomavirus types in advanced squamous cell carcinoma among individuals of African vs non-African ancestry. JAMA Netw Open. 2021;4:e216481.
42. Piana A, Sotgiu G, Cocuzza C, et al. High HPV-51 prevalence in invasive cervical cancers: results of a pre-immunization survey in North Sardinia, Italy. PLoS One. 2013;8:e63395.
43. Gonzalez-Bosquet E, Gibert M, Serra M, et al. Candidate HPV genotypes not included in the 9-valent vaccine for prevention of CIN 2–3. Int J Gynecol Cancer. 2020;30:954–958.
44. Gallegos-Bolaños J, Rivera-Domínguez JA, Presno-Bernal JM, et al. High prevalence of co-infection between human papillomavirus (HPV) 51 and 52 in Mexican population. BMC Cancer. 2017;17:531.
45. Huang YQ, Li JJ, Rush MG, et al. HPV-16-related DNA sequences in Kaposi’s sarcoma. Lancet. 1992;339:515–518.
46. Perez-Reyes N, Halbert CL, Smith PP, et al. Immortalization of primary human smooth muscle cells. Proc Natl Acad Sci U S A. 1992;89:1224–1228.
47. Zhang Y, Guo X, Zhang Y, et al. Studies on transformation of human esophageal fibroblasts mediated by a retroviral vector containing the E6E7 ORFs of human papillomavirus type 16. Zhonghua Shi Yan He Lin Chuang Bing Du Xue Za Zhi. 2004;18:223–226.
48. Arends MJ, Wyllie AH, Bird CC. Human papillomavirus type 18 is associated with less apoptosis in fibroblast tumours than human papillomavirus type 16. Br J Cancer. 1995;72:646–649.
49. Cornelissen MTE, Smits HL, Briet MA, et al. Uniformity of the splicing pattern of the E6/E7 transcripts in human papillomavirus type 16-transformed human fibroblasts, human cervical premalignant lesions and carcinomas. J Gen Virol. 1990;71:1243–1246.
50. Nishida M, Miyamoto S, Kato H, et al. Transcriptional repression of smooth-muscle α-actin gene associated with human papillomavirus type 16 E7 expression. Mol Carcinog. 1995;13:157–165.
51. Lawson JS, Glenn WK, Tran DD, et al. Identification of human papilloma viruses in atheromatous coronary artery disease. Front Cardiovasc Med. 2015;2:17.
52. Bonin LR, Madden K, Shera K, et al. Generation and characterization of human smooth muscle cell lines derived from atherosclerotic plaque. Arterioscler Thromb Vasc Biol. 1999;19:575–587.
53. Luo Y, Yao Q, Ren C, et al. HPV infection status in cervical metaplastic carcinomas. J Clin Pathol. 2015;68:170–172.
54. Grayson W, Taylor LF, Cooper K. Carcinosarcoma of the uterine cervix: a report of eight cases with immunohistochemical analysis and evaluation of human papillomavirus status. Am J Surg Pathol. 2001;25:338–347.
55. Ikoma S, Nicolas M, Jagirdar J, et al. Chondrosarcoma-like metastasis from a poorly differentiated uterine cervical squamous cell carcinoma. A unique morphology and diagnostic pitfall in cytology. Diagn Cytopathol. 2017;45:750–753.
56. Hickman RA, Bradshaw AD, Cassai N, et al. A rare case of anal carcinosarcoma with human papilloma virus infection in both biphasic tumor elements: an immunohistochemical, molecular and ultrastructural study. Papillomavirus Res. 2016;2:164–166.
57. Poblet E, Pascual A, Godínez JM, et al. Human papillomavirus-associated penile sarcomatoid carcinoma. J Cutan Pathol. 2008;35:559–565.
58. Velazquez EF, Melamed J, Barreto JE, et al. Sarcomatoid carcinoma of the penis: a clinicopathologic study of 15 cases. Am J Surg Pathol. 2005;29:1152–1158.
59. Bishop JA, Andreasen S, Hang J-F, et al. HPV-related multiphenotypic sinonasal carcinoma HPV-related carcinoma with adenoid cystic carcinoma-like features. Am J Surg Pathol. 2017;41:1690–1701.
60. Di Gilio AR, Cormio G, Resta L, et al. Rapid growth of myxoid leiomyosarcoma of the vulva during pregnancy: a case report. Int J Gynecol Cancer. 2004;14:172–175.
61. Tjalma WAA, Colpaert CGA. Myxoid leiomyosarcoma of the vulva. Gynecol Oncol. 2005;96:548–551.
62. Mowers EL, Shank JJ, Frisch N, et al. Myxoid leiomyosarcoma of the Bartholin gland. Obstet Gynecol. 2014;124:433–435.
63. Guo X, Jo VY, Mills AM, et al. Clinically relevant molecular subtypes in leiomyosarcoma. Clin Cancer Res. 2015;21:3501–3511.
64. Simanaviciene V, Popendikyte V, Gudleviciene Z, et al. Different DNA methylation pattern of HPV16, HPV18 and HPV51 genomes in asymptomatic HPV infection as compared to cervical neoplasia. Virology. 2015;484:227–233.
Keywords:

leiomyosarcoma; vagina; comprehensive genomic profiling; HPV; ATRX

Copyright © 2022 The Author(s). Published by Wolters Kluwer Health, Inc.