The EN-TEx resource of multi-tissue personal epigenomes & variant-impact models.

Understanding how genetic variants impact molecular phenotypes is a key goal of functional genomics, currently hindered by reliance on a single haploid reference genome. Here, we present the EN-TEx resource of 1,635 open-access datasets from four donors (∼30 tissues × âˆ¼15 assays). The datasets are mapped to matched, diploid genomes with long-read phasing and structural variants, instantiating a catalog of >1 million allele-specific loci. These loci exhibit coordinated activity along haplotypes and are less conserved than corresponding, non-allele-specific ones. Surprisingly, a deep-learning transformer model can predict the allele-specific activity based only on local nucleotide-sequence context, highlighting the importance of transcription-factor-binding motifs particularly sensitive to variants. Furthermore, combining EN-TEx with existing genome annotations reveals strong associations between allele-specific and GWAS loci. It also enables models for transferring known eQTLs to difficult-to-profile tissues (e.g., from skin to heart). Overall, EN-TEx provides rich data and generalizable models for more accurate personal functional genomics.

Trans Effects on Gene Expression Can Drive Omnigenic Inheritance.

Early genome-wide association studies (GWASs) led to the surprising discovery that, for typical complex traits, most of the heritability is due to huge numbers of common variants with tiny effect sizes. Previously, we argued that new models are needed to understand these patterns. Here, we provide a formal model in which genetic contributions to complex traits are partitioned into direct effects from core genes and indirect effects from peripheral genes acting in trans. We propose that most heritability is driven by weak trans-eQTL SNPs, whose effects are mediated through peripheral genes to impact the expression of core genes. In particular, if the core genes for a trait tend to be co-regulated, then the effects of peripheral variation can be amplified such that nearly all of the genetic variance is driven by weak trans effects. Thus, our model proposes a framework for understanding key features of the architecture of complex traits.

Impact of Genetic Polymorphisms on Human Immune Cell Gene Expression.

While many genetic variants have been associated with risk for human diseases, how these variants affect gene expression in various cell types remains largely unknown. To address this gap, the DICE (database of immune cell expression, expression quantitative trait loci [eQTLs], and epigenomics) project was established. Considering all human immune cell types and conditions studied, we identified cis-eQTLs for a total of 12,254 unique genes, which represent 61% of all protein-coding genes expressed in these cell types. Strikingly, a large fraction (41%) of these genes showed a strong cis-association with genotype only in a single cell type. We also found that biological sex is associated with major differences in immune cell gene expression in a highly cell-specific manner. These datasets will help reveal the effects of disease risk-associated genetic polymorphisms on specific immune cell types, providing mechanistic insights into how they might influence pathogenesis (https://dice-database.org).

Partitioning gene-mediated disease heritability without eQTLs.

Unknown SNP-to-gene regulatory architecture complicates efforts to link noncoding GWAS associations with genes implicated by sequencing or functional studies. eQTLs are often used to link SNPs to genes, but expression in bulk tissue explains a small fraction of disease heritability. A simple but successful approach has been to link SNPs with nearby genes via base pair windows, but genes may often be regulated by SNPs outside their window. We propose the abstract mediation model (AMM) to estimate (1) the fraction of heritability mediated by the closest or kth-closest gene to each SNP and (2) the mediated heritability enrichment of a gene set (e.g., genes with rare-variant associations). AMM jointly estimates these quantities by matching the decay in SNP enrichment with distance from genes in the gene set. Across 47 complex traits and diseases, we estimate that the closest gene to each SNP mediates 27% (SE: 6%) of heritability and that a substantial fraction is mediated by genes outside the ten closest. Mendelian disease genes are strongly enriched for common-variant heritability; for example, just 21 dyslipidemia genes mediate 25% of LDL heritability (211× enrichment, p = 0.01). Among brain-related traits, genes involved in neurodevelopmental disorders are only about 4× enriched, but gene expression patterns are highly informative, as they have detectable differences in per-gene heritability even among weakly brain-expressed genes.

Pleiotropy with sex-specific traits reveals genetic aspects of sex differences in Parkinson's disease.

Parkinson's disease is an age-related neurodegenerative disorder with a higher incidence in males than females. The causes for this sex difference are unknown. Genome-wide association studies (GWAS) have identified 90 Parkinson's disease risk loci, but the genetic studies have not found sex-specific differences in allele frequency on autosomal chromosomes or sex chromosomes. Genetic variants, however, could exert sex-specific effects on gene function and regulation of gene expression. To identify genetic loci that might have sex-specific effects, we studied pleiotropy between Parkinson's disease and sex-specific traits. Summary statistics from GWASs were acquired from large-scale consortia for Parkinson's disease (n cases = 13 708; n controls = 95 282), age at menarche (n = 368 888 females) and age at menopause (n = 69 360 females). We applied the conditional/conjunctional false discovery rate (FDR) method to identify shared loci between Parkinson's disease and these sex-specific traits. Next, we investigated sex-specific gene expression differences in the superior frontal cortex of both neuropathologically healthy individuals and Parkinson's disease patients (n cases = 61; n controls = 23). To provide biological insights to the genetic pleiotropy, we performed sex-specific expression quantitative trait locus (eQTL) analysis and sex-specific age-related differential expression analysis for genes mapped to Parkinson's disease risk loci. Through conditional/conjunctional FDR analysis we found 11 loci shared between Parkinson's disease and the sex-specific traits age at menarche and age at menopause. Gene-set and pathway analysis of the genes mapped to these loci highlighted the importance of the immune response in determining an increased disease incidence in the male population. Moreover, we highlighted a total of nine genes whose expression or age-related expression in the human brain is influenced by genetic variants in a sex-specific manner. With these analyses we demonstrated that the lack of clear sex-specific differences in allele frequencies for Parkinson's disease loci does not exclude a genetic contribution to differences in disease incidence. Moreover, further studies are needed to elucidate the role that the candidate genes identified here could have in determining a higher incidence of Parkinson's disease in the male population.

Cis-eQTLs in seven duck tissues identify novel candidate genes for growth and carcass traits.

BACKGROUND: Expression quantitative trait loci (eQTL) studies aim to understand the influence of genetic variants on gene expression. The colocalization of eQTL mapping and GWAS strategy could help identify essential candidate genes and causal DNA variants vital to complex traits in human and many farm animals. However, eQTL mapping has not been conducted in ducks. It is desirable to know whether eQTLs within GWAS signals contributed to duck economic traits. RESULTS: In this study, we conducted an eQTL analysis using publicly available RNA sequencing data from 820 samples, focusing on liver, muscle, blood, adipose, ovary, spleen, and lung tissues. We identified 113,374 cis-eQTLs for 12,266 genes, a substantial fraction 39.1% of which were discovered in at least two tissues. The cis-eQTLs of blood were less conserved across tissues, while cis-eQTLs from any tissue exhibit a strong sharing pattern to liver tissue. Colocalization between cis-eQTLs and genome-wide association studies (GWAS) of 50 traits uncovered new associations between gene expression and potential loci influencing growth and carcass traits. SRSF4, GSS, and IGF2BP1 in liver, NDUFC2 in muscle, ELF3 in adipose, and RUNDC1 in blood could serve as the candidate genes for duck growth and carcass traits. CONCLUSIONS: Our findings highlight substantial differences in genetic regulation of gene expression across duck primary tissues, shedding light on potential mechanisms through which candidate genes may impact growth and carcass traits. Furthermore, this availability of eQTL data offers a valuable resource for deciphering further genetic association signals that may arise from ongoing extensive endeavors aimed at enhancing duck production traits.

Stratification of risk of progression to colectomy in ulcerative colitis via measured and predicted gene expression.

An important goal of clinical genomics is to be able to estimate the risk of adverse disease outcomes. Between 5% and 10% of individuals with ulcerative colitis (UC) require colectomy within 5 years of diagnosis, but polygenic risk scores (PRSs) utilizing findings from genome-wide association studies (GWASs) are unable to provide meaningful prediction of this adverse status. By contrast, in Crohn disease, gene expression profiling of GWAS-significant genes does provide some stratification of risk of progression to complicated disease in the form of a transcriptional risk score (TRS). Here, we demonstrate that a measured TRS based on bulk rectal gene expression in the PROTECT inception cohort study has a positive predictive value approaching 50% for colectomy. Single-cell profiling demonstrates that the genes are active in multiple diverse cell types from both the epithelial and immune compartments. Expression quantitative trait locus (QTL) analysis identifies genes with differential effects at baseline and week 52 follow-up, but for the most part, differential expression associated with colectomy risk is independent of local genetic regulation. Nevertheless, a predicted polygenic transcriptional risk score (PPTRS) derived by summation of transcriptome-wide association study (TWAS) effects identifies UC-affected individuals at 5-fold elevated risk of colectomy with data from the UK Biobank population cohort studies, independently replicated in an NIDDK-IBDGC dataset. Prediction of gene expression from relatively small transcriptome datasets can thus be used in conjunction with TWASs for stratification of risk of disease complications.

Genetic susceptibility of urolithiasis: comprehensive results from genome-wide analysis.

BACKGROUND: The pathogenesis of urolithiasis is multi-factorial and genetic factors have been shown to play a significant role in the development of urolithiasis. We tried to apply genome-wide Mendelian randomization (MR) analysis and figure out reliable gene susceptibility of urolithiasis from the largest samples to date in two independent genome-wide association studies (GWAS) database of European ancestry. METHODS: We extracted summary statistics of expression quantitative trait locus (eQTL) from eQTLGen consortium. Urolithiasis phenotype information was obtained from both FinnGen Biobank and UK Biobank. Multiple two-sample MR analysis with a Bonferroni-corrected P threshold (P < 2.5e-06) was conducted. The primary endpoint was the causal effect calculated by random-effect inverse variance weighted (IVW) method. Sensitivity analysis, volcano plots, scatter plots, and regional plots were also performed and visualized. RESULTS: After multiple MR tests between 19942 eQTLs and urolithiasis phenotype from both cohorts, 30 common eQTLs with consistent effect size direction were found to be causally associated with urolithiasis risk. Finally only one gene (LMAN2) was simultaneously identified among all top significant eQTLs from both FinnGen Biobank (beta = 0.6758, se = 0.0327, P = 6.775e-95) and UK Biobank (beta = 0.0044, se = 0.0009, P = 2.417e-06). We also found that LMAN2 was with the largest beta effect size on urolithiasis phenotype from the two cohorts. CONCLUSION: We for the first time implemented genome-wide MR analysis to investigate the genetic susceptibility of urolithiasis in general population of European ancestry. Our results provided novel insights into common genetic variants of urinary stone disease, which was of great help to subsequent researches.

Exploration of the causality of frailty index on psoriasis: A Mendelian randomization study.

BACKGROUND: Frailty is associated with a variety of diseases, but the relationship between frailty and psoriasis remains unclear. METHODS: First, we conducted a two-sample Mendelian randomization based on genome-wide association studies (GWAS) to investigate genetic causality between frailty index and common diseases in dermatology. Inverse variance weighted was used to estimate causality. Second, expression quantitative trait locus (eQTLs) analysis was conducted to identify the genes affected by Single nucleotide polymorphisms (SNPs). Third, we performed function and pathway enrichment, transcriptome-wide association studies (TWAS) analysis based on eQTLs. RESULTS: It was shown that the rise of frailty index could increase the risk of psoriasis (IVW, beta = 0.916, OR = 2.500, 95%CI:1.418-4.408, p = 0.002) through Mendelian randomization (MR), and there was no heterogeneity and pleiotropy. There was no causality between the frailty index and other common diseases in dermatology. We found 31 eQTLs based on strongly correlated SNPs in the causality. TWAS analysis found that the expressions of four genes were closely related to psoriasis, including HLA-DQA1, HLA-DQA2, HLA-DRB1 and HLA-DQB1. CONCLUSION: It suggested that the frailty index had a significant positive causality on the risk of psoriasis, which was well documented by combined genomic, transcriptome, and proteome analyses.

Identifying novel regulatory effects for clinically relevant genes through the study of the Greek population.

BACKGROUND: Expression quantitative trait loci (eQTL) studies provide insights into regulatory mechanisms underlying disease risk. Expanding studies of gene regulation to underexplored populations and to medically relevant tissues offers potential to reveal yet unknown regulatory variants and to better understand disease mechanisms. Here, we performed eQTL mapping in subcutaneous (S) and visceral (V) adipose tissue from 106 Greek individuals (Greek Metabolic study, GM) and compared our findings to those from the Genotype-Tissue Expression (GTEx) resource. RESULTS: We identified 1,930 and 1,515 eGenes in S and V respectively, over 13% of which are not observed in GTEx adipose tissue, and that do not arise due to different ancestry. We report additional context-specific regulatory effects in genes of clinical interest (e.g. oncogene ST7) and in genes regulating responses to environmental stimuli (e.g. MIR21, SNX33). We suggest that a fraction of the reported differences across populations is due to environmental effects on gene expression, driving context-specific eQTLs, and suggest that environmental effects can determine the penetrance of disease variants thus shaping disease risk. We report that over half of GM eQTLs colocalize with GWAS SNPs and of these colocalizations 41% are not detected in GTEx. We also highlight the clinical relevance of S adipose tissue by revealing that inflammatory processes are upregulated in individuals with obesity, not only in V, but also in S tissue. CONCLUSIONS: By focusing on an understudied population, our results provide further candidate genes for investigation regarding their role in adipose tissue biology and their contribution to disease risk and pathogenesis.

Evaluation of noninvasive biospecimens for transcriptome studies.

Transcriptome studies disentangle functional mechanisms of gene expression regulation and may elucidate the underlying biology of disease processes. However, the types of tissues currently collected typically assay a single post-mortem timepoint or are limited to investigating cell types found in blood. Noninvasive tissues may improve disease-relevant discovery by enabling more complex longitudinal study designs, by capturing different and potentially more applicable cell types, and by increasing sample sizes due to reduced collection costs and possible higher enrollment from vulnerable populations. Here, we develop methods for sampling noninvasive biospecimens, investigate their performance across commercial and in-house library preparations, characterize their biology, and assess the feasibility of using noninvasive tissues in a multitude of transcriptomic applications. We collected buccal swabs, hair follicles, saliva, and urine cell pellets from 19 individuals over three to four timepoints, for a total of 300 unique biological samples, which we then prepared with replicates across three library preparations, for a final tally of 472 transcriptomes. Of the four tissues we studied, we found hair follicles and urine cell pellets to be most promising due to the consistency of sample quality, the cell types and expression profiles we observed, and their performance in disease-relevant applications. This is the first study to thoroughly delineate biological and technical features of noninvasive samples and demonstrate their use in a wide array of transcriptomic and clinical analyses. We anticipate future use of these biospecimens will facilitate discovery and development of clinical applications.

Enhancer Domains Predict Gene Pathogenicity and Inform Gene Discovery in Complex Disease.

Non-coding transcriptional regulatory elements are critical for controlling the spatiotemporal expression of genes. Here, we demonstrate that the sizes and number of enhancers linked to a gene reflect its disease pathogenicity. Moreover, genes with redundant enhancer domains are depleted of cis-acting genetic variants that disrupt gene expression, and they are buffered against the effects of disruptive non-coding mutations. Our results demonstrate that dosage-sensitive genes have evolved a robustness to the disruptive effects of genetic variation by expanding their regulatory domains. This solves a puzzle about why genes associated with human disease are depleted of cis-eQTLs (cis-expression quantitative trait loci), suggesting that this relationship might complicate gene identification in causal genome-wide association studies (GWASs) using eQTL information, and establishes a framework for identifying non-coding regulatory variation with phenotypic consequences.

Genetic Architecture of Gene Expression in European and African Americans: An eQTL Mapping Study in GENOA.

Most existing expression quantitative trait locus (eQTL) mapping studies have been focused on individuals of European ancestry and are underrepresented in other populations including populations with African ancestry. Lack of large-scale well-powered eQTL mapping studies in populations with African ancestry can both impede the dissemination of eQTL mapping results that would otherwise benefit individuals with African ancestry and hinder the comparable analysis for understanding how gene regulation is shaped through evolution. We fill this critical knowledge gap by performing a large-scale in-depth eQTL mapping study on 1,032 African Americans (AA) and 801 European Americans (EA) in the GENOA cohort. We identified a total of 354,931 eSNPs in AA and 371,309 eSNPs in EA, with 112,316 eSNPs overlapped between the two. We found that eQTL harboring genes (eGenes) are enriched in metabolic pathways and tend to have higher SNP heritability compared to non-eGenes. We found that eGenes that are common in the two populations tend to be less conserved than eGenes that are unique to one population, which are less conserved than non-eGenes. Through conditional analysis, we found that eGenes in AA tend to harbor more independent eQTLs than eGenes in EA, suggesting potentially diverse genetic architecture underlying expression variation in the two populations. Finally, the large sample sizes in GENOA allow us to construct accurate expression prediction models in both AA and EA, facilitating powerful transcriptome-wide association studies. Overall, our results represent an important step toward revealing the genetic architecture underlying expression variation in African Americans.

A Multi-tissue Transcriptome Analysis of Human Metabolites Guides Interpretability of Associations Based on Multi-SNP Models for Gene Expression.

There is particular interest in transcriptome-wide association studies (TWAS) gene-level tests based on multi-SNP predictive models of gene expression-for identifying causal genes at loci associated with complex traits. However, interpretation of TWAS associations may be complicated by divergent effects of model SNPs on phenotype and gene expression. We developed an iterative modeling scheme for obtaining multi-SNP models of gene expression and applied this framework to generate expression models for 43 human tissues from the Genotype-Tissue Expression (GTEx) Project. We characterized the performance of single- and multi-SNP models for identifying causal genes in GWAS data for 46 circulating metabolites. We show that: (A) multi-SNP models captured more variation in expression than did the top cis-eQTL (median 2-fold improvement); (B) predicted expression based on multi-SNP models was associated (false discovery rate < 0.01) with metabolite levels for 826 unique gene-metabolite pairs, but, after stepwise conditional analyses, 90% were dominated by a single eQTL SNP; (C) among the 35% of associations where a SNP in the expression model was a significant cis-eQTL and metabolomic-QTL (met-QTL), 92% demonstrated colocalization between these signals, but interpretation was often complicated by incomplete overlap of QTLs in multi-SNP models; and (D) using a "truth" set of causal genes at 61 met-QTLs, the sensitivity was high (67%), but the positive predictive value was low, as only 8% of TWAS associations (19% when restricted to colocalized associations at met-QTLs) involved true causal genes. These results guide the interpretation of TWAS and highlight the need for corroborative data to provide confident assignment of causality.

Expression quantitative trait loci mapping identified PtrXB38 as a key hub gene in adventitious root development in Populus.

Plant establishment requires the formation and development of an extensive root system with architecture modulated by complex genetic networks. Here, we report the identification of the PtrXB38 gene as an expression quantitative trait loci (eQTL) hotspot, mapped using 390 leaf and 444 xylem Populus trichocarpa transcriptomes. Among predicted targets of this trans-eQTL were genes involved in plant hormone responses and root development. Overexpression of PtrXB38 in Populus led to significant increases in callusing and formation of both stem-born roots and base-born adventitious roots. Omics studies revealed that genes and proteins controlling auxin transport and signaling were involved in PtrXB38-mediated adventitious root formation. Protein-protein interaction assays indicated that PtrXB38 interacts with components of endosomal sorting complexes required for transport machinery, implying that PtrXB38-regulated root development may be mediated by regulating endocytosis pathway. Taken together, this work identified a crucial root development regulator and sheds light on the discovery of other plant developmental regulators through combining eQTL mapping and omics approaches.

Causality assessment of circulating Vitamin D level on venous thromboembolism: A Mendelian randomization study.

BACKGROUND AND AIMS: The associations of vitamin D level with venous thromboembolism (VTE) reported in observational studies, whereas these causal associations were uncertain in European population. Therefore, we used Mendelian randomization (MR) method to explore the causal associations between 25-hydroxyvitamin D (25(OH)D) concentrations and the risk of VTE and its subtypes [including deep vein thrombosis (DVT) and pulmonary embolism (PE)]. METHODS AND RESULTS: We used three kinds of genetic instruments to proxy the exposure of 25(OH)D, including genetic variants significantly associated with 25(OH)D, expression quantitative trait loci of 25(OH)D target genes, and genetic variants within or nearby 25(OH)D target genes. MR analyses did not provide any evidence for the associations of 25(OH)D levels with VTE and its subtypes (p > 0.05). The summary-data-based MR (SMR) analyses indicated that elevated expression of VDR was associated with decreased risk of VTE (OR = 0.81; 95% CI, 0.65-0.998; p = 0.047) and PE (OR = 0.67; 95% CI, 0.50-0.91; p = 0.011), and expression of AMDHD1 was associated with PE (OR = 0.93; 95% CI, 0.88-0.99; p = 0.027). MR analysis provided a significant causal effect of 25(OH)D level mediated by gene AMDHD1 on PE risk (OR = 0.09; 95% CI, 0.01-0.60; p = 0.012). CONCLUSION: Our MR analysis did not support causal association of 25(OH)D level with the risk of VTE and its subtypes. In addition, the expression of VDR and AMDHD1 involved in vitamin D metabolism showed a strong association with VTE or PE and might represent targets for these conditions.

Genetic variants associated mRNA stability in lung.

BACKGROUND: Expression quantitative trait loci (eQTLs) analyses have been widely used to identify genetic variants associated with gene expression levels to understand what molecular mechanisms underlie genetic traits. The resultant eQTLs might affect the expression of associated genes through transcriptional or post-transcriptional regulation. In this study, we attempt to distinguish these two types of regulation by identifying genetic variants associated with mRNA stability of genes (stQTLs). RESULTS: Here, we presented a computational framework that takes advantage of recently developed methods to infer the mRNA stability of genes based on RNA-seq data and performed association analysis to identify stQTLs. Using the Genotype-Tissue Expression (GTEx) lung RNA-Seq data, we identified a total of 142,801 stQTLs for 3942 genes and 186,132 eQTLs for 4751 genes from 15,122,700 genetic variants for 13,476 genes on the autosomes, respectively. Interestingly, our results indicated that stQTLs were enriched in the CDS and 3'UTR regions, while eQTLs are enriched in the CDS, 3'UTR, 5'UTR, and upstream regions. We also found that stQTLs are more likely than eQTLs to overlap with RNA binding protein (RBP) and microRNA (miRNA) binding sites. Our analyses demonstrate that simultaneous identification of stQTLs and eQTLs can provide more mechanistic insight on the association between genetic variants and gene expression levels.

Multivariate Genome-wide Association Analysis of a Cytokine Network Reveals Variants with Widespread Immune, Haematological, and Cardiometabolic Pleiotropy.

Cytokines are essential regulatory components of the immune system, and their aberrant levels have been linked to many disease states. Despite increasing evidence that cytokines operate in concert, many of the physiological interactions between cytokines, and the shared genetic architecture that underlies them, remain unknown. Here, we aimed to identify and characterize genetic variants with pleiotropic effects on cytokines. Using three population-based cohorts (n = 9,263), we performed multivariate genome-wide association studies (GWAS) for a correlation network of 11 circulating cytokines, then combined our results in meta-analysis. We identified a total of eight loci significantly associated with the cytokine network, of which two (PDGFRB and ABO) had not been detected previously. In addition, conditional analyses revealed a further four secondary signals at three known cytokine loci. Integration, through the use of Bayesian colocalization analysis, of publicly available GWAS summary statistics with the cytokine network associations revealed shared causal variants between the eight cytokine loci and other traits; in particular, cytokine network variants at the ABO, SERPINE2, and ZFPM2 loci showed pleiotropic effects on the production of immune-related proteins, on metabolic traits such as lipoprotein and lipid levels, on blood-cell-related traits such as platelet count, and on disease traits such as coronary artery disease and type 2 diabetes.

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

Interspecific 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.

eQTLs play critical roles in regulating gene expression and identifying key regulators in rice.

The regulation of gene expression plays an essential role in both the phenotype and adaptation of plants. Transcriptome sequencing enables simultaneous identification of exonic variants and quantification of gene expression. Here, we sequenced the leaf transcriptomes of 287 rice accessions from around the world and obtained a total of 177 853 high-quality single nucleotide polymorphisms after filtering. Genome-wide association study identified 44 354 expression quantitative trait loci (eQTLs), which regulate the expression of 13 201 genes, as well as 17 local eQTL hotspots and 96 distant eQTL hotspots. Furthermore, a transcriptome-wide association study screened 21 candidate genes for starch content in the flag leaves at the heading stage. HS002 was identified as a significant distant eQTL hotspot with five downstream genes enriched for diterpene antitoxin synthesis. Co-expression analysis, eQTL analysis, and linkage mapping together demonstrated that bHLH026 acts as a key regulator to activate the expression of downstream genes. The transgenic assay revealed that bHLH026 is an important regulator of diterpenoid antitoxin synthesis and enhances the disease resistance of rice. These findings improve our knowledge of the regulatory mechanisms of gene expression variation and complex regulatory networks of the rice genome and will facilitate genetic improvement of cultivated rice varieties.