Population-scale tissue transcriptomics maps long non-coding RNAs to complex disease.

Long non-coding RNA (lncRNA) genes have well-established and important impacts on molecular and cellular functions. However, among the thousands of lncRNA genes, it is still a major challenge to identify the subset with disease or trait relevance. To systematically characterize these lncRNA genes, we used Genotype Tissue Expression (GTEx) project v8 genetic and multi-tissue transcriptomic data to profile the expression, genetic regulation, cellular contexts, and trait associations of 14,100 lncRNA genes across 49 tissues for 101 distinct complex genetic traits. Using these approaches, we identified 1,432 lncRNA gene-trait associations, 800 of which were not explained by stronger effects of neighboring protein-coding genes. This included associations between lncRNA quantitative trait loci and inflammatory bowel disease, type 1 and type 2 diabetes, and coronary artery disease, as well as rare variant associations to body mass index.

Organ aging signatures in the plasma proteome track health and disease.

Animal studies show aging varies between individuals as well as between organs within an individual1-4, but whether this is true in humans and its effect on age-related diseases is unknown. We utilized levels of human blood plasma proteins originating from specific organs to measure organ-specific aging differences in living individuals. Using machine learning models, we analysed aging in 11 major organs and estimated organ age reproducibly in five independent cohorts encompassing 5,676 adults across the human lifespan. We discovered nearly 20% of the population show strongly accelerated age in one organ and 1.7% are multi-organ agers. Accelerated organ aging confers 20-50% higher mortality risk, and organ-specific diseases relate to faster aging of those organs. We find individuals with accelerated heart aging have a 250% increased heart failure risk and accelerated brain and vascular aging predict Alzheimer's disease (AD) progression independently from and as strongly as plasma pTau-181 (ref. 5), the current best blood-based biomarker for AD. Our models link vascular calcification, extracellular matrix alterations and synaptic protein shedding to early cognitive decline. We introduce a simple and interpretable method to study organ aging using plasma proteomics data, predicting diseases and aging effects.

Interaction molecular QTL mapping discovers cellular and environmental modifiers of genetic regulatory effects.

Bulk-tissue molecular quantitative trait loci (QTLs) have been the starting point for interpreting disease-associated variants, and context-specific QTLs show particular relevance for disease. Here, we present the results of mapping interaction QTLs (iQTLs) for cell type, age, and other phenotypic variables in multi-omic, longitudinal data from the blood of individuals of diverse ancestries. By modeling the interaction between genotype and estimated cell-type proportions, we demonstrate that cell-type iQTLs could be considered as proxies for cell-type-specific QTL effects, particularly for the most abundant cell type in the tissue. The interpretation of age iQTLs, however, warrants caution because the moderation effect of age on the genotype and molecular phenotype association could be mediated by changes in cell-type composition. Finally, we show that cell-type iQTLs contribute to cell-type-specific enrichment of diseases that, in combination with additional functional data, could guide future functional studies. Overall, this study highlights the use of iQTLs to gain insights into the context specificity of regulatory effects.

Clonal haematopoiesis of indeterminate potential and atrial fibrillation: an east Asian cohort study.

BACKGROUND AND AIMS: Both clonal haematopoiesis of indeterminate potential (CHIP) and atrial fibrillation (AF) are age-related conditions. This study investigated the potential role of CHIP in the development and progression of AF. METHODS: Deep-targeted sequencing of 24 CHIP mutations (a mean depth of coverage = 1000×) was performed in 1004 patients with AF and 3341 non-AF healthy subjects. Variant allele fraction ≥ 2.0% indicated the presence of CHIP mutations. The association between CHIP and AF was evaluated by the comparison of (i) the prevalence of CHIP mutations between AF and non-AF subjects and (ii) clinical characteristics discriminated by CHIP mutations within AF patients. Furthermore, the risk of clinical outcomes-the composite of heart failure, ischaemic stroke, or death-according to the presence of CHIP mutations in AF was investigated from the UK Biobank cohort. RESULTS: The mean age was 67.6 ± 6.9 vs. 58.5 ± 6.5 years in AF (paroxysmal, 39.0%; persistent, 61.0%) and non-AF cohorts, respectively. CHIP mutations with a variant allele fraction of ≥2.0% were found in 237 (23.6%) AF patients (DNMT3A, 13.5%; TET2, 6.6%; and ASXL1, 1.5%) and were more prevalent than non-AF subjects [356 (10.7%); P < .001] across the age. After multivariable adjustment (age, sex, smoking, body mass index, diabetes, and hypertension), CHIP mutations were 1.4-fold higher in AF [adjusted odds ratio (OR) 1.38; 95% confidence interval 1.10-1.74, P < .01]. The ORs of CHIP mutations were the highest in the long-standing persistent AF (adjusted OR 1.50; 95% confidence interval 1.14-1.99, P = .004) followed by persistent (adjusted OR 1.44) and paroxysmal (adjusted OR 1.33) AF. In gene-specific analyses, TET2 somatic mutation presented the highest association with AF (adjusted OR 1.65; 95% confidence interval 1.05-2.60, P = .030). AF patients with CHIP mutations were older and had a higher prevalence of diabetes, a longer AF duration, a higher E/E', and a more severely enlarged left atrium than those without CHIP mutations (all P < .05). In UK Biobank analysis of 21 286 AF subjects (1297 with CHIP and 19 989 without CHIP), the CHIP mutation in AF is associated with a 1.32-fold higher risk of a composite clinical event (heart failure, ischaemic stroke, or death). CONCLUSIONS: CHIP mutations, primarily DNMT3A or TET2, are more prevalent in patients with AF than non-AF subjects whilst their presence is associated with a more progressive nature of AF and unfavourable clinical outcomes.

Nonsense-mediated decay is highly stable across individuals and tissues.

Precise interpretation of the effects of rare protein-truncating variants (PTVs) is important for accurate determination of variant impact. Current methods for assessing the ability of PTVs to induce nonsense-mediated decay (NMD) focus primarily on the position of the variant in the transcript. We used RNA sequencing of the Genotype Tissue Expression v.8 cohort to compute the efficiency of NMD using allelic imbalance for 2,320 rare (genome aggregation database minor allele frequency ≤ 1%) PTVs across 809 individuals in 49 tissues. We created an interpretable predictive model using penalized logistic regression in order to evaluate the comprehensive influence of variant annotation, tissue, and inter-individual variation on NMD. We found that variant position, allele frequency, the inclusion of ultra-rare and singleton variants, and conservation were predictive of allelic imbalance. Furthermore, we found that NMD effects were highly concordant across tissues and individuals. Due to this high consistency, we demonstrate in silico that utilizing peripheral tissues or cell lines provides accurate prediction of NMD for PTVs.

An integrated approach to identify environmental modulators of genetic risk factors for complex traits.

Complex traits and diseases can be influenced by both genetics and environment. However, given the large number of environmental stimuli and power challenges for gene-by-environment testing, it remains a critical challenge to identify and prioritize specific disease-relevant environmental exposures. We propose a framework for leveraging signals from transcriptional responses to environmental perturbations to identify disease-relevant perturbations that can modulate genetic risk for complex traits and inform the functions of genetic variants associated with complex traits. We perturbed human skeletal-muscle-, fat-, and liver-relevant cell lines with 21 perturbations affecting insulin resistance, glucose homeostasis, and metabolic regulation in humans and identified thousands of environmentally responsive genes. By combining these data with GWASs from 31 distinct polygenic traits, we show that the heritability of multiple traits is enriched in regions surrounding genes responsive to specific perturbations and, further, that environmentally responsive genes are enriched for associations with specific diseases and phenotypes from the GWAS Catalog. Overall, we demonstrate the advantages of large-scale characterization of transcriptional changes in diversely stimulated and pathologically relevant cells to identify disease-relevant perturbations.

Transcriptional profiling of isogenic Friedreich ataxia neurons and effect of an HDAC inhibitor on disease signatures.

Friedreich ataxia (FRDA) is a neurodegenerative disorder caused by transcriptional silencing of the frataxin (FXN) gene, resulting in loss of the essential mitochondrial protein frataxin. Based on the knowledge that a GAA·TTC repeat expansion in the first intron of FXN induces heterochromatin, we previously showed that 2-aminobenzamide-type histone deacetylase inhibitors (HDACi) increase FXN mRNA levels in induced pluripotent stem cell (iPSC)-derived FRDA neurons and in circulating lymphocytes from patients after HDACi oral administration. How the reduced expression of frataxin leads to neurological and other systemic symptoms in FRDA patients remains unclear. Similar to other triplet-repeat disorders, it is unknown why FRDA affects only specific cell types, primarily the large sensory neurons of the dorsal root ganglia and cardiomyocytes. The combination of iPSC technology and genome-editing techniques offers the unique possibility to address these questions in a relevant cell model of FRDA, obviating confounding effects of variable genetic backgrounds. Here, using "scarless" gene-editing methods, we created isogenic iPSC lines that differ only in the length of the GAA·TTC repeats. To uncover the gene expression signatures due to the GAA·TTC repeat expansion in FRDA neuronal cells and the effect of HDACi on these changes, we performed RNA-seq-based transcriptomic analysis of iPSC-derived central nervous system (CNS) and isogenic sensory neurons. We found that cellular pathways related to neuronal function, regulation of transcription, extracellular matrix organization, and apoptosis are affected by frataxin loss in neurons of the CNS and peripheral nervous system and that these changes are partially restored by HDACi treatment.

Peripheral blood gene expression reveals an inflammatory transcriptomic signature in Friedreich's ataxia patients.

Transcriptional changes in Friedreich's ataxia (FRDA), a rare and debilitating recessive Mendelian neurodegenerative disorder, have been studied in affected but inaccessible tissues-such as dorsal root ganglia, sensory neurons and cerebellum-in animal models or small patient series. However, transcriptional changes induced by FRDA in peripheral blood, a readily accessible tissue, have not been characterized in a large sample. We used differential expression, association with disability stage, network analysis and enrichment analysis to characterize the peripheral blood transcriptome and identify genes that were differentially expressed in FRDA patients (n = 418) compared with both heterozygous expansion carriers (n = 228) and controls (n = 93 739 individuals in total), or were associated with disease progression, resulting in a disease signature for FRDA. We identified a transcriptional signature strongly enriched for an inflammatory innate immune response. Future studies should seek to further characterize the role of peripheral inflammation in FRDA pathology and determine its relevance to overall disease progression.

regionalpcs: improved discovery of DNA methylation associations with complex traits.

We have developed the regional principal components (rPCs) method, a novel approach for summarizing gene-level methylation. rPCs address the challenge of deciphering complex epigenetic mechanisms in diseases like Alzheimer's disease (AD). In contrast to traditional averaging, rPCs leverage principal components analysis to capture complex methylation patterns across gene regions. Our method demonstrated a 54% improvement in sensitivity over averaging in simulations, offering a robust framework for identifying subtle epigenetic variations. Applying rPCs to the AD brain methylation data in ROSMAP, combined with cell type deconvolution, we uncovered 838 differentially methylated genes associated with neuritic plaque burden-significantly outperforming conventional methods. Integrating methylation quantitative trait loci (meQTL) with genome-wide association studies (GWAS) identified 17 genes with potential causal roles in AD, including MS4A4A and PICALM. Our approach is available in the Bioconductor package regionalpcs, opening avenues for research and facilitating a deeper understanding of the epigenetic landscape in complex diseases.

PRC2-AgeIndex as a universal biomarker of aging and rejuvenation.

DNA methylation (DNAm) is one of the most reliable biomarkers of aging across mammalian tissues. While the age-dependent global loss of DNAm has been well characterized, DNAm gain is less characterized. Studies have demonstrated that CpGs which gain methylation with age are enriched in Polycomb Repressive Complex 2 (PRC2) targets. However, whole-genome examination of all PRC2 targets as well as determination of the pan-tissue or tissue-specific nature of these associations is lacking. Here, we show that low-methylated regions (LMRs) which are highly bound by PRC2 in embryonic stem cells (PRC2 LMRs) gain methylation with age in all examined somatic mitotic cells. We estimated that this epigenetic change represents around 90% of the age-dependent DNAm gain genome-wide. Therefore, we propose the "PRC2-AgeIndex," defined as the average DNAm in PRC2 LMRs, as a universal biomarker of cellular aging in somatic cells which can distinguish the effect of different anti-aging interventions.

Epigenetic and proteomic signatures associate with clonal hematopoiesis expansion rate.

Clonal hematopoiesis of indeterminate potential (CHIP), whereby somatic mutations in hematopoietic stem cells confer a selective advantage and drive clonal expansion, not only correlates with age but also confers increased risk of morbidity and mortality. Here, we leverage genetically predicted traits to identify factors that determine CHIP clonal expansion rate. We used the passenger-approximated clonal expansion rate method to quantify the clonal expansion rate for 4,370 individuals in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) cohort and calculated polygenic risk scores for DNA methylation aging, inflammation-related measures and circulating protein levels. Clonal expansion rate was significantly associated with both genetically predicted and measured epigenetic clocks. No associations were identified with inflammation-related lab values or diseases and CHIP expansion rate overall. A proteome-wide search identified predicted circulating levels of myeloid zinc finger 1 and anti-Müllerian hormone as associated with an increased CHIP clonal expansion rate and tissue inhibitor of metalloproteinase 1 and glycine N-methyltransferase as associated with decreased CHIP clonal expansion rate. Together, our findings identify epigenetic and proteomic patterns associated with the rate of hematopoietic clonal expansion.

Interaction molecular QTL mapping discovers cellular and environmental modifiers of genetic regulatory effects.

Bulk tissue molecular quantitative trait loci (QTLs) have been the starting point for interpreting disease-associated variants, while context-specific QTLs show particular relevance for disease. Here, we present the results of mapping interaction QTLs (iQTLs) for cell type, age, and other phenotypic variables in multi-omic, longitudinal data from blood of individuals of diverse ancestries. By modeling the interaction between genotype and estimated cell type proportions, we demonstrate that cell type iQTLs could be considered as proxies for cell type-specific QTL effects. The interpretation of age iQTLs, however, warrants caution as the moderation effect of age on the genotype and molecular phenotype association may be mediated by changes in cell type composition. Finally, we show that cell type iQTLs contribute to cell type-specific enrichment of diseases that, in combination with additional functional data, may guide future functional studies. Overall, this study highlights iQTLs to gain insights into the context-specificity of regulatory effects.

Multiset correlation and factor analysis enables exploration of multi-omics data.

Multi-omics datasets are becoming more common, necessitating better integration methods to realize their revolutionary potential. Here, we introduce multi-set correlation and factor analysis (MCFA), an unsupervised integration method tailored to the unique challenges of high-dimensional genomics data that enables fast inference of shared and private factors. We used MCFA to integrate methylation markers, protein expression, RNA expression, and metabolite levels in 614 diverse samples from the Trans-Omics for Precision Medicine/Multi-Ethnic Study of Atherosclerosis multi-omics pilot. Samples cluster strongly by ancestry in the shared space, even in the absence of genetic information, while private spaces frequently capture dataset-specific technical variation. Finally, we integrated genetic data by conducting a genome-wide association study (GWAS) of our inferred factors, observing that several factors are enriched for GWAS hits and trans-expression quantitative trait loci. Two of these factors appear to be related to metabolic disease. Our study provides a foundation and framework for further integrative analysis of ever larger multi-modal genomic datasets.

Clonal hematopoiesis is associated with protection from Alzheimer's disease.

Clonal hematopoiesis of indeterminate potential (CHIP) is a premalignant expansion of mutated hematopoietic stem cells. As CHIP-associated mutations are known to alter the development and function of myeloid cells, we hypothesized that CHIP may also be associated with the risk of Alzheimer's disease (AD), a disease in which brain-resident myeloid cells are thought to have a major role. To perform association tests between CHIP and AD dementia, we analyzed blood DNA sequencing data from 1,362 individuals with AD and 4,368 individuals without AD. Individuals with CHIP had a lower risk of AD dementia (meta-analysis odds ratio (OR) = 0.64, P = 3.8 × 10-5), and Mendelian randomization analyses supported a potential causal association. We observed that the same mutations found in blood were also detected in microglia-enriched fraction of the brain in seven of eight CHIP carriers. Single-nucleus chromatin accessibility profiling of brain-derived nuclei in six CHIP carriers revealed that the mutated cells comprised a large proportion of the microglial pool in the samples examined. While additional studies are required to validate the mechanistic findings, these results suggest that CHIP may have a role in attenuating the risk of AD.

Transcriptomics and chromatin accessibility in multiple African population samples.

Mapping the functional human genome and impact of genetic variants is often limited to European-descendent population samples. To aid in overcoming this limitation, we measured gene expression using RNA sequencing in lymphoblastoid cell lines (LCLs) from 599 individuals from six African populations to identify novel transcripts including those not represented in the hg38 reference genome. We used whole genomes from the 1000 Genomes Project and 164 Maasai individuals to identify 8,881 expression and 6,949 splicing quantitative trait loci (eQTLs/sQTLs), and 2,611 structural variants associated with gene expression (SV-eQTLs). We further profiled chromatin accessibility using ATAC-Seq in a subset of 100 representative individuals, to identity chromatin accessibility quantitative trait loci (caQTLs) and allele-specific chromatin accessibility, and provide predictions for the functional effect of 78.9 million variants on chromatin accessibility. Using this map of eQTLs and caQTLs we fine-mapped GWAS signals for a range of complex diseases. Combined, this work expands global functional genomic data to identify novel transcripts, functional elements and variants, understand population genetic history of molecular quantitative trait loci, and further resolve the genetic basis of multiple human traits and disease.

Integration of genetic colocalizations with physiological and pharmacological perturbations identifies cardiometabolic disease genes.

BACKGROUND: Identification of causal genes for polygenic human diseases has been extremely challenging, and our understanding of how physiological and pharmacological stimuli modulate genetic risk at disease-associated loci is limited. Specifically, insulin resistance (IR), a common feature of cardiometabolic disease, including type 2 diabetes, obesity, and dyslipidemia, lacks well-powered genome-wide association studies (GWAS), and therefore, few associated loci and causal genes have been identified. METHODS: Here, we perform and integrate linkage disequilibrium (LD)-adjusted colocalization analyses across nine cardiometabolic traits (fasting insulin, fasting glucose, insulin sensitivity, insulin sensitivity index, type 2 diabetes, triglycerides, high-density lipoprotein, body mass index, and waist-hip ratio) combined with expression and splicing quantitative trait loci (eQTLs and sQTLs) from five metabolically relevant human tissues (subcutaneous and visceral adipose, skeletal muscle, liver, and pancreas). To elucidate the upstream regulators and functional mechanisms for these genes, we integrate their transcriptional responses to 21 relevant physiological and pharmacological perturbations in human adipocytes, hepatocytes, and skeletal muscle cells and map their protein-protein interactions. RESULTS: We identify 470 colocalized loci and prioritize 207 loci with a single colocalized gene. Patterns of shared colocalizations across traits and tissues highlight different potential roles for colocalized genes in cardiometabolic disease and distinguish several genes involved in pancreatic ß-cell function from others with a more direct role in skeletal muscle, liver, and adipose tissues. At the loci with a single colocalized gene, 42 of these genes were regulated by insulin and 35 by glucose in perturbation experiments, including 17 regulated by both. Other metabolic perturbations regulated the expression of 30 more genes not regulated by glucose or insulin, pointing to other potential upstream regulators of candidate causal genes. CONCLUSIONS: Our use of transcriptional responses under metabolic perturbations to contextualize genetic associations from our custom colocalization approach provides a list of likely causal genes and their upstream regulators in the context of IR-associated cardiometabolic risk.

Methylation differences in Alzheimer's disease neuropathologic change in the aged human brain.

Alzheimer's disease (AD) is the most common cause of dementia with advancing age as its strongest risk factor. AD neuropathologic change (ADNC) is known to be associated with numerous DNA methylation changes in the human brain, but the oldest old (> 90 years) have so far been underrepresented in epigenetic studies of ADNC. Our study participants were individuals aged over 90 years (n = 47) from The 90+ Study. We analyzed DNA methylation from bulk samples in eight precisely dissected regions of the human brain: middle frontal gyrus, cingulate gyrus, entorhinal cortex, dentate gyrus, CA1, substantia nigra, locus coeruleus and cerebellar cortex. We deconvolved our bulk data into cell-type-specific (CTS) signals using computational methods. CTS methylation differences were analyzed across different levels of ADNC. The highest amount of ADNC related methylation differences was found in the dentate gyrus, a region that has so far been underrepresented in large scale multi-omic studies. In neurons of the dentate gyrus, DNA methylation significantly differed with increased burden of amyloid beta (Aß) plaques at 5897 promoter regions of protein-coding genes. Amongst these, higher Aß plaque burden was associated with promoter hypomethylation of the Presenilin enhancer 2 (PEN-2) gene, one of the rate limiting genes in the formation of gamma-secretase, a multicomponent complex that is responsible in part for the endoproteolytic cleavage of amyloid precursor protein into Aß peptides. In addition to novel ADNC related DNA methylation changes, we present the most detailed array-based methylation survey of the old aged human brain to date. Our open-sourced dataset can serve as a brain region reference panel for future studies and help advance research in aging and neurodegenerative diseases.

Clonal hematopoiesis of indeterminate potential, DNA methylation, and risk for coronary artery disease.

Age-related changes to the genome-wide DNA methylation (DNAm) pattern observed in blood are well-documented. Clonal hematopoiesis of indeterminate potential (CHIP), characterized by the age-related acquisition and expansion of leukemogenic mutations in hematopoietic stem cells (HSCs), is associated with blood cancer and coronary artery disease (CAD). Epigenetic regulators DNMT3A and TET2 are the two most frequently mutated CHIP genes. Here, we present results from an epigenome-wide association study for CHIP in 582 Cardiovascular Health Study (CHS) participants, with replication in 2655 Atherosclerosis Risk in Communities (ARIC) Study participants. We show that DNMT3A and TET2 CHIP have distinct and directionally opposing genome-wide DNAm association patterns consistent with their regulatory roles, albeit both promoting self-renewal of HSCs. Mendelian randomization analyses indicate that a subset of DNAm alterations associated with these two leading CHIP genes may promote the risk for CAD.

GADD45A is a protective modifier of neurogenic skeletal muscle atrophy.

Neurogenic muscle atrophy is the loss of skeletal muscle mass and function that occurs with nerve injury and in denervating diseases, such as amyotrophic lateral sclerosis. Aside from prompt restoration of innervation and exercise where feasible, there are currently no effective strategies for maintaining skeletal muscle mass in the setting of denervation. We conducted a longitudinal analysis of gene expression changes occurring in atrophying skeletal muscle and identified growth arrest and DNA damage-inducible A (Gadd45a) as a gene that shows one of the earliest and most sustained increases in expression in skeletal muscle after denervation. We evaluated the role of this induction using genetic mouse models and found that mice lacking GADD45A showed accelerated and exacerbated neurogenic muscle atrophy, as well as loss of fiber type identity. Our genetic analyses demonstrate that, rather than directly contributing to muscle atrophy as proposed in earlier studies, GADD45A induction likely represents a protective negative feedback response to denervation. Establishing the downstream effectors that mediate this protective effect and the pathways they participate in may yield new opportunities to modify the course of muscle atrophy.

Transcriptomic profiling of whole blood in 22q11.2 reciprocal copy number variants reveals that cell proportion highly impacts gene expression.

22q11.2 reciprocal copy number variants (CNVs) offer a powerful quasi-experimental "reverse-genetics" paradigm to elucidate how gene dosage (i.e., deletions and duplications) disrupts the transcriptome to cause further downstream effects. Clinical profiles of 22q11.2 CNV carriers indicate that disrupted gene expression causes alterations in neuroanatomy, cognitive function, and psychiatric disease risk. However, interpreting transcriptomic signal in bulk tissue requires careful consideration of potential changes in cell composition. We first characterized transcriptomic dysregulation in peripheral blood from reciprocal 22q11.2 CNV carriers using differential expression analysis and weighted gene co-expression network analysis (WGCNA) to identify modules of co-expressed genes. We also assessed for group differences in cell composition and re-characterized transcriptomic differences after accounting for cell type proportions and medication usage. Finally, to explore whether CNV-related transcriptomic changes relate to downstream phenotypes associated with 22q11.2 CNVs, we tested for associations of gene expression with neuroimaging measures and behavioral traits, including IQ and psychosis or ASD diagnosis. 22q11.2 deletion carriers (22qDel) showed widespread expression changes at the individual gene as well as module eigengene level compared to 22q11.2 duplication carriers (22qDup) and controls. 22qDup showed increased expression of 5 genes within the 22q11.2 locus, and CDH6 located outside of the locus. Downregulated modules in 22qDel implicated altered immune and inflammatory processes. Celltype deconvolution analyses revealed significant differences between CNV and control groups in T-cell, mast cell, and macrophage proportions; differential expression of individual genes between groups was substantially attenuated after adjusting for cell composition. Individual gene, module eigengene, and cell proportions were not significantly associated with psychiatric or neuroanatomic traits. Our findings suggest broad immune-related dysfunction in 22qDel and highlight the importance of understanding differences in cell composition when interpreting transcriptomic changes in clinical populations. Results also suggest novel directions for future investigation to test whether 22q11.2 CNV effects on macrophages have implications for brain-related microglial function that may contribute to psychiatric phenotypes in 22q11.2 CNV carriers.