Single-Cell and Spatial Transcriptome Analysis of Dermal Fibroblast Development in Perinatal Mouse Skin: Dynamic Lineage Differentiation and Key Driver Genes

Several single-cell RNA studies of developing mouse skin have elucidated the molecular and cellular processes involved in skin development. However, they have primarily focused on either the fetal or early postnatal period, leaving a gap in our understanding of skin development. In this study, we conducted a comprehensive time-series analysis by combining single-cell RNA-sequencing datasets collected at different stages of development (embryonic days 13.5, 14.5, and 16.5 and postnatal days 0, 2, and 4) and validated our findings through multipanel in situ spatial transcriptomics. Our analysis indicated that embryonic fibroblasts exhibit heterogeneity from a very early stage and that the rapid determination of each lineage occurs within days after birth. The expression of putative key driver genes, including Hey1, Ebf1, Runx3, and Sox11 for the dermal papilla trajectory; Lrrc15 for the dermal sheath trajectory; Zfp536 and Nrn1 for the papillary fibroblast trajectory; and Lrrn4cl and Mfap5 for the fascia fibroblast trajectory, was detected in the corresponding, spatially identified cell types. Finally, cell-to-cell interaction analysis indicated that the dermal papilla lineage is the primary source of the noncanonical Wnt pathway during skin development. Together, our study provides a transcriptomic reference for future research in the field of skin development and regeneration.

RNA expression of 6 genes from metastatic mucosal gastric cancer serves as the global prognostic marker for gastric cancer with functional validation

BACKGROUND

Molecular analysis of advanced tumors can increase tumor heterogeneity and selection bias. We developed a robust prognostic signature for gastric cancer by comparing RNA expression between very rare early gastric cancers invading only mucosal layer (mEGCs) with lymph node metastasis (Npos) and those without metastasis (Nneg).

METHODS

Out of 1003 mEGCs, all Npos were matched to Nneg using propensity scores. Machine learning approach comparing Npos and Nneg was used to develop prognostic signature. The function and robustness of prognostic signature was validated using cell lines and external datasets.

RESULTS

Extensive machine learning with cross-validation identified the prognostic classifier consisting of four overexpressed genes (HDAC5, NPM1, DTX3, and PPP3R1) and two downregulated genes (MED12 and TP53), and enabled us to develop the risk score predicting poor prognosis. Cell lines engineered to high-risk score showed increased invasion, migration, and resistance to 5-FU and Oxaliplatin but maintained sensitivity to an HDAC inhibitor. Mouse models after tail vein injection of cell lines with high-risk score revealed increased metastasis. In three external cohorts, our risk score was identified as the independent prognostic factor for overall and recurrence-free survival.

CONCLUSION

The risk score from the 6-gene classifier can successfully predict the prognosis of gastric cancer.

Exploring the role of apolipoprotein E gene promoter polymorphisms in susceptibility to normal-tension glaucoma in a Korean population

Glaucoma, particularly primary open-angle glaucoma (POAG), poses a significant global health concern. Distinguished by intraocular pressure (IOP), POAG encompasses high-tension glaucoma (HTG) and normal-tension glaucoma (NTG). Apolipoprotein E (APOE) is a multifaceted protein with roles in lipid metabolism, neurobiology, and neurodegenerative diseases. However, controversies persist regarding the impact of APOE single-nucleotide polymorphisms (SNPs) on open-angle glaucoma and NTG. This study aimed to identify APOE-specific SNPs influencing NTG risk. A cohort of 178 patients with NTG recruited from Uijeongbu St. Mary's Hospital and 32,858 individuals from the Korean Genome and Epidemiology Study (KoGES) cohort were included in the analysis. Genotype and haplotype analyses were performed on three promoter SNPs (rs449647, rs769446, and rs405509) and two exonic SNPs (rs429358 and rs7412) located on chromosome 19. Among the five SNPs, rs769446 genotypes exhibited significant differences between cases and controls. The minor allele C of rs769446 emerged as a protective factor against NTG. Furthermore, haplotype analysis of the five SNPs revealed that the A-T-G-T-T haplotype was a statistically significant risk factor for NTG. This study indicated an association between APOE promoter SNPs and NTG in the Korean population. Further studies are required to understand how APOE promoter SNPs contribute to NTG pathogenesis.

Placental expression quantitative trait loci in an East Asian population

Expression quantitative trait loci (eQTL) analysis measures the contribution of genetic variation in gene expression on complex traits. Although this methodology has been used to examine gene regulation in numerous human tissues, eQTL research in solid tissues is relatively lacking. We conducted eQTL analysis on placentas collected from an East Asian population in an effort to identify gene regulatory mechanisms in this tissue. Placentas (n = 102) were collected at the time of cesarean delivery. mRNA was extracted, sequenced with NGS, and compared with matched maternal and fetal DNA arrays performed using maternal and neonatal cord blood. Linear regression modeling was performed using tensorQTL. Fine-mapping along with epigenomic annotation was used to select putative functional variants. We identified 2,703 coding genes that contained at least one eQTL with statistical significance (false discovery rate <0.05). After fine-mapping, we found 108 previously unreported eQTL variants with posterior inclusion probability >0.1. Of these, 19% were located in genomic regions with evidence from public placental epigenome suggesting that they may be functionally relevant. For example, variant rs28379289 located in the placenta-specific regulatory region changes the binding affinity of transcription factor leading to higher expression of LGALS3, which is known to affect placental function. This study expands the knowledge base of regulatory elements within the human placenta and identifies 108 previously unreported placenta eQTL signals, which are listed in our publicly available GMI eQTL database. Further studies are needed to identify and characterize genetic regulatory mechanisms that affect placental function in normal pregnancy and placenta-related diseases.

Identification of New Pathogenic Variants of Hereditary Diffuse Gastric Cancer

Purpose

Hereditary diffuse gastric cancer (HDGC) presents a significant genetic predisposition, notably linked to mutations in the CDH1 and CTNNA1. However, the genetic basis for over half of HDGC cases remains unidentified. The aim of this study is to identify novel pathogenic variants in HDGC and evaluate their protein expression.

Materials and Methods

Among 20 qualifying families, two were selected based on available pedigree and DNA. Whole genome sequencing (WGS) on DNA extracted from blood and whole exome sequencing (WES) on DNA from formalin-fixed paraffin-embedded tissues were performed to find potential pathogenic variants in HDGC. After selection of a candidate variant, functional validation and enrichment analysis were performed.

Results

As a result of WGS, three candidate germline mutations (EPHA5, MCOA2, and RHOA) were identified in one family. After literature review and in silico analyses, the RHOA mutation (R129W) was selected as a candidate. This mutation was found in two gastric cancer patients within the family. In functional validation, it showed RhoA overexpression and a higher GTP-bound state in the RhoaR129W mutant. Decreased phosphorylation at Ser127/397 suggested altered YAP1 regulation in the Rho-ROCK pathway. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses linked RhoAR129W overexpression to changed migration/adhesion in MKN1 cell line. However, this RHOA mutation (R129W) was not found in index patients in other families.

Conclusion

The RHOA mutation (R129W) emerges as a potential causative gene for HDGC, but only in one family, indicating a need for further studies to understand its role in HDGC pathogenesis fully.

Central neurocytoma exhibits radial glial cell signatures with FGFR3 hypomethylation and overexpression

We explored the genomic events underlying central neurocytoma (CN), a rare neoplasm of the central nervous system, via multiomics approaches, including whole-exome sequencing, bulk and single-nuclei RNA sequencing, and methylation sequencing. We identified FGFR3 hypomethylation leading to FGFR3 overexpression as a major event in the ontogeny of CN that affects crucial downstream events, such as aberrant PI3K-AKT activity and neuronal development pathways. Furthermore, we found similarities between CN and radial glial cells based on analyses of gene markers and CN tumor cells and postulate that CN tumorigenesis is due to dysregulation of radial glial cell differentiation into neurons. Our data demonstrate the potential role of FGFR3 as one of the leading drivers of tumorigenesis in CN.

Case Series of Soft Tissue Sarcoma Patients with Brain Metastasis with Implications from Genomic and Transcriptomic Analysis

PURPOSE

Brain metastasis rarely occurs in soft tissue sarcoma (STS). Here, we present five cases of STS with brain metastases with genetic profiles.

MATERIALS AND METHODS

We included five patients from Seoul National University Hospital who were diagnosed with STS with metastasis to the brain. Tissue from the brain metastasis along with that from the primary site or other metastases were used for DNA and RNA sequencing to identify genetic profiles. Gene expression profiles were compared with sarcoma samples from The Cancer Genome Atlas.

RESULTS

The overall survival after diagnosis of brain metastasis ranged from 2.2 to 34.3 months. Comparison of mutational profiles between brain metastases and matched primary or other metastatic samples showed similar profiles. In two patients, copy number variation profiles between brain metastasis and other tumors showed several differences including MYCL, JUN, MYC, and DDR2 amplification. Gene ontology analysis showed that the group of genes significantly highly expressed in the brain metastasis samples was enriched in the G-protein coupled receptor activity, structural constituent of chromatin, protein heterodimerization activity, and binding of DNA, RNA, and protein. Gene set enrichment analysis showed enrichment in the pathway of neuroactive ligand-receptor interaction and systemic lupus erythematosus.

CONCLUSION

The five patients had variable ranges of clinical courses and outcomes. Genomic and transcriptomic analysis of STS with brain metastasis implicates possible involvement of complex expression modification and epigenetic changes rather than the addition of single driver gene alteration.

Brain hypothyroidism silences the immune response of microglia in Alzheimer's disease animal model

Thyroid hormone (TH) imbalance is linked to the pathophysiology of reversible dementia and Alzheimer's disease (AD). It is unclear whether tissue hypothyroidism occurs in the AD brain and how it affects on AD pathology. We find that decreased iodothyronine deiodinase 2 is correlated with hippocampal hypothyroidism in early AD model mice before TH alterations in the blood. TH deficiency leads to spontaneous activation of microglia in wild-type mice under nonstimulated conditions, resulting in lowered innate immune responses of microglia in response to inflammatory stimuli or amyloid-β. In AD model mice, TH deficiency aggravates AD pathology by reducing the disease-associated microglia population and microglial phagocytosis. We find that TH deficiency reduces microglial ecto-5'-nucleotidase (CD73) and inhibition of CD73 leads to impaired innate immune responses in microglia. Our findings reveal that TH shapes microglial responses to inflammatory stimuli including amyloid-β, and brain hypothyroidism in early AD model mice aggravates AD pathology by microglial dysfunction.

Potential Perturbations of Critical Cancer-regulatory Genes in Triple-Negative Breast Cancer Cells Within the Humanized Microenvironment of Patient-derived Xenograft Models

PURPOSE

In this study, we aimed to establish humanized patient-derived xenograft (PDX) models for triple-negative breast cancer (TNBC) using cord blood (CB) hematopoietic stem cells (HSCs). Additionally, we attempted to characterize the immune microenvironment of the humanized PDX model to understand the potential implications of altered tumor-immune interactions in the humanized PDX model on the behavior of TNBC cells.

METHODS

To establish a humanized mouse model, high-purity CD34+ HSCs from CB were transplanted into immunodeficient NOD scid γ mice. Peripheral and intratumoral immune cell compositions of humanized and non-humanized mice were compared. Additionally, RNA sequencing of the tumor tissues was performed to characterize the gene expression features associated with humanization.

RESULTS

After transplanting the CD34+ HSCs, CD45+ human immune cells appeared within five weeks. A humanized mouse model showed viable human immune cells in the peripheral blood, lymphoid organs, and in the tumor microenvironment. Humanized TNBC PDX models showed varying rates of tumor growth compared to that of non-humanized mice. RNA sequencing of the tumor tissue showed significant alterations in tumor tissues from the humanized models. tumor necrosis factor receptor superfamily member 11B (TNFRSF11B) is a shared downregulated gene in tumor tissues from humanized models. Silencing of TNFRSF11B in TNBC cell lines significantly reduced cell proliferation, migration, and invasion in vitro. Additionally, TNFRSF11B silenced cells showed decreased tumorigenicity and metastatic capacity in vivo.

CONCLUSION

Humanized PDX models successfully recreated tumor-immune interactions in TNBC. TNFRSF11B, a commonly downregulated gene in humanized PDX models, may play a key role in tumor growth and metastasis. Differential tumor growth rates and gene expression patterns highlighted the complexities of the immune response in the tumor microenvironment of humanized PDX models.

Targeting the m6A RNA methyltransferase METTL3 attenuates the development of kidney fibrosis

Kidney fibrosis is a major mechanism underlying chronic kidney disease (CKD). N6-methyladenosine (m6A) RNA methylation is associated with organ fibrosis. We investigated m6A profile alterations and the inhibitory effect of RNA methylation in kidney fibrosis in vitro (TGF-β-treated HK-2 cells) and in vivo (unilateral ureteral obstruction [UUO] mouse model). METTL3-mediated signaling was inhibited using siRNA in vitro or the METTL3-specific inhibitor STM2457 in vivo and in vitro. In HK-2 cells, METTL3 protein levels increased in a dose- and time-dependent manner along with an increase in the cellular m6A levels. In the UUO model, METTL3 expression and m6A levels were significantly increased. Transcriptomic and m6A profiling demonstrated that epithelial-to-mesenchymal transition- and inflammation-related pathways were significantly associated with RNA m6A methylation. Genetic and pharmacologic inhibition of METTL3 in HK-2 cells decreased TGF-β-induced fibrotic marker expression. STM2457-induced inhibition of METTL3 attenuated the degree of kidney fibrosis in vivo. Furthermore, METTL3 protein expression was significantly increased in the tissues of CKD patients with diabetic or IgA nephropathy. Therefore, targeting alterations in RNA methylation could be a potential therapeutic strategy for treating kidney fibrosis.

Genetic Studies of Actinic Keratosis Development: Where Are We Now?

Actinic keratosis (AK) is a common precancerous skin lesion that can develop into cutaneous squamous cell carcinoma (CSCC). AK is characterized by atypical keratinocytes in the skin's outer layer and is commonly found in sun-exposed areas. Like many precancerous lesions, the development of AK is closely associated with genetic mutations. The molecular biology and transcriptional mechanisms underlying AK development are not well understood. Ultraviolet (UV) light exposure, especially UVA and UVB radiation, is a significant risk factor for AK, causing DNA damage and mutagenic effects. Besides UV exposure, comorbidities like diabetes, rheumatoid arthritis, and psoriasis may also influence AK development. AK patients have shown associations with various internal malignancies, indicating potential vulnerability in cancer-associated genes. Treatment for AK includes cryosurgery, electrodesiccation and curettage, chemotherapeutic creams, photodynamic therapy, or topical immune-modulators. Genomic studies have identified genetic aberrations in AK, with common mutations found in genes like TP53, NOTCH1, and NOTCH2. The progression from AK to CSCC involves chromosomal aberrations and alterations in oncogenes and tumor-suppressor genes. The functional relationships among these genes are not fully understood, but network analysis provides insights into their potential mechanisms. Further research is needed to enhance our understanding of AK's pathogenesis and develop novel therapeutic approaches.

Profiling of Microbial Landscape in Lung of Chronic Obstructive Pulmonary Disease Patients Using RNA Sequencing

Purpose

The aim of the study was to use RNA sequencing (RNA-seq) data of lung from chronic obstructive pulmonary disease (COPD) patients to identify the bacteria that are most commonly detected. Additionally, the study sought to investigate the differences in these infections between normal lung tissues and those affected by COPD.

Patients and Methods

We re-analyzed RNA-seq data of lung from 99 COPD patients and 93 non-COPD smokers to determine the extent to which the metagenomes differed between the two groups and to assess the reliability of the metagenomes. We used unmapped reads in the RNA-seq data that were not aligned to the human reference genome to identify more common infections in COPD patients.

Results

We identified 18 bacteria that exhibited significant differences between the COPD and non-COPD smoker groups. Among these, Yersinia enterocolitica was found to be more than 30% more abundant in COPD. Additionally, we observed difference in detection rate based on smoking history. To ensure the accuracy of our findings and distinguish them from false positives, we double-check the metagenomic profile using Basic Local Alignment Search Tool (BLAST). We were able to identify and remove specific species that might have been misclassified as other species in Kraken2 but were actually Staphylococcus aureus, as identified by BLAST analysis.

Conclusion

This study highlighted the method of using unmapped reads, which were not typically used in sequencing data, to identify microorganisms present in patients with lung diseases such as COPD. This method expanded our understanding of the microbial landscape in COPD and provided insights into the potential role of microorganisms in disease development and progression.

Multidimensional fragmentomic profiling of cell-free DNA released from patient-derived organoids

BACKGROUND

Fragmentomics, the investigation of fragmentation patterns of cell-free DNA (cfDNA), has emerged as a promising strategy for the early detection of multiple cancers in the field of liquid biopsy. However, the clinical application of this approach has been hindered by a limited understanding of cfDNA biology. Furthermore, the prevalence of hematopoietic cell-derived cfDNA in plasma complicates the in vivo investigation of tissue-specific cfDNA other than that of hematopoietic origin. While conventional two-dimensional cell lines have contributed to research on cfDNA biology, their limited representation of in vivo tissue contexts underscores the need for more robust models. In this study, we propose three-dimensional organoids as a novel in vitro model for studying cfDNA biology, focusing on multifaceted fragmentomic analyses.

RESULTS

We established nine patient-derived organoid lines from normal lung airway, normal gastric, and gastric cancer tissues. We then extracted cfDNA from the culture medium of these organoids in both proliferative and apoptotic states. Using whole-genome sequencing data from cfDNA, we analyzed various fragmentomic features, including fragment size, footprints, end motifs, and repeat types at the end. The distribution of cfDNA fragment sizes in organoids, especially in apoptosis samples, was similar to that found in plasma, implying occupancy by mononucleosomes. The footprints determined by sequencing depth exhibited distinct patterns depending on fragment sizes, reflecting occupancy by a variety of DNA-binding proteins. Notably, we discovered that short fragments (< 118 bp) were exclusively enriched in the proliferative state and exhibited distinct fragmentomic profiles, characterized by 3 bp palindromic end motifs and specific repeats.

CONCLUSIONS

In conclusion, our results highlight the utility of in vitro organoid models as a valuable tool for studying cfDNA biology and its associated fragmentation patterns. This, in turn, will pave the way for further enhancements in noninvasive cancer detection methodologies based on fragmentomics.

Development of Novel Focal Irradiation Tool for High-Precision Irradiation Using Clinical Brachytherapy System

PURPOSE/OBJECTIVE(S)

Several small animals, including mice, are used to conduct research on state-of-the-art radiation therapy techniques or treatment-related toxicity. However, it is difficult to conduct the focal irradiation to a shallow depth on small animals, because irradiation using LINAC has limitations in energy and field size. The purpose of this paper was to develop a focal irradiation tool for high-precision irradiation and to evaluate beam characteristics.

MATERIALS/METHODS

We designed the collimator of 1 mm diameter consisting of tungsten material for high-precision irradiation applied to the clinical brachytherapy system and the percent depth dose and horizontal profile were measured. We compared the depth dose and horizontal profile with 4 mm diameter SRS cone for 6 MV in LINAC. We measured the PDD and horizontal profile using EBT3 film for high-precision irradiation of 1 mm diameter using Ir-192 source. In case of 4 mm diameter, the beam was measured using edge detector. In addition, all measurements were compared with the results of planning tool simulation.

RESULTS

In case of focal irradiation tool, the maximum dose showed at the surface for both measurement and simulation, and 26% and 32% doses at 1 mm depth, respectively. In addition, FWHM at a 1 mm depth showed that high-precision irradiation was possible with measurement and simulation results of 1.86 and 1.28 mm. In case of LINAC, the maximum dose was showed at a depth of 1 cm and 0.8 cm in the measurement and simulation, respectively. Even if the smallest cone is used, the FWHM at a dmax depth was 4.0 mm in both simulation and measurement.

CONCLUSION

We overcame the limitation for energy and field size through the focal irradiation tool for high-precision irradiation. The focal irradiation tool enables high dose delivery to the shallow depth. In addition, small FWHM reduced dose delivery to the periphery at a specific depth and enabled accurate dose delivery. These results mean that the focal irradiation tool can be useful in small animal experiments that require accurate doses near the shallow depth.

Tagmentation-based analysis reveals the clonal behavior of CAR-T cells in association with lentivector integration sites

Integration site (IS) analysis is essential in ensuring safety and efficacy of gene therapies when integrating vectors are used. Although clinical trials of gene therapy are rapidly increasing, current methods have limited use in clinical settings because of their lengthy protocols. Here, we describe a novel genome-wide IS analysis method, "detection of the integration sites in a time-efficient manner, quantifying clonal size using tagmentation sequencing" (DIStinct-seq). In DIStinct-seq, a bead-linked Tn5 transposome is used, allowing the sequencing library to be prepared within a single day. We validated the quantification performance of DIStinct-seq for measuring clonal size with clones of known IS. Using ex vivo chimeric antigen receptor (CAR)-T cells, we revealed the characteristics of lentiviral IS. We then applied it to CAR-T cells collected at various times from tumor-engrafted mice, detecting 1,034-6,233 IS. Notably, we observed that the highly expanded clones had a higher integration frequency in the transcription units and vice versa in genomic safe harbors (GSH). Also, in GSH, persistent clones had more frequent IS. Together with these findings, the new IS analysis method will help to improve the safety and efficacy of gene therapies.

Association of Phosphorylated Pyruvate Dehydrogenase with Pyruvate Kinase M2 Promotes PKM2 Stability in Response to Insulin

Insulin is a crucial signalling molecule that primarily functions to reduce blood glucose levels through cellular uptake of glucose. In addition to its role in glucose homeostasis, insulin has been shown to regulate cell proliferation. Specifically, insulin enhances the phosphorylation of pyruvate dehydrogenase E1α (PDHA1) at the Ser293 residue and promotes the proliferation of HepG2 hepatocellular carcinoma cells. Furthermore, we previously observed that p-Ser293 PDHA1 bound with pyruvate kinase M2 (PKM2) as confirmed by coimmunoprecipitation. In this study, we used an in silico analysis to predict the structural conformation of the two binding proteins. However, the function of the protein complex remained unclear. To investigate further, we treated cells with si-PDHA1 and si-PKM2, which led to a reduction in PKM2 and p-Ser293 PDHA1 levels, respectively. Additionally, we found that the PDHA S293A dephospho-mimic reduced PKM2 levels and its associated enzyme activity. Treatment with MG132 and leupeptin impeded the PDHA1 S293A-mediated PKM2 reduction. These results suggest that the association between p-PDHA1 and PKM2 promotes their stability and protects them from protein degradation. Of interest, we observed that p-PDHA1 and PKM2 were localized in the nucleus in liver cancer patients. Under insulin stimulation, the knockdown of both PDHA1 and PKM2 led to a reduction in the expression of common genes, including KDMB1. These findings suggest that p-PDHA1 and PKM2 play a regulatory role in these proteins' expression and induce tumorigenesis in response to insulin.

A genome-wide by PM10 exposure interaction study for blood pressure in Korean adults

Blood pressure (BP) is a typical complex trait, and the genetic susceptibility of individuals to changes in BP induced by air pollution exposure is different. Although interactions of exposure to air pollutants with several candidate genes have been identified, genome-wide interaction studies (GWISs) are needed to understand the association between them with BP. Therefore, we aimed to discover the unique genetic loci for BP that interact with exposure to air pollutants in Korean adults. We ultimately included 1868 participants in the discovery step and classified them into groups of those with low-to-moderate exposure and high exposure to average annual concentration of particulate matter with an aerodynamic diameter ≤ 10 μm (PM10). Because none of the single nucleotide polymorphisms (SNPs) achieved a genome-wide level of significance of pint < 5 × 10-8 for either systolic BP (SBP) or diastolic BP (DBP), we considered the top 10 ranking SNPs for each BP trait. To validate these suggestive SNPs, we finally selected six genetic variants for SBP and five variants for DBP, respectively. In a replication result for SBP, only one SNP (rs12914147) located in an intergenic region of the NR2F2 showed a significant interaction. We also identified several genetic susceptibility loci (e.g., CHST11, TEK, and ITGA1) implicated in candidate mechanisms such as inflammation and oxidative stress in the discovery step, although their interaction effects were not replicated. Our study reports the first GWIS finding to our knowledge, and the association between exposure to PM10 and BP levels may be determined in part by several newly discovered genetic suggestive loci, including NR2F2.

A whole-genome reference panel of 14,393 individuals for East Asian populations accelerates discovery of rare functional variants

Underrepresentation of non-European (EUR) populations hinders growth of global precision medicine. Resources such as imputation reference panels that match the study population are necessary to find low-frequency variants with substantial effects. We created a reference panel consisting of 14,393 whole-genome sequences including more than 11,000 Asian individuals. Genome-wide association studies were conducted using the reference panel and a population-specific genotype array of 72,298 subjects for eight phenotypes. This panel yields improved imputation accuracy of rare and low-frequency variants within East Asian populations compared with the largest reference panel. Thirty-nine previously unidentified associations were found, and more than half of the variants were East Asian specific. We discovered genes with rare protein-altering variants, including LTBP1 for height and GPR75 for body mass index, as well as putative regulatory mechanisms for rare noncoding variants with cell type-specific effects. We suggest that this dataset will add to the potential value of Asian precision medicine.

A single-cell atlas of in vitro multiculture systems uncovers the in vivo lineage trajectory and cell state in the human lung

We present an in-depth single-cell atlas of in vitro multiculture systems on human primary airway epithelium derived from normal and diseased lungs of 27 individual donors. Our large-scale single-cell profiling identified new cell states and differentiation trajectories of rare airway epithelial cell types in human distal lungs. By integrating single-cell datasets of human lung tissues, we discovered immune-primed subsets enriched in lungs and organoids derived from patients with chronic respiratory disease. To demonstrate the full potential of our platform, we further illustrate transcriptomic responses to various respiratory virus infections in vitro airway models. Our work constitutes a single-cell roadmap for the cellular and molecular characteristics of human primary lung cells in vitro and their relevance to human tissues in vivo.

Earlier Age at Type 2 Diabetes Diagnosis Is Associated With Increased Genetic Risk of Cardiovascular Disease

OBJECTIVE

We investigated genetic risk of cardiovascular disease (CVD) by age at type 2 diabetes (T2D) diagnosis.

RESEARCH DESIGN AND METHODS

We compared incident CVD events by age at T2D diagnosis using UK Biobank (N = 12,321) and the Seoul National University Hospital (SNUH) cohort (N = 1,165). Genetic risk was quantified using polygenic risk score (PRS).

RESULTS

Individuals with earlier T2D diagnosis had higher CVD risk. In UK Biobank, the effect size of coronary artery disease (CAD) PRS on incident CAD was largest in individuals diagnosed with T2D at ages 30-39 years (hazard ratio 2.25; 95% CI 1.56-3.26) and decreased as age at diagnosis increased: ages 40-49 (1.51; 1.30-1.75), 50-59 (1.36; 1.24-1.50), and 60-69 years (1.30; 1.14-1.48) (Pinteraction = 0.0031). A similar trend was observed in the SNUH cohort. This increased genetic risk associated with earlier T2D diagnosis was largely mitigated by a healthy lifestyle.

CONCLUSIONS

Individuals with an earlier T2D diagnosis have a higher genetic risk of CAD, and this information could be used to tailor lifestyle interventions.

Triple-negative breast cancer-derived extracellular vesicles promote a hepatic pre-metastatic niche via a cascade of microenvironment remodeling

Triple negative breast cancer (TNBC) patients often develop metastases in visceral organs including the liver but the detailed molecular mechanisms of TNBC liver metastasis is not clearly understood. In this study, we tried to dissect the process of pre-metastatic niche formation in liver by using patient-derived xenograft (PDX) models of TNBC with different metastatic propensity. RNA sequencing of TNBC PDX models that successfully metastasized to liver showed up-regulation of Cx3cr1 gene in the liver microenvironment. In syngeneic breast cancer models, the Cx3cr1 up-regulation in liver preceded the development of cancer cell metastasis and was the results of recruitment of CX3CR1-expressing macrophages. The recruitment was induced by the CX3CL1 production from the liver endothelial cells and this CX3CL1-CX3CR1 signaling in the pre-metastatic niche resulted in up-regulation of MMP9 that promoted macrophage migration and cancer cell invasion. Additionally, our data suggest that the extra-cellular vesicles derived from the breast cancer cells induced the TNF-alpha expression in liver which leads to the CX3CL1 up-regulation. Lastly, the plasma CX3CL1 levels in 155 breast cancer patients were significantly associated with development of liver metastasis. Implications: Our data provides previously unknown cascades regarding the molecular education of pre-metastatic niche in liver for TNBC.

Abstract P2-21-06: The relationship between the gene expression difference of fibroblasts in the tumor microenvironment and metastasis of triple-negative breast cancer

Background Triple negative breast cancer (TNBC) is a breast cancer subtype with increased risk of distant metastasis, especially to lung and liver. In addition to the intrinsic characteristics of tumor cells, the molecular features of the tumor microenvironment cells can also regulate the metastasis process of TNBC. In this study, we tried to characterize the transcriptomic features of microenvironment fibroblasts associated with the distant metastasis of TNBC. Method We established PDX models by obtaining tumor tissues from 26 triple-negative breast cancer patients and transplanting them into the mammary fat pads of immunodeficient mice. When the PDX tumors grew sufficiently in the fat pads, we excised them and saved the mice to assess whether each PDX tumor was capable of developing spontaneous distant metastasis. We used bulk RNA sequencing and single cell RNA sequencing to determine the transcriptomic features associated with the distant metastasis. Results We classified the 26 TNBC PDX models based on their in vivo metastasis capacity and the clinical outcomes of the corresponding patients. Non-metastatic PDX models were defined as the lack of developing distant metastasis in either mouse experiment or clinical follow-up data (n=5). PDX models were defined as metastatic models when they developed distant metastasis in both mouse and human (n=4). We performed single cell RNA sequencing using the nine PDX tumors and observed a substantial gene expression differences of fibroblasts in the tumor microenvironment between the metastatic and non-metastatic PDX models. The fibroblasts of the metastatic PDX models showed up-regulation of genes involved in the response to hypoxia, inflammatory response, neutrophil chemotaxis, cellular response to IL-1, positive regulation of cell proliferation, and negative regulation of apoptotic process. Moreover, this difference of gene expression between metastatic and non-metastatic TNBC PDX models were not identified by the bulk RNA sequencing data. Conclusion By using the single cell RNA sequencing using the TNBC PDX models, we identified key features of fibroblasts in the tumor microenvironment that contributes to the metastasis capacity of TNBC. Further studies are on-going to elucidate the underlying mechanisms.

Citation Format: Woo Hang Heo, Woochan Lee, SONGBIN LI, Jong-Il Kim, Hyeong-Gon Moon. The relationship between the gene expression difference of fibroblasts in the tumor microenvironment and metastasis of triple-negative breast cancer [abstract]. In: Proceedings of the 2022 San Antonio Breast Cancer Symposium; 2022 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2023;83(5 Suppl):Abstract nr P2-21-06.

An integrative approach for exploring the nature of fibroepithelial neoplasms

BACKGROUND

Malignant phyllodes tumour (MPT) is a rare breast malignancy with epithelial and mesenchymal features. Currently, there are no appropriate research models or effective targeted therapeutic approaches for MPT.

METHODS

We collected fresh frozen tissues from nine patients with MPT and performed whole-exome and RNA sequencing. Additionally, we established patient-derived xenograft (PDX) models from patients with MPT and tested the efficacy of targeting dysregulated pathways in MPT using the PDX model from one MPT.

RESULTS

MPT has unique molecular characteristics when compared to breast cancers of epithelial origin and can be classified into two groups. The PDX model derived from one patient with MPT showed that the mouse epithelial component increased during tumour growth. Moreover, targeted inhibition of platelet-derived growth factor receptor (PDGFR) and phosphoinositide 3-kinase (PI3K)/mammalian target of rapamycin (mTOR) by imatinib mesylate and PKI-587 showed in vivo tumour suppression effects.

CONCLUSIONS

This study revealed the molecular profiles of MPT that can lead to molecular classification and potential targeted therapy, and suggested that the MPT PDX model can be a useful tool for studying the pathogenesis of fibroepithelial neoplasms and for preclinical drug screening to find new therapeutic strategies for MPT.

Abstract B019: Identify features of fibroblast and macrophage which regulate metastasis of triple-negative breast cancer

The tumor microenvironment (TME), which the cancer is known to remodel, can interact with cancer in various ways, such as promoting or hindering metastasis. The challenge with comprehending TME is that each patient's individual variance makes it difficult to carry out a consistent analysis. We used the PDX model to analyze the alterations entirely brought about by cancer cells of each patient in the same environment exclusively. In 121 breast cancer patients, we created a PDX model and performed WES and RNA-seq. We repeatedly performed metastatic induction tests on the models of 26 TNBC patients among the PDX models we generated. The replicated experiments confirmed that cancer cells with frequent metastasis usually spread well, while cancer cells with seldom metastases rarely metastasize. For each of the five samples that had the greatest and worst rates of metastasis, the pre-metastatic xenograft tumor was analyzed using scRNA-seq. Metastasis could not be predicted by only cancer cells because there was no difference between the two groups. Still, fibroblasts and macrophages had distinct gene markers before metastasis which can predict metastasis. To explore TME that develops pre-metastatic niches (PMN) in organs that would be metastasis targets, we developed a metastatic model using 4T1 cell lines. Daily lung macrophage RNA-seq, lung WES, and lung scRNA-seq procedures were carried out for four weeks to examine the precise timing of metastasis and the daily changes in the transcriptome. Surprisingly, a few days before the presence of any cancer cells, the gene expression of lung macrophages had already undergone significant modification. Additionally, only injections of tumor-conditioned media that mimicked metastasis could activate PMN in lung macrophages.

Citation Format: Woochan Lee, Hyewon Chung, Woohang Heo, Hyeong-Gon Moon, Seung Hyeok Seok, Jong-Il Kim. Identify features of fibroblast and macrophage which regulate metastasis of triple-negative breast cancer [abstract]. In: Proceedings of the AACR Special Conference: Cancer Metastasis; 2022 Nov 14-17; Portland, OR. Philadelphia (PA): AACR; Cancer Res 2022;83(2 Suppl_2):Abstract nr B019.

AKR1C2 Promotes Metastasis and Regulates the Molecular Features of Luminal Androgen Receptor Subtype in Triple Negative Breast Cancer Cells

PURPOSE

Patients with triple-negative breast cancer (TNBC) have an increased risk of distant metastasis compared to those with other subtypes. In this study, we aimed to identify the genes associated with distant metastasis in TNBC and their underlying mechanisms.

METHODS

We established patient-derived xenograft (PDX) models using surgically resected breast cancer tissues from 31 patients with TNBC. Among these, 15 patients subsequently developed distant metastases. Candidate metastasis-associated genes were identified using RNA sequencing. In vitro wound healing, proliferation, migration, and invasion assays and in vivo tumor xenograft and metastasis assays were performed to determine the functional importance of aldo-keto reductase family 1 member C2 (AKR1C2). Additionally, we used the METABRIC dataset to investigate the potential role of AKR1C2 in regulating TNBC subtypes and their downstream signaling activities.

RESULTS

RNA sequencing of primary and PDX tumors showed that genes involved in steroid hormone biosynthesis, including AKR1C2, were significantly upregulated in patients who subsequently developed metastasis. In vitro and in vivo assays showed that silencing of AKR1C2 resulted in reduced cell proliferation, migration, invasion, tumor growth, and incidence of lung metastasis. AKR1C2 was upregulated in the luminal androgen receptor (LAR) subtype of TNBC in the METABRIC dataset, and AKR1C2 silencing resulted in the downregulation of LAR classifier genes in TNBC cell lines. The androgen receptor (AR) gene was a downstream mediator of AKR1C2-associated phenotypes in TNBC cells. AKR1C2 expression was associated with gene expression pathways that regulate AR expression, including JAK-STAT signaling or interleukin 6 (IL-6). The levels of phospho-signal transducer and activator of transcription and IL-6, along with secreted IL-6, were significantly downregulated in AKR1C2-silenced TNBC cells.

CONCLUSION

Our data indicate that AKR1C2 is an important regulator of cancer growth and metastasis in TNBC and may be a critical determinant of LAR subtype features.

Whole genomic approach in mutation discovery of infantile spasms patients

Infantile spasms (IS) are a clinically and genetically heterogeneous group of epilepsy disorders in early infancy. The genetic backgrounds of IS have been gradually unraveled along with the increased application of next-generation sequencing (NGS). However, to date, only selected genomic regions have been sequenced using a targeted approach in most cases of IS, and the genetic etiologies of the majority of patients remain unknown. We conducted a proof-of-concept study using whole-genome sequencing (WGS) for the genetic diagnosis of IS. We included 16 patients with IS for this study, and WGS was applied as a first-tier test for genetic diagnosis. In total, we sequenced the whole genomes of 28 participants, including the genomes of six patients, which were sequenced with those of their parents. Among variants identified, we focused on those located in epilepsy or seizure-associated genes. We used two different methods to call relevant large deletions from WGS results. We found pathogenic or likely pathogenic variants in four patients (25.0%); a de novo variant in HDAC4, compound heterozygous variants in GRM7, and heterozygous variants in CACNA1E and KMT2E. We also selected two more candidate variants in SOX5 and SHROOM4 intronic regions. Although there are currently several difficulties in applying WGS for genetic diagnosis, especially in clinical interpretation of non-coding variants, we believe that developing sequencing technologies would overcome these hurdles in the near future. Considering the vast genetic heterogeneity and the substantial portion of patients with unknown etiologies, further studies using whole genomic approaches are necessary for patients with IS.

Whole-Genome Sequencing and Genomic Variant Analysis of Kazakh Individuals

Kazakhstan, the ninth-largest country in the world, is located along the Great Silk Road and connects Europe with Asia. Historically, its territory has been inhabited by nomadic tribes, and modern-day Kazakhstan is a multiethnic country with a dominant Kazakh population. We sequenced and analyzed the genomes of five ethnic Kazakhs at high coverage using the Illumina HiSeq2000 next-generation sequencing platform. The five Kazakhs yielded a total number of base pairs ranging from 87,308,581,400 to 107,526,741,301. On average, 99.06% were properly mapped. Based on the Het/Hom and Ti/Tv ratios, the quality of the genomic data ranged from 1.35 to 1.49 and from 2.07 to 2.08, respectively. Genetic variants were identified and annotated. Functional analysis of the genetic variants identified several variants that were associated with higher risks of metabolic and neurogenerative diseases. The present study showed high levels of genetic admixture of Kazakhs that were comparable to those of other Central Asians. These whole-genome sequence data of healthy Kazakhs could contribute significantly to biomedical studies of common diseases as their findings could allow better insight into the genotype–phenotype relations at the population level.

Northeastern Asian and Jomon-related genetic structure in the Three Kingdoms period of Gimhae, Korea

The genetic history of prehistoric and protohistoric Korean populations is not well understood because only a small number of ancient genomes are available. Here, we report the first paleogenomic data from the Korean Three Kingdoms period, a crucial point in the cultural and historic formation of Korea. These data comprise eight shotgun-sequenced genomes from ancient Korea (0.7×-6.1× coverage). They were derived from two archeological sites in Gimhae: the Yuha-ri shell mound and the Daesung-dong tumuli, the latter being the most important funerary complex of the Gaya confederacy. All individuals are from between the 4th and 5th century CE and are best modeled as an admixture between a northern China Bronze Age genetic source and a source of Jomon-related ancestry that shares similarities with the present-day genomes from Japan. The observed substructure and proportion of Jomon-related ancestry suggest the presence of two genetic groups within the population and diversity among the Gaya population. We could not correlate the genomic differences between these two groups with either social status or sex. All the ancient individuals' genomic profiles, including phenotypically relevant SNPs associated with hair and eye color, facial morphology, and myopia, imply strong genetic and phenotypic continuity with modern Koreans for the last 1,700 years.

Abstract 5726: Splicing-based classifier for gastric cancer identifies epithelial-mesenchymal transition subtypes associated with survival

Alternatively spliced RNA isoforms are a hallmark of tumors, but their nature, prevalence, and clinical implications in gastric cancer are unknown. We systematically profiled the splicing landscape of 83 gastric tumors and matched normal mucosa. We identified and experimentally validated eight splicing events that can classify all gastric cancers into three subtypes: Epithelial-splicing, Mesenchymal-splicing, and Hybrid-splicing. These subtypes were associated with distinct molecular signatures and epithelial-mesenchymal transition markers. Subtype-specific splicing events were enriched in motifs for splicing factors RBM24 and ESRP1, which were upregulated in Mesenchymal-splicing and Epithelial-splicing tumors, respectively. A simple classifier based only on RNA levels of RBM24 and ESRP1, and which is thus readily implementable in the clinic, is sufficient to distinguish gastric cancer subtypes and predict patient survival in multiple independent patient cohorts

Citation Format: Yukyung Jun, Yun-Suhk Suh, SungHee Park, Jieun Lee, Jong-Il Kim, Sanghyuk Lee, Wan-Ping Lee, Olga Anczuków, Han-Kwang Yang, Charles Lee. Splicing-based classifier for gastric cancer identifies epithelial-mesenchymal transition subtypes associated with survival [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 5726.

First Report of Pseudomonas cichorii Causing Bacterial vein necrosis on Perilla plants [Perilla frutescens (L.) Britton.] in South Korea

Perilla (Perilla frutescens L.) is the second most important upland crop and the third largest edible oil crop in Korea (Shin and Kim 1994). During a disease survey in Busan, Korea in September 2021, symptoms of vein necrosis were observed in perilla plants, with incidences of approximately 30% and 50% in two fields. Symptoms of spots on the perilla appeared as leaf dryness and spots with water-soaked blotches largely concentrated on the mid-veins of leaves. The lesions were initiated with water-soaked spots on the leaf or stem and gradually turned black or brown. Necrosis was also observed in the stems. A bacterium was isolated on Luria-Bertani (LB) agar from diseased leaf tissues that were surface-disinfected with 70% ethyl alcohol for 3-5 min and then washed with sterile water three times. Three pieces of sterilized leaf tissue (size: 0.5 × 0.5 cm) were mixed with 500 µL sterile water for 30 min, and then the suspension was serially diluted and spread on LB agar. Subsequently, isolates were cultivated on LB agar and King's Medium B agar (KMB) (Schaad et al. 2001), and they were predominantly cream-colored and circular bacterial colonies with undulated margins. The bacterial colonies on KMB displayed fluorescence under 365 nm UV light. The isolates were analyzed with the GEN III MicroPlate (Biolog, Hayward, CA, USA), and all isolates were identified as Pseudomonas cichorii, a devastating plant bacterium that damages a wide range of host plants worldwide, including in South Korea (Hikichi et al. 2013; Ramkumar et al. 2015). To identify the species of the bacterial pathogen, genomic DNA of four isolates (BS4922, BS4167, BS4345, and BS4560) was extracted, and the 16S rRNA gene and hrcRST gene were amplified with universal primers, 27F/1492R and Hcr1/Hcr2, and sequencing was then done (Patel et al. 2019). In the BLAST analysis, the 16S rRNA sequences (GenBank OM060656, OM275434, OM275435, OM275436) showed a 100% and 99% similarity to P. cichorii strains MAFF 302698 (AB724286) and P. cichorii strain Pc-Gd-4 (KU923373), respectively. Further, hrcRST gene sequences (GenBank OM143596, OM268864, OM268865, and OM268866) showed high similarity (>99%) with P. cichorii strain P16-51 (MG518230). A pathogenicity test of the four isolates was performed on 3 - 4 weeks old perilla plants by creating wounds with a needle on the lower leaves and stems, and then the plants were inoculated by spraying inoculum (108 CFU/ml). The plants that served as the negative control were wounded and sprayed with unsterilized water. The inoculated perilla plants were placed in a greenhouse at 28 ± 2oC , 80-85% relative humidity, and a natural photoperiod. The inoculation site began to show symptoms of water-soaked brown lesions. Disease symptoms such as leaf dryness, water-soaked blotches on the mid-vein of leaves, and necrosis on plant stems were observed in the inoculated plants 7-10 days after inoculation, whereas the plants of the negative control group did not show any symptoms. The bacteria were re-isolated from the diseased tissues of the plants, and DNA sequence analysis identified them as P. cichorii. Additionally, all isolates induced hypersensitivity reactions in tobacco and tomato leaves within 24 h after inoculation. To our knowledge, this is the first report of P. cichorii infecting perilla in South Korea. The findings in this study will provide the basic information for the development of diagnostic tools and management measures against P. cichorii in perilla.

Cross-Talk between Wnt Signaling and Src Tyrosine Kinase

Src, a non-receptor tyrosine kinase, was first discovered as a prototype oncogene and has been shown to critical for cancer progression for a variety of tissues. Src activity is regulated by a number of post-translational modifications in response to various stimuli. Phosphorylations of Src Tyr419 (human; 416 in chicken) and Src Tyr530 (human; 527 in chicken) have been known to be critical for activation and inactivation of Src, respectively. Wnt signaling regulates a variety of cellular functions including for development and cell proliferation, and has a role in certain diseases such as cancer. Wnt signaling is carried out through two pathways: β-catenin-dependent canonical and β-catenin-independent non-canonical pathways as Wnt ligands bind to their receptors, Frizzled, LRP5/6, and ROR1/2. In addition, many signaling components including Axin, APC, Damm, Dishevelled, JNK kinase and Rho GTPases contribute to these canonical and non-canonical Wnt pathways. However, the communication between Wnt signaling and Src tyrosine kinase has not been well reviewed as Src regulates Wnt signaling through LRP6 tyrosine phosphorylation. GSK-3β phosphorylated by Wnt also regulates Src activity. As Wnt signaling and Src mutually regulate each other, it is noted that aberrant regulation of these components give rise to various diseases including typically cancer, and as such, merit a closer look.

Expression-based species deconvolution and realignment removes misalignment error in multispecies single-cell data

Background

Although single-cell RNA sequencing of xenograft samples has been widely used, no comprehensive bioinformatics pipeline is available for human and mouse mixed single-cell analyses. Considering the numerous homologous genes across the human and mouse genomes, misalignment errors should be evaluated, and a new algorithm is required. We assessed the extents and effects of misalignment errors and exonic multi-mapping events when using human and mouse combined reference data and developed a new bioinformatics pipeline with expression-based species deconvolution to minimize errors. We also evaluated false-positive signals presumed to originate from ambient RNA of the other species and address the importance to computationally remove them.

Result

Error when using combined reference account for an average of 0.78% of total reads, but such reads were concentrated to few genes that were greatly affected. Human and mouse mixed single-cell data, analyzed using our pipeline, clustered well with unmixed data and showed higher k-nearest-neighbor batch effect test and Local Inverse Simpson’s Index scores than those derived from Cell Ranger (10 × Genomics). We also applied our pipeline to multispecies multisample single-cell library containing breast cancer xenograft tissue and successfully identified all samples using genomic array and expression. Moreover, diverse cell types in the tumor microenvironment were well captured.

Conclusion

We present our bioinformatics pipeline for mixed human and mouse single-cell data, which can also be applied to pooled libraries to obtain cost-effective single-cell data. We also address misalignment, multi-mapping error, and ambient RNA as a major consideration points when analyzing multispecies single-cell data.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04676-0.

Deconvolution of Adult T-Cell Leukemia/Lymphoma With Single-Cell RNA-Seq Using Frozen Archived Skin Tissue Reveals New Subset of Cancer-Associated Fibroblast

Adult T-cell Leukemia/Lymphoma (ATLL) is a rare aggressive T-cell malignancy caused by human T-cell leukemia virus type 1 (HTLV-1) infection. However, little is known about the underlying activated molecular pathways at the single cell level. Moreover, the intercellular communications between the tumor microenvironment (TME) and tumor cells in this malignancy are currently unknown. Difficulties in harvesting fresh tissue in a clinical setting have hampered our deeper understanding of this malignancy. Herein, we examined ATLL using archived fresh frozen tissue after biopsy using single-cell RNA sequencing (scRNA-seq) with T-cell receptor (TCR) clonal analysis. Highly clonal tumor cells showed multiple activating pathways, suggesting dynamic evolution of the malignancy. By dissecting diverse cell types comprising the TME, we identified a novel subset of cancer-associated fibroblast, which showed enriched epidermal growth factor receptor (EGFR)-related transcripts including early growth response 1 and 2 (EGR1 and EGR2). Cancer associated fibroblasts (CAFs) of ATLL play an important role for CD4 T-cell proliferation via FGF7-FGF1 and PDGFA-PDGFRA/B signaling, and CAFs, particularly EGR-enriched, are also associated with CD8 and NKT expansion by EGFR. These findings suggest a potential targeted therapeutic pathway to better treat this neoplasm.

Experimental development of the epigenomic library construction method to elucidate the epigenetic diversity and causal relationship between epigenome and transcriptome at a single-cell level

The method of single-cell RNA sequencing has been rapidly developed, and numerous experiments have been conducted over the past decade. Their results allow us to recognize various subpopulations and rare cell states in tissues, tumors, and immune systems that are previously unidentified, and guide us to understand fundamental biological processes that determine cell identity based on single-cell gene expression profiles. However, it is still challenging to understand the principle of comprehensive gene regulation that determines the cell fate only with transcriptome, a consequential output of the gene expression program. To elucidate the mechanisms related to the origin and maintenance of comprehensive single-cell transcriptome, we require a corresponding single-cell epigenome, which is a differentiated information of each cell with an identical genome. This review deals with the current development of single-cell epigenomic library construction methods, including multi-omics tools with crucial factors and additional requirements in the future focusing on DNA methylation, chromatin accessibility, and histone post-translational modifications. The study of cellular differentiation and the disease occurrence at a single-cell level has taken the first step with single-cell transcriptome and is now taking the next step with single-cell epigenome.

Comprehensive Analysis of Alternative Splicing in Gastric Cancer Identifies Epithelial-Mesenchymal Transition Subtypes Associated with Survival

Alternatively spliced RNA isoforms are a hallmark of tumors, but their nature, prevalence, and clinical implications in gastric cancer have not been comprehensively characterized. We systematically profiled the splicing landscape of 83 gastric tumors and matched normal mucosa, identifying and experimentally validating eight splicing events that can classify all gastric cancers into three subtypes: epithelial-splicing (EpiS), mesenchymal-splicing (MesS), and hybrid-splicing. These subtypes were associated with distinct molecular signatures and epithelial-mesenchymal transition markers. Subtype-specific splicing events were enriched in motifs for splicing factors RBM24 and ESRP1, which were upregulated in MesS and EpiS tumors, respectively. A simple classifier based only on RNA levels of RBM24 and ESRP1, which can be readily implemented in the clinic, was sufficient to distinguish gastric cancer subtypes and predict patient survival in multiple independent patient cohorts. Overall, this study provides insights into alternative splicing in gastric cancer and the potential clinical utility of splicing-based patient classification.

SIGNIFICANCE

This study presents a comprehensive analysis of alternative splicing in the context of patient classification, molecular mechanisms, and prognosis in gastric cancer.

Image Correlation-Based Method to Assess Ciliary Beat Frequency in Human Airway Organoids

Ciliary movements within the human airway are essential for maintaining a clean lung environment. Motile cilia have a characteristic ciliary beat frequency (CBF). However, CBF measurement with current video microscopic techniques can be error-prone due to the use of the single-point Fourier transformation, which is often biased for ciliary measurements. Herein, we describe a new video microscopy technique that harnesses a metric of motion-contrast imaging and image correlation for CBF analysis. It can provide objective and selective CBF measurements for individual motile cilia and generate CBF maps for the imaged area. The measurement performance of our methodology was validated with in vitro human airway organoid models that simulated an actual human airway epithelium. The CBF determined for the region of interest (ROI) was equal to that obtained with manual counting. The signal redundancy problem of conventional methods was not observed. Moreover, the obtained CBF measurements were robust to optical focal shifts, and exhibited spatial heterogeneity and temperature dependence. This technique can be used to evaluate ciliary movement in respiratory tracts and determine whether it is non-synchronous or aperiodic in patients. Therefore, our observations suggest that the proposed method can be clinically adapted as a screening tool to diagnose ciliopathies.

Glucose metabolic profiles evaluated by PET associated with molecular characteristic landscape of gastric cancer

BACKGROUND

Although FDG-PET is widely used in cancer, its role in gastric cancer (GC) is still controversial due to variable [¹⁸F]fluorodeoxyglucose ([¹⁸F]FDG) uptake. Here, we sought to develop a genetic signature to predict high FDG-avid GC to plan individualized PET and investigate the molecular landscape of GC and its association with glucose metabolic profiles noninvasively evaluated by [¹⁸F]FDG-PET.

METHODS

Based on a genetic signature, PETscore, representing [¹⁸F]FDG avidity, was developed by imaging data acquired from thirty patient-derived xenografts (PDX). The PETscore was validated by [¹⁸F]FDG-PET data and gene expression data of human GC. The PETscore was associated with genomic and transcriptomic profiles of GC using The Cancer Genome Atlas.

RESULTS

Five genes, PLS1, PYY, HBQ1, SLC6A5, and NAT16, were identified for the predictive model for [¹⁸F]FDG uptake of GC. The PETscore was validated in independent PET data of human GC with qRT-PCR and RNA-sequencing. By applying PETscore on TCGA, a significant association between glucose uptake and tumor mutational burden as well as genomic alterations were identified.

CONCLUSION

Our findings suggest that molecular characteristics are underlying the diverse metabolic profiles of GC. Diverse glucose metabolic profiles may apply to precise diagnostic and therapeutic approaches for GC.

A Genome-Wide Association Study on Abdominal Adiposity-Related Traits in Adult Korean Men

INTRODUCTION

Although previous genome-wide association studies (GWASs) have identified genetic susceptibility loci for abdominal adiposity, GWASs on Asian samples remain scarce. Therefore, we performed a GWAS for abdominal adipose tissue depots in a Korean population.

METHODS

A total of 1,937 Korean men were included in the study. Areas of abdominal fat were quantified by computed tomography. We performed a GWAS analysis under an additive model, and a replication study was conducted on 480 additional Korean adult men.

RESULTS

In the discovery step, we identified a total of 10 single-nucleotide polymorphisms (SNPs) associated with adiposity indicators (p < 1 × 10-5). The top SNP, rs1028014, for visceral adipose tissue (VAT) was located in the ZMAT4 gene and remained significant after adjustment for body mass index (BMI). Three additional SNPs were also associated with VAT-adj-BMI and located within the SLC26A10, FAM155A, and COL4A1-COL4A2 genes, respectively. In addition, we identified a SNP (rs4668224) of the MYO3B gene for visceral-to-subcutaneous fat ratio. For subcutaneous adipose tissue and total adipose tissue, two (rs6585735 and rs363527) and three SNPs (rs1487892, rs9357565, and rs1985358) were found, respectively. Overall, eight SNPs were used in the replication study; however, none of the SNPs reached our level of significance for replication (p < 0.0063). Nevertheless, rs4773144 of COL4A1-COL4A2 for VAT-adj-BMI was the most interesting SNP identified in previous GWASs for coronary artery disease (based on the same risk allele "G"), along with functional effects.

CONCLUSION

This study suggests for the first time that an SNP (rs4773144) of COL4A1-COL4A2 may contribute to the increase in VAT level, especially in adult Korean men.

Calsequestrin 2 overexpression in breast cancer increases tumorigenesis and metastasis by modulating the tumor microenvironment

The spatial tumor shape is determined by the complex interactions between tumor cells and their microenvironment. Here, we investigated the role of a newly identified breast cancer-related gene, calsequestrin 2 (CASQ2), in tumor-microenvironment interactions during tumor growth and metastasis. We analyzed gene expression and three-dimensional tumor shape data from the breast cancer dataset of The Cancer Genome Atlas (TCGA) and identified CASQ2 as a potential regulator of tumor-microenvironment interaction. In TCGA breast cancer cases containing information of three-dimensional tumor shapes, CASQ2 mRNA showed the highest correlation with the spatial tumor shapes. Furthermore, we investigated the expression pattern of CASQ2 in human breast cancer tissues. CASQ2 was not detected in breast cancer cell lines in vitro but was induced in the xenograft tumors and human breast cancer tissues. To evaluate the role of CASQ2, we established CASQ2-overexpressing breast cancer cell lines for in vitro and in vivo experiments. CASQ2 overexpression in breast cancer cells resulted in a more aggressive phenotype and altered epithelial-mesenchymal transition (EMT) markers in vitro. CASQ2 overexpression induced cancer-associated fibroblast characteristics along with increased hypoxia-inducible factor 1α (HIF1α) expression in stromal fibroblasts. CASQ2 overexpression accelerated tumorigenesis, induced collagen structure remodeling, and increased distant metastasis in vivo. CASQ2 conferred more metaplastic features to triple-negative breast cancer cells. Our data suggest that CASQ2 is a key regulator of breast cancer tumorigenesis and metastasis by modulating diverse aspects of tumor-microenvironment interactions.

Feasibility of sodium long chain polyphosphate as a potential growth promoter in broilers

The objective of this study was to evaluate in vitro antimicrobial and anti-biofilm activity of sodium long chain polyphosphate (SLCPP) and effect of dietary supplementation of SLCPP on growth performance, organ characteristics, blood metabolites, and intestinal microflora of broilers. Antimicrobial activities of SLCPP were observed against Escherichia coli O157:H7, Listeria monocytogenes, Salmonella enterica ser. Pullorum, Shigella sonnei, Klebsiella pneumonia, Pseudomonas aeruginosa in agar well diffusion assay. In addition, SLCPP demonstrated good anti-biofilm activity against K. pneumonia and P. aeruginosa. Furthermore, to investigate the dietary effect of SLCPP, a total of 480 1-day-old male Ross 308 broiler chicks were randomly allotted to three dietary treatment groups (4 replicates per group, 40 birds in each replicate): an antibiotic-free corn-soybean meal basal diet (NC); basal diet + enramycin 0.01% (PC); and basal diet + 0.1% SLCPP (SPP). The experiment lasted for 35 days. Results showed that birds fed with SLCPP had higher body weight (BW) and average daily gain (ADG), and lower feed conversion ratio (FCR) during the grower phase (days 7 to 21) (p < 0.05). Except for blood urea nitrogen, all other blood biochemical parameters remained unaffected by the dietary supplementation of SLCPP. Compared to the control group, lengths of the duodenum and ileum in the SPP group were significantly shorter (p < 0.05). Moreover, counts of lactic acid bacteria (LAB), total aerobes, and Streptococcus spp. in jejunum as well as LAB in cecum were increased in the SPP group than in the PC group (p < 0.05). These results suggest that dietary supplementation of SLCPP might promote the growth of broilers in their early growth phase.

Thyroid nodules in childhood-onset Hashimoto's thyroiditis: Frequency, risk factors, follow-up course and genetic alterations of thyroid cancer

OBJECTIVE

We evaluated the frequency, risk factors and the follow-up outcomes of thyroid nodules, and genetic alterations in thyroid cancer, in youth with childhood-onset Hashimoto thyroiditis (HT) residing in an iodine-sufficient country.

DESIGN

A retrospective cohort study.

PATIENTS AND MEASUREMENTS

A total of 213 patients (194 females, mean age 10.6 years at the time of HT diagnosis) were ultrasonographically evaluated. Thyroid nodules were categorized using the Korean Thyroid Imaging Reporting and Data System (K-TIRADS) and the American College of Radiology Thyroid Imaging Reporting and Data System (ACR-TI-RADS).

RESULTS

Thyroid nodules were detected in 40 (18.8%) patients over a median follow-up period of 3.4 years, usually after the onset of puberty. A family history of thyroid disease (hazard ratio 2.1, p = .031) was predictive of thyroid nodule detection. Papillary thyroid carcinoma (PTC) was diagnosed in 9 (4.2% of all and 22.5% of nodule-positive patients). The malignant nodules had a higher K-TIRADS or ACR-TI-RADS risk level compared with benign nodules (p < .01 for both). Genetic alterations were revealed in 7 (BRAF^V600E in 6 and RET-ERC1 fusion in 1) of the eight available tumour tissue samples. None showed evidence of disease over a median follow-up period of 3.4 years.

CONCLUSIONS

The nodule detection rate was 18.8%, with a 22.5% risk of malignancy among the detected nodules in childhood-onset HT patients, showing increased risk in those with a family history. Additional large-scale studies are required to evaluate the usefulness of K-TIRADS or ACR-TI-RADS risk level for the differentiation of paediatric thyroid nodules.

Radiation Response Prediction Model based on Integrated Clinical and Genomic Data Analysis

Purpose

The value of the genomic profiling by targeted gene-sequencing on radiation therapy response prediction was evaluated through integrated analysis including clinical information. Radiation response prediction model was constructed based on the analyzed findings.

Materials and Methods

Patients who had the tumor sequenced using institutional cancer panel after informed consent and received radiotherapy for the measurable disease served as the target cohort. Patients with irradiated tumor locally controlled for more than 6 months after radiotherapy were defined as the durable local control (DLC) group, otherwise, non-durable local control (NDLC) group. Significant genomic factors and domain knowledge were used to develop the Bayesian Network model to predict radiotherapy response.

Results

Altogether, 88 patients were collected for analysis. Of those, 41 (43.6%) and 47 (54.4%) patients were classified as the NDLC and DLC group, respectively. Somatic mutations of NOTCH2 and BCL were enriched in the NDLC group, whereas, mutations of CHEK2, MSH2, and NOTCH1 were more frequently found in the DLC group. Altered DNA repair pathway was associated with better local failure-free survival (HR 0.40, 95%CI 0.19-0.86, p=0.014). Smoking somatic signature was found more frequently in the DLC group. AUC of the Bayesian Network model predicting probability of 6-month local control was 0.83.

Conclusion

Durable radiation response was associated with alterations of DNA repair pathway and smoking somatic signature. Bayesian network model could provide helpful insights for high precision radiotherapy. However, these findings should be verified in prospective cohort for further individualization.

Sequencing Cell-free Fetal DNA in Pregnant Women With GCK-MODY

CONTEXT

Individuals with monogenic diabetes due to inactivating glucokinase (GCK) variants typically do not require treatment, except potentially during pregnancy. In pregnancy, fetal GCK genotype determines whether treatment is indicated, but noninvasive methods are not clinically available.

OBJECTIVE

This work aims to develop a method to determine fetal GCK genotype noninvasively using maternal cell-free fetal DNA.

METHODS

This was a proof-of-concept study involving 3 pregnant women with a causal GCK variant that used information from 1) massive parallel sequencing of maternal plasma cell-free DNA, 2) direct haplotype sequences of maternal genomic DNA, and 3) the paternal genotypes to estimate relative haplotype dosage of the pathogenic variant-linked haplotype. Statistical testing of variant inheritance was performed using a sequential probability ratio test (SPRT).

RESULTS

In each of the 3 cases, plasma cell-free DNA was extracted once between gestational weeks 24 and 36. The fetal fraction of cell-free DNA ranged from 21.8% to 23.0%. Paternal homozygous alleles that were identical to the maternal GCK variant-linked allele were not overrepresented in the cell-free DNA. Paternal homozygous alleles that were identical to the maternal wild-type-linked allele were significantly overrepresented. Based on the SPRT, we predicted that all 3 cases did not inherit the GCK variant. Postnatal infant genotyping confirmed our prediction in each case.

CONCLUSION

We have successfully implemented a noninvasive method to predict fetal GCK genotype using cell-free DNA in 3 pregnant women carrying an inactivating GCK variant. This method could guide tailoring of hyperglycemia treatment in pregnancies of women with GCK monogenic diabetes.

Predictive biomarkers for 5-fluorouracil and oxaliplatin-based chemotherapy in gastric cancers via profiling of patient-derived xenografts

Gastric cancer (GC) is commonly treated by chemotherapy using 5-fluorouracil (5-FU) derivatives and platinum combination, but predictive biomarker remains lacking. We develop patient-derived xenografts (PDXs) from 31 GC patients and treat with a combination of 5-FU and oxaliplatin, to determine biomarkers associated with responsiveness. When the PDXs are defined as either responders or non-responders according to tumor volume change after treatment, the responsiveness of PDXs is significantly consistent with the respective clinical outcomes of the patients. An integrative genomic and transcriptomic analysis of PDXs reveals that pathways associated with cell-to-cell and cell-to-extracellular matrix interactions enriched among the non-responders in both cancer cells and the tumor microenvironment (TME). We develop a 30-gene prediction model to determine the responsiveness to 5-FU and oxaliplatin-based chemotherapy and confirm the significant poor survival outcomes among cases classified as non-responder-like in three independent GC cohorts. Our study may inform clinical decision-making when designing treatment strategies.

The Origin and Composition of Korean Ethnicity Analyzed by Ancient and Present-Day Genome Sequences

Koreans are thought to be an ethnic group of admixed northern and southern subgroups. However, the exact genetic origins of these two remain unclear. In addition, the past admixture is presumed to have taken place on the Korean peninsula, but there is no genomic scale analysis exploring the origin, composition, admixture, or the past migration of Koreans. Here, 88 Korean genomes compared with 91 other present-day populations showed two major genetic components of East Siberia and Southeast Asia. Additional paleogenomic analysis with 115 ancient genomes from Pleistocene hunter-gatherers to Iron Age farmers showed a gradual admixture of Tianyuan (40 ka) and Devil’s gate (8 ka) ancestries throughout East Asia and East Siberia up until the Neolithic era. Afterward, the current genetic foundation of Koreans may have been established through a rapid admixture with ancient Southern Chinese populations associated with Iron Age Cambodians. We speculate that this admixing trend initially occurred mostly outside the Korean peninsula followed by continuous spread and localization in Korea, corresponding to the general admixture trend of East Asia. Over 70% of extant Korean genetic diversity is explained to be derived from such a recent population expansion and admixture from the South.

NTRK- and RET-fusion-directed therapy in pediatric thyroid cancer yields a tumor response and radioiodine uptake

BACKGROUND

Molecular characterization in pediatric papillary thyroid cancer (PTC), distinct from adult PTC, is important for developing molecular targeted therapies for progressive 131I-refractory PTC.

METHODS

PTC samples from 106 pediatric patients (age: 4.3-19.8 years; 21 boys) who attended Seoul National University Hospital (January 1983-March 2020) were available for genomic profiling. Previous transcriptome data from 125 adult PTCs were used for comparison.

RESULTS

Genetic drivers were found in 80 tumors; 31 with fusion oncogenes (RET in 21, ALK in 6, and NTRK1/3 in 4), 47 with point mutations (BRAFV600E in 41, TERTC228T in 2, and DICER1 variants in 5), and 2 with amplifications. Fusion-oncogene PTCs, predominantly detected in younger patients, presented with a more advanced stage and showed more recurrent or persistent disease than BRAFV600E PTCs, which were detected mostly in adolescents. Pediatric fusion PTCs (in those aged < 10 years) showed lower expression of thyroid differentiation genes, including SLC5A5, than adult fusion PTCs. Two girls with progressive 131I-refractory lung metastases harboring a TPR-NTRK1 or CCDC6-RET fusion received fusion-targeted therapy; larotrectinib and selpercatinib decreased the tumor extent and restored radioiodine uptake. The girl with the CCDC6-RET fusion received 131I therapy combined with selpercatinib, leading to a tumor response. In vitro 125I uptake and 131I clonogenic assays showed that larotrectinib inhibited growth and restored radioiodine avidity.

CONCLUSIONS

In pediatric fusion-oncogene PTC cases with 131I-refractory advanced disease, selective fusion-directed therapy may restore radioiodine avidity and lead to a dramatic tumor response, underscoring the importance of molecular testing in pediatric PTC patients.

FUNDING

The Ministry of Science, ICT & Future Planning (grant number NRF-2016R1A2B4012417 91 and 2019R1A2C2084332), the Ministry of Health & Welfare, Republic of Korea (grant number 92 H14C1277), the Ministry of Education (grant number 2020R1A6A1A03047972), and the Seoul 93 National University Hospital Research Fund (grant number 04-2015-0830).