Blockade of STING activation alleviates microglial dysfunction and a broad spectrum of Alzheimer's disease pathologies.

Abnormal glial activation promotes neurodegeneration in Alzheimer's disease (AD), the most common cause of dementia. Stimulation of the cGAS-STING pathway induces microglial dysfunction and sterile inflammation, which exacerbates AD. We showed that inhibiting STING activation can control microglia and ameliorate a wide spectrum of AD symptoms. The cGAS-STING pathway is required for the detection of ectopic DNA and the subsequent immune response. Amyloid-ß (Aß) and tau induce mitochondrial stress, which causes DNA to be released into the cytoplasm of microglia. cGAS and STING are highly expressed in Aß plaque-associated microglia, and neuronal STING is upregulated in the brains of AD model animals. The presence of the APOE ε4 allele, an AD risk factor, also upregulated both proteins. STING activation was necessary for microglial NLRP3 activation, proinflammatory responses, and type-I-interferon responses. Pharmacological STING inhibition reduced a wide range of AD pathogenic features in AppNL-G-F/hTau double-knock-in mice. An unanticipated transcriptome shift in microglia reduced gliosis and cerebral inflammation. Significant reductions in the Aß load, tau phosphorylation, and microglial synapse engulfment prevented memory loss. To summarize, our study describes the pathogenic mechanism of STING activation as well as its potential as a therapeutic target in AD.

Genome-wide association study for metabolic syndrome reveals APOA5 single nucleotide polymorphisms with multilayered effects in Koreans.

BACKGROUND AND PURPOSE: Genome-wide association studies (GWAS) of metabolic syndrome (MetS) have predominantly focused on non-Asian populations, with limited representation from East Asian cohorts. Moreover, previous GWAS analyses have primarily emphasized the significance of top single nucleotide polymorphisms (SNPs), poorly explaining other SNP signals in linkage disequilibrium. This study aimed to reveal the interaction between rs651821 and rs2266788, the principal variants of apolipoprotein A5 (APOA5), within the most significant loci identified through GWAS on MetS. METHODS: GWAS on MetS and its components was conducted using the data from the Korean Genome and Epidemiology Study (KoGES) city cohort comprising 58,600 individuals with available biochemical, demographic, lifestyle factors, and the most significant APOA5 locus was analyzed further in depth. RESULTS: According to GWAS of MetS and its diagnostic components, a significant association between the APOA5 SNPs rs651821/rs2266788 and MetS/triglycerides/high-density lipoprotein phenotypes was revealed. However, a conditional analysis employing rs651821 unveiled a reversal in the odds ratio for rs2266788. Therefore, rs651821 and rs2266788 emerged as independent and opposing signals in the extended GWAS analysis, i.e., the multilayered effects. Further gene-environment interaction analyses regarding lifestyle factors such as smoking, alcohol consumption, and physical activity underscored these multilayered effects. CONCLUSION: This study unveils the intricate interplay between rs651821 and rs2266788 derived from MetS GWAS. Removing the influence of lead SNP reveals an independent protective signal associated with rs2266788, suggesting a multilayered effect between these SNPs. These findings underline the need for novel perspectives in future MetS GWAS.

Harnessing Institutionally Developed Clinical Targeted Sequencing to Improve Patient Survival in Breast Cancer: A Seven-Year Experience.

Purpose: Considering the high disease burden and unique features of Asian patients with breast cancer (BC), it is essential to have a comprehensive view of genetic characteristics in this population. An institutional targeted sequencing platform was developed through the Korea Research-Driven Hospitals project and was incorporated into clinical practice. This study explores the use of targeted next-generation sequencing (NGS) and its outcomes in patients with advanced/metastatic BC in the real world. Materials and Methods: We reviewed the results of NGS tests administered to BC patients using a customized sequencing platform - FiRST Cancer Panel (FCP) - over seven years. We systematically described clinical translation of FCP for precise diagnostics, personalized therapeutic strategies, and unraveling disease pathogenesis. Results: NGS tests were conducted on 548 samples from 522 patients with BC. 97.6% of tested samples harbored at least one pathogenic alteration. The common alterations included mutations in TP53(56.2%), PIK3CA(31.2%), GATA3(13.8%), BRCA2(10.2%), and amplifications of CCND1(10.8%), FGF19(10.0%), and ERBB2(9.5%). NGS analysis of ERBB2 amplification correlated well with HER2 immunohistochemistry and in situ hybridization. RNA panel analyses found potentially actionable and prognostic fusion genes. FCP effectively screened for potentially germline pathogenic/likely pathogenic mutation. 10.3% of BC patients received matched therapy guided by NGS, resulting in a significant overall survival advantage (p=0.022), especially for metastatic BCs. . Conclusion: Clinical NGS provided multifaceted benefits, deepening our understanding of the disease, improving diagnostic precision, and paving the way for targeted therapies. The concrete advantages of FCP highlight the importance of multi-gene testing for BC, especially for metastatic conditions.

Microglia Gravitate toward Amyloid Plaques Surrounded by Externalized Phosphatidylserine via TREM2.

Microglia play a crucial role in synaptic elimination by engulfing dystrophic neurons via triggering receptors expressed on myeloid cells 2 (TREM2). They are also involved in the clearance of beta-amyloid (Aß) plaques in Alzheimer's disease (AD); nonetheless, the driving force behind TREM2-mediated phagocytosis of beta-amyloid (Aß) plaques remains unknown. Here, using advanced 2D/3D/4D co-culture systems with loss-of-function mutations in TREM2 (a frameshift mutation engineered in exon 2) brain organoids/microglia/assembloids, it is identified that the clearance of Aß via TREM2 is accelerated by externalized phosphatidylserine (ePtdSer) generated from dystrophic neurons surrounding the Aß plaques. Moreover, it is investigated whether microglia from both sporadic (CRISPR-Cas9-based APOE4 lines) and familial (APPNL-G-F/MAPT double knock-in mice) AD models show reduced levels of TREM2 and lack of phagocytic activity toward ePtdSer-positive Aß plaques. Herein new insight is provided into TREM2-dependent microglial phagocytosis of Aß plaques in the context of the presence of ePtdSer during AD progression.

Image-based deep learning model using DNA methylation data predicts the origin of cancer of unknown primary.

Cancer of unknown primary (CUP) is a rare type of metastatic cancer in which the origin of the tumor is unknown. Since the treatment strategy for patients with metastatic tumors depends on knowing the primary site, accurate identification of the origin site is important. Here, we developed an image-based deep-learning model that utilizes a vision transformer algorithm for predicting the origin of CUP. Using DNA methylation dataset of 8,233 primary tumors from The Cancer Genome Atlas (TCGA), we categorized 29 cancer types into 18 organ classes and extracted 2,312 differentially methylated CpG sites (DMCs) from non-squamous cancer group and 420 DMCs from squamous cell cancer group. Using these DMCs, we created organ-specific DNA methylation images and used them for model training and testing. Model performance was evaluated using 394 metastatic cancer samples from TCGA (TCGA-meta) and 995 samples (693 primary and 302 metastatic cancers) obtained from 20 independent external studies. We identified that the DNA methylation image reveals a distinct pattern based on the origin of cancer. Our model achieved an overall accuracy of 96.95 % in the TCGA-meta dataset. In the external validation datasets, our classifier achieved overall accuracies of 96.39 % and 94.37 % in primary and metastatic tumors, respectively. Especially, the overall accuracies for both primary and metastatic samples of non-squamous cell cancer were exceptionally high, with 96.79 % and 96.85 %, respectively.

Higher Genetic Risk for Type 2 Diabetes Is Associated With a Faster Decline of ß-Cell Function in an East Asian Population.

OBJECTIVE: While most genetic variants of type 2 diabetes (T2D) are suggested to be associated with ß-cell dysfunction cross sectionally, their association with the longitudinal change of ß-cell function remains largely unknown. RESEARCH DESIGN AND METHODS: We analyzed data from 6,311 participants without T2D at baseline (mean [SD] age 51.6 [8.7] years) from a community-based prospective cohort in Korea. Participants underwent biennial 2-h 75-g oral glucose tolerance tests (OGTTs) during 14 years of follow-up, and the OGTT-derived disposition index (DI) was used as a marker for ß-cell function. Genetic risk was quantified using the genome-wide polygenic risk score (PRS) and was stratified into low (1st quintile), intermediate (2nd-4th quintiles), and high (5th quintile) genetic risk. Lifestyle was assessed according to Life's Essential 8. RESULTS: During a mean follow-up of 10.9 years, 374 (29.6%), 851 (22.5%), and 188 (14.9%) participants developed T2D in the high, intermediate, and low genetic risk groups, respectively. Compared with the low genetic risk group, participants in the high genetic risk group had a 25% lower DI at baseline. Furthermore, in longitudinal analysis, we observed a 1.83-fold faster decline in log2-transformed DI per year (-0.034 vs. -0.019, P = 2.1 × 10-3; per 1-SD increase in T2D PRS, P = 1.2 × 10-4). Healthy lifestyle attenuated the rate of decline in DI across all genetic risk groups. CONCLUSIONS: Individuals with a higher genetic risk for T2D exhibited not only a lower OGTT-derived ß-cell function at baseline but also a notably more rapid decline during follow-up. This information could be used to enable a focused precision prevention with lifestyle intervention.

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.

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.

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 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.

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.

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.

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.

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.

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.

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.

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.