Epigenetic variation impacts individual differences in the transcriptional response to influenza infection.

Humans display remarkable interindividual variation in their immune response to identical challenges. Yet, our understanding of the genetic and epigenetic factors contributing to such variation remains limited. Here we performed in-depth genetic, epigenetic and transcriptional profiling on primary macrophages derived from individuals of European and African ancestry before and after infection with influenza A virus. We show that baseline epigenetic profiles are strongly predictive of the transcriptional response to influenza A virus across individuals. Quantitative trait locus (QTL) mapping revealed highly coordinated genetic effects on gene regulation, with many cis-acting genetic variants impacting concomitantly gene expression and multiple epigenetic marks. These data reveal that ancestry-associated differences in the epigenetic landscape can be genetically controlled, even more than gene expression. Lastly, among QTL variants that colocalized with immune-disease loci, only 7% were gene expression QTL, while the remaining genetic variants impact epigenetic marks, stressing the importance of considering molecular phenotypes beyond gene expression in disease-focused studies.

Adjusting for genetic confounders in transcriptome-wide association studies improves discovery of risk genes of complex traits.

Many methods have been developed to leverage expression quantitative trait loci (eQTL) data to nominate candidate genes from genome-wide association studies. These methods, including colocalization, transcriptome-wide association studies (TWAS) and Mendelian randomization-based methods; however, all suffer from a key problem-when assessing the role of a gene in a trait using its eQTLs, nearby variants and genetic components of other genes' expression may be correlated with these eQTLs and have direct effects on the trait, acting as potential confounders. Our extensive simulations showed that existing methods fail to account for these 'genetic confounders', resulting in severe inflation of false positives. Our new method, causal-TWAS (cTWAS), borrows ideas from statistical fine-mapping and allows us to adjust all genetic confounders. cTWAS showed calibrated false discovery rates in simulations, and its application on several common traits discovered new candidate genes. In conclusion, cTWAS provides a robust statistical framework for gene discovery.

GoM DE: interpreting structure in sequence count data with differential expression analysis allowing for grades of membership.

Parts-based representations, such as non-negative matrix factorization and topic modeling, have been used to identify structure from single-cell sequencing data sets, in particular structure that is not as well captured by clustering or other dimensionality reduction methods. However, interpreting the individual parts remains a challenge. To address this challenge, we extend methods for differential expression analysis by allowing cells to have partial membership to multiple groups. We call this grade of membership differential expression (GoM DE). We illustrate the benefits of GoM DE for annotating topics identified in several single-cell RNA-seq and ATAC-seq data sets.

A new Bayesian factor analysis method improves detection of genes and biological processes affected by perturbations in single-cell CRISPR screening.

Clustered regularly interspaced short palindromic repeats (CRISPR) screening coupled with single-cell RNA sequencing has emerged as a powerful tool to characterize the effects of genetic perturbations on the whole transcriptome at a single-cell level. However, due to its sparsity and complex structure, analysis of single-cell CRISPR screening data is challenging. In particular, standard differential expression analysis methods are often underpowered to detect genes affected by CRISPR perturbations. We developed a statistical method for such data, called guided sparse factor analysis (GSFA). GSFA infers latent factors that represent coregulated genes or gene modules; by borrowing information from these factors, it infers the effects of genetic perturbations on individual genes. We demonstrated through extensive simulation studies that GSFA detects perturbation effects with much higher power than state-of-the-art methods. Using single-cell CRISPR data from human CD8+ T cells and neural progenitor cells, we showed that GSFA identified biologically relevant gene modules and specific genes affected by CRISPR perturbations, many of which were missed by existing methods, providing new insights into the functions of genes involved in T cell activation and neurodevelopment.

Single-cell genomics improves the discovery of risk variants and genes of atrial fibrillation.

Genome-wide association studies (GWAS) have linked hundreds of loci to cardiac diseases. However, in most loci the causal variants and their target genes remain unknown. We developed a combined experimental and analytical approach that integrates single cell epigenomics with GWAS to prioritize risk variants and genes. We profiled accessible chromatin in single cells obtained from human hearts and leveraged the data to study genetics of Atrial Fibrillation (AF), the most common cardiac arrhythmia. Enrichment analysis of AF risk variants using cell-type-resolved open chromatin regions (OCRs) implicated cardiomyocytes as the main mediator of AF risk. We then performed statistical fine-mapping, leveraging the information in OCRs, and identified putative causal variants in 122 AF-associated loci. Taking advantage of the fine-mapping results, our novel statistical procedure for gene discovery prioritized 46 high-confidence risk genes, highlighting transcription factors and signal transduction pathways important for heart development. In summary, our analysis provides a comprehensive map of AF risk variants and genes, and a general framework to integrate single-cell genomics with genetic studies of complex traits.

Asymmetric total synthesis and anti-hepatocellular carcinoma profile of enantiopure euphopilolide and jolkinolide E.

A flexible asymmetric synthesis of both enantiomers of euphopilolide (1) and jolkinolide E (2) [(+)-and (-)-1, (+)-and (-)-2] has been accomplished. This synthesis features an intramolecular oxa-Pauson-Khand reaction (o-PKR) to expeditiously construct the challenging tetracyclic [6.6.6.5] abietane-type diterpene framework, elegantly showcasing the complexity-generating features of o-PKR synthetic methodology leveraging on a judiciously chosen suitable chiral pool scaffold. Furthermore, the anti-hepatocellular carcinoma (HCC) activity of synthetic (-)-euphopilolide (1), (-)-jolkinolide E (2) and their analogues was evaluated. We found that (-)-euphopilolide (1) and (-)-jolkinolide E (2) inhibited the proliferation and induced apoptosis in HCC cells. These findings lay a good foundation for further pharmacology studies of abietane lactone derivatives and provide valuable insight for the development of anti-HCC small molecule drug of natural product origin.

Dimeric oxyberberine CT4-1 targets LINC02331 to induce cytotoxicity and inhibit chemoresistance via suppressing Wnt/ß-catenin signaling in hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) is a type of cancer characterized by high recurrence rates. Overcoming chemoresistance can reduce HCC recurrence and improve patients' prognosis. This work aimed to identify HCC chemoresistance-associated long non-coding RNA (lncRNA) and find an effective drug targeting the identified lncRNA for ameliorating the chemoresistance. In this investigation, bioinformatics analysis based on The Cancer Genome Atlas revealed a new chemoresistance index and suggested LINC02331 as an HCC chemoresistance and patients' prognosis-associated lncRNA that served as an independent prognostic indicator. Moreover, LINC02331 promoted DNA damage repair, DNA replication, and epithelial-mesenchymal transition as well as attenuated cell cycle arrest and apoptosis through regulating Wnt/ß-catenin signaling, thus stimulating HCC resistance to cisplatin cytotoxicity, proliferation, and metastasis. Interestingly, we developed a novel oxidative coupling approach to synthesize a dimeric oxyberberine CT4-1, which exerted superior anti-HCC activities without obvious side effects measured by in vivo mice model and could downregulate LINC02331 mice model and could downregulate LINC02331 to mitigate LINC02331-induced HCC progression by suppressing Wnt/ß-catenin signaling. RNA sequencing analyses verified the involvement of CT4-1-affected differential expression genes in dysregulated pathways and processes, including Wnt, DNA damage repair, cell cycle, DNA replication, apoptosis, and cell adhesion molecules. Furthermore, CT4-1 was demonstrated to be an effective cytotoxic drug in ameliorating HCC patients' prognosis with a prediction model constructed based on RNA-sequencing data from CT4-1-treated cancer cells and public cancer database. In summary, HCC chemoresistance-associated LINC02331 independently predicted poor patients' prognosis and enhanced HCC progression by promoting resistance to cisplatin cytotoxicity, proliferation, and metastasis. Targeting LINC02331 by the dimeric oxyberberine CT4-1 that exhibited synergistic cytotoxicity with cisplatin could alleviate HCC progression and improve patients' prognosis. Our study identified LINC02331 as an alternative target and suggested CT4-1 as an effective cytotoxic drug in HCC treatment.

GoM DE: interpreting structure in sequence count data with differential expression analysis allowing for grades of membership.

Parts-based representations, such as non-negative matrix factorization and topic modeling, have been used to identify structure from single-cell sequencing data sets, in particular structure that is not as well captured by clustering or other dimensionality reduction methods. However, interpreting the individual parts remains a challenge. To address this challenge, we extend methods for differential expression analysis by allowing cells to have partial membership to multiple groups. We call this grade of membership differential expression (GoM DE). We illustrate the benefits of GoM DE for annotating topics identified in several single-cell RNA-seq and ATAC-seq data sets.

Analysis of the water state and microfractal dimension of tilapia fillets during freezing and thawing.

To study the effects of freezing and thawing times and freezing temperatures on the water state and microstructure of tilapia fillets, experiments on tilapia fillets were carried out at -4 and -18°C with one to four freezing and thawing cycles (FTCs). Low-field nuclear magnetic resonance (LF-NMR) and nuclear magnetic resonance imaging were used to observe the water state after different treatments, and scanning electron microscopy (SEM) and frozen sections were used to observe the microstructure changes. Fractal dimension (FD) was used to quantitatively characterize the microstructure of the fish tissue, and the correlation between FD and fish fillet quality parameters was studied by principal component analysis (PCA). The findings showed that with the increase of FTCs, the thawing loss increased, and the water holding capacity (WHC) fell. FTCs cause a decrease in immobilized water and an increase in free water in the fillet. This indicates the migration of immobilized water to free water. SEM and frozen slice images showed that the growth of ice crystals led to the destruction of myogenic fibers. A decrease in freezing temperature inhibited ice crystal growth. The FD value dropped in accordance with an increase in FTCs. PCA demonstrated that the WHC, NMR data, and FD value had a strong correlation with the quality changes in the tilapia fillets. Therefore, FD and water state can reflect the quality characteristics of tilapia fillets. PRACTICAL APPLICATION: The water migration in tilapia fillets is detected with LF-NMR, and the microscopic image may be quantified using the FD value. Both approaches can offer fresh perspectives on how to assess the quality of tilapia fillets and reflect changes in their quality.

Avellis syndrome with ipsilateral prosopalgia, glossopharyngeal neuralgia, and central post-stroke pain: A case report and literature review.

RATIONALE: Avellis syndrome is a rare bulbar syndrome. The main lesions may involve nucleus ambiguus and the lateral spinothalamic tract. The typical reported clinical manifestations are hoarseness, dysphagia, pain, and temperature disturbance of contralateral body. The manifestations, however, may vary. We aim to report new manifestations of Avellis syndrome in this report. PATIENT CONCERNS: A 47-year-old Chinese peasant woman who felt sudden dizziness, nausea when she was doing the laundry was referred to our department from other hospital. She vomited the stomach contents once and complained numbness of the left trunk and limbs as well as coughing while drinking. The patient presented with palatopharyngeal paralysis, Horner syndrome, and diminished pain as well as temperature sensation in the contralateral face, trunk, and limbs. She also presented with ipsilateral prosopalgia, glossopharyngeal neuralgia, and central poststroke pain. DIAGNOSES: T2-weighted MRI images demonstrated a high-signal intensity lesion in the right medulla oblongata which indicated a banded infarction site. The patient was diagnosed with medulla oblongata infarction, Avellis syndrome, Horner syndrome, dysphagia, hemiparesthesia, ipsilateral prosopalgia, glossopharyngeal neuralgia, and central poststroke pain. INTERVENTIONS: The patient was administrated aspirin to prevent the aggregation of platelet and rosuvastatin tablets to regulate lipids as well as to stabilize vascular plaque. She was injected with butylphthalide sodium chloride to improve nerve nutritional status and carbamazepine was prescribed to deal with prosopalgia and glossopharyngeal neuralgia. Gabapentin and pregabalin were administrated to deal with the central poststroke pain. OUTCOMES: The symptoms of prosopalgia as well as glossopharyngeal neuralgia were gone, and dizziness, dysphagia, and Horner syndrome were significantly alleviated when she was discharged from the hospital while the interventions showed little effect on central poststroke pain. LESSONS: We reported a case of Avellis syndrome who manifested as the typical reported manifestations. The patient, what's more, presented with ipsilateral trigeminal, glossopharyngeal neuralgia, and central poststroke pain which were described for the first time. It is of great significance for clinicians to recognize the typical as well as other manifestations which helps to make a clear diagnosis.

Focused Ultrasound Promotes the Delivery of Gastrodin and Enhances the Protective Effect on Dopaminergic Neurons in a Mouse Model of Parkinson's Disease.

Parkinson's disease (PD) is the second most common chronic neurodegenerative disease globally; however, it lacks effective treatment at present. Focused ultrasound (FUS) combined with microbubbles could increase the efficacy of drug delivery to specific brain regions and is becoming a promising technology for the treatment of central nervous system diseases. In this study, we explored the therapeutic potential of FUS-mediated blood-brain barrier (BBB) opening of the left striatum to deliver gastrodin (GAS) in a subacute PD mouse model induced by 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP). The concentration of GAS in the left hemisphere was detected by ultra-high performance liquid chromatography electrospray Q-Orbitrap mass spectrometry (UHPLC/ESI Q-Orbitrap) and the distribution of tyrosine hydroxylase (TH) neurons was detected by immunohistochemical staining. The expression of TH, Dopamine transporter (DAT), cleaved-caspase-3, B-cell lymphoma 2 (Bcl-2), brain-derived neurotrophic factor (BDNF), postsynaptic density protein 95 (PSD-95), and synaptophysin (SYN) protein were detected by western blotting. Analysis showed that the concentration of GAS in the left hemisphere of PD mice increased by approximately 1.8-fold after the BBB was opened. FUS-mediated GAS delivery provided optimal neuroprotective effects and was superior to the GAS or FUS control group. In addition, FUS enhanced GAS delivery significantly increased the expression of Bcl-2, BDNF, PSD-95, and SYN protein in the left striatum (P < 0.05) and reduced the levels of cleaved-caspase-3 remarkably (P = 0.001). In conclusion, the enhanced delivery by FUS effectively strengthened the protective effect of GAS on dopaminergic neurons which may be related to the reinforcement of the anti-apoptotic activity and the expression of synaptic-related proteins in the striatum. Data suggests that FUS-enhanced GAS delivery may represent a new strategy for PD treatment.

Profiling the quantitative occupancy of myriad transcription factors across conditions by modeling chromatin accessibility data.

Over a thousand different transcription factors (TFs) bind with varying occupancy across the human genome. Chromatin immunoprecipitation (ChIP) can assay occupancy genome-wide, but only one TF at a time, limiting our ability to comprehensively observe the TF occupancy landscape, let alone quantify how it changes across conditions. We developed TF occupancy profiler (TOP), a Bayesian hierarchical regression framework, to profile genome-wide quantitative occupancy of numerous TFs using data from a single chromatin accessibility experiment (DNase- or ATAC-seq). TOP is supervised, and its hierarchical structure allows it to predict the occupancy of any sequence-specific TF, even those never assayed with ChIP. We used TOP to profile the quantitative occupancy of hundreds of sequence-specific TFs at sites throughout the genome and examined how their occupancies changed in multiple contexts: in approximately 200 human cell types, through 12 h of exposure to different hormones, and across the genetic backgrounds of 70 individuals. TOP enables cost-effective exploration of quantitative changes in the landscape of TF binding.

Low-Intensity Focused Ultrasound Alleviates Spasticity and Increases Expression of the Neuronal K-Cl Cotransporter in the L4-L5 Sections of Rats Following Spinal Cord Injury.

Low-intensity focused ultrasound (LIFU) has been shown to provide effective activation of the spinal cord neurocircuits. The aim of this study was to investigate the effects of LIFU in order to alleviate spasticity following spinal cord injury (SCI) by activating the spinal neurocircuits and increasing the expression of the neuronal K-Cl cotransporter KCC2. Adult male Sprague Dawley (SD) rats (220-300 g) were randomly divided into a sham control group, a LIFU- group, and a LIFU+ group. The mechanical threshold hold (g) was used to evaluate the behavioral characteristics of spasm. Electromyography (EMG) was used to assess activation of the spinal cord neurocircuits and muscle spontaneous contraction. Spasticity was assessed by frequency-dependent depression (FDD). The expression of KCC2 of the lumbar spinal cord was determined via western blot (WB) and immunofluorescence (IF) staining. The spinal cord neurocircuits were activated by LIFU simulation, which significantly reduced the mechanical threshold (g), FDD, and EMG recordings (s) after 4 weeks of treatment. WB and IF staining both demonstrated that the expression of KCC2 was reduced in the LIFU- group (P < 0.05). After 4 weeks of LIFU stimulation, expression of KCC2 had significantly increased (P < 0.05) in the LIFU+ group compared with the LIFU- group. Thus, we hypothesized that LIFU treatment can alleviate spasticity effectively and upregulate the expression of KCC2 in the L4-L5 section of SCI rats.

Treatment Combining Focused Ultrasound with Gastrodin Alleviates Memory Deficit and Neuropathology in an Alzheimer's Disease-Like Experimental Mouse Model.

Alzheimer's disease (AD) is the most common type of dementia but lacks effective treatment at present. Gastrodin (GAS) is a phenolic glycoside extracted from the traditional Chinese herb-Gastrodia elata-and has been reported as a potential therapeutic agent for AD. However, its efficiency is reduced for AD patients due to its limited BBB permeability. Studies have demonstrated the feasibility of opening the blood-brain barrier (BBB) via focused ultrasound (FUS) to overcome the obstacles preventing medicines from blood flow into the brain tissue. We explored the therapeutic potential of FUS-mediated BBB opening combined with GAS in an AD-like mouse model induced by unilateral intracerebroventricular (ICV) injection of Aß 1-42. Mice were divided into 5 groups: control, untreated, GAS, FUS and FUS+GAS. Combined treatment (FUS+GAS) rather than single intervention (GAS or FUS) alleviated memory deficit and neuropathology of AD-like mice. The time that mice spent in the novel arm was prolonged in the Y-maze test after 15-day intervention, and the waste-cleaning effect was remarkably increased. Contents of Aß, tau, and P-tau in the observed (also the targeted) hippocampus were reduced. BDNF, synaptophysin (SYN), and PSD-95 were upregulated in the combined group. Overall, our results demonstrate that FUS-mediated BBB opening combined with GAS injection exerts the potential to alleviate memory deficit and neuropathology in the AD-like experimental mouse model, which may be a novel strategy for AD treatment.

Impact of the magnetic field-assisted freezing on the moisture content, water migration degree, microstructure, fractal dimension, and the quality of the frozen tilapia.

In this study, we determined the effect of a magnetic field applied during refrigeration in improving the quality of frozen tilapia. Alternating magnetic fields of 10 G, 20 G, 30 G, 40 G, and 50 G were applied during a low-temperature freezing treatment on the back, abdomen, and tail of tilapia. The control group was set at 0 G. A correlation analysis for the fish films after treating with different magnetic field strengths was carried out. The results showed that when the magnetic field was applied to assist freezing, the frozen quality of the tilapia was significantly improved, and the water separation and residual damage were reduced. The felled muscle tissue decreased, the fractal dimension value increased, the hardness decreased, and the elasticity increased. However, the impact of the magnetic field on the quality of the frozen tilapia did not change with an increase in the magnetic field strength. The effect on the back samples was more prominent when the fish were exposed to the magnetic field strength of 40 or 50 G. A magnetic field strength of 50 G was the most effective for the abdominal and tail samples. However, no significant difference was observed in the groups exposed to 10 and 20 G of magnetic fields.

Effects of Noninvasive Low-Intensity Focus Ultrasound Neuromodulation on Spinal Cord Neurocircuits In Vivo.

Although neurocircuits can be activated by focused ultrasound stimulation, it is unclear whether this is also true for spinal cord neurocircuits. In this study, we used low-intensity focused ultrasound (LIFU) to stimulate lumbar 4-lumbar 5 (L4-L5) segments of the spinal cord of normal Sprague Dawley rats with a clapper. The activation of the spinal cord neurocircuits enhanced soleus muscle contraction as measured by electromyography (EMG). Neuronal activation and injury were assessed by EMG, western blotting (WB), immunofluorescence, hematoxylin and eosin (H&E) staining, Nissl staining, enzyme-linked immunosorbent assay (ELISA), immunohistochemistry (IHC), somatosensory evoked potentials (SEPs), motor evoked potentials (MEPs), and the Basso-Beattie-Bresnahan locomotor rating scale. When the LIFU intensity was more than 0.5 MPa, LIFU stimulation induced soleus muscle contraction and increased the EMG amplitudes (P < 0.05) and the number of c-fos- and GAD65-positive cells (P < 0.05). When the LIFU intensity was 3.0 MPa, the LIFU stimulation led to spinal cord damage and decreased SEP amplitudes for electrophysiological assessment (P < 0.05); this resulted in coagulation necrosis, structural destruction, neuronal loss in the dorsal horn by H&E and Nissl staining, and increased expression of GFAP, IL-1ß, TNF-α, and caspase-3 by IHC, ELISA, and WB (P < 0.05). These results show that LIFU can activate spinal cord neurocircuits and that LIFU stimulation with an irradiation intensity ≤1.5 MPa is a safe neurostimulation method for the spinal cord.

Transcranial Ultrasound Stimulation of the Anterior Cingulate Cortex Reduces Neuropathic Pain in Mice.

Focused ultrasound (FUS) is a potential tool for treating chronic pain by modulating the central nervous system. Herein, we aimed to determine whether transcranial FUS stimulation of the anterior cingulate cortex (ACC) effectively improved chronic pain in the chronic compress injury mice model at different stages of neuropathic pain. The mechanical threshold of pain was recorded in the nociceptive tests. We found FUS stimulation elevated the mechanical threshold of pain in both short-term (p < 0.01) and long-term (p < 0.05) experiments. Furthermore, we determined protein expression differences in ACC between the control group, the intervention group, and the Sham group to analyze the underlying mechanism of FUS stimulation in improving neuropathic pain. Additionally, the results showed FUS stimulation led to alterations in differential proteins in long-term experiments, including cellular processes, cellular signaling, and information storage and processing. Our findings indicate FUS may effectively alleviate mechanical neuropathic pain via the ACC's stimulation, especially in the chronic state.

Altered transcriptional and chromatin responses to rhinovirus in bronchial epithelial cells from adults with asthma.

There is a life-long relationship between rhinovirus (RV) infection and the development and clinical manifestations of asthma. In this study we demonstrate that cultured primary bronchial epithelial cells from adults with asthma (n = 9) show different transcriptional and chromatin responses to RV infection compared to those without asthma (n = 9). Both the number and magnitude of transcriptional and chromatin responses to RV were muted in cells from asthma cases compared to controls. Pathway analysis of the transcriptionally responsive genes revealed enrichments of apoptotic pathways in controls but inflammatory pathways in asthma cases. Using promoter capture Hi-C we tethered regions of RV-responsive chromatin to RV-responsive genes and showed enrichment of these regions and genes at asthma GWAS loci. Taken together, our studies indicate a delayed or prolonged inflammatory state in cells from asthma cases and highlight genes that may contribute to genetic risk for asthma.

Allele-specific open chromatin in human iPSC neurons elucidates functional disease variants.

Most neuropsychiatric disease risk variants are in noncoding sequences and lack functional interpretation. Because regulatory sequences often reside in open chromatin, we reasoned that neuropsychiatric disease risk variants may affect chromatin accessibility during neurodevelopment. Using human induced pluripotent stem cell (iPSC)-derived neurons that model developing brains, we identified thousands of genetic variants exhibiting allele-specific open chromatin (ASoC). These neuronal ASoCs were partially driven by altered transcription factor binding, overrepresented in brain gene enhancers and expression quantitative trait loci, and frequently associated with distal genes through chromatin contacts. ASoCs were enriched for genetic variants associated with brain disorders, enabling identification of functional schizophrenia risk variants and their cis-target genes. This study highlights ASoC as a functional mechanism of noncoding neuropsychiatric risk variants, providing a powerful framework for identifying disease causal variants and genes.

Genetic analyses support the contribution of mRNA N6-methyladenosine (m6A) modification to human disease heritability.

N6-methyladenosine (m6A) plays important roles in regulating messenger RNA processing. Despite rapid progress in this field, little is known about the genetic determinants of m6A modification and their role in common diseases. In this study, we mapped the quantitative trait loci (QTLs) of m6A peaks in 60 Yoruba (YRI) lymphoblastoid cell lines. We found that m6A QTLs are largely independent of expression and splicing QTLs and are enriched with binding sites of RNA-binding proteins, RNA structure-changing variants and transcriptional features. Joint analysis of the QTLs of m6A and related molecular traits suggests that the downstream effects of m6A are heterogeneous and context dependent. We identified proteins that mediate m6A effects on translation. Through integration with data from genome-wide association studies, we show that m6A QTLs contribute to the heritability of various immune and blood-related traits at levels comparable to splicing QTLs and roughly half of expression QTLs. By leveraging m6A QTLs in a transcriptome-wide association study framework, we identified putative risk genes of these traits.