Abstract:Integrative approaches that simultaneously model multi-omics data have gained increasing popularity because they provide holistic system biology views of multiple or all components in a biological system of interest. Canonical correlation analysis (CCA) is a correlation-based integrative method. It was initially designed to extract latent features shared between two assays by finding the linear combinations of features – referred to as canonical vectors (CVs) – within each assay that achieve maximal across-assay correlation. Sparse multiple CCA (SMCCA), a widely-used derivative of CCA, allows more than two assays but can result in non-orthogonal CVs when applied to high-dimensional data. Here, we incorporated a variation of the Gram-Schmidt (GS) algorithm with SMCCA to improve orthogonality among CVs. Applying our SMCCA-GS method to proteomics and methylomics data from the Multi-Ethnic Study of Atherosclerosis (MESA) and Jackson Heart Study (JHS), we identified strong associations between blood cell counts and protein abundance. This finding suggests that adjustment of blood cell composition should be considered in protein-based association studies. Importantly, CVs obtained from two independent cohorts demonstrate transferability across the cohorts. For example, proteomic CVs learned from JHS explain similar amounts of blood cell count phenotypic variance in MESA, explaining 39.0% ~ 50.0% variation in JHS and 38.9% ~ 49.1% in MESA, similar transferability was observed for other omics-CV-trait pairs. This suggests that biologically meaningful and cohort-agnostic variation is captured by CVs. We further developed Sparse Supervised Multiple CCA (SSMCCA) to allow supervised integration analysis for more than two assays. We anticipate that applying our SMCCA-GS and SSMCCA on various cohorts would help identify cohort-agnostic biologically meaningful relationships between multi-omics data and phenotypic traits. Author Summary Comprehensive understanding of human complex traits may benefit from incorporation of molecular features from multiple biological layers such as genome, epigenome, transcriptome, proteome, and metabolome. CCA is a correlation-based method for multi-omics data which reduces the dimension of each omic assay to several orthogonal components – commonly referred to as canonical vectors (CVs). The widely-used SMCCA method allows effective dimension reduction and integration of multi-omics data, but suffers from potentially highly correlated CVs when applied to high-dimensional omics data. Here, we improve the statistical independence among the CVs by adopting a variation of the GS algorithm. We applied our SMCCA-GS method to proteomic and methylomic data from two cohort studies, MESA and JHS. Our results reveal a pronounced effect of blood cell counts on protein abundance, strongly suggesting blood cell composition adjustment in protein-based association studies may be necessary. Finally, we present SSMCCA which allows supervised CCA analysis for the association between one phenotype of interest and more than two assays. We anticipate that SMCCA-GS would help reveal meaningful system-level factors from biological processes involving features from multiple assays; and SSMCCA would further empower interrogation of these factors for phenotypic traits related to health and diseases. ### Competing Interest Statement LMR is a consultant for the TOPMed Administrative Coordinating Center (through WeStat).

MOCHA: advanced statistical modeling of scATAC-seq data enables functional genomic inference in large human disease cohorts

MOCHA's advanced statistical modeling of scATAC-seq data enables functional genomic inference in large human cohorts

Functional Inference of Gene Regulation Using Single-Cell Multi-Omics

Single-cell multi-ome regression models identify functional and disease-associated enhancers and enable chromatin potential analysis

An ultra high-throughput method for single-cell joint analysis of open chromatin and transcriptome

Profiling the genetic determinants of chromatin accessibility with scalable single-cell CRISPR screens

scAMACE: Model-based approach to the joint analysis of single-cell data on chromatin accessibility, gene expression and methylation

SnapATAC: A Comprehensive Analysis Package for Single Cell ATAC-seq

Uncovering Cross-Cohort Molecular Features with Multi-Omics Integration Analysis

Massively parallel single-cell chromatin landscapes of human immune cell development and intratumoral T cell exhaustion

Mitochondrial single-cell ATAC-seq for high-throughput multi-omic detection of mitochondrial genotypes and chromatin accessibility

Efficient Generation of Paired Single‐Cell Multiomics Profiles by Deep Learning

Stabilized marker gene identification and functional annotation from single-cell transcriptomic data

SCALA: A complete solution for multimodal analysis of single-cell Next Generation Sequencing data

Single-cell CUT&Tag analysis of chromatin modifications in differentiation and tumor progression

A component overlapping attribute clustering (COAC) algorithm for single-cell RNA sequencing data analysis and potential pathobiological implications

CrossMP: Enabling Cross-Modality Translation between Single-Cell RNA-Seq and Single-Cell ATAC-Seq through Web-Based Portal

Abstract LB234: Identifying cancer vulnerabilities associated with changes in chromatin accessibility by simultaneously profiling hundreds of cancer cell lines with ATAC-seq

SMGR: a joint statistical method for integrative analysis of single-cell multi-omics data

Abstract 4411: Enhancing scalability and consistency in clinical multiomics via an optimized fixed cell ATAC-seq method

Enhancement and Imputation of Peak Signal Enables Accurate Cell-Type Classification in scATAC-seq