B cell super-enhancers and regulatory clusters recruit AID tumorigenic activity.

The antibody gene mutator activation-induced cytidine deaminase (AID) promiscuously damages oncogenes, leading to chromosomal translocations and tumorigenesis. Why nonimmunoglobulin loci are susceptible to AID activity is unknown. Here, we study AID-mediated lesions in the context of nuclear architecture and the B cell regulome. We show that AID targets are not randomly distributed across the genome but are predominantly grouped within super-enhancers and regulatory clusters. Unexpectedly, in these domains, AID deaminates active promoters and eRNA(+) enhancers interconnected in some instances over megabases of linear chromatin. Using genome editing, we demonstrate that 3D-linked targets cooperate to recruit AID-mediated breaks. Furthermore, a comparison of hypermutation in mouse B cells, AID-induced kataegis in human lymphomas, and translocations in MEFs reveals that AID damages different genes in different cell types. Yet, in all cases, the targets are predominantly associated with topological complex, highly transcribed super-enhancers, demonstrating that these compartments are key mediators of AID recruitment.

Genomic basis for RNA alterations in cancer.

Transcript alterations often result from somatic changes in cancer genomes1. Various forms of RNA alterations have been described in cancer, including overexpression2, altered splicing3 and gene fusions4; however, it is difficult to attribute these to underlying genomic changes owing to heterogeneity among patients and tumour types, and the relatively small cohorts of patients for whom samples have been analysed by both transcriptome and whole-genome sequencing. Here we present, to our knowledge, the most comprehensive catalogue of cancer-associated gene alterations to date, obtained by characterizing tumour transcriptomes from 1,188 donors of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA)5. Using matched whole-genome sequencing data, we associated several categories of RNA alterations with germline and somatic DNA alterations, and identified probable genetic mechanisms. Somatic copy-number alterations were the major drivers of variations in total gene and allele-specific expression. We identified 649 associations of somatic single-nucleotide variants with gene expression in cis, of which 68.4% involved associations with flanking non-coding regions of the gene. We found 1,900 splicing alterations associated with somatic mutations, including the formation of exons within introns in proximity to Alu elements. In addition, 82% of gene fusions were associated with structural variants, including 75 of a new class, termed 'bridged' fusions, in which a third genomic location bridges two genes. We observed transcriptomic alteration signatures that differ between cancer types and have associations with variations in DNA mutational signatures. This compendium of RNA alterations in the genomic context provides a rich resource for identifying genes and mechanisms that are functionally implicated in cancer.

Comprehensive evaluation of BRCA1/2 variant interpretation ability among laboratories in China.

High-quality interpretation of BRCA1/2 variants plays a critical role in the clinical practice of precision medicine. However, a comprehensive system to evaluate the quality and accuracy of variant interpretation has yet to be established. This study investigates the performance of an interpretation system in evaluating the capacities of BRCA1/2 interpretation among distinct laboratories in China. The evaluation system is based on a reference database that contains 750 different variants in BRCA1/2 Evaluation was performed among 41 laboratories in China. We classified their performance into five levels. Only level A was considered qualified. This level allows for a 0.3% error rate for clinical decision-related misinterpretation; 26 of 41 laboratories (63%) met the qualified standard, while 7 laboratories were at levels D and E, which indicated egregious mistakes and systemic problems in variant interpretation. Due to strict quality demands, the interpretation of several variants was amended, which largely influenced the quality rate. The number of qualified laboratories would decrease from 26 to 17 if those incorrect recommended interpretations were not corrected. This evaluation system provides a potential approach for standardisation of variant interpretation and lowers the discordance of variant interpretation between different laboratories. A well-designed interpretation ability evaluation is essential to evaluate the interpretation level of laboratories before they provide service in real-world clinical settings.

Genome sequencing of 320 Chinese children with epilepsy: a clinical and molecular study.

The aim of this study is to evaluate the diagnostic value of genome sequencing in children with epilepsy, and to provide genome sequencing-based insights into the molecular genetic mechanisms of epilepsy to help establish accurate diagnoses, design appropriate treatments and assist in genetic counselling. We performed genome sequencing on 320 Chinese children with epilepsy, and interpreted single-nucleotide variants and copy number variants of all samples. The complete pedigree and clinical data of the probands were established and followed up. The clinical phenotypes, treatments, prognoses and genotypes of the patients were analysed. Age at seizure onset ranged from 1 day to 17 years, with a median of 4.3 years. Pathogenic/likely pathogenic variants were found in 117 of the 320 children (36.6%), of whom 93 (29.1%) had single-nucleotide variants, 22 (6.9%) had copy number variants and two had both single-nucleotide variants and copy number variants. Single-nucleotide variants were most frequently found in SCN1A (10/95, 10.5%), which is associated with Dravet syndrome, followed by PRRT2 (8/95, 8.4%), which is associated with benign familial infantile epilepsy, and TSC2 (7/95, 7.4%), which is associated with tuberous sclerosis. Among the copy number variants, there were three with a length <25 kilobases. The most common recurrent copy number variants were 17p13.3 deletions (5/24, 20.8%), 16p11.2 deletions (4/24, 16.7%), and 7q11.23 duplications (2/24, 8.3%), which are associated with epilepsy, developmental retardation and congenital abnormalities. Four particular 16p11.2 deletions and two 15q11.2 deletions were considered to be susceptibility factors contributing to neurodevelopmental disorders associated with epilepsy. The diagnostic yield was 75.0% in patients with seizure onset during the first postnatal month, and gradually decreased in patients with seizure onset at a later age. Forty-two patients (13.1%) were found to be specifically treatable for the underlying genetic cause identified by genome sequencing. Three of them received corresponding targeted therapies and demonstrated favourable prognoses. Genome sequencing provides complete genetic diagnosis, thus enabling individualized treatment and genetic counselling for the parents of the patients. Genome sequencing is expected to become the first choice of methods for genetic testing of patients with epilepsy.

Genetic characteristics of human papillomavirus type 16, 18, 52 and 58 in southern China.

Persistent infections of high-risk human papillomaviruses (HPVs) are the leading cause of cervical cancers. We collected cervical exfoliated cell samples from females in Changsha city, Hunan Province and obtained 338 viral genomes of four major HPV types, including HPV 16 (n = 82), 18 (n = 35), 52 (n = 121) and 58 (n = 100). The lineage/sublineage distribution of the four HPVs confirmed previous epidemiological reports, with the predominant prevailing sublineage as A4 (50%), A1 (37%) and A3 (13%) for HPV16, A1 (83%) for HPV18, B2 (86%) for HPV52 and A1 (65%), A3 (19%) and A2 (12%) for HPV58. We also identified two potentially novel HPV18 sublineages, i.e. A6 and A7. Virus mutation analysis further revealed the presence of HPV16 and HPV58 sublineages associated with potentially high oncogenicity. These findings expanded our knowledge of the HPV genetic diversity in China, providing valuable evidence to facilitate HPV DNA screening, vaccine effectiveness evaluation and control strategy development.

Heterogeneous immunogenomic features and distinct escape mechanisms in multifocal hepatocellular carcinoma.

BACKGROUND & AIMS: The presence of multifocal tumors, developed either from intrahepatic metastasis (IM) or multicentric occurrence (MO), is a distinct feature of hepatocellular carcinoma (HCC). Immunogenomic characterization of multifocal HCC is important for understanding immune escape in different lesions and developing immunotherapy. METHODS: We combined whole-exome/transcriptome sequencing, multiplex immunostaining, immunopeptidomes, T cell receptor (TCR) sequencing and bioinformatic analyses of 47 tumors from 15 patients with HCC and multifocal lesions. RESULTS: IM and MO demonstrated distinct clonal architecture, mutational spectrum and genetic susceptibility. The immune microenvironment also displayed spatiotemporal heterogeneity, such as less T cell and more M2 macrophage infiltration in IM and higher expression of inhibitory immune checkpoints in MO. Similar to mutational profiles, shared neoantigens and TCR repertoires among tumors from the same patients were abundant in IM but scarce in MO. Combining neoantigen prediction and immunopeptidomes identified T cell-specific neoepitopes and achieved a high verification rate in vitro. Immunoediting mainly occurred in MO but not IM, due to the relatively low immune infiltration. Loss of heterozygosity of human leukocyte antigen (HLA) alleles, identified in 17% of multifocal HCC, hampered the ability of major histocompatibility complex to present neoantigens, especially in IM. An integrated analysis of Immunoscore, immunoediting, TCR clonality and HLA loss of heterozygosity in each tumor could stratify patients into 2 groups based on whether they have a high or low risk of recurrence (p = 0.038). CONCLUSION: Our study comprehensively characterized the genetic structure, neoepitope landscape, T cell profile and immunoediting status that collectively shape tumor evolution and could be used to optimize personalized immunotherapies for multifocal HCC. LAY SUMMARY: Immunogenomic features of multifocal hepatocellular carcinoma (HCC) are important for understanding immune-escape mechanisms and developing more effective immunotherapy. Herein, comprehensive immunogenomic characterization showed that diverse genomic structures within multifocal HCC would leave footprints on the immune landscape. Only a few tumors were under the control of immunosurveillance, while others evaded the immune system through multiple mechanisms that led to poor prognosis. Our study revealed heterogeneous immunogenomic landscapes and immune-constrained tumor evolution, the understanding of which could be used to optimize personalized immunotherapies for multifocal HCC.

Genetic landscape of hepatitis B virus-associated diffuse large B-cell lymphoma.

Hepatitis B virus (HBV) infection is endemic in some parts of Asia, Africa, and South America and remains to be a significant public health problem in these areas. It is known as a leading risk factor for the development of hepatocellular carcinoma, but epidemiological studies have also shown that the infection may increase the incidence of several types of B-cell lymphoma. Here, by characterizing altogether 275 Chinese diffuse large B-cell lymphoma (DLBCL) patients, we showed that patients with concomitant HBV infection (surface antigen positive [HBsAg+]) are characterized by a younger age, a more advanced disease stage at diagnosis, and reduced overall survival. Furthermore, by whole-genome/exome sequencing of 96 tumors and the respective peripheral blood samples and targeted sequencing of 179 tumors from these patients, we observed an enhanced rate of mutagenesis and a distinct set of mutation targets in HBsAg+ DLBCL genomes, which could be partially explained by the activities of APOBEC and activation-induced cytidine deaminase. By transcriptome analysis, we further showed that the HBV-associated gene expression signature is contributed by the enrichment of genes regulated by BCL6, FOXO1, and ZFP36L1. Finally, by analysis of immunoglobulin heavy chain gene sequences, we showed that an antigen-independent mechanism, rather than a chronic antigenic simulation model, is favored in HBV-related lymphomagenesis. Taken together, we present the first comprehensive genomic and transcriptomic study that suggests a link between HBV infection and B-cell malignancy. The genetic alterations identified in this study may also provide opportunities for development of novel therapeutic strategies.

Comprehensive genomic variation profiling of cervical intraepithelial neoplasia and cervical cancer identifies potential targets for cervical cancer early warning.

BACKGROUND: To better understand the pathogenesis of cervical cancer (CC), we systematically analysed the genomic variation and human papillomavirus (HPV) integration profiles of cervical intraepithelial neoplasia (CIN) and CC. METHODS: We performed whole-genome sequencing or whole-exome sequencing of 102 tumour-normal pairs and human papillomavirus probe capture sequencing of 45 CCs, 44 CIN samples and 25 normal cervical samples, and constructed strict integrated workflow of genomic analysis. RESULTS: Mutational analysis identified eight significantly mutated genes in CC including four genes (FAT1, MLL3, MLL2 and FADD), which have not previously been reported in CC. Targetable alterations were identified in 55.9% of patients. In addition, HPV integration breakpoints occurred in 97.8% of the CC samples, 70.5% of the CIN samples and 42.8% of the normal cervical samples with HPV infection. Integrations of high-risk HPV strains in CCs, including HPV16, 18, 33 and 58, also occurred in the CIN samples. Moreover, gene mutations were detected in 52% of the CIN specimens, and 54.8% of these mutations occurred in genes that also mutated in CCs. CONCLUSION: Our results lay the foundation for a deep understanding of the molecular mechanisms and finding new diagnostic and therapeutic targets of CC.

Nationwide population genetic screening improves outcomes of newborn screening for hearing loss in China.

PURPOSE: The benefits of concurrent newborn hearing and genetic screening have not been statistically proven due to limited sample sizes and outcome data. To fill this gap, we analyzed outcomes of newborns with genetic screening results. METHODS: Newborns in China were screened for 20 hearing-loss-related genetic variants from 2012 to 2017. Genetic results were categorized as positive, at-risk, inconclusive, or negative. Hearing screening results, risk factors, and up-to-date hearing status were followed up via phone interviews. RESULTS: Following up 12,778 of 1.2 million genetically screened newborns revealed a higher rate of hearing loss by three months of age among referrals from the initial hearing screening (60% vs. 5.0%, P < 0.001) and a lower rate of lost-to-follow-up/documentation (5% vs. 22%, P < 0.001) in the positive group than in the inconclusive group. Importantly, genetic screening detected 13% more hearing-impaired infants than hearing screening alone and identified 2,638 (0.23% of total) newborns predisposed to preventable ototoxicity undetectable by hearing screening. CONCLUSION: Incorporating genetic screening improves the effectiveness of newborn hearing screening programs by elucidating etiologies, discerning high-risk subgroups for vigilant management, identifying additional children who may benefit from early intervention, and informing at-risk newborns and their maternal relatives of increased susceptibility to ototoxicity.

Genetic basis of PD-L1 overexpression in diffuse large B-cell lymphomas.

Diffuse large B-cell lymphoma (DLBCL) is one of the most common and aggressive types of B-cell lymphoma. Deregulation of proto-oncogene expression after a translocation, most notably to the immunoglobulin heavy-chain locus (IGH), is one of the hallmarks of DLBCL. Using whole-genome sequencing analysis, we have identified the PD-L1/PD-L2 locus as a recurrent translocation partner for IGH in DLBCL. PIM1 and TP63 were also identified as novel translocation partners for PD-L1/PD-L2 Fluorescence in situ hybridization was furthermore used to rapidly screen an expanded DLBCL cohort. Collectively, a subset of samples was found to be affected by gains (12%), amplifications (3%), and translocations (4%) of the PD-L1/PD-L2 locus. RNA sequencing data coupled with immunohistochemistry revealed that these cytogenetic alterations correlated with increased expression of PD-L1 but not of PD-L2 Moreover, cytogenetic alterations affecting the PD-L1/PD-L2 locus were more frequently observed in the non-germinal center B cell-like (non-GCB) subtype of DLBCL. These findings demonstrate the genetic basis of PD-L1 overexpression in DLBCL and suggest that treatments targeting the PD-1-PD-L1/PD-L2 axis might benefit DLBCL patients, especially those belonging to the more aggressive non-GCB subtype.

Clonal Characteristics of Circulating B Lymphocyte Repertoire in Primary Biliary Cholangitis.

Primary biliary cholangitis (PBC) is an autoimmune liver disease characterized by elevated serum anti-mitochondrial Ab and lymphocyte-mediated bile duct damage. This study was designed to reveal the clonal characteristics of B lymphocyte repertoire in patients with PBC to facilitate better understanding of its pathogenesis and better management of these patients. Using high-throughput sequencing of Ig genes, we analyzed the repertoire of circulating B lymphocytes in 43 patients with PBC, and 34 age- and gender-matched healthy controls. Compared with healthy controls, PBC patients showed 1) a gain of 14 new clones and a loss of 8 clones; 2) a significant clonal expansion and increased relative IgM abundance, which corresponded with the elevated serum IgM level; 3) a significant reduction of clonal diversity and somatic hypermutations in class-switched sequences, which suggested a general immunocompromised status; 4) the reduction of clonal diversity and enhancement of clonal expansion were more obvious at the cirrhotic stage; and 5) treatment with ursodeoxycholic acid could increase the clonal diversity and reduce clonal expansion of the IgM repertoire, with no obvious effect on the somatic hypermutation level. Our data suggest that PBC is a complex autoimmune disease process with evidence of B lymphocyte clonal gains and losses, Ag-dependent ogligoclonal expansion, and a generally compromised immune reserve. This new insight into the pathogenesis of PBC opens up the prospect of studying disease-relevant B cells to better diagnose and treat this devastating disease.

The stepwise evolution of the exome during acquisition of docetaxel resistance in breast cancer cells.

BACKGROUND: Resistance to taxane-based therapy in breast cancer patients is a major clinical problem that may be addressed through insight of the genomic alterations leading to taxane resistance in breast cancer cells. In the current study we used whole exome sequencing to discover somatic genomic alterations, evolving across evolutionary stages during the acquisition of docetaxel resistance in breast cancer cell lines. RESULTS: Two human breast cancer in vitro models (MCF-7 and MDA-MB-231) of the step-wise acquisition of docetaxel resistance were developed by exposing cells to 18 gradually increasing concentrations of docetaxel. Whole exome sequencing performed at five successive stages during this process was used to identify single point mutational events, insertions/deletions and copy number alterations associated with the acquisition of docetaxel resistance. Acquired coding variation undergoing positive selection and harboring characteristics likely to be functional were further prioritized using network-based approaches. A number of genomic changes were found to be undergoing evolutionary selection, some of which were likely to be functional. Of the five stages of progression toward resistance, most resistance relevant genomic variation appeared to arise midway towards fully resistant cells corresponding to passage 31 (5 nM docetaxel) for MDA-MB-231 and passage 16 (1.2 nM docetaxel) for MCF-7, and where the cells also exhibited a period of reduced growth rate or arrest, respectively. MCF-7 cell acquired several copy number gains on chromosome 7, including ABC transporter genes, including ABCB1 and ABCB4, as well as DMTF1, CLDN12, CROT, and SRI. For MDA-MB-231 numerous copy number losses on chromosome X involving more than 30 genes was observed. Of these genes, CASK, POLA1, PRDX4, MED14 and PIGA were highly prioritized by the applied network-based gene ranking approach. At higher docetaxel concentration MCF-7 subclones exhibited a copy number loss in E2F4, and the gene encoding this important transcription factor was down-regulated in MCF-7 resistant cells. CONCLUSIONS: Our study of the evolution of acquired docetaxel resistance identified several genomic changes that might explain development of docetaxel resistance. Interestingly, the most relevant resistance-associated changes appeared to originate midway through the evolution towards fully resistant cell lines. Our data suggest that no single genomic event sufficiently predicts resistance to docetaxel, but require genomic alterations affecting multiple pathways that in concert establish the final resistance stage.

The potential role of Alu Y in the development of resistance to SN38 (Irinotecan) or oxaliplatin in colorectal cancer.

BACKGROUND: Irinotecan (SN38) and oxaliplatin are chemotherapeutic agents used in the treatment of colorectal cancer. However, the frequent development of resistance to these drugs represents a considerable challenge in the clinic. Alus as retrotransposons comprise 11% of the human genome. Genomic toxicity induced by carcinogens or drugs can reactivate Alus by altering DNA methylation. Whether or not reactivation of Alus occurs in SN38 and oxaliplatin resistance remains unknown. RESULTS: We applied reduced representation bisulfite sequencing (RRBS) to investigate the DNA methylome in SN38 or oxaliplatin resistant colorectal cancer cell line models. Moreover, we extended the RRBS analysis to tumor tissue from 14 patients with colorectal cancer who either did or did not benefit from capecitabine + oxaliplatin treatment. For the clinical samples, we applied a concept of 'DNA methylation entropy' to estimate the diversity of DNA methylation states of the identified resistance phenotype-associated methylation loci observed in the cell line models. We identified different loci being characteristic for the different resistant cell lines. Interestingly, 53% of the identified loci were Alu sequences- especially the Alu Y subfamily. Furthermore, we identified an enrichment of Alu Y sequences that likely results from increased integration of new copies of Alu Y sequence in the drug-resistant cell lines. In the clinical samples, SOX1 and other SOX gene family members were shown to display variable DNA methylation states in their gene regions. The Alu Y sequences showed remarkable variation in DNA methylation states across the clinical samples. CONCLUSION: Our findings imply a crucial role of Alu Y in colorectal cancer drug resistance. Our study underscores the complexity of colorectal cancer aggravated by mobility of Alu elements and stresses the importance of personalized strategies, using a systematic and dynamic view, for effective cancer therapy.

Genomic and epigenomic analyses of monozygotic twins discordant for congenital renal agenesis.

Monozygotic twins have been widely studied to distinguish genetic and environmental factors in the pathogenesis of human diseases. For renal agenesis, the one-sided absence of renal tissue, the relative contributions of genetic and environmental factors to its pathogenesis are still unclear. In this study of a pair of monozygotic twins discordant for congenital renal agenesis, the genomic profile was analyzed from a set of blood samples using high-throughput exome-capture sequencing to detect single-nucleotide polymorphisms (SNPs), copy number variations (CNVs), and insertions and deletions (indels). Also, an epigenomic analysis used reduced-representation bisulfite sequencing to detect differentially methylated regions (DMRs). No discordant SNPs, CNVs, or indels were confirmed, but 514 DMRs were detected. KEGG analysis indicated the DMRs localized to 10 signaling pathways and 25 genes, including the mitogen-activated protein kinase pathway and 6 genes (FGF18, FGF12, PDGFRA, MAPK11, AMH, CTBP1) involved in organ development. Although methylation results from our adult patient and her sister may not represent the pattern that was present during kidney development, we could at least confirm a lack of obvious differences at the genome level, which suggests that nongenetic factors may be involved in the pathogenesis of renal agenesis.

IL-17A deficiency inhibits lung cancer-induced osteoclastogenesis by promoting apoptosis of osteoclast precursor cells.

Osteoclasts are crucial in the events leading to bone metastasis of lung cancer. Interleukin-17A (IL-17A) affects osteogenesis by regulating the survival of osteoclast precursors (OCPs) and is enriched in lung cancer cells. However, how factors derived from tumor cells that metastasize to bone affect osteoclastogenesis remains poorly understood. We examined whether IL-17A derived from lung cancer cells affects osteoclast differentiation by regulating OCP apoptosis. IL-17A expression was inhibited in A549 non-small cell lung cancer cells using RNA interference. Compared with conditioned medium (CM) from A549 cells (A549-CM), CM from IL-17A-deficient A549 cells (A549-si-CM) suppressed osteoclastogenesis. The mRNA expression of osteoclast-specific genes was downregulated following A549-si-CM treatment. Furthermore, A549-si-CM promoted osteoclast precursor apoptosis at an early stage of osteoclastogenesis, which was related to the promotion of caspase-3 expression by A549-si-CM during osteoclast differentiation. In vivo experiments also showed that inhibition of IL-17A expression in A549 cells reduced osteoclast activation and bone tissue destruction. Collectively, our results indicate that IL-17A deficiency inhibits lung cancer-induced osteoclast differentiation by promoting apoptosis of osteoclast precursors in the early stage of osteoclast formation and that IL-17A is a potential therapeutic target for cancer-associated bone resorption in patients with lung cancer.

Germline Mutational Landscape and Novel Targetable RAD51D Variant in Chinese Patients With Ovarian Cancer.

PURPOSE: Genetic variants of ovarian cancer (OV) show ethnic differences, but data from the Chinese population are still insufficient. Here, we elucidate the inheritance landscape in Chinese patients with OV and examine the functional implications of a Chinese-enriched RAD51D variant. METHODS: Between 2015 and 2018, 373 consecutive patients with OV were prospectively enrolled. Variants of BRCA1/2, other homologous recombination repair (HRR) genes, and DNA mismatch repair (MMR) genes were analyzed using next-generation sequencing. An enriched RAD51D variant was identified, and its functional effects were examined using Cell Counting Kit-8, colony formation, transwell migration, and drug sensitivity assays. RESULTS: Overall, 31.1% (116/373) of patients had at least one pathogenic or likely pathogenic germline variant. BRCA1 and BRCA2 accounted for 16.09% and 5.36%, respectively, with one patient having both variants. In addition, 32 (8.58%) patients carried other HRR gene variants, whereas three (0.8%) patients had MMR gene variants. The RAD51D variant ranked third (8/373, 2.1%), and its rate was much higher than that in other populations. Remarkably, all eight patients harbored the RAD51D K91fs variant (c.270_271dup, p.Lys91Ilefs*13) and demonstrated satisfactory platinum response and favorable prognosis. This variant confers enhanced sensitivity to poly (ADP-ribose) polymerase inhibitors in OV cells. However, the effects on platinum sensitivity were inconsistent across different cell lines. Against the background of the TP53 variant, RAD51D K91fs variant showed increased sensitivity to cisplatin. CONCLUSION: Our study revealed the inheritance landscape of OV and identified an enriched RAD51D variant in Chinese patients with OV. This can serve as an important reference for OV management and a potential therapeutic target.

Circulating cell-free DNA methylation-based multi-omics analysis allows early diagnosis of pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is a highly aggressive cancer with a 5-year survival rate of 7.2% in China. However, effective approaches for diagnosis of PDAC are limited. Tumor-originating genomic and epigenomic aberration in circulating free DNA (cfDNA) have potential as liquid biopsy biomarkers for cancer diagnosis. Our study aims to assess the feasibility of cfDNA-based liquid biopsy assay for PDAC diagnosis. In this study, we performed parallel genomic and epigenomic profiling of plasma cfDNA from Chinese PDAC patients and healthy individuals. Diagnostic models were built to distinguish PDAC patients from healthy individuals. Cancer-specific changes in cfDNA methylation landscape were identified, and a diagnostic model based on six methylation markers achieved high sensitivity (88.7% for overall cases and 78.0% for stage I patients) and specificity (96.8%), outperforming the mutation-based model significantly. Moreover, the combination of the methylation-based model with carbohydrate antigen 19-9 (CA19-9) levels further improved the performance (sensitivity: 95.7% for overall cases and 95.5% for stage I patients; specificity: 93.3%). In conclusion, our findings suggest that both methylation-based and integrated liquid biopsy assays hold promise as non-invasive tools for detection of PDAC.

Consistent signatures in the human gut microbiome of old- and young-onset colorectal cancer.

The incidence of young-onset colorectal cancer (yCRC) has been increasing in recent decades, but little is known about the gut microbiome of these patients. Most studies have focused on old-onset CRC (oCRC), and it remains unclear whether CRC signatures derived from old patients are valid in young patients. To address this, we assembled the largest yCRC gut metagenomes to date from two independent cohorts and found that the CRC microbiome had limited association with age across adulthood. Differential analysis revealed that well-known CRC-associated taxa, such as Clostridium symbiosum, Peptostreptococcus stomatis, Parvimonas micra and Hungatella hathewayi were significantly enriched (false discovery rate <0.05) in both old- and young-onset patients. Similar strain-level patterns of Fusobacterium nucleatum, Bacteroides fragilis and Escherichia coli were observed for oCRC and yCRC. Almost all oCRC-associated metagenomic pathways had directionally concordant changes in young patients. Importantly, CRC-associated virulence factors (fadA, bft) were enriched in both oCRC and yCRC compared to their respective controls. Moreover, the microbiome-based classification model had similar predication accuracy for CRC status in old- and young-onset patients, underscoring the consistency of microbial signatures across different age groups.