Study / S5QHT8DXE2024-02-14

Integrated multiplexed assays of variant effect reveal determinants of catechol-O-methyltransferase gene expression

Ian Hoskins, Shilpa Rao, Charisma Tante, Can Cenik

About this study

Multiplexed assays of variant effect are powerful methods to profile the consequences of rare variants on gene expression and organismal fitness. Yet, few studies have integrated several multiplexed assays to map variant effects on gene expression in coding sequences. Here, we pioneered a multiplexed assay based on polysome profiling to measure variant effects on translation at scale, uncovering single-nucleotide variants that increase or decrease ribosome load. By combining high-throughput ribosome load data with multiplexed mRNA and protein abundance readouts, we mapped the cis-regulatory landscape of thousands of catechol-O-methyltransferase (COMT) variants from RNA to protein and found numerous coding variants that alter COMT expression. Finally, we trained machine learning models to map signatures of variant effects on COMT gene expression and uncovered both directional and divergent impacts across expression layers. Our analyses reveal expression phenotypes for thousands of variants in COMT and highlight variant effects on both single and multiple layers of expression. Our findings prompt future studies that integrate several multiplexed assays for the readout of gene expression.

Full author list & citation

Ian Hoskins, Shilpa Rao, Charisma Tante, Can Cenik. Integrated multiplexed assays of variant effect reveal determinants of catechol-O-methyltransferase gene expression. 2024-02-14. https://doi.org/10.1038/s44320-024-00018-9

Experiments 4

E006ZS8NV

COMT heavy-polysome ribosome-load MAVE (F4)

A codon-saturation COMT library covering MB-COMT codons 40–74 and 136–158 was stably integrated into the HEK293T landing-pad cell line. Four biological replicates were profiled by amplicon sequencing, and this table reports the published ALDEx2 variant effects for polysome metafraction 4 (heavy polysome) relative to total RNA.

Targeted Genomic Integration MPRAHumanGRCh38
Explore data
E11UT7DTO

COMT protein abundance flow-cytometry MAVE

A codon-saturation COMT library covering MB-COMT codons 40–74 and 136–158 was stably integrated into the HEK293T landing-pad cell line. Four biological replicates were sorted by flow cytometry and profiled by amplicon sequencing; this table reports the published ALDEx2 effects for the low-protein flow population P3 relative to total genomic DNA.

Targeted Genomic Integration MPRAHumanGRCh38
Explore data
E6VXT3E2N

COMT total RNA abundance MAVE

A codon-saturation COMT library covering MB-COMT codons 40–74 and 136–158 was stably integrated into the HEK293T landing-pad cell line. Four biological replicates were profiled by amplicon sequencing, and this table reports the published ALDEx2 variant effects for total RNA abundance relative to total genomic DNA.

Targeted Genomic Integration MPRAHumanGRCh38
Explore data
E6XI1YFXO

COMT mid-polysome ribosome-load MAVE (F3)

A codon-saturation COMT library covering MB-COMT codons 40–74 and 136–158 was stably integrated into the HEK293T landing-pad cell line. Four biological replicates were profiled by amplicon sequencing, and this table reports the published ALDEx2 variant effects for polysome metafraction 3 (mid polysome) relative to total RNA.

Targeted Genomic Integration MPRAHumanGRCh38
Explore data

Raw source data 25 files

Original supplemental and deposited inputs retained for this study. Download files individually or together as a ZIP; nested folders are preserved. Source reuse terms apply, and sequencing reads may be omitted.

Download all 25 files (ZIP)mavedb/experiment-set.jsonmavedb/score-set-a-1.jsonmavedb/score-set-b-1.jsonmavedb/score-set-c-1.jsonmavedb/score-set-d-1.jsonmavedb/scores-a-1.csvmavedb/scores-b-1.csvmavedb/scores-c-1.csvmavedb/scores-d-1.csvREADME.mdreferences/pDEST_HC_Rec_Bxb_v2_UTR5_Flag_COMT_moxGFP_LPS.fareferences/pDEST_HC_Rec_Bxb_v2_UTR5_Flag_COMT_moxGFP_LPS.gffreferences/pDEST_HC_Rec_Bxb_v2_UTR5_Flag_COMT_moxGFP_LPS_coding.bedreferences/pDEST_HC_Rec_Bxb_v2_UTR5_Flag_COMT_moxGFP_LPS_insert.bedreferences/pDEST_HC_Rec_Bxb_v2_UTR5_Flag_COMT_moxGFP_LPS_ROI_primers_v2.bedreferences/pDEST_HC_Rec_Bxb_v2_UTR5_Flag_COMT_moxGFP_LPS_target_ROI.bedsupplementary/Dataset_EV1.xlssupplementary/Dataset_EV2.xlssupplementary/Dataset_EV3.xlssupplementary/Dataset_EV4.xlssupplementary/Dataset_EV5.xlszenodo/comt_library_yields.txtzenodo/COMT_log10_freq.txtzenodo/COMT_log_freq.txtzenodo/IQLRs.txt

Cite OpenMPRA

Cite the OpenMPRA database. Include your access date because the collection changes over time.

Please also cite the source studies when using their data.