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.

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.

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.

Integrating Gene Expression with Summary Association Statistics to Identify Genes Associated with 30 Complex Traits.

Although genome-wide association studies (GWASs) have identified thousands of risk loci for many complex traits and diseases, the causal variants and genes at these loci remain largely unknown. Here, we introduce a method for estimating the local genetic correlation between gene expression and a complex trait and utilize it to estimate the genetic correlation due to predicted expression between pairs of traits. We integrated gene expression measurements from 45 expression panels with summary GWAS data to perform 30 multi-tissue transcriptome-wide association studies (TWASs). We identified 1,196 genes whose expression is associated with these traits; of these, 168 reside more than 0.5 Mb away from any previously reported GWAS significant variant. We then used our approach to find 43 pairs of traits with significant genetic correlation at the level of predicted expression; of these, eight were not found through genetic correlation at the SNP level. Finally, we used bi-directional regression to find evidence that BMI causally influences triglyceride levels and that triglyceride levels causally influence low-density lipoprotein. Together, our results provide insight into the role of gene expression in the susceptibility of complex traits and diseases.

Genome-wide expression quantitative trait locus analysis in a recombinant inbred line population for trait dissection in peanut.

The transcriptome connects genome to the gene function and ultimate phenome in biology. So far, transcriptomic approach was not used in peanut for performing trait mapping in bi-parental populations. In this research, we sequenced the whole transcriptome in immature seeds in a peanut recombinant inbred line (RIL) population and explored thoroughly the landscape of transcriptomic variations and its genetic basis. The comprehensive analysis identified total 49 691 genes in RIL population, of which 92 genes followed a paramutation-like expression pattern. Expression quantitative trait locus (eQTL) analysis identified 1207 local eQTLs and 15 837 distant eQTLs contributing to the whole-genome transcriptomic variation in peanut. There were 94 eQTL hot spot regions detected across the genome with the dominance of distant eQTL. By integrating transcriptomic profile and annotation analyses, we unveiled a putative candidate gene and developed a linked marker InDel02 underlying a major QTL responsible for purple testa colour in peanut. Our result provided a first understanding of genetic basis of whole-genome transcriptomic variation in peanut and illustrates the potential of the transcriptome-aid approach in dissecting important traits in non-model plants.

A cis-eQTL allele regulating reduced expression of CHI3L1 is associated with late-onset adult asthma in Japanese cohorts.

BACKGROUND: The chitinase-like protein YKL-40 plays a major role in inhibiting the inflammasome. Deregulation of inflammasome activation is emerging as a key modulator of pathologic airway inflammation in patients with asthma. We determined whether cis-expression quantitative trait loci (eQTLs) of the gene that encodes YKL-40, chitinase 3-like 1 (CHI3L1), are involved in the onset of asthma or in specific asthma phenotypes. METHODS: This case-control study, which was conducted at the University of Tsukuba, Japan, included a total of 2709 adults from the Tsukuba genome-wide association study (GWAS) cohort (734 healthy volunteers and 237 asthma patients), the Tsukuba replication cohort (375 healthy adult volunteers and 381 adult asthma patients), and the Hokkaido replication cohort (554 healthy adult volunteers and 428 adult asthma patients). Among 34 cis-eQTLs in CHI3L1 in the lung, rs946261 was associated with adult asthma in these Japanese cohorts. The genetic impact of rs946261 on asthma was also examined according to the age at onset and adult asthma clusters. RESULTS: In the Tsukuba GWAS cohort, the C allele at rs946261 was significantly associated with reduced expression of CHI3L1 mRNA in the lung and with development of asthma (odds ratio (OR) 1.27; P = 0.036). The association was also observed following analysis of the three Japanese cohorts (OR 1.16; P = 0.013). A stronger association was found with late-onset asthma that developed at 41 years of age or later (OR 1.24; 95% confidence interval (CI) 1.07-1.45; P = 0.0058) and with a specific asthma phenotype characterized by late onset, less atopy, and mild airflow obstruction (OR 1.29; 95% CI 1.03-1.61; P = 0.027). CONCLUSIONS: The genotype consisting of the cis-eQTL allele that reduces expression of CHI3L1 was specifically associated with late-onset adult asthma. Given the important role of YKL-40 in many pathophysiological processes, including cell growth, migration, chemotaxis, reorganization, and tissue remodeling, it may be involved in an important pathogenic role in the establishment of inflammation and remodeling in asthmatic airways. Our findings may indicate the presence of a specific endotype related to exaggerated activation of YKL-40 in the pathogenesis of late-onset adult asthma.

Genome-wide association studies and expression-based quantitative trait loci analyses reveal roles of HCT2 in caffeoylquinic acid biosynthesis and its regulation by defense-responsive transcription factors in Populus.

3-O-caffeoylquinic acid, also known as chlorogenic acid (CGA), functions as an intermediate in lignin biosynthesis in the phenylpropanoid pathway. It is widely distributed among numerous plant species and acts as an antioxidant in both plants and animals. Using GC-MS, we discovered consistent and extreme variation in CGA content across a population of 739 4-yr-old Populus trichocarpa accessions. We performed genome-wide association studies (GWAS) from 917 P. trichocarpa accessions and expression-based quantitative trait loci (eQTL) analyses to identify key regulators. The GWAS and eQTL analyses resolved an overlapped interval encompassing a hydroxycinnamoyl-CoA:shikimate hydroxycinnamoyl transferase 2 (PtHCT2) that was significantly associated with CGA and partially characterized metabolite abundances. PtHCT2 leaf expression was significantly correlated with CGA abundance and it was regulated by cis-eQTLs containing W-box for WRKY binding. Among all nine PtHCT homologs, PtHCT2 is the only one that responds to infection by the fungal pathogen Sphaerulina musiva (a Populus pathogen). Validation using protoplast-based transient expression system suggests that PtHCT2 is regulated by the defense-responsive WRKY. These results are consistent with reports of CGA functioning as an antioxidant in response to biotic stress. This study provides insights into data-driven and omics-based inference of gene function in woody species.

Single cell eQTL analysis identifies cell type-specific genetic control of gene expression in fibroblasts and reprogrammed induced pluripotent stem cells.

BACKGROUND: The discovery that somatic cells can be reprogrammed to induced pluripotent stem cells (iPSCs) has provided a foundation for in vitro human disease modelling, drug development and population genetics studies. Gene expression plays a critical role in complex disease risk and therapeutic response. However, while the genetic background of reprogrammed cell lines has been shown to strongly influence gene expression, the effect has not been evaluated at the level of individual cells which would provide significant resolution. By integrating single cell RNA-sequencing (scRNA-seq) and population genetics, we apply a framework in which to evaluate cell type-specific effects of genetic variation on gene expression. RESULTS: Here, we perform scRNA-seq on 64,018 fibroblasts from 79 donors and map expression quantitative trait loci (eQTLs) at the level of individual cell types. We demonstrate that the majority of eQTLs detected in fibroblasts are specific to an individual cell subtype. To address if the allelic effects on gene expression are maintained following cell reprogramming, we generate scRNA-seq data in 19,967 iPSCs from 31 reprogramed donor lines. We again identify highly cell type-specific eQTLs in iPSCs and show that the eQTLs in fibroblasts almost entirely disappear during reprogramming. CONCLUSIONS: This work provides an atlas of how genetic variation influences gene expression across cell subtypes and provides evidence for patterns of genetic architecture that lead to cell type-specific eQTL effects.

Integrative genomic deconvolution of rheumatoid arthritis GWAS loci into gene and cell type associations.

BACKGROUND: Although genome-wide association studies (GWAS) have identified over 100 genetic loci associated with rheumatoid arthritis (RA), our ability to translate these results into disease understanding and novel therapeutics is limited. Most RA GWAS loci reside outside of protein-coding regions and likely affect distal transcriptional enhancers. Furthermore, GWAS do not identify the cell types where the associated causal gene functions. Thus, mapping the transcriptional regulatory roles of GWAS hits and the relevant cell types will lead to better understanding of RA pathogenesis. RESULTS: We combine the whole-genome sequences and blood transcription profiles of 377 RA patients and identify over 6000 unique genes with expression quantitative trait loci (eQTLs). We demonstrate the quality of the identified eQTLs through comparison to non-RA individuals. We integrate the eQTLs with immune cell epigenome maps, RA GWAS risk loci, and adjustment for linkage disequilibrium to propose target genes of immune cell enhancers that overlap RA risk loci. We examine 20 immune cell epigenomes and perform a focused analysis on primary monocytes, B cells, and T cells. CONCLUSIONS: We highlight cell-specific gene associations with relevance to RA pathogenesis including the identification of FCGR2B in B cells as possessing both intragenic and enhancer regulatory GWAS hits. We show that our RA patient cohort derived eQTL network is more informative for studying RA than that from a healthy cohort. While not experimentally validated here, the reported eQTLs and cell type-specific RA risk associations can prioritize future experiments with the goal of elucidating the regulatory mechanisms behind genetic risk associations.