Analysis of blood methylation quantitative trait loci in East Asians reveals ancestry-specific impacts on complex traits.

Methylation quantitative trait loci (mQTLs) are essential for understanding the role of DNA methylation changes in genetic predisposition, yet they have not been fully characterized in East Asians (EAs). Here we identified mQTLs in whole blood from 3,523 Chinese individuals and replicated them in additional 1,858 Chinese individuals from two cohorts. Over 9% of mQTLs displayed specificity to EAs, facilitating the fine-mapping of EA-specific genetic associations, as shown for variants associated with height. Trans-mQTL hotspots revealed biological pathways contributing to EA-specific genetic associations, including an ERG-mediated 233 trans-mCpG network, implicated in hematopoietic cell differentiation, which likely reflects binding efficiency modulation of the ERG protein complex. More than 90% of mQTLs were shared between different blood cell lineages, with a smaller fraction of lineage-specific mQTLs displaying preferential hypomethylation in the respective lineages. Our study provides new insights into the mQTL landscape across genetic ancestries and their downstream effects on cellular processes and diseases/traits.

Reliable biological and multi-omics research through biometrology.

Metrology is the science of measurement and its applications, whereas biometrology is the science of biological measurement and its applications. Biometrology aims to achieve accuracy and consistency of biological measurements by focusing on the development of metrological traceability, biological reference measurement procedures, and reference materials. Irreproducibility of biological and multi-omics research results from different laboratories, platforms, and analysis methods is hampering the translation of research into clinical uses and can often be attributed to the lack of biologists' attention to the general principles of metrology. In this paper, the progresses of biometrology including metrology on nucleic acid, protein, and cell measurements and its impacts on the improvement of reliability and comparability in biological research are reviewed. Challenges in obtaining more reliable biological and multi-omics measurements due to the lack of primary reference measurement procedures and new standards for biological reference materials faced by biometrology are discussed. In the future, in addition to establishing reliable reference measurement procedures, developing reference materials from single or multiple parameters to multi-omics scale should be emphasized. Thinking in way of biometrology is warranted for facilitating the translation of high-throughput omics research into clinical practices.

High Intensity Focused Ultrasound-Driven Nanomotor for Effective Ferroptosis-Immunotherapy of TNBC.

The heterogeneity of triple-negative breast cancers (TNBC) remains challenging for various treatments. Ferroptosis, a recently identified form of cell death resulting from the unrestrained peroxidation of phospholipids, represents a potential vulnerability in TNBC. In this study, a high intensity focused ultrasound (HIFU)-driven nanomotor is developed for effective therapy of TNBC through induction of ferroptosis. Through bioinformatics analysis of typical ferroptosis-associated genes in the FUSCCTNBC dataset, gambogic acid is identified as a promising ferroptosis drug and loaded it into the nanomotor. It is found that the rapid motion of nanomotors propelled by HIFU significantly enhanced tumor accumulation and penetration. More importantly, HIFU not only actuated nanomotors to trigger effective ferroptosis of TNBC cells, but also drove nanomotors to activate ferroptosis-mediated antitumor immunity in primary and metastatic TNBC models, resulting in effective tumor regression and prevention of metastases. Overall, HIFU-driven nanomotors show great potential for ferroptosis-immunotherapy of TNBC.

Integrated multiomic profiling of breast cancer in the Chinese population reveals patient stratification and therapeutic vulnerabilities.

Molecular profiling guides precision treatment of breast cancer; however, Asian patients are underrepresented in publicly available large-scale studies. We established a comprehensive multiomics cohort of 773 Chinese patients with breast cancer and systematically analyzed their genomic, transcriptomic, proteomic, metabolomic, radiomic and digital pathology characteristics. Here we show that compared to breast cancers in white individuals, Asian individuals had more targetable AKT1 mutations. Integrated analysis revealed a higher proportion of HER2-enriched subtype and correspondingly more frequent ERBB2 amplification and higher HER2 protein abundance in the Chinese HR+HER2+ cohort, stressing anti-HER2 therapy for these individuals. Furthermore, comprehensive metabolomic and proteomic analyses revealed ferroptosis as a potential therapeutic target for basal-like tumors. The integration of clinical, transcriptomic, metabolomic, radiomic and pathological features allowed for efficient stratification of patients into groups with varying recurrence risks. Our study provides a public resource and new insights into the biology and ancestry specificity of breast cancer in the Asian population, offering potential for further precision treatment approaches.

Quartet metabolite reference materials for inter-laboratory proficiency test and data integration of metabolomics profiling.

BACKGROUND: Various laboratory-developed metabolomic methods lead to big challenges in inter-laboratory comparability and effective integration of diverse datasets. RESULTS: As part of the Quartet Project, we establish a publicly available suite of four metabolite reference materials derived from B lymphoblastoid cell lines from a family of parents and monozygotic twin daughters. We generate comprehensive LC-MS-based metabolomic data from the Quartet reference materials using targeted and untargeted strategies in different laboratories. The Quartet multi-sample-based signal-to-noise ratio enables objective assessment of the reliability of intra-batch and cross-batch metabolomics profiling in detecting intrinsic biological differences among the four groups of samples. Significant variations in the reliability of the metabolomics profiling are identified across laboratories. Importantly, ratio-based metabolomics profiling, by scaling the absolute values of a study sample relative to those of a common reference sample, enables cross-laboratory quantitative data integration. Thus, we construct the ratio-based high-confidence reference datasets between two reference samples, providing "ground truth" for inter-laboratory accuracy assessment, which enables objective evaluation of quantitative metabolomics profiling using various instruments and protocols. CONCLUSIONS: Our study provides the community with rich resources and best practices for inter-laboratory proficiency tests and data integration, ensuring reliability of large-scale and longitudinal metabolomic studies.

Haplotype-resolved assemblies and variant benchmark of a Chinese Quartet.

BACKGROUND: Recent state-of-the-art sequencing technologies enable the investigation of challenging regions in the human genome and expand the scope of variant benchmarking datasets. Herein, we sequence a Chinese Quartet, comprising two monozygotic twin daughters and their biological parents, using four short and long sequencing platforms (Illumina, BGI, PacBio, and Oxford Nanopore Technology). RESULTS: The long reads from the monozygotic twin daughters are phased into paternal and maternal haplotypes using the parent-child genetic map and for each haplotype. We also use long reads to generate haplotype-resolved whole-genome assemblies with completeness and continuity exceeding that of GRCh38. Using this Quartet, we comprehensively catalogue the human variant landscape, generating a dataset of 3,962,453 SNVs, 886,648 indels (< 50 bp), 9726 large deletions (≥ 50 bp), 15,600 large insertions (≥ 50 bp), 40 inversions, 31 complex structural variants, and 68 de novo mutations which are shared between the monozygotic twin daughters. Variants underrepresented in previous benchmarks owing to their complexity-including those located at long repeat regions, complex structural variants, and de novo mutations-are systematically examined in this study. CONCLUSIONS: In summary, this study provides high-quality haplotype-resolved assemblies and a comprehensive set of benchmarking resources for two Chinese monozygotic twin samples which, relative to existing benchmarks, offers expanded genomic coverage and insight into complex variant categories.

Quartet DNA reference materials and datasets for comprehensively evaluating germline variant calling performance.

BACKGROUND: Genomic DNA reference materials are widely recognized as essential for ensuring data quality in omics research. However, relying solely on reference datasets to evaluate the accuracy of variant calling results is incomplete, as they are limited to benchmark regions. Therefore, it is important to develop DNA reference materials that enable the assessment of variant detection performance across the entire genome. RESULTS: We established a DNA reference material suite from four immortalized cell lines derived from a family of parents and monozygotic twins. Comprehensive reference datasets of 4.2 million small variants and 15,000 structural variants were integrated and certified for evaluating the reliability of germline variant calls inside the benchmark regions. Importantly, the genetic built-in-truth of the Quartet family design enables estimation of the precision of variant calls outside the benchmark regions. Using the Quartet reference materials along with study samples, batch effects are objectively monitored and alleviated by training a machine learning model with the Quartet reference datasets to remove potential artifact calls. Moreover, the matched RNA and protein reference materials and datasets from the Quartet project enables cross-omics validation of variant calls from multiomics data. CONCLUSIONS: The Quartet DNA reference materials and reference datasets provide a unique resource for objectively assessing the quality of germline variant calls throughout the whole-genome regions and improving the reliability of large-scale genomic profiling.

The Quartet Data Portal: integration of community-wide resources for multiomics quality control.

The Quartet Data Portal facilitates community access to well-characterized reference materials, reference datasets, and related resources established based on a family of four individuals with identical twins from the Quartet Project. Users can request DNA, RNA, protein, and metabolite reference materials, as well as datasets generated across omics, platforms, labs, protocols, and batches. Reproducible analysis tools allow for objective performance assessment of user-submitted data, while interactive visualization tools support rapid exploration of reference datasets. A closed-loop "distribution-collection-evaluation-integration" workflow enables updates and integration of community-contributed multiomics data. Ultimately, this portal helps promote the advancement of reference datasets and multiomics quality control.

Molecular classification of hormone receptor-positive HER2-negative breast cancer.

Hormone receptor-positive (HR+)/human epidermal growth factor receptor 2-negative (HER2-) breast cancer is the most prevalent type of breast cancer, in which endocrine therapy resistance and distant relapse remain unmet challenges. Accurate molecular classification is urgently required for guiding precision treatment. We established a large-scale multi-omics cohort of 579 patients with HR+/HER2- breast cancer and identified the following four molecular subtypes: canonical luminal, immunogenic, proliferative and receptor tyrosine kinase (RTK)-driven. Tumors of these four subtypes showed distinct biological and clinical features, suggesting subtype-specific therapeutic strategies. The RTK-driven subtype was characterized by the activation of the RTK pathways and associated with poor outcomes. The immunogenic subtype had enriched immune cells and could benefit from immune checkpoint therapy. In addition, we developed convolutional neural network models to discriminate these subtypes based on digital pathology for potential clinical translation. The molecular classification provides insights into molecular heterogeneity and highlights the potential for precision treatment of HR+/HER2- breast cancer.

Quartet RNA reference materials improve the quality of transcriptomic data through ratio-based profiling.

Certified RNA reference materials are indispensable for assessing the reliability of RNA sequencing to detect intrinsically small biological differences in clinical settings, such as molecular subtyping of diseases. As part of the Quartet Project for quality control and data integration of multi-omics profiling, we established four RNA reference materials derived from immortalized B-lymphoblastoid cell lines from four members of a monozygotic twin family. Additionally, we constructed ratio-based transcriptome-wide reference datasets between two samples, providing cross-platform and cross-laboratory 'ground truth'. Investigation of the intrinsically subtle biological differences among the Quartet samples enables sensitive assessment of cross-batch integration of transcriptomic measurements at the ratio level. The Quartet RNA reference materials, combined with the ratio-based reference datasets, can serve as unique resources for assessing and improving the quality of transcriptomic data in clinical and biological settings.

Prediction of base editor off-targets by deep learning.

Due to the tolerance of mismatches between gRNA and targeting sequence, base editors frequently induce unwanted Cas9-dependent off-target mutations. Here, to develop models to predict such off-targets, we design gRNA-off- target pairs for adenine base editors (ABEs) and cytosine base editors (CBEs) and stably integrate them into the human cells. After five days of editing, we obtain valid efficiency datasets of 54,663 and 55,727 off-targets for ABEs and CBEs, respectively. We use the datasets to train deep learning models, resulting in ABEdeepoff and CBEdeepoff, which can predict off-target sites. We use these tools to predict off-targets for a panel of endogenous loci and achieve Spearman correlation values varying from 0.710 to 0.859. Finally, we develop an integrated tool that is freely accessible via an online web server http://www.deephf.com/#/bedeep/bedeepoff . These tools could facilitate minimizing the off-target effects of base editing.

Correcting batch effects in large-scale multiomics studies using a reference-material-based ratio method.

BACKGROUND: Batch effects are notoriously common technical variations in multiomics data and may result in misleading outcomes if uncorrected or over-corrected. A plethora of batch-effect correction algorithms are proposed to facilitate data integration. However, their respective advantages and limitations are not adequately assessed in terms of omics types, the performance metrics, and the application scenarios. RESULTS: As part of the Quartet Project for quality control and data integration of multiomics profiling, we comprehensively assess the performance of seven batch effect correction algorithms based on different performance metrics of clinical relevance, i.e., the accuracy of identifying differentially expressed features, the robustness of predictive models, and the ability of accurately clustering cross-batch samples into their own donors. The ratio-based method, i.e., by scaling absolute feature values of study samples relative to those of concurrently profiled reference material(s), is found to be much more effective and broadly applicable than others, especially when batch effects are completely confounded with biological factors of study interests. We further provide practical guidelines for implementing the ratio based approach in increasingly large-scale multiomics studies. CONCLUSIONS: Multiomics measurements are prone to batch effects, which can be effectively corrected using ratio-based scaling of the multiomics data. Our study lays the foundation for eliminating batch effects at a ratio scale.

High-throughput autoantibody profiling of different stages of Schistosomiasis japonica.

Infection by the Schistosoma japonicum can result in acute, chronic and late-stage manifestations. The latter often presents with severe organ failures and premature death. Importantly, infection can also produce autoimmune phenomena reflected by the development of autoantibodies. We wished to explore and profile the presence of autoantibodies in sera of patients with different stages of S. japonicum infection with the added aim of providing a reference assisting diagnosis. Blood samples from 55 patients with chronic and 20 patients with late-stage schistosomiasis japonica together, with a control group of 50 healthy people were randomly investigated against a microarray of 121 different autoantigens. In addition, the frequency of antibodies against Schistosoma egg antigen (SEA) was examined. In the sera from patients with chronic schistosomiasis japonica, 14 different highly expressed autoantibodies were detected, while patients with late-stage schistosomiasis were found to express as many as 43 autoantibody specificities together with a significantly higher frequency of antibodies against SEA compared to the control group. The findings presented suggest that autoantibody-based classification of schistosomiasis japonica represents a promising approach for the elucidation of subtypes of the disease. This approach may reflect differential disease mechanisms, which could ultimately lead to better treatment.

Overexpression of CD99 is associated with tumor adaptiveness and indicates the tumor recurrence and therapeutic responses in gliomas.

Glioma undergoes adaptive changes, leading to poor prognosis and resistance to treatment. CD99 influences the migration and invasion of glioma cells and plays an oncogene role. However, whether CD99 can affect the adaptiveness of gliomas is still lacking in research, making its clinical value underestimated. Here, we enrolled our in-house and public multiomics datasets for bioinformatic analysis and conducted immunohistochemistry staining to investigate the role of CD99 in glioma adaptive response and its clinical implications. CD99 is expressed in more adaptative glioma subtypes and cell states. Under hypoxic conditions, CD99 is upregulated in glioma cells and is associated with angiogenesis and metabolic adaptations. Gliomas with over-expressed CD99 also increased the immunosuppressive tumor-associated macrophages. The relevance with tumor adaptiveness of CD99 presented clinical significance. We discovered that CD99 overexpression is associated with short-time recurrence and validated its prognostic value. Additionally, Glioma patients with high expression of CD99 were resistant to chemotherapy and radiotherapy. The CD99 expression was also related to anti-angiogenic and immune checkpoint inhibitor therapy response. Inhibitors of the PI3K-AKT pathway have therapeutic potential against CD99-overexpressing gliomas. Our study identified CD99 as a biomarker characterizing the adaptive response in glioma. Gliomas with high CD99 expression are highly tolerant to stress conditions such as hypoxia and antitumor immunity, making treatment responses dimmer and tumor progression. Therefore, for patients with CD99-overexpressing gliomas, tumor adaptiveness should be fully considered during treatment to avoid drug resistance, and closer clinical monitoring should be carried out to improve the prognosis.

Differences of molecular events driving pathological and radiological progression of lung adenocarcinoma.

BACKGROUND: Ground-glass opacity (GGO)-like lung adenocarcinoma (LUAD) has been detected increasingly in the clinic and its inert property and superior survival indicate unique biological characteristics. However, we do not know much about them, which hampers identification of key reasons for the inert property of GGO-like LUAD. METHODS: Using whole-exome sequencing and RNA sequencing, taking into account both radiological and pathological classifications of the same 197 patients concomitantly, we systematically interrogate genes driving the progression from GGO to solid nodule and potential reasons for the inertia of GGO. Using flow cytometry and IHC, we validated the abundance of immune cells and activity of cell proliferation. FINDINGS: Identifying the differences between GGO and solid nodule, we found adenocarcinoma in situ/minimally invasive adenocarcinoma (AIS/MIA) and GGO-like LUAD exhibited lower TP53 mutation frequency and less active cell proliferation-related pathways than solid nodule in LUAD. Identifying the differences in GGO between AIS/MIA and LUAD, we noticed that EGFR mutation frequency and CNV load were significantly higher in LUAD than in AIS/MIA. Regulatory T cell was also higher in LUAD, while CD8+ T cell decreased from AIS/MIA to LUAD. Finally, we constructed a transcriptomic signature to quantify the development from GGO to solid nodule, which was an independent predictor of patients' prognosis in 11 external LUAD datasets. INTERPRETATION: Our results provide deeper insights into the indolent nature of GGO and provide a molecular basis for the treatment of GGO-like LUAD. FUNDING: This study was supported in part by the National Natural Science Foundation of China (32170657), the National Natural Science Foundation of China (82203037), and Shanghai Sailing Program (22YF1408900).

De novo assembly of a chromosome-level reference genome of the ornamental butterfly Sericinus montelus based on nanopore sequencing and Hi-C analysis.

Sericinus montelus (Lepidoptera, Papilionidae, Parnassiinae) is a high-value ornamental swallowtail butterfly species widely distributed in Northern and Central China, Japan, Korea, and Russia. The larval stage of this species feeds exclusively on Aristolochia plants. The Aristolochia species is well known for its high levels of aristolochic acids (AAs), which have been found to be carcinogenic for numerous animals. The swallowtail butterfly is among the few that can feed on these toxic host plants. However, the genetic adaptation of S. montelus to confer new abilities for AA tolerance has not yet been well explored, largely due to the limited genomic resources of this species. This study aimed to present a chromosome-level reference genome for S. montelus using the Oxford Nanopore long-read sequencing, Illumina short-read sequencing, and Hi-C technology. The final assembly was composed of 581.44 Mb with an expected genome size of 619.27 Mb. Further, 99.98% of the bases could be anchored onto 30 chromosomes. The N50 of contigs and scaffolds was 5.74 and 19.12 Mb, respectively. Approximately 48.86% of the assembled genome was suggested to be repeat elements, and 13,720 protein-coding genes were predicted in the current assembly. The phylogenetic analysis indicated that S. montelus diverged from the common ancestor of swallowtails about 58.57-80.46 million years ago. Compared with related species, S. montelus showed a significant expansion of P450 gene family members, and positive selections on eloa, heatr1, and aph1a resulted in the AA tolerance for S. montelus larva. The de novo assembly of a high-quality reference genome for S. montelus provided a fundamental genomic tool for future research on evolution, genome genetics, and toxicology of the swallowtail butterflies.

PAnno: A pharmacogenomics annotation tool for clinical genomic testing.

Introduction: Next-generation sequencing (NGS) technologies have been widely used in clinical genomic testing for drug response phenotypes. However, the inherent limitations of short reads make accurate inference of diplotypes still challenging, which may reduce the effectiveness of genotype-guided drug therapy. Methods: An automated Pharmacogenomics Annotation tool (PAnno) was implemented, which reports prescribing recommendations and phenotypes by parsing the germline variant call format (VCF) file from NGS and the population to which the individual belongs. Results: A ranking model dedicated to inferring diplotypes, developed based on the allele (haplotype) definition and population allele frequency, was introduced in PAnno. The predictive performance was validated in comparison with four similar tools using the consensus diplotype data of the Genetic Testing Reference Materials Coordination Program (GeT-RM) as ground truth. An annotation method was proposed to summarize prescribing recommendations and classify drugs into avoid use, use with caution, and routine use, following the recommendations of the Clinical Pharmacogenetics Implementation Consortium (CPIC), etc. It further predicts phenotypes of specific drugs in terms of toxicity, dosage, efficacy, and metabolism by integrating the high-confidence clinical annotations in the Pharmacogenomics Knowledgebase (PharmGKB). PAnno is available at https://github.com/PreMedKB/PAnno. Discussion: PAnno provides an end-to-end clinical pharmacogenomics decision support solution by resolving, annotating, and reporting germline variants.

Patient Assessment and Therapy Planning Based on Homologous Recombination Repair Deficiency.

Defects in genes involved in the DNA damage response cause homologous recombination repair deficiency (HRD). HRD is found in a subgroup of cancer patients for several tumor types, and it has a clinical relevance to cancer prevention and therapies. Accumulating evidence has identified HRD as a biomarker for assessing the therapeutic response of tumor cells to poly(ADP-ribose) polymerase inhibitors and platinum-based chemotherapies. Nevertheless, the biology of HRD is complex, and its applications and the benefits of different HRD biomarker assays are controversial. This is primarily due to inconsistencies in HRD assessments and definitions (gene-level tests, genomic scars, mutational signatures, or a combination of these methods) and difficulties in assessing the contribution of each genomic event. Therefore, we aim to review the biological rationale and clinical evidence of HRD as a biomarker. This review provides a blueprint for the standardization and harmonization of HRD assessments.

Tackling the translational challenges of multi-omics research in the realm of European personalised medicine: A workshop report.

Personalised medicine (PM) presents a great opportunity to improve the future of individualised healthcare. Recent advances in -omics technologies have led to unprecedented efforts characterising the biology and molecular mechanisms that underlie the development and progression of a wide array of complex human diseases, supporting further development of PM. This article reflects the outcome of the 2021 EATRIS-Plus Multi-omics Stakeholder Group workshop organised to 1) outline a global overview of common promises and challenges that key European stakeholders are facing in the field of multi-omics research, 2) assess the potential of new technologies, such as artificial intelligence (AI), and 3) establish an initial dialogue between key initiatives in this space. Our focus is on the alignment of agendas of European initiatives in multi-omics research and the centrality of patients in designing solutions that have the potential to advance PM in long-term healthcare strategies.