Expression Quantitative Trait Loci Analysis Research Articles

Fine mapping spatiotemporal mechanisms of genetic variants underlying cardiac traits and disease

The causal variants and genes underlying thousands of cardiac GWAS signals have yet to be identified. Here, we leverage spatiotemporal information on 966 RNA-seq cardiac samples and perform an expression quantitative trait locus (eQTL) analysis detecting eQTLs considering both eGenes and eIsoforms. We identify 2,578 eQTLs associated with a specific developmental stage-, tissue- and/or cell type. Colocalization between eQTL and GWAS signals of five cardiac traits identified variants with high posterior probabilities for being causal in 210 GWAS loci. Pulse pressure GWAS loci are enriched for colocalization with fetal- and smooth muscle- eQTLs; pulse rate with adult- and cardiac muscle- eQTLs; and atrial fibrillation with cardiac muscle- eQTLs. Fine mapping identifies 79 credible sets with five or fewer SNPs, of which 15 were associated with spatiotemporal eQTLs. Our study shows that many cardiac GWAS variants impact traits and disease in a developmental stage-, tissue- and/or cell type-specific fashion.

Distilling functional variations for human UGT2B4 upstream region based on selection signals and implications for phenotypes of Neanderthal and Denisovan

Our previous work identified one region upstream human UGT2B4 (UDP glucuronosyltransferase family 2 member B4) which is associated with breast cancer and under balancing selection. However, the distribution, functional variation and molecular mechanism underlying breast cancer and balancing selection remain unclear. In current study, the two haplotypes with deep divergence are described by analyzing 1000 genomes project data and observed to be with high frequencies in all human populations. Through population genetics analysis and genome annotation, the potential functional region is identified and verified by reporter gene assay. Further mutagenesis indicates that the functional mutations are rs66862535 and rs68096061. Both SNPs can alter the interaction efficiency of transcription factor POU2F1 (POU class 2 homeobox 1). Through chromosome conformation capture, it is identified that the enhancer containing these two SNPs can interact with UGT2B4 promoter. Expression quantitative trait loci analysis indicates that UGT2B4 expression is dependent on the genotype of this locus. The common haplotype in human is lost in four genomes of archaic hominins, which suggests that Neanderthal and Denisovan should present relatively lower UGT2B4 expression and further higher steroid hormone level. This study provides new insight into the contribution of ancient population structure to human phenotypes.

Brain expression quantitative trait locus and network analyses reveal downstream effects and putative drivers for brain-related diseases

Identification of therapeutic targets from genome-wide association studies (GWAS) requires insights into downstream functional consequences. We harmonized 8,613 RNA-sequencing samples from 14 brain datasets to create the MetaBrain resource and performed cis- and trans-expression quantitative trait locus (eQTL) meta-analyses in multiple brain region- and ancestry-specific datasets (n ≤ 2,759). Many of the 16,169 cortex cis-eQTLs were tissue-dependent when compared with blood cis-eQTLs. We inferred brain cell types for 3,549 cis-eQTLs by interaction analysis. We prioritized 186 cis-eQTLs for 31 brain-related traits using Mendelian randomization and co-localization including 40 cis-eQTLs with an inferred cell type, such as a neuron-specific cis-eQTL (CYP24A1) for multiple sclerosis. We further describe 737 trans-eQTLs for 526 unique variants and 108 unique genes. We used brain-specific gene-co-regulation networks to link GWAS loci and prioritize additional genes for five central nervous system diseases. This study represents a valuable resource for post-GWAS research on central nervous system diseases.

Integrative Functional Genomics Identifies Systemic Lupus Erythematosus Causal Genetic Variant in the IRF5 Risk Locus.

IRF5 plays a crucial role in the development of lupus. Genome-wide association studies have identified several systemic lupus erythematosus (SLE) risk single-nucleotide polymorphisms (SNPs) enriched in the IRF5 locus. However, no comprehensive genome editing-based functional analysis exists to establish a direct link between these variants and altered IRF5 expression, particularly for enhancer variants. This study was undertaken to dissect the regulatory function and mechanisms of SLE IRF5 enhancer risk variants and to explore the utilization of clustered regularly interspaced short palindromic repeat interference (CRISPRi) to regulate the expression of disease risk gene to intervene in the disease. Epigenomic profiles and expression quantitative trait locus analysis were applied to prioritize putative functional variants in the IRF5 locus. CRISPR-mediated deletion, activation, and interference were performed to investigate the genetic function of rs4728142. Allele-specific chromatin immunoprecipitation-quantitative polymerase chain reaction and allele-specific formaldehyde-assisted isolation of regulatory element-quantitative polymerase chain reaction were used to decipher the mechanism of alleles differentially regulating IRF5 expression. The CRISPRi approach was used to evaluate the intervention effect in monocytes from SLE patients. SLE risk SNP rs4728142 was located in an enhancer region, indicating a disease-related regulatory function, and risk allele rs4728142-A was closely associated with increased IRF5 expression. We demonstrated that an rs4728142-containing region could act as an enhancer to regulate the expression of IRF5. Moreover, rs4728142 affected the binding affinity of zinc finger and BTB domain-containing protein 3 (ZBTB3), a transcription factor involved in regulation. Furthermore, in monocytes from SLE patients, CRISPR-based interference with the regulation of this enhancer attenuated the production of disease-associated cytokines. These results demonstrate that the rs4728142-A allele increases the SLE risk by affecting ZBTB3 binding, chromatin status, and regulating IRF5 expression, establishing a biologic link between genetic variation and lupus pathogenesis.

Pleiotropy analysis between lobar intracerebral hemorrhage and CSF β-amyloid highlights new and established associations.

Combining biologically related traits in genome-wide association studies (GWAS) increases the power for genetic discovery. Given the established relationship between lobar intracerebral hemorrhage (ICH) and cerebral amyloid angiopathy (CAA), and between the latter and levels of cerebrospinal fluid amyloid-β 42 (CSF-Aβ42), we leveraged genetic predisposition for lower CSF-Aβ42 levels as a proxy phenotype for CAA to identify new genes associated with lobar ICH. We used publicly available GWAS data for CSF-Aβ42 levels (n = 3146) and for lobar ICH (n = 2094). First, we evaluated the association between lobar ICH risk and CSF-Aβ42 in lobar ICH patients using a polygenic risk score (PRS) for CSF-Aβ42. Next, we conducted multi-trait analysis of GWAS (MTAG) for pleiotropy analysis of lobar ICH and CSF-Aβ42. MTAG results were further tested using Expression Quantitative Trait Locus and Differential Gene Expression Analyses. CSF-Aβ42 PRS was associated with lobar ICH risk (p = 0.04). MTAG analysis identified a novel association within CDH9 (rs1007589; minor allele frequency = 0.09; MTAG p = 5.4 × 10-8; lobar ICH odds ratio = 1.4 and p = 2.4 × 10-3; CSF-Aβ42 β = -0.03 and p = 4.5 × 10-6). rs1007589 was significantly associated with the expression levels of CDH9 in temporal and occipital cortices, regions known to preferentially accumulate microhemorrhages in CAA. Our pleiotropy analysis suggested a variant possibly implicated with lobar ICH driven by amyloid-related mechanisms in CDH9 and associated with differential expression in brain regions characteristically affected by CAA. CDH9 is one subtype of the cadherin superfamily, which regulates intercellular adhesion, is involved in blood-brain barrier integrity, and is elevated in Alzheimer's disease patients. Further analyses are warranted to understand the effects of the variant on the pathogenesis of ICH and its clinical significance.

Genetic associations between bipolar disorder and brain structural phenotypes.

Patients with bipolar disorder (BD) and their first-degree relatives exhibit alterations in brain volume and cortical structure, whereas the underlying genetic mechanisms remain unclear. In this study, based on the published genome-wide association studies (GWAS), the extent of polygenic overlap between BD and 15 brain structural phenotypes was investigated using linkage disequilibrium score regression and MiXeR tool, and the shared genomic loci were discovered by conjunctional false discovery rate (conjFDR) and expression quantitative trait loci (eQTL) analyses. MiXeR estimated the overall measure of polygenic overlap between BD and brain structural phenotypes as 4-53% on a 0-100% scale (as quantified by the Dice coefficient). Subsequent conjFDR analyses identified 54 independent loci (71 risk single-nucleotide polymorphisms) jointly associated with BD and brain structural phenotypes with a conjFDR < 0.05, among which 33 were novel that had not been reported in the previous BD GWAS. Follow-up eQTL analyses in respective brain regions both confirmed well-known risk genes (e.g. CACNA1C, NEK4, GNL3, MAPK3) and discovered novel risk genes (e.g. LIMK2 and CAMK2N2). This study indicates a substantial shared genetic basis between BD and brain structural phenotypes, and provides novel insights into the developmental origin of BD and related biological mechanisms.

Novel Genetic Variants in TP37, PIK3R1, CALM1, and PLCG2 of the Neurotrophin Signaling Pathway Are Associated with the Progression from Mild Cognitive Impairment to Alzheimer's Disease.

Alzheimer's disease (AD) is a common neurodegenerative disease and mild cognitive impairment (MCI) is considered as the prodromal stage of AD. Previous studies showed that changes in the neurotrophin signaling pathway could lead to cognitive decline in AD. However, the association of single nucleotide polymorphisms (SNPs) in genes that are involved in this pathway with AD progression from MCI remains unclear. We investigated the associations between SNPs involved in the neurotrophin signaling pathway with AD progression. We performed single-locus analysis to identify neurotrophin-signaling-related SNPs associated with the AD progression using 767 patients from the Alzheimer's Disease Neuroimaging Initiative study and 1,373 patients from the National Alzheimer's Coordinating Center study. We constructed polygenic risk scores (PRSs) using the identified independent non-APOE SNPs and evaluated its prediction performance on AD progression. We identified 25 SNPs significantly associated with AD progression with Bayesian false-discovery probability ≤0.8. Based on the linkage disequilibrium clumping and expression quantitative trait loci analysis, we found 6 potentially functional SNPs that were associated with AD progression independently. The PRS analysis quantified the combined effects of these SNPs on longitudinal cognitive assessments and biomarkers from cerebrospinal fluid and neuroimaging. The addition of PRSs to the prediction model for 3-year progression to AD from MCI significantly increased the predictive accuracy. Genetic variants in the specific genes of the neurotrophin signaling pathway are predictors of AD progression. eQTL analysis supports that these SNPs regulate expression of key genes involved in the neurotrophin signaling pathway.

Genome-Wide Association Study on Longitudinal Change in Fasting Plasma Glucose in Korean Population.

Genome-wide association studies (GWAS) on type 2 diabetes mellitus (T2DM) have identified more than 400 distinct genetic loci associated with diabetes and nearly 120 loci for fasting plasma glucose (FPG) and fasting insulin level to date. However, genetic risk factors for the longitudinal deterioration of FPG have not been thoroughly evaluated. We aimed to identify genetic variants associated with longitudinal change of FPG over time. We used two prospective cohorts in Korean population, which included a total of 10,528 individuals without T2DM. GWAS of repeated measure of FPG using linear mixed model was performed to investigate the interaction of genetic variants and time, and meta-analysis was conducted. Genome-wide complex trait analysis was used for heritability calculation. In addition, expression quantitative trait loci (eQTL) analysis was performed using the Genotype-Tissue Expression project. A small portion (4%) of the genome-wide single nucleotide polymorphism (SNP) interaction with time explained the total phenotypic variance of longitudinal change in FPG. A total of four known genetic variants of FPG were associated with repeated measure of FPG levels. One SNP (rs11187850) showed a genome-wide significant association for genetic interaction with time. The variant is an eQTL for NOC3 like DNA replication regulator (NOC3L) gene in pancreas and adipose tissue. Furthermore, NOC3L is also differentially expressed in pancreatic β-cells between subjects with or without T2DM. However, this variant was not associated with increased risk of T2DM nor elevated FPG level. We identified rs11187850, which is an eQTL of NOC3L, to be associated with longitudinal change of FPG in Korean population.

Rs77283072 influences breast cancer susceptibility by regulating CDKN2A expression.

Breast cancer is the cancer type with the highest morbidity rates in women, and previous genome-wide association studies (GWASs) have suggested that the single nucleotide polymorphism (SNP) rs1011970 is significantly associated with this disease. An analysis of data from the 1000 Genomes Project demonstrated that there is an SNP, rs77283072, in almost complete linkage disequilibrium with rs1011970, which should therefore present the same signal in a GWAS. However, the actual causal SNP and its associated underlying mechanism have yet to be elucidated. Therefore, the present study evaluated the role of rs77283072 in terms of its association with breast cancer. A dual-luciferase assay was performed, which demonstrated that the two alleles of rs1011970 did not exhibit significantly different reporter gene activity. However, the A allele of rs77283072 exhibited a significant increase in relative luciferase activity compared with the G allele, which suggested that rs77283072 was the causal SNP for breast cancer. Chromosome conformation capture demonstrated that the enhancer containing rs77283072 interacted with the promoter of cyclin-dependent kinase inhibitor 2A (CDKN2A). Furthermore, expression quantitative trait locus analysis demonstrated that the expression of CDKN2A was dependent on the genotype of rs77283072. Taken together, the findings of the present study provided novel insights into the mechanism underlying how the genetic variation in this locus was able to influence breast cancer susceptibility and further the treatment for this disease.

Whole-exome sequence analysis of anthropometric traits illustrates challenges in identifying effects of rare genetic variants.

Anthropometric traits, measuring body size and shape, are highly heritable and significant clinical risk factors for cardiometabolic disorders. These traits have been extensively studied in genome-wide association studies (GWASs), with hundreds of genome-wide significant loci identified. We performed a whole-exome sequence analysis of the genetics of height, body mass index (BMI) and waist/hip ratio (WHR). We meta-analyzed single-variant and gene-based associations of whole-exome sequence variation with height, BMI, and WHR in up to 22,004 individuals, and we assessed replication of our findings in up to 16,418 individuals from 10 independent cohorts from Trans-Omics for Precision Medicine (TOPMed). We identified four trait associations with single-nucleotide variants (SNVs; two for height and two for BMI) and replicated the LECT2 gene association with height. Our expression quantitative trait locus (eQTL) analysis within previously reported GWAS loci implicated CEP63 and RFT1 as potential functional genes for known height loci. We further assessed enrichment of SNVs, which were monogenic or syndromic variants within loci associated with our three traits. This led to the significant enrichment results for height, whereas we observed no Bonferroni-corrected significance for all SNVs. With a sample size of ∼20,000 whole-exome sequences in our discovery dataset, our findings demonstrate the importance of genomic sequencing in genetic association studies, yet they also illustrate the challenges in identifying effects of rare genetic variants.

The circular RNA landscape in multiple sclerosis: Disease-specific associated variants and exon methylation shape circular RNA expression profile.

Circular RNAs (circRNAs) are a class of non-coding RNAs increasingly emerging as crucial actors in the pathogenesis of human diseases, including autoimmune and neurological disorders as multiple sclerosis (MS). Despite several efforts, the mechanisms regulating circRNAs expression are still largely unknown and the circRNA profile and regulation in MS-relevant cell models has not been completely investigated. In this work, we aimed at exploring the global landscape of circRNA expression in MS patients, also evaluating a possible correlation with their genetic and epigenetic background. We performed RNA-seq experiments on circRNA-enriched samples, derived from peripheral blood mononuclear cells (PBMCs) of 10 MS patients and 10 matched controls and performed differential circRNA expression. The genetic background was evaluated using array genotyping, and an expression quantitative trait loci (eQTL) analysis was carried out. Expression analysis revealed 166 differentially expressed circRNAs in MS patients, 125 of which are downregulated. One of the top dysregulated circRNAs, hsa_circ_0007990, derives from the PGAP3 gene, encoding a protein relevant for the control of autoimmune responses. The downregulation of this circRNA was confirmed in two independent replication cohorts, suggesting its implementation as a possible RNA-based biomarker. The eQTL analysis evidenced a significant association between 89 MS-associated loci and the expression of at least one circRNA, suggesting that MS-associated variants could impact on disease pathogenesis by altering circRNA profiles. Finally, we found a significant correlation between exon methylation and circRNA expression levels, supporting the hypothesis that epigenetic features may play an important role in the definition of the cell circRNA pool. We described the circRNA expression profile of PBMCs in MS patients, suggesting that MS-associated variants may tune the expression levels of circRNAs acting as "circ-QTLs", and proposing a role for exon-based DNA methylation in regulating circRNA expression.

Association of RNA m7G Modification Gene Polymorphisms with Pediatric Glioma Risk.

Glioma stemming from glial cells of the central nervous system (CNS) is one of the leading causes of cancer death in childhood. The genetic predisposition of glioma is not fully understood. METTL1-WDR4 methyltransferase complex is implicated in tumorigenesis by catalyzing N7-methylguanosine (m7G) modification of RNA. This study is aimed at determining the association of glioma risk with three polymorphisms (rs2291617, rs10877013, and rs10877012) in METTL1 and five polymorphisms (rs2156315 rs2156316, rs6586250, rs15736, and rs2248490) in WDR4 gene in children of Chinese Han. We enrolled 314 cases and 380 controls from three independent hospitals. Genotypes of these polymorphisms were determined using the TaqMan assay. We found the WDR4 gene rs15736 was significantly associated with reduced glioma risk (GA/AA vs. GG: adjusted odds ratio = 0.63, 95%confidence interval = 0.42 - 0.94, P = 0.023) out of the eight studied polymorphisms. Stratified analyses showed that the association of rs15736 with the risk of glioma remained significant in children aged 60 months or older, girls, the subgroups with astrocytic tumors, or grade I + II glioma. We also found the combined effects of five WDR4 gene polymorphisms on glioma risk. Finally, expression quantitative trait locus (eQTL) analyses elucidated that the rs15736 polymorphism was related to the expression level of WDR4 and neighboring gene cystathionine-beta-synthase (CBS). Our finding provided evidence of a causal association between WDR4 gene polymorphisms and glioma susceptibility in Chinese Han children.

Potentially Functional Genetic Variants in the NRF2 Signaling Pathway Genes are Associated With HBV-related Hepatocellular Carcinoma Survival.

The nuclear factor E2-related factor 2 (NRF2) signaling pathway is one of the most important cell defense pathways. However, it is unclear whether genetic variants in NRF2 signaling pathway genes are associated with the survival of hepatitis B virus (HBV)-related hepatocellular carcinoma (HCC). In the present study, we utilized a new hypothesis-driven approach based on biological pathways to investigate the associations between 17919 single nucleotide polymorphisms (SNPs) in 137 NRF2 signaling pathway genes and the overall survival (OS) of 866 patients with HBV-related HCC. As a result, two independent SNPs with potential biological function were identified to be significantly associated with HBV-related HCC OS: [SLC2A9 rs28643326 T>C: hazard ratio (HR) = 0.74, 95% confidence interval (95% CI) = 0.62-0.89, P < 0.001 and SLC5A10 rs2472711 G>T: HR = 0.81, 95% CI = 0.71-0.93, P = 0.003, respectively]. The expression quantitative trait loci (eQTL) analysis further revealed that the rs28643326 C allele was significantly associated with increased levels of SLC2A9 mRNA expression (P < 0.001), and higher mRNA expression levels of SLC2A9 in adjacent normal liver tissues were associated with better survival. Although the association between the rs2472711 T allele and the mRNA expression of SLC5A10 was not statistically significant (P = 0.200), the fact that rs2472711 is located at the DNase I hypersensitivity site and is a marker for promoter and enhancer histones also suggests that it may have the function of regulating its corresponding gene expression. In conclusion, genetic variants of NRF2 signaling pathway genes may serve as potential prognostic biomarkers for HBV-related HCC and also provide a solid basis for further mechanistic exploration.

Integration Analysis of Epigenetic-related m6A-SNPs Associated with Atrial Fibrillation

Objective: Numerous single nucleotide polymorphisms (SNPs) have been identified as genetic contributors to atrial fibrillation (AF). The aim of this study was to investigate the effects of genome-wide N6-methyladenosine (m6A)-SNPs on AF. Method: m6A-SNPs were identified by analysis of raw data from published AF GWAS datasets and the list of m6A-SNPs from the m6AVar database. Expression quantitative trait loci (eQTL) analysis was conducted to evaluate the effects of m6A-SNPs on gene expression. The expression of linked genes was validated in three independent AF-associated gene expression datasets (GSE14975, GSE108660 and GSE2240). Results: A total of 1429 (6.2%) unique m6A-SNPs that were significantly associated with AF were identified. Seventeen m6A-SNPs in 14 genes reached genome-wide significance. Eight m6A-SNPs demonstrated eQTL signals. Four m6A-SNPs (rs383692, rs3211105, rs1061259 and rs1152582) exhibited strong cis-eQTL signals associated with the gene expression levels of SMIM8, JMJD1C and SYNE2. SYNE2 and SMIM8 had differential gene expression levels between the AF and sinus rhythm groups. In addition, SYNE2 expression was uniformly downregulated in AF samples compared with normal control samples in the three datasets. Conclusions: Our results provide the first demonstration that m6A-SNPs are strongly associated with AF, and extend understanding of m6A modification as a potential biological pathway underlying AF.

S-22-6: RENAL EXPRESSION OF MICRORNAS AND THEIR ASSOCIATION WITH BLOOD PRESSURE REGULATION - INSIGHTS FROM SMALL RNA-SEQUENCING OF 493 HUMAN KIDNEYS

The kidney is a critical regulator of blood pressure (BP) and a contributor to cardiovascular health. Recently, microRNAs (miRNAs) have emerged as potential mediators of genetic predisposition to cardiovascular and kidney disorders, yet their overall role in blood pressure control has not been fully characterised. We sought to identify and catalogue all renal miRNAs, investigate which of them are genetically regulated and assess whether their expression-modifying single nucleotide polymorphisms (SNPs) overlap with genetic variants associated with BP in previous genome-wide association studies (GWAS). We first characterised the entire renal miRNA transcriptome through de novo small RNA-sequencing of 493 human kidney samples. We then integrated these expression profiles with genotype information to perform expression quantitative trait loci (eQTL) analysis - the largest of its kind - and identify genetically regulated miRNAs and their regulatory SNPs. We then examined the extent to which miRNA-regulating variants overlapped with the curated list of BP-associated GWAS SNPs. We identified 1,459 mature miRNAs expressed in the human kidney, and the top 50 most highly expressed miRNAs accounted for 24.2% of the measured expression. We also discovered that 78% of miRNA genes expressed in the kidney were intragenic to other genes and 58% mapped to an intron of those genes. After correction for multiple testing, our eQTL analysis had identified 140 genetically regulated miRNAs in the kidney. They were controlled by 19,901 unique SNPs (eSNPs) with a total of 21,805 miRNA-eSNP pairs. Functional annotation of miRNA eSNPs demonstrated 2.54-fold enrichment for kidney enhancer regions and 1.69-fold enrichment for kidney promoters when compared to randomly selected SNPs. Through overlap with GWAS SNPs for systolic, diastolic and pulse pressures, we identified 8 kidney BP miRNA eSNPs responsible for the regulation of 9 mature miRNAs. Collectively, we generated a catalogue of miRNAs expressed in the human kidney, annotated novel kidney relevant miRNAs and their regulatory variants and discovered 9 renal microRNAs that may act as mediators of genetic predisposition to high BP.

THE TREND TOWARDS PHENOME-WIDE ASSOCIATION STUDIES (PheWASs) IN COVID-19 RESEARCH

Background: Coronavirus Disease-2019 (COVID-19) appears in individuals asymptomatically and in various symptomatic forms. Symptomatic diversity can result in diagnosis failures, hospitalization, admission to intensive care, multi-organ failure, and death. The causes and risk factors of the severity of disease symptoms are uncertain. This uncertainty can only be resolved by elucidating the effects of host genes and genetic variations on different phenotypes. Aim: This review aimed to emphasize the importance of large-scale genotype-phenotype correlation studies in elucidating the phenotypic diversity in COVID-19 disease. Methods: All publications related to Phenome-Wide Association Study (PheWAS) in the PubMed database were searched. PheWAS studies applied to COVID-19 patients have been identified. In addition, studies applied to the genome-wide association study (GWAS)- Electronic health records (EHRs) data and additionally matched to the gene expression data were systematically reviewed. The latest PheWAS methodology and its importance in Large-scale genotype-phenotype correlations are discussed within the context of published COVID-19 studies. Results: According to our PubMed search data, there are few PheWAS studies on COVID-19 disease. This review explains the use of PheWAS studies applied to health records and GWAS data, and colocalization studies applied to expression quantitative trait locus (eQTL) analysis to understand the phenotypic variability of COVID-19. Discussion; Although there is a very limited number of PheWAS studies on COVID-19 diseases, these studies have obtained important data. At the current stage, there is a need for such studies in COVID-19 research. Conclusions: PheWAS is an ideal method for large-scale genotype-phenotype correlation studies that can reveal genetic diversity and phenotypic diversity in the pathophysiology of the disease.

Applying multi-omics data to study the genetic background of bovine respiratory disease infection in feedlot crossbred cattle.

Bovine respiratory disease (BRD) is the most common and costly infectious disease affecting the wellbeing and productivity of beef cattle in North America. BRD is a complex disease whose development is dependent on environmental factors and host genetics. Due to the polymicrobial nature of BRD, our understanding of the genetic and molecular mechanisms underlying the disease is still limited. This knowledge would augment the development of better genetic/genomic selection strategies and more accurate diagnostic tools to reduce BRD prevalence. Therefore, this study aimed to utilize multi-omics data (genomics, transcriptomics, and metabolomics) analyses to study the genetic and molecular mechanisms of BRD infection. Blood samples of 143 cattle (80 BRD; 63 non-BRD animals) were collected for genotyping, RNA sequencing, and metabolite profiling. Firstly, a genome-wide association study (GWAS) was performed for BRD susceptibility using 207,038 SNPs. Two SNPs (Chr5:25858264 and BovineHD1800016801) were identified as associated (p-value <1 × 10-5) with BRD susceptibility. Secondly, differential gene expression between BRD and non-BRD animals was studied. At the significance threshold used (log2FC>2, logCPM>2, and FDR<0.01), 101 differentially expressed (DE) genes were identified. These DE genes significantly (p-value <0.05) enriched several immune responses related functions such as inflammatory response. Additionally, we performed expression quantitative trait loci (eQTL) analysis and identified 420 cis-eQTLs and 144 trans-eQTLs significantly (FDR <0.05) associated with the expression of DE genes. Interestingly, eQTL results indicated the most significant SNP (Chr5:25858264) identified via GWAS was a cis-eQTL for DE gene GPR84. This analysis also demonstrated that an important SNP (rs209419196) located in the promoter region of the DE gene BPI significantly influenced the expression of this gene. Finally, the abundance of 31 metabolites was significantly (FDR <0.05) different between BRD and non-BRD animals, and 17 of them showed correlations with multiple DE genes, which shed light on the interactions between immune response and metabolism. This study identified associations between genome, transcriptome, metabolome, and BRD phenotype of feedlot crossbred cattle. The findings may be useful for the development of genomic selection strategies for BRD susceptibility, and for the development of new diagnostic and therapeutic tools.

A novel BCL11A polymorphism influences gene expression, therapeutic response and epilepsy risk: A multicenter study.

Genetic factors have been found to be associated with the efficacy and adverse reactions of antiseizure medications. BCL11A is an important regulator of the development of neuronal networks. However, the role of BCL11A in epilepsy remains unclear. This study aimed to evaluate the genetic association of BCL11A with the susceptibility to develop epileptic seizures and therapeutic response of patients with epilepsy in Han Chinese. We matched 450 epilepsy cases with 550 healthy controls and 131 drug-resistant epilepsy patients with 319 drug-responsive epilepsy patients from two different centers. Genetic association analysis, genetic interaction analysis, expression quantitative trait loci analysis and protein-protein interaction analysis were conducted. Our results showed that rs2556375 not only increases susceptibility to develop epileptic seizures (OR = 2.700, 95% = 1.366-5.338, p = 0.004 and OR = 2.984, 95% = 1.401-6.356, p = 0.005, respectively), but also increases the risk of drug resistance(OR = 21.336, 95%CI =2.489-183.402, p = 0.005). The interaction between rs2556375 and rs12477097 results in increased risk for pharma coresistant. In addition, rs2556375 regulated BCL11A expression in human brain tissues (p = 0.0096 and p = 0.033, respectively). Furthermore, the protein encoded by BCL11A interacted with targets of approved antiepileptic drugs. BCL11A may be a potential therapeutic target for epilepsy. Rs2556375 may increase the risks of epilepsy and drug resistance by regulating BCL11A expression in human brain tissues. Moreover, the interaction between rs2556375 and rs12477097 results in increased risk for drug resistance.

BOT-3 GENOME-WIDE ASSOCIATION STUDY OF INTRACRANIAL GERM CELL TUMORS: A COMMON DELETION ATBAK1 ATTENUATES THE ENHANCER ACTIVITY AND CONFERS RISK FOR THE BRAIN TUMORS IN CHILDREN ADOLESCENTS AND YOUNG ADULTS

Abstract Intracranial germ cell tumors (IGCTs) are rare brain tumors that mainly occur in children, adolescents and young adults with a particularly high incidence in East Asian populations. The biological basis of these tumors is still largely unknown. We conduct a genome-wide association study (GWAS) of 133 patients with IGCTs and 762 controls of Japanese ancestry by using Infinium Asian Screening Array and other molecular biology methods. A common 4-bp deletion polymorphism in an enhancer adjacent to BAK1 is significantly associated with the disease risk (rs3831846; P = 2.4 × 10−9, odds ratio = 2.46 [95% CI: 1.83-3.31], minor allele frequency = 0.43). Rs3831846 is in strong linkage disequilibrium with a testicular GCTs susceptibility variant rs210138. Expression quantitative trait locus (eQTL) analysis using the GTEx dataset reveals that the risk allele of rs3831846 has a down-regulating effect on BAK1 expression in a wide range of tissues. Further in-vitro reporter assays reveal rs3831846 to be a functional variant attenuating the enhancer activity. BAK1 is a pro-apoptotic member of the BCL-2 family, thus our results suggested that the risk allele may contribute to IGCTs predisposition through down-regulating BAK1 expression. Risk alleles of testicular GCTs derived from the European GWAS show significant positive correlations in the effect sizes with the Japanese IGCTs GWAS (P = 1.3 × 10-4, Spearman’s ρ = 0.48). Of the 57 loci, 11 exhibit significant association with IGCTs and these loci were implicated in a broad range of biological pathways, including KIT/KITLG signaling, apoptosis regulation, and telomerase activity. The risk allele frequency of rs3831846 is higher in east Asia than Europe (0.49 vs. 0.20), which may provide a partial explanation for the ethnic difference in incidence. Nevertheless, our results suggest the shared genetic susceptibility of GCTs beyond ethnicity and primary sites.

Identification of Genetic Loci and Sex Differences Associated with Amyloid Burden Using the AD‐BXD Mouse Population

AbstractBackgroundIn light of the rising prevalence of both diet‐induced obesity and Alzheimer’s Disease (AD) in the developed world, the influence of a high‐fat diet (HFD) on AD progression was investigated. Beta‐amyloid 1‐42 (Aβ42), in particular, was examined as little is known about the extent to which an individual’s genetics, sex, and diet determine the onset and progression of this classic AD neuropathology.MethodAD‐BXD mice were generated by crossing female 5XFAD mice to male members of the BXD reference panel (Neuner, et al. 2019). Cortical Aβ42 levels were quantified (ELISA) from four hundred and twenty‐two mice (235 female, 187 male) across 52 AD‐BXD strains fed either a 6kcal% extruded control diet (n = 227) or 45kcal% HFD (n = 195) at 10 weeks of age until 3, 6, or 14 months old.ResultWhile a 45kcal% HFD did not alter cortical Aβ42 levels in 5XFAD mice relative to control diet, there was a main effect of genetic background (strain) and sex. Specifically, wide variation in cortical Aβ42 was observed across the AD‐BXD mouse population that was highly heritable (h2≈0.9), and females displayed greater Aβ42 levels than males. Quantitative trait locus (QTL) mapping revealed significant genetic loci associated with increased amyloid burden with aging, indicating sequence differences on chromosomes (Chr) 19 and X that may underlie the observed sex differences when age was included as covariates. Furthermore, we identified several missense mutations within our Chr 19 QTL that could explain differences in Aβ42 phenotype, and expression QTL analysis is underway to identify putative regulatory variants of cortical amyloid burden.ConclusionWe found that genetic factors, strain and sex, are key determinants of cortical amyloid burden – but not diet (environmental factor) suggesting involvement of other regulators/pathways in eliciting cognitive decline in AD‐BXD mouse population (Dunn et. al., AAIC 2021) Resolving genetic modifiers of amyloid burden located on Chr 19 and X may reveal novel targets for anti‐amyloid therapeutics and provide a possible mechanistic explanation for why females (XX) are at a higher risk for developing AD than males (XY).