Development of Spondyloarthritis After COVID-19 in HLA-B27-Positive Monozygotic Twins: Case Reports With Single Cell Transcriptome Profiling.

A subset of spondyloarthritis (SpA) called 'reactive arthritis' is triggered by causal pathogens, usually bacteria related to venereal disease or gastrointestinal infection. During the outbreak of coronavirus disease 2019 (COVID-19), there have been case reports about SpA after COVID-19, but the causality is still elusive. We described cases of 23-year-old monozygotic twins both diagnosed with SpA after COVID-19. The probable linkage between SpA and COVID-19 was elaborated with our cases as well as literature reviews. Of note, shared genetic traits by monozygotic twins, particularly HLA-B27 positivity, might have contributed to their susceptibility to COVID-19-induced SpA. Moreover, single-cell transcriptome analysis revealed that the transcriptomic profile of peripheral compartment of SpA after COVID-19 was distinctive from that of typical radiographic axial SpA as shown by differential expression of ribosomal protein S26 (RPS26) and small nucleolar RNA host gene 5 (SNHG5) in nearly all subsets of peripheral blood mononuclear cells.

Differential effects of periodontal microbiome on the rheumatoid factor induction during rheumatoid arthritis pathogenesis.

Association between exposure to periodontal bacteria and development of autoantibodies related to rheumatoid arthritis (RA) has been widely accepted; however, direct causal relationship between periodontal bacteria and rheumatoid factor (RF) is currently not fully understood. We investigated whether periodontal bacteria could affect RF status. Patients with preclinical, new-onset, or chronic RA underwent periodontal examination, and investigation of subgingival microbiome via 16S rRNA sequencing. Degree of arthritis and RF induction was examined in collagen-induced arthritis (CIA) mice that were orally inoculated with different periodontal bacteria species. Subsequently, single-cell RNA sequencing analysis of the mouse spleen cells was performed. Patients with preclinical RA showed an increased abundance of the Porphyromonadacae family in the subgingival microbiome compared to those with new-onset or chronic RA, despite comparable periodontitis severity among them. Notably, a distinct subgingival microbial community was found between patients with high-positive RF and those with negative or low-positive RF (p=0.022). Oral infections with the periodontal pathogens P. gingivalis and Treponema denticola in CIA mice similarly enhanced arthritis score, but resulted in different levels of RF induction. Genes related to B cell receptor signaling, B cell proliferation, activation, and differentiation, and CD4+ T cell costimulation and cytokine production were involved in the differential induction of RF in mice exposed to different bacteria. In summary, periodontal microbiome might shape RF status by affecting the humoral immune response during RA pathogenesis.

An efficient and tunable parameter to improve variant calling for whole genome and exome sequencing data.

Next generation sequencing (NGS) has traditionally been performed in various fields including agricultural to clinical and there are so many sequencing platforms available in order to obtain accurate and consistent results. However, these platforms showed amplification bias when facilitating variant calls in personal genomes. Here, we sequenced whole genomes and whole exomes from ten Korean individuals using Illumina and Ion Proton, respectively to find the vulnerability and accuracy of NGS platform in the GC rich/poor area. Overall, a total of 1013 Gb reads from Illumina and ~39.1 Gb reads from Ion Proton were analyzed using BWA-GATK variant calling pipeline. Furthermore, conjunction with the VQSR tool and detailed filtering strategies, we achieved high-quality variants. Finally, each of the ten variants from Illumina only, Ion Proton only, and intersection was selected for Sanger validation. The validation results revealed that Illumina platform showed higher accuracy than Ion Proton. The described filtering methods are advantageous for large population-based whole genome studies designed to identify common and rare variations associated with complex diseases.

Genomic analysis reveals secondary glioblastoma after radiotherapy in a subset of recurrent medulloblastomas.

Despite great advances in understanding of molecular pathogenesis and achievement of a high cure rate in medulloblastoma, recurrent medulloblastomas are still dismal. Additionally, misidentification of secondary malignancies due to histological ambiguity leads to misdiagnosis and eventually to inappropriate treatment. Nevertheless, the genomic characteristics of recurrent medulloblastomas are poorly understood, largely due to a lack of matched primary and recurrent tumor tissues. We performed a genomic analysis of recurrent tumors from 17 pediatric medulloblastoma patients. Whole transcriptome sequencing revealed that a subset of recurrent tumors initially diagnosed as locally recurrent medulloblastomas are secondary glioblastomas after radiotherapy, showing high similarity to the non-G-CIMP proneural subtype of glioblastoma. Further analysis, including whole exome sequencing, revealed missense mutations or complex gene fusion events in PDGFRA with augmented expression in the secondary glioblastomas after radiotherapy, implicating PDGFRA as a putative driver in the development of secondary glioblastomas after treatment exposure. This result provides insight into the possible application of PDGFRA-targeted therapy in these second malignancies. Furthermore, genomic alterations of TP53 including 17p loss or germline/somatic mutations were also found in most of the secondary glioblastomas after radiotherapy, indicating a crucial role of TP53 alteration in the process. On the other hand, analysis of recurrent medulloblastomas revealed that the most prevalent alterations are the loss of 17p region including TP53 and gain of 7q region containing EZH2 which already exist in primary tumors. The 7q gain events are frequently accompanied by high expression levels of EZH2 in both primary and recurrent medulloblastomas, which provides a clue to a new therapeutic target to prevent recurrence. Considering the fact that it is often challenging to differentiate between recurrent medulloblastomas and secondary glioblastomas after radiotherapy, our findings have major clinical implications both for correct diagnosis and for potential therapeutic interventions in these devastating diseases.

Genetic Progression of High Grade Prostatic Intraepithelial Neoplasia to Prostate Cancer.

BACKGROUND: Although high grade prostatic intraepithelial neoplasia (HGPIN) is considered a neoplastic lesion that precedes prostate cancer (PCA), the genomic structures of HGPIN remain unknown. OBJECTIVE: Identification of the genomic landscape of HGPIN and the genomic differences between HGPIN and PCA that may drive the progression to PCA. DESIGN, SETTINGS, AND PARTICIPANTS: We analyzed 20 regions of paired HGPIN and PCA from six patients using whole-exome sequencing and array-comparative genomic hybridization. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Somatic mutation and copy number alteration (CNA) profiles of paired HGPIN and PCA were measured and compared. RESULTS AND LIMITATIONS: The number of total mutations and CNAs of HGPINs were significantly fewer than those of PCAs. Mutations in FOXA1 and CNAs (1q and 8q gains) were detected in both HGPIN and PCA ('common'), suggesting their roles in early PCA development. Mutations in SPOP, KDM6A, and KMT2D were 'PCA-specific', suggesting their roles in HGPIN progression to PCA. The 8p loss was either 'common' or 'PCA-specific'. In-silico estimation of evolutionary ages predicted that HGPIN genomes were much younger than PCA genomes. Our data show that PCAs are direct descendants of HGPINs in most cases that require more genomic alterations to progress to PCA. The nature of heterogeneous HGPIN population that might attenuate genomic signals should further be studied. CONCLUSIONS: HGPIN genomes harbor relatively fewer mutations and CNAs than PCA but require additional hits for the progression. PATIENT SUMMARY: In this study, we suggest a systemic diagram from high grade prostatic intraepithelial neoplasia (HGPIN) to prostate cancer (PCA). Our results provide a clue to explain the long latency from HGPIN to PCA and provide useful information for the genetic diagnosis of HGPIN and PCA.

Genomic landscape of endometrial stromal sarcoma of uterus.

Although recurrent gene fusions such as JAZF1-JJAZ1 are considered driver events for endometrial stromal sarcoma (ESS) development, other genomic alterations remain largely unknown. In this study, we performed whole-exome sequencing, transcriptome sequencing and copy number profiling for five ESSs (three low-grade ESS (LG-ESS) and two undifferentiated uterine sarcomas (UUSs)). All three LG-ESSs exhibited either one of JAZF1-SUZ12, JAZF1-PHF1 and MEAF6-PHF1 fusions, whereas the two UUSs did not. All ESSs except one LG-ESS exhibited copy number alterations (CNAs), many of which encompassed cancer-related genes. In UUSs, five CNAs encompassing cancer-related genes (EZR, CDH1, RB1, TP53 and PRKAR1A) accompanied their expressional changes, suggesting that they might stimulate UUS development. We found 81 non-silent mutations (35 from LG-ESSs and 46 from UUSs) that included 15 putative cancer genes catalogued in cancer-related databases, including PPARG and IRF4 mutations. However, they were non-recurrent and did not include any well-known mutations, indicating that point mutations may not be a major driver for ESS development. Our data show that gene fusions and CNAs are the principal drivers for LG-ESS and USS, respectively, but both may require additional genomic alterations including point mutations. These differences may explain the different biologic behaviors between LG-ESS and UUS. Our findings suggest that ESS development requires point mutations and CNAs as well as the gene fusions.

Clonal origins and parallel evolution of regionally synchronous colorectal adenoma and carcinoma.

Although the colorectal adenoma-to-carcinoma sequence represents a classical cancer progression model, the evolution of the mutational landscape underlying this model is not fully understood. In this study, we analyzed eight synchronous pairs of colorectal high-grade adenomas and carcinomas, four microsatellite-unstable (MSU) and four-stable (MSS) pairs, using whole-exome sequencing. In the MSU adenoma-carcinoma pairs, we observed no subclonal mutations in adenomas that became fixed in paired carcinomas, suggesting a 'parallel' evolution of synchronous adenoma-to-carcinoma, rather than a 'stepwise' evolution. The abundance of indel (in MSU and MSS pairs) and microsatellite instability (in MSU pairs) was noted in the later adenoma- or carcinoma-specific mutations, indicating that the mutational processes and functional constraints operative in early and late colorectal carcinogenesis are different. All MSU cases exhibited clonal, truncating mutations in ACVR2A, TGFBR2, and DNA mismatch repair genes, but none were present in APC or KRAS. In three MSS pairs, both APC and KRAS mutations were identified as both early and clonal events, often accompanying clonal copy number changes. An MSS case uniquely exhibited clonal ERBB2 amplification, followed by APC and TP53 mutations as carcinoma-specific events. Along with the previously unrecognized clonal origins of synchronous colorectal adenoma-carcinoma pairs, our study revealed that the preferred sequence of mutational events during colorectal carcinogenesis can be context-dependent.

Subclonal Genomic Architectures of Primary and Metastatic Colorectal Cancer Based on Intratumoral Genetic Heterogeneity.

PURPOSE: The intratumoral heterogeneity (ITH) and the evolution of genomic architectures associated with the development of distant metastases are not well understood in colorectal cancers. EXPERIMENTAL DESIGN: We performed multiregion biopsies of primary and liver metastatic regions from five colorectal cancers with whole-exome sequencing and copy number profiling. RESULTS: In addition to a substantial level of genetic ITH, multiregion genetic profiling identifies the subclonal mutational architecture, leading to the region-based or spatial categorization of somatic mutations and the inference of intratumoral evolutionary history of cancers. The universal mutations (those observed in all the regional biopsies) are enriched in known cancer genes such as APC and TP53 with distinct mutational spectra compared with biopsy- or region-specific mutations, suggesting that major operative mutational mechanisms and their selective pressures are not constant across the metastatic progression. The phylogenies inferred from genomic data show branching evolutionary patterns where some primary biopsies are often segregated with metastastic lesions. Our analyses also revealed that copy number changes such as the chromosomal gains of c-MYC and chromothripsis can be region specific and the potential source of genetic ITH. CONCLUSIONS: Our data show that the genetic ITH is prevalent in colorectal cancer serving as a potential driving force to generate metastasis-initiating clones and also as a means to infer the intratumoral evolutionary history of cancers. The paucity of recurrent metastasis-clonal events suggests that colorectal cancer distant metastases may not follow a uniform course of genomic evolution, which should be considered in the genetic diagnosis and the selection of therapeutic targets for the advanced colorectal cancer.

Genomic differences between pure ductal carcinoma in situ and synchronous ductal carcinoma in situ with invasive breast cancer.

Although ductal carcinoma in situ (DCIS) precedes invasive ductal carcinoma (IDC), the related genomic alterations remain unknown. To identify the genomic landscape of DCIS and better understand the mechanisms behind progression to IDC, we performed whole-exome sequencing and copy number profiling for six cases of pure DCIS and five pairs of synchronous DCIS and IDC. Pure DCIS harbored well-known mutations (e.g., TP53, PIK3CA and AKT1), copy number alterations (CNAs) and chromothripses, but had significantly fewer driver genes and co-occurrence of mutation/CNAs than synchronous DCIS-IDC. We found neither recurrent nor significantly mutated genes with synchronous DCIS-IDC compared to pure DCIS, indicating that there may not be a single determinant for pure DCIS progression to IDC. Of note, synchronous DCIS genomes were closer to IDC than pure DCIS. Among the clinicopathologic parameters, progesterone receptor (PR)-negative status was associated with increased mutations, CNAs, co-occurrence of mutations/CNAs and driver mutations. Our results indicate that although pure DCIS has already acquired some drivers, more changes are needed to progress to IDC. In addition, IDC-associated DCIS is more aggressive than pure DCIS at genomic level and should really be considered IDC. Finally, the data suggest that PR-negativity could be used to predict aggressive breast cancer genotypes.

Progression of naive intraepithelial neoplasia genome to aggressive squamous cell carcinoma genome of uterine cervix.

Although cervical intraepithelial neoplasia (CIN) is considered a neoplasia, its genomic alterations remain unknown. For this, we performed whole-exome sequencing and copy number profiling of three CINs, a microinvasive carcinoma (MIC) and four cervical squamous cell carcinomas (CSCC). Both total mutation and driver mutation numbers of the CINs were significantly fewer than those of the MIC/CSCCs (P = 0.036 and P = 0.018, respectively). Importantly, PIK3CA was altered in all MIC/CSCCs by either mutation or amplification, but not in CINs. The CINs harbored significantly lower numbers of copy number alterations (CNAs) than the MIC/CSCCs as well (P = 0.036). Pathway analysis predicted that the MIC/CSCCs were enriched with cancer-related signalings such as cell adhesion, mTOR signaling pathway and cell migration that were depleted in the CINs. The mutation-based estimation of evolutionary ages identified that CIN genomes were younger than MIC/CSCC genomes. The data indicate that CIN genomes harbor unfixed mutations in addition to human papilloma virus infection but require additional driver hits such as PIK3CA, TP53, STK11 and MAPK1 mutations for CSCC progression. Taken together, our data may explain the long latency from CIN to CSCC progression and provide useful information for molecular diagnosis of CIN and CSCC.

The mutational burdens and evolutionary ages of early gastric cancers are comparable to those of advanced gastric cancers.

Early gastric cancers (EGCs) precede advanced gastric cancers (AGCs), with a favourable prognosis compared to AGC. To understand the progression mechanism of EGC to AGC, it is required to disclose EGC and AGC genomes in mutational and evolutionary perspectives. We performed whole-exome sequencing and copy number profiling of nine microsatellite (MS)-unstable (MSI-H) (five EGCs and four AGCs) and eight MS-stable (MSS) gastric cancers (four EGCs and four AGCs). In the cancers, we observed well-known driver mutations (TP53, APC, PIK3CA, ARID1A, and KRAS) that were enriched in cancer-related pathways, including chromatin remodelling and tyrosine kinase activity. The MSI-H genomes harboured ten times more mutations, but were largely depleted of copy number alterations (CNAs) compared to the MSS cancers. Interestingly, EGC genomes showed a comparable level of mutations to AGC in terms of the number, sequence composition, and functional consequences (potential driver mutations and affected pathways) of mutations. Furthermore, the CNAs between EGC and AGC genomes were not significantly different in either MSI-H and MSS. Evolutionary analyses using somatic mutations and MSI as molecular clocks further identified that EGC genomes were as old as AGC genomes in both MSS and MSI-H cancers. Our results suggest that the genetic makeup for gastric cancer may already be achieved in EGC genomes and that the time required for transition to AGC may be relatively short. Also, the data suggest a possibility that the mutational profiles obtained from early biopsies may be useful in the clinical settings for the molecular diagnosis and therapeutics of gastric cancer patients.

Regional biases in mutation screening due to intratumoural heterogeneity of prostate cancer.

Intratumoural heterogeneity (ITH) leads to regional biases of the mutational landscape in a single tumour and may influence the single biopsy-based clinical diagnosis and treatment decision. To evaluate the extent of ITH in unifocal prostate cancers (PCAs), we analysed multiple regional biopsies from three PCAs, using whole-exome sequencing, DNA copy number and gene expression profiling analyses. A substantial level of ITH was identified, in that 0-61% and 18-71% of somatic variants were common or private, respectively, within a given cancer. The enhanced mutation detection rate in the combined sequencing dataset across intratumoural biopsies was demonstrated with respect to the total number of mutations identified in a given tumour. Allele frequencies of the mutations were positively correlated with the levels of intratumoural recurrence (private < shared < common), but some common mutations showed low allele frequency, suggesting that not all were clonally fixed. Regional biases in the presentation of a well-known TMPRSS2-ERG fusion was noted in one PCA and the somatic mutation- and copy number-based phylogenetic relationships between intratumoural biopsies were largely concordant. Genes showing intratumoural expression variability were commonly enriched in the molecular function of eicosanoid metabolism and PCA-relevant clinical markers. Taken together, our analyses identified a substantial level of genetic ITH in unifocal PCAs at the mutation, copy number and expression levels, which should be taken into account for the identification of biomarkers in the clinical setting.

MAP: Mutation Arranger for Defining Phenotype-Related Single-Nucleotide Variant.

Next-generation sequencing (NGS) is widely used to identify the causative mutations underlying diverse human diseases, including cancers, which can be useful for discovering the diagnostic and therapeutic targets. Currently, a number of single-nucleotide variant (SNV)-calling algorithms are available; however, there is no tool for visualizing the recurrent and phenotype-specific mutations for general researchers. In this study, in order to support defining the recurrent mutations or phenotype-specific mutations from NGS data of a group of cancers with diverse phenotypes, we aimed to develop a user-friendly tool, named mutation arranger for defining phenotype-related SNV (MAP). MAP is a user-friendly program with multiple functions that supports the determination of recurrent or phenotype-specific mutations and provides graphic illustration images to the users. Its operation environment, the Microsoft Windows environment, enables more researchers who cannot operate Linux to define clinically meaningful mutations with NGS data from cancer cohorts.