Gene Expression Quantitative Trait Loci Research Articles

#2347 Multi-omics analysis identifies traditional Chinese medicine targets for chronic kidney disease

Abstract Background and Aims As a major public health problem, chronic kidney disease (CKD) has substantial comorbidities and disease burden. It is necessary to discover new therapeutic targets considering its irreversibly. Traditional Chinese medicine (TCM) has become indispensable for the integrated treatment of CKD. The current study aims to explore potential new TCM targets in high-throughput experiment- and reference-guided database (HERB) that are related to renal phenotypes utilizing multi-omics. Method We conducted comprehensive Mendelian Randomization and Bayesian colocalization on more than 2,000 HERB target integrating genomes and transcriptomes, using summary statistics of genome-wide association study and cis gene expression quantitative trait loci. The results were validated through differential gene expression in public database, previous literature reports and clinical kidney biopsy specimens. Results In all, we discovered 34 targets associated with renal phenotypes among the 2,421 HERB TCM targets in the data mining part. Among them, most of the targets can be validated by differential gene expression analysis in GEO and Nephroseq, and have been reported by basic research, thus confirming the reliability of this method. In addition, PRKCI, KNG1, PKN3, GGT7, GBA, APEH and GSTA2 are expected to be potential new therapeutic targets for CKD. Immunohistochemical staining using clinical kidney biopsy specimens of PKN3 and KNG1, two targets with the highest evidence level, further confirm the results. Conclusion Integration of multi-omics can extensively and efficiently screen potential drug targets related to CKD. The new targets discovered in the current study, especially PKN3 and KNG1, need to be further validated in basic research.

Determinants of gastric cancer immune escape identified from non-coding immune-landscape quantitative trait loci

The landscape of non-coding mutations in cancer progression and immune evasion is largely unexplored. Here, we identify transcrptome-wide somatic and germline 3′ untranslated region (3′-UTR) variants from 375 gastric cancer patients from The Cancer Genome Atlas. By performing gene expression quantitative trait loci (eQTL) and immune landscape QTL (ilQTL) analysis, we discover 3′-UTR variants with cis effects on expression and immune landscape phenotypes, such as immune cell infiltration and T cell receptor diversity. Using a massively parallel reporter assay, we distinguish between causal and correlative effects of 3′-UTR eQTLs in immune-related genes. Our approach identifies numerous 3′-UTR eQTLs and ilQTLs, providing a unique resource for the identification of immunotherapeutic targets and biomarkers. A prioritized ilQTL variant signature predicts response to immunotherapy better than standard-of-care PD-L1 expression in independent patient cohorts, showcasing the untapped potential of non-coding mutations in cancer.

Identification and characterization of whole blood gene expression and splicing quantitative trait loci during early to mid-lactation of dairy cattle

BackgroundCharacterization of regulatory variants (e.g., gene expression quantitative trait loci, eQTL; gene splicing QTL, sQTL) is crucial for biologically interpreting molecular mechanisms underlying loci associated with complex traits. However, regulatory variants in dairy cattle, particularly in specific biological contexts (e.g., distinct lactation stages), remain largely unknown. In this study, we explored regulatory variants in whole blood samples collected during early to mid-lactation (22–150 days after calving) of 101 Holstein cows and analyzed them to decipher the regulatory mechanisms underlying complex traits in dairy cattle.ResultsWe identified 14,303 genes and 227,705 intron clusters expressed in the white blood cells of 101 cattle. The average heritability of gene expression and intron excision ratio explained by cis-SNPs is 0.28 ± 0.13 and 0.25 ± 0.13, respectively. We identified 23,485 SNP-gene expression pairs and 18,166 SNP-intron cluster pairs in dairy cattle during early to mid-lactation. Compared with the 2,380,457 cis-eQTLs reported to be present in blood in the Cattle Genotype-Tissue Expression atlas (CattleGTEx), only 6,114 cis-eQTLs (P < 0.05) were detected in the present study. By conducting colocalization analysis between cis-e/sQTL and the results of genome-wide association studies (GWAS) from four traits, we identified a cis-e/sQTL (rs109421300) of the DGAT1 gene that might be a key marker in early to mid-lactation for milk yield, fat yield, protein yield, and somatic cell score (PP4 > 0.6). Finally, transcriptome-wide association studies (TWAS) revealed certain genes (e.g., FAM83H and TBC1D17) whose expression in white blood cells was significantly (P < 0.05) associated with complex traits.ConclusionsThis study investigated the genetic regulation of gene expression and alternative splicing in dairy cows during early to mid-lactation and provided new insights into the regulatory mechanisms underlying complex traits of economic importance.

Integrated analysis of single cell-RNA sequencing and Mendelian randomization identifies lactate dehydrogenase B as a target of melatonin in ischemic stroke.

Despite the success of single-cell RNA sequencing in identifying cellular heterogeneity in ischemic stroke, clarifying the mechanisms underlying these associations of differently expressed genes remains challenging. Several studies that integrate gene expression and gene expression quantitative trait loci (eQTLs) with genome wide-association study (GWAS) data to determine their causal role have been proposed. Here, we combined Mendelian randomization (MR) framework and single cell (sc) RNA sequencing to study how differently expressed genes (DEGs) mediating the effect of gene expression on ischemic stroke. The hub gene was further validated in the invitro model. We identified 2339 DEGs in 10 cell clusters. Among these DEGs, 58 genes were associated with the risk of ischemic stroke. After external validation with eQTL dataset, lactate dehydrogenase B (LDHB) is identified to be positively associated with ischemic stroke. The expression of LDHB has also been validated in sc RNA-seq with dominant expression in microglia and astrocytes, and melatonin is able to reduce the LDHB expression and activity invitro ischemic models. Our study identifies LDHB as a novel biomarker for ischemic stroke via combining the sc RNA-seq and MR analysis.

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.

Circadian clock-related genome-wide mendelian randomization identifies putatively genes for ulcerative colitis and its comorbidity

BackgroundCircadian rhythm is crucial to the function of the immune system. Disorders of the circadian rhythm can contribute to inflammatory diseases such as Ulcerative colitis (UC). This Mendelian Randomization (MR) analysis applies genetic tools to represent the aggregated statistical results of exposure to circadian rhythm disorders and UC and its comorbidities, allowing for causal inferences.MethodsSummary statistics of protein, DNA methylation and gene expression quantitative trait loci in individuals of European ancestry (pQTL, mQTL, and eQTL, respectively) were used. Genetic variants located within or near 152 circadian clock-related genes and closely related to circadian rhythm disorders were selected as instrumental variables. Causal relationships with UC and its comorbidities were then estimated through employed Summary data-based Mendelian Randomization (SMR) and Inverse-Variance-Weighted MR (IVW-MR).ResultsThrough preliminary SMR analysis, we identified a potential causal relationship between circadian clock-related genes and UC along with its comorbidities, which was further confirmed by IVW-MR analysis. Our study identified strong evidence of positive correlation involving seven overlapping genes (CSNK1E, OPRL1, PIWIL2, RORC, MAX, PPP5C, and AANAT) through MWAS and TWAS in UC, four overlapping genes (OPRL1, CHRNB2, FBXL17, and SIRT1) in UC with PSC, and three overlapping genes (ARNTL, USP7, and KRAS) in UC with arthropathy.ConclusionsThis SMR study demonstrates the causal effect of circadian rhythm disorders in UC and its comorbidities. Furthermore, our investigation pinpointed candidate genes that could potentially serve as drug targets.

Liver regulatory mechanisms of noncoding variants at lipid and metabolic trait loci

Genome-wide association studies (GWASs) have identified hundreds of risk loci for liver disease and lipid-related metabolic traits, although identifying their target genes and molecular mechanisms remains challenging. We predicted target genes at GWAS signals by integrating them with molecular quantitative trait loci for liver gene expression (eQTL) and liver chromatin accessibility QTL (caQTL). We predicted specific regulatory caQTL variants at four GWAS signals located near EFHD1, LITAF, ZNF329, and GPR180. Using transcriptional reporter assays, we determined that caQTL variants rs13395911, rs11644920, rs34003091, and rs9556404 exhibit allelic differences in regulatory activity. We also performed a protein binding assay for rs13395911 and found that FOXA2 differentially interacts with the alleles of rs13395911. For variants rs13395911 and rs11644920 in putative enhancer regulatory elements, we used CRISPRi to demonstrate that repression of the enhancers altered the expression of the predicted target and/or nearby genes. Repression of the element at rs13395911 reduced the expression of EFHD1, and repression of the element at rs11644920 reduced the expression of LITAF, SNN, and TXNDC11. Finally, we showed that EFHD1 is a metabolically active gene in HepG2 cells. Together, these results provide key steps to connect genetic variants with cellular mechanisms and help elucidate the causes of liver disease.

Osteoporosis and osteoarthritis: a bi-directional Mendelian randomization study

ObjectiveTo investigate the causal relationship between low bone mineral density (BMD) and osteoarthritis (OA) using Mendelian randomization (MR) design.MethodsTwo-sample bi-directional MR analyses were performed using summary-level information on OA traits from UK Biobank and arcOGEN. Sensitivity analyses including MR-Egger, simple median, weighted median, MR pleiotropy residual sum, and outlier approaches were utilized in conjunction with inverse variance weighting (IVW). Gene ontology (GO) enrichment analyses and expression quantitative trait locus (eQTL) colocalization analyses were used to investigate the potential mechanism and shared genes between osteoporosis (OP) and OA.ResultsThe IVW method revealed that genetically predicted low femoral neck BMD was significantly linked with hip (β = 0.105, 95% CI: 0.023–0.188) and knee OA (β = 0.117, 95% CI: 0.049–0.184), but not with other site-specific OA. Genetically predicted low lumber spine BMD was significantly associated with OA at any sites (β = 0.048, 95% CI: 0.011–0.085), knee OA (β = 0.101, 95% CI: 0.045–0.156), and hip OA (β = 0.150, 95% CI: 0.077–0.224). Only hip OA was significantly linked with genetically predicted reduced total bone BMD (β = 0.092, 95% CI: 0.010–0.174). In the reverse MR analyses, no evidence for a causal effect of OA on BMD was found. GO enrichment analysis and eQTL analysis illustrated that DDN and SMAD-3 were the most prominent co-located genes.ConclusionsThese findings suggested that OP may be causally linked to an increased risk of OA, indicating that measures to raise BMD may be effective in preventing OA. More research is required to determine the underlying processes via which OP causes OA.

Identifying actionable druggable targets for breast cancer: Mendelian randomization and population-based analyses

Drug repurposing provides a cost-effective approach to address the need for breast cancer prevention and therapeutics. We aimed to identify actionable druggable targets using Mendelian randomization (MR) and then validate the candidate drugs using population-based analyses. We identified genetic instruments for 1406 actionable targets of approved non-oncological drugs based on gene expression, DNA methylation, and protein expression quantitative trait loci (eQTL, mQTL, and pQTL, respectively). Genome-wide association study (GWAS) summary statistics were obtained from the Breast Cancer Association Consortium (122,977 cases, 105,974 controls). We further conducted a nested case-control study using data retrieved from Swedish registers to validate the candidate drugs that were identified from MR analyses. We identified six significant MR associations with gene expression levels (TUBB, MDM2, CSK, ULK3, MC1R and KCNN4) and two significant associations with gene methylation levels across 21 CpG islands (RPS23 and MAPT). Results from the nested case-control study showed that the use of raloxifene (targeting MAPT) was associated with 35% reduced breast cancer risk (odds ratio, OR, 0.65; 95% confidence interval, CI, 0.51-0.83). However, usage of estradiol, tolterodine, and nitrofurantoin (also targeting MAPT) was associated with increased breast cancer risk, with adjusted ORs and 95% CI of 1.10 (1.07-1.13), 1.16 (1.09-1.24), and 1.09 (1.05-1.13), respectively. The effect of raloxifene and nitrofurantoin lost significance in further validation analyses using active-comparator and new-user design. This large-scale MR analysis, combined with population-based validation, identified eight druggable target genes for breast cancer and suggested that raloxifene is an effective chemoprevention against breast cancer. Swedish Research Council, Cancerfonden, Crafoordska Stiftelsen, Allmänna Sjukhusets i Malmö Stiftelsen för bekämpande av cancer, 111 Project and MAS cancer.

Causality between depression and ankylosing spondylitis in a European population: Results from a Mendelian randomization analysis.

The aim of this study was to explore the application of Mendelian randomization (MR) Egger and inverse variance weighted (IVW) in a causal effect on depression and ankylosing spondylitis (AS). Instrumental variables (IVs) were determined using genome-wide association studies. The 2-sample MR analysis was conducted by MR Egger to test the causal effect between depression and AS. The pleiotropy of potential instrumental variables was evaluated. The results of MR Egger and IVW were further compared. A total of 3 single nucleotide polymorphisms as the construct IVs were included. IVW results showed a significant causal effect between depression and AS (P < .001). Depression could promote the risk of AS (odds ratio = 1.060, 95% confidence interval: 1.026-1.094). However, the MR Egger showed no causal effect (P = .311). Heterogeneity statistics suggested that no heterogeneity was existed (P > .05). It was also suggested that there was no horizontal pleiotropy in IVs (MR Egger intercept: -0.0004, P = .471). Reverse MR analysis suggested that there was no causal effect between AS and depression (P > .05). Gene expression quantitative trait locus (QTLs) suggested that rs2517601 and RNF39 were positively correlated (beta = 1.066, P < .001). Depression may be one of the causes of AS by MR analysis in a European population. We can estimate the causal effect based on IVW when horizontal pleiotropy is very tiny.

EQTL Catalogue 2023: New datasets, X chromosome QTLs, and improved detection and visualisation of transcript-level QTLs.

The eQTL Catalogue is an open database of uniformly processed human molecular quantitative trait loci (QTLs). We are continuously updating the resource to further increase its utility for interpreting genetic associations with complex traits. Over the past two years, we have increased the number of uniformly processed studies from 21 to 31 and added X chromosome QTLs for 19 compatible studies. We have also implemented Leafcutter to directly identify splice-junction usage QTLs in all RNA sequencing datasets. Finally, to improve the interpretability of transcript-level QTLs, we have developed static QTL coverage plots that visualise the association between the genotype and average RNA sequencing read coverage in the region for all 1.7 million fine mapped associations. To illustrate the utility of these updates to the eQTL Catalogue, we performed colocalisation analysis between vitamin D levels in the UK Biobank and all molecular QTLs in the eQTL Catalogue. Although most GWAS loci colocalised both with eQTLs and transcript-level QTLs, we found that visual inspection could sometimes be used to distinguish primary splicing QTLs from those that appear to be secondary consequences of large-effect gene expression QTLs. While these visually confirmed primary splicing QTLs explain just 6/53 of the colocalising signals, they are significantly less pleiotropic than eQTLs and identify a prioritised causal gene in 4/6 cases.

Overexpression of REDUCED WALL ACETYLATION C increases xylan acetylation and biomass recalcitrance in Populus.

Plant lignocellulosic biomass, i.e. secondary cell walls of plants, is a vital alternative source for bioenergy. However, the acetylation of xylan in secondary cell walls impedes the conversion of biomass to biofuels. Previous studies have shown that REDUCED WALL ACETYLATION (RWA) proteins are directly involved in the acetylation of xylan but the regulatory mechanism of RWAs is not fully understood. In this study, we demonstrate that overexpression of a Populus trichocarpa PtRWA-C gene increases the level of xylan acetylation and increases the lignin content and S/G ratio, ultimately yielding poplar woody biomass with reduced saccharification efficiency. Furthermore, through gene coexpression network and expression quantitative trait loci (eQTL) analysis, we found that PtRWA-C was regulated not only by the secondary cell wall hierarchical regulatory network but also by an AP2 family transcription factor HARDY (HRD). Specifically, HRD activates PtRWA-C expression by directly binding to the PtRWA-C promoter, which is also the cis-eQTL for PtRWA-C. Taken together, our findings provide insights into the functional roles of PtRWA-C in xylan acetylation and consequently saccharification and shed light on synthetic biology approaches to manipulate this gene and alter cell wall properties. These findings have substantial implications for genetic engineering of woody species, which could be used as a sustainable source of biofuels, valuable biochemicals, and biomaterials.

A computational method for cell type-specific expression quantitative trait loci mapping using bulk RNA-seq data

Mapping cell type-specific gene expression quantitative trait loci (ct-eQTLs) is a powerful way to investigate the genetic basis of complex traits. A popular method for ct-eQTL mapping is to assess the interaction between the genotype of a genetic locus and the abundance of a specific cell type using a linear model. However, this approach requires transforming RNA-seq count data, which distorts the relation between gene expression and cell type proportions and results in reduced power and/or inflated type I error. To address this issue, we have developed a statistical method called CSeQTL that allows for ct-eQTL mapping using bulk RNA-seq count data while taking advantage of allele-specific expression. We validated the results of CSeQTL through simulations and real data analysis, comparing CSeQTL results to those obtained from purified bulk RNA-seq data or single cell RNA-seq data. Using our ct-eQTL findings, we were able to identify cell types relevant to 21 categories of human traits.

Mycobacterium tuberculosis Gene Expression Associated With Fluoroquinolone Resistance and Efflux Pump Inhibition.

We evaluated the relationship between response to efflux pump inhibition in fluoroquinolone-resistant Mycobacterium tuberculosis (Mtb) isolates and differences in gene expression and expression quantitative trait loci (eQTL). We determined ofloxacin minimum inhibitory concentration (MIC) for ofloxacin-resistant and -susceptible Mtb isolates without and with the efflux pump inhibitor verapamil. We performed RNA sequencing (RNA-seq), whole genome sequencing (WGS), and eQTL analysis, focusing on efflux pump, transport, and secretion-associated genes. Of 42 ofloxacin-resistant Mtb isolates, 27 had adequate WGS coverage and acceptable RNA-seq quality. Of these 27, 7 had >2-fold reduction in ofloxacin MIC with verapamil; 6 had 2-fold reduction, and 14 had <2-fold reduction. Five genes (including Rv0191) had significantly increased expression in the MIC fold change >2 compared to <2 groups. Among regulated genes, 31 eQTLs (without ofloxacin) and 35 eQTLs (with ofloxacin) had significant allele frequency differences between MIC fold change >2 and <2 groups. Of these, Rv1410c, Rv2459, and Rv3756c (without ofloxacin) and Rv0191 and Rv3756c (with ofloxacin) have previously been associated with antituberculosis drug resistance. In this first reported eQTL analysis in Mtb, Rv0191 had increased gene expression and significance in eQTL analysis, making it a candidate for functional evaluation of efflux-mediated fluoroquinolone resistance in Mtb.

Deletion mapping of regulatory elements for GATA3 in T cells reveals a distal enhancer involved in allergic diseases

GATA3 is essential for Tcell differentiation and is surrounded by genome-wide association study (GWAS) hits for immune traits. Interpretation of these GWAS hits is challenging because gene expression quantitative trait locus (eQTL) studies lack power to detect variants with small effects on gene expression in specific cell types and the genome region containing GATA3 contains dozens of potential regulatory sequences. To map regulatory sequences for GATA3, we performed a high-throughput tiling deletion screen of a 2 Mb genome region in Jurkat Tcells. This revealed 23 candidate regulatory sequences, all but one of which is within the same topological-associating domain (TAD) as GATA3. We then performed a lower-throughput deletion screen to precisely map regulatory sequences in primary T helper 2 (Th2) cells. We tested 25 sequences with ∼100bp deletions and validated five of the strongest hits with independent deletion experiments. Additionally, we fine-mapped GWAS hits for allergic diseases in a distal regulatory element, 1 Mb downstream of GATA3, and identified 14 candidate causal variants. Small deletions spanning the candidate variant rs725861 decreased GATA3 levels in Th2 cells, and luciferase reporter assays showed regulatory differences between its two alleles, suggesting a causal mechanism for this variant in allergic diseases. Our study demonstrates the power of integrating GWAS signals with deletion mapping and identifies critical regulatory sequences for GATA3.

Exploiting the mediating role of the metabolome to unravel transcript-to-phenotype associations

Despite the success of genome-wide association studies (GWASs) in identifying genetic variants associated with complex traits, understanding the mechanisms behind these statistical associations remains challenging. Several methods that integrate methylation, gene expression, and protein quantitative trait loci (QTLs) with GWAS data to determine their causal role in the path from genotype to phenotype have been proposed. Here, we developed and applied a multi-omics Mendelian randomization (MR) framework to study how metabolites mediate the effect of gene expression on complex traits. We identified 216 transcript-metabolite-trait causal triplets involving 26 medically relevant phenotypes. Among these associations, 58% were missed by classical transcriptome-wide MR, which only uses gene expression and GWAS data. This allowed the identification of biologically relevant pathways, such as between ANKH and calcium levels mediated by citrate levels and SLC6A12 and serum creatinine through modulation of the levels of the renal osmolyte betaine. We show that the signals missed by transcriptome-wide MR are found, thanks to the increase in power conferred by integrating multiple omics layer. Simulation analyses show that with larger molecular QTL studies and in case of mediated effects, our multi-omics MR framework outperforms classical MR approaches designed to detect causal relationships between single molecular traits and complex phenotypes.

Germline Cancer Gene Expression Quantitative Trait Loci Are Associated with Local and Global Tumor Mutations.

Analysis of associations between the germline genetic background and somatic mutations in patients with cancer suggests that germline variants can influence local and global tumor mutations by altering expression of cancer-related genes. See related commentary by Kar, p. 1165.

Mitochondrial related genome-wide Mendelian randomization identifies putatively causal genes for multiple cancer types.

Mitochondrial dysfunction is a hallmark of cancer. However, it is unclear whether it is a cause of cancer. This two-sample Mendelian randomization (MR) analyses, uses genetic instruments to proxy the exposure of mitochondrial dysfunction and cancer summary statistics as outcomes, allowing for causal inferences. Summary statistics from 18 common cancers (2107-491,974 participants), gene expression, DNA methylation and protein expression quantitative trait loci (eQTL, mQTL and pQTL, respectively, 1000-31,684 participants) on individuals of European ancestry, were included. Genetic variants located within or close to the 1136 mitochondrial-related genes (in cis) and robustly associated with the mitochondrial molecular alterations were used as instrumental variables, and their causal associations with cancers were examined using summary-data-based MR (SMR) analyses. An additional five MR methods were used as sensitivity analyses to confirm the casual associations. A Bayesian test for colocalization between mitochondrial molecular QTLs and cancer risk loci was performed to provide insights into the potential regulatory mechanisms of risk variants on cancers. We identified potential causal relationships between mitochondrial-related genes and breast, prostate, gastric, lung cancer and melanoma by primary SMR analyses. The sensitivity and the colocalization analyses further refined four genes that have causal effects on three types of cancer. We found strong evidence of positive association of FDPS expression level with breast cancer risk (OR per SD, 0.66; 95% CI, 0.49-0.83; P=9.77×10-7), NSUN4 expression level with both breast cancer risk (OR per SD, 1.05; 95% CI, 1.03-1.07; P=5.24×10-6) and prostate cancer risk (OR per SD, 1.06; 95% CI, 1.03-1.09; P=1.01×10-5), NSUN4 methylation level with both breast and prostate cancer risk, and VARS2 methylation level with lung cancer risk. This data-driven MR study demonstrated the causal role of mitochondrial dysfunction in multiple cancers. Furthermore, this study identified candidate genes that can be the targets of potential pharmacological agents for cancer prevention. This work was supported by Styrelsen för Allmänna Sjukhusets i Malmö Stiftelse för bekämpande av cancer (20211025).

Integrative Analyses Identify KCNJ15 as a Candidate Gene in Patients with Epilepsy.

Although there is accumulating evidence that genetic factors play a vital role in the pathogenesis of epilepsy, few epilepsy-associated genes have been identified. Additionally, the role of KCNJ15 in epilepsy has not been evaluated so far. Here, we performed differentially expressed gene analysis, expression quantitative trait loci analysis, gene co-expression analysis, and protein-protein interaction analysis to evaluate the role of KCNJ15 in epilepsy. Analysis of gene expression and expression quantitative trait loci data revealed that KCNJ15 was significantly downregulated in patients with epilepsy (adjusted P = 0.0146 and log2 Fold change = - 1.0025), and an epilepsy-associated polymorphism (rs2833098) was linked to altered KCNJ15 expression level in human temporal lobe brain tissue (P = 0.0036). Gene co-expression analysis revealed that KCNJ15 was co-expressed with genes that have been reported to be associated with epilepsy in human brain tissue. Furthermore, protein-protein interaction analysis revealed strong supportive evidence for the role of KCNJ15 in epilepsy. Our results show that KCNJ15 may be a candidate target for epilepsy. Functional analysis of KCNJ15 may provide novel insights for epilepsy treatment.

Expressed genes and their new alleles identification during fibre elongation reveal the genetic factors underlying improvements of fibre length in cotton.

SummaryInterspecific breeding in cotton takes advantage of genetic recombination among desirable genes from different parental lines. However, the expression new alleles (ENAs) from crossovers within genic regions and their significance in fibre length (FL) improvement are currently not understood. Here, we generated resequencing genomes of 191 interspecific backcross inbred lines derived from CRI36 (Gossypium hirsutum) × Hai7124 (Gossypium barbadense) and 277 dynamic fibre transcriptomes to identify the ENAs and extremely expressed genes (eGenes) potentially influencing FL, and uncovered the dynamic regulatory network of fibre elongation. Of 35 420 eGenes in developing fibres, 10 366 ENAs were identified and preferentially distributed in chromosomes subtelomeric regions. In total, 1056–1255 ENAs showed transgressive expression in fibres at 5–15 dpa (days post‐anthesis) of some BILs, 520 of which were located in FL‐quantitative trait locus (QTLs) and GhFLA9 (recombination allele) was identified with a larger effect for FL than GhFLA9 of CRI36 allele. Using ENAs as a type of markers, we identified three novel FL‐QTLs. Additionally, 456 extremely eGenes were identified that were preferentially distributed in recombination hotspots. Importantly, 34 of them were significantly associated with FL. Gene expression quantitative trait locus analysis identified 1286, 1089 and 1059 eGenes that were colocalized with the FL trait at 5, 10 and 15 dpa, respectively. Finally, we verified the Ghir_D10G011050 gene linked to fibre elongation by the CRISPR‐cas9 system. This study provides the first glimpse into the occurrence, distribution and expression of the developing fibres genes (especially ENAs) in an introgression population, and their possible biological significance in FL.