Expression Genome-wide Association Studies Research Articles

HASTY-mediated miRNA dynamics modulate nitrogen starvation-induced leaf senescence in Arabidopsis

Nitrogen (N) deficiency responses are essential for plant survival and reproduction. Here, via an expression genome-wide association study (eGWAS), we reveal a mechanism that regulates microRNA (miRNA) dynamics necessary for N deficiency responses in Arabidopsis. Differential expression levels of three NAC transcription factor (TF) genes involved in leaf N deficiency responses among Arabidopsis accessions are most significantly associated with polymorphisms in HASTY (HST), which encodes an importin/exportin family protein responsible for the generation of mature miRNAs. HST acts as a negative regulator of N deficiency-induced leaf senescence, and the disruption and overexpression of HST differently modifies miRNA dynamics in response to N deficiency, altering levels of miRNAs targeting transcripts. Interestingly, N deficiency prevents the interaction of HST with HST-interacting proteins, DCL1 and RAN1, and some miRNAs. This suggests that HST-mediated regulation of miRNA dynamics collectively controls regulations mediated by multiple N deficiency response-associated NAC TFs, thereby being central to the N deficiency response network.

Expression genome-wide association study identifies key regulatory variants enriched with metabolic and immune functions in four porcine tissues

BackgroundIntegration of high throughput DNA genotyping and RNA-sequencing data enables the discovery of genomic regions that regulate gene expression, known as expression quantitative trait loci (eQTL). In pigs, efforts to date have been mainly focused on purebred lines for traits with commercial relevance as such growth and meat quality. However, little is known on genetic variants and mechanisms associated with the robustness of an animal, thus its overall health status. Here, the liver, lung, spleen, and muscle transcriptomes of 100 three-way crossbred female finishers were studied, with the aim of identifying novel eQTL regulatory regions and transcription factors (TFs) associated with regulation of porcine metabolism and health-related traits.ResultsAn expression genome-wide association study with 535,896 genotypes and the expression of 12,680 genes in liver, 13,310 genes in lung, 12,650 genes in spleen, and 12,595 genes in muscle resulted in 4,293, 10,630, 4,533, and 6,871 eQTL regions for each of these tissues, respectively. Although only a small fraction of the eQTLs were annotated as cis-eQTLs, these presented a higher number of polymorphisms per region and significantly stronger associations with their target gene compared to trans-eQTLs. Between 20 and 115 eQTL hotspots were identified across the four tissues. Interestingly, these were all enriched for immune-related biological processes. In spleen, two TFs were identified: ERF and ZNF45, with key roles in regulation of gene expression.ConclusionsThis study provides a comprehensive analysis with more than 26,000 eQTL regions identified that are now publicly available. The genomic regions and their variants were mostly associated with tissue-specific regulatory roles. However, some shared regions provide new insights into the complex regulation of genes and their interactions that are involved with important traits related to metabolism and immunity.

Variation in TaSPL6-D confers salinity tolerance in bread wheat by activating TaHKT1;5-D while preserving yield-related traits.

Na+ exclusion from above-ground tissues via the Na+-selective transporter HKT1;5 is a major salt-tolerance mechanism in crops. Using the expression genome-wide association study and yeast-one-hybrid screening, we identified TaSPL6-D, a transcriptional suppressor of TaHKT1;5-D in bread wheat. SPL6 also targeted HKT1;5 in rice and Brachypodium. A 47-bp insertion in the first exon of TaSPL6-D resulted in a truncated peptide, TaSPL6-DIn, disrupting TaHKT1;5-D repression exhibited by TaSPL6-DDel. Overexpressing TaSPL6-DDel, but not TaSPL6-DIn, led to inhibited TaHKT1;5-D expression and increased salt sensitivity. Knockout of TaSPL6-DDel in two wheat genotypes enhanced salinity tolerance, which was attenuated by a further TaHKT1;5-D knockdown. Spike development was preserved in Taspl6-dd mutants but not in Taspl6-aabbdd mutants. TaSPL6-DIn was mainly present in landraces, and molecular-assisted introduction of TaSPL6-DIn from a landrace into a leading wheat cultivar successfully improved yield on saline soils. The SPL6-HKT1;5 module offers a target for the molecular breeding of salt-tolerant crops.

Large-scale transcriptomic and genomic analyses reveal a novel functional gene SERPINB6 for chicken carcass traits

BackgroundCarcass traits are crucial indicators of meat production efficiency. However, the molecular regulatory mechanisms associated with these traits remain unclear.ResultsIn this study, we conducted comprehensive transcriptomic and genomic analyses on 399 Tiannong partridge chickens to identify key genes and variants associated with carcass traits and to elucidate the underlying regulatory mechanisms. Based on association analyses with the elastic net (EN) model, we identified 12 candidate genes (AMY1A, AP3B2, CEBPG, EEF2, EIF4EBP1, FGFR1, FOXD3, GOLM1, LOC107052698, PABPC1, SERPINB6 and TBC1D16) for 4 carcass-related traits, namely live weight, dressed weight, eviscerated weight, and breast muscle weight. SERPINB6 was identified as the only overlapping gene by 3 analyses, EN model analysis, weighted gene co-expression network analysis and differential expression analysis. Cell-level experiments confirmed that SERPINB6 promotes the proliferation of chicken DF1 cells and primary myoblasts. Further expression genome-wide association study and association analysis indicated that rs317934171 is the critical site that enhances SERPINB6 expression. Furthermore, a dual-luciferase reporter assay proved that gga-miR-1615 targets the 3′UTR of SERPINB6.ConclusionsCollectively, our findings reveal that SERPINB6 serves as a novel gene for chicken carcass traits by promoting fibroblast and myoblast proliferation. Additionally, the downstream variant rs317934171 regulates SERPINB6 expression. These results identify a new target gene and molecular marker for the molecular mechanisms of chicken carcass traits.

A maize semi-dwarf mutant reveals a GRAS transcription factor involved in brassinosteroid signaling.

Brassinosteroids (BR) and gibberellins (GA) regulate plant height and leaf angle in maize (Zea mays). Mutants with defects in BR or GA biosynthesis or signaling identify components of these pathways and enhance our knowledge about plant growth and development. In this study, we characterized three recessive mutant alleles of GRAS transcription factor 42 (gras42) in maize, a GRAS transcription factor gene orthologous to the DWARF AND LOW TILLERING (DLT) gene of rice (Oryza sativa). These maize mutants exhibited semi-dwarf stature, shorter and wider leaves, and more upright leaf angle. Transcriptome analysis revealed a role for GRAS42 as a determinant of BR signaling. Analysis of the expression consequences from loss of GRAS42 in the gras42-mu1021149 mutant indicated a weak loss of BR signaling in the mutant, consistent with its previously demonstrated role in BR signaling in rice. Loss of BR signaling was also evident by the enhancement of weak BR biosynthetic mutant alleles in double mutants of nana plant1-1 and gras42-mu1021149. The gras42-mu1021149 mutant had little effect on GA-regulated gene expression, suggesting that GRAS42 is not a regulator of core GA signaling genes in maize. Single-cell expression data identified gras42 expressed among cells in the G2/M phase of the cell cycle consistent with its previously demonstrated role in cell cycle gene expression in Arabidopsis (Arabidopsis thaliana). Cis-acting natural variation controlling GRAS42 transcript accumulation was identified by expression genome-wide association study (eGWAS) in maize. Our results demonstrate a conserved role for GRAS42/SCARECROW-LIKE 28 (SCL28)/DLT in BR signaling, clarify the role of this gene in GA signaling, and suggest mechanisms of tillering and leaf angle control by BR.

Identification of genomic regions, genetic variants and gene networks regulating candidate genes for lipid metabolism in pig muscle

The intramuscular fat content and fatty acid composition of porcine meat have a significant impact on its quality and nutritional value. This research aimed to investigate the expression of 45 genes involved in lipid metabolism in the longissimus dorsi muscle of three experimental pig backcrosses, with a 25% of Iberian background. To achieve this objective, we conducted an expression Genome-Wide Association Study (eGWAS) using gene expression levels in muscle measured by high-throughput real-time qPCR for 45 target genes and genotypes from the PorcineSNP60 BeadChip or Axiom Porcine Genotyping Array and 65 single nucleotide polymorphisms (SNPs) located in 20 genes genotyped by a custom-designed Taqman OpenArray in a cohort of 354 animals. The eGWAS analysis identified 301 eSNPs associated with 18 candidate genes (ANK2, APOE, ARNT, CIITA, CPT1A, EGF, ELOVL6, ELOVL7, FADS3, FASN, GPAT3, NR1D2, NR1H2, PLIN1, PPAP2A, RORA, RXRA and UCP3). Three cis-eQTL (expression quantitative trait loci) were identified for GPAT3, RXRA, and UCP3 genes, which indicates that a genetic polymorphism proximal to the same gene is affecting its expression. Furthermore, 24 trans-eQTLs were detected, and eight candidate regulatory genes were located in these genomic regions. Additionally, two trans-regulatory hotspots in Sus scrofa chromosomes 13 and 15 were identified. Moreover, a co-expression analysis performed on 89 candidate genes and the fatty acid composition revealed the regulatory role of four genes (FABP5, PPARG, SCD, and SREBF1). These genes modulate the levels of α-linolenic, arachidonic, and oleic acids, as well as regulating the expression of other candidate genes associated with lipid metabolism. The findings of this study offer novel insights into the functional regulatory mechanism of genes involved in lipid metabolism, thereby enhancing our understanding of this complex biological process.

A role for heritable transcriptomic variation in maize adaptation to temperate environments

BackgroundTranscription bridges genetic information and phenotypes. Here, we evaluated how changes in transcriptional regulation enable maize (Zea mays), a crop originally domesticated in the tropics, to adapt to temperate environments.ResultWe generated 572 unique RNA-seq datasets from the roots of 340 maize genotypes. Genes involved in core processes such as cell division, chromosome organization and cytoskeleton organization showed lower heritability of gene expression, while genes involved in anti-oxidation activity exhibited higher expression heritability. An expression genome-wide association study (eGWAS) identified 19,602 expression quantitative trait loci (eQTLs) associated with the expression of 11,444 genes. A GWAS for alternative splicing identified 49,897 splicing QTLs (sQTLs) for 7614 genes. Genes harboring both cis-eQTLs and cis-sQTLs in linkage disequilibrium were disproportionately likely to encode transcription factors or were annotated as responding to one or more stresses. Independent component analysis of gene expression data identified loci regulating co-expression modules involved in oxidation reduction, response to water deprivation, plastid biogenesis, protein biogenesis, and plant-pathogen interaction. Several genes involved in cell proliferation, flower development, DNA replication, and gene silencing showed lower gene expression variation explained by genetic factors between temperate and tropical maize lines. A GWAS of 27 previously published phenotypes identified several candidate genes overlapping with genomic intervals showing signatures of selection during adaptation to temperate environments.ConclusionOur results illustrate how maize transcriptional regulatory networks enable changes in transcriptional regulation to adapt to temperate regions.

Polymorphisms in cis-elements confer SAUR26 gene expression difference for thermo-response natural variation in Arabidopsis.

The SAUR26 subfamily genes play an important role in conferring variations of thermo-responsiveness of growth architecture among natural accessions of Arabidopsis thaliana. The expression variations are critical for their activity variations, but how expression variations are generated is unknown. We identified genetic loci for gene expression variations through expression genome-wide association study (eGWAS) and investigated their mechanisms through molecular analyses. We found that cis elements are the major determinants for expression variations of SAUR26, SAUR27, and SAUR28. Polymorphisms in the promoter region likely impact PIF4 regulation while those at the 3'UTR affect mRNA stability to generate variations in SAUR26 expression levels. These polymorphisms also differentially affect the mRNA stability of SAUR26 at two temperatures. This study reveals two mechanisms involving cis elements in generating gene expression diversity, which is likely important for local adaptations in Arabidopsis natural accessions.

Expression genome-wide association study identifies that phosphatidylinositol-derived signalling regulates ALUMINIUM SENSITIVE3 expression under aluminium stress in the shoots of Arabidopsis thaliana

To identify unknown regulatory mechanisms leading to aluminium (Al)-induction of the Al tolerance gene ALS3, we conducted an expression genome-wide association study (eGWAS) for ALS3 in the shoots of 95 Arabidopsis thaliana accessions in the presence of Al. The eGWAS was conducted using a mixed linear model with 145,940 genome-wide single nucleotide polymorphisms (SNPs) and the association results were validated using reverse genetics. We found that many SNPs from the eGWAS were associated with genes related to phosphatidylinositol metabolism as well as stress signal transduction, including Ca2+signals, inter-connected in a co-expression network. Of these, PLC9, CDPK32, ANAC071, DIR1, and a hypothetical protein (AT4G10470) possessed amino acid sequence/ gene expression level polymorphisms that were significantly associated with ALS3 expression level variation. Furthermore, T-DNA insertion mutants of PLC9, CDPK32, and ANAC071 suppressed shoot ALS3 expression in the presence of Al. This study clarified the regulatory mechanisms of ALS3 expression in the shoot and provided genetic evidence of the involvement of phosphatidylinositol-derived signal transduction under Al stress.

Identification of eQTLs associated with lipid metabolism in Longissimus dorsi muscle of pigs with different genetic backgrounds

Intramuscular fat content and its fatty acid composition affect porcine meat quality and its nutritional value. The present work aimed to identify genomic variants regulating the expression in the porcine muscle (Longissimus dorsi) of 45 candidate genes for lipid metabolism and fatty acid composition in three experimental backcrosses based on the Iberian breed. Expression genome-wide association studies (eGWAS) were performed between the muscle gene expression values, measured by real-time quantitative PCR, and the genotypes of 38,426 SNPs distributed along all chromosomes. The eGWAS identified 186 eSNPs located in ten Sus scrofa regions and associated with the expression of ACSM5, ACSS2, ATF3, DGAT2, FOS and IGF2 (FDR < 0.05) genes. Two expression quantitative trait loci (eQTLs) for IGF2 and ACSM5 were classified as cis-acting eQTLs, suggesting a mutation in the same gene affecting its expression. Conversely, ten eQTLs showed trans-regulatory effects on gene expression. When the eGWAS was performed for each backcross independently, only three common trans-eQTL regions were observed, indicating different regulatory mechanisms or allelic frequencies among the breeds. In addition, hotspot regions regulating the expression of several genes were detected. Our results provide new data to better understand the functional regulatory mechanisms of lipid metabolism genes in muscle.

Alzheimer's Disease Risk Variant rs2373115 Regulates GAB2 and NARS2 Expression in Human Brain Tissues.

Genetic association studies have identified significant association between the GAB2 rs2373115 variant and Alzheimer's disease (AD). However, it is unknown whether rs2373115 affects the regulation of nearby genes. Here, we evaluate the potential effect of rs2373115 on gene expression using multiple eQTL (expression quantitative trait loci) datasets from human brain tissues from the Mayo Clinic brain expression genome-wide association study (eGWAS), the UK Brain Expression Consortium (UKBEC), the Genotype-Tissue Expression (GTEx) project, and the Brain xQTL Serve. Our findings indicate that the rs2373115 C allele is associated with increased NARS2 expression, and both reduced and increased GAB2 expression in human tissues. Using a large-scale AD case-control expression dataset, we found increased GAB2 expression and reduced NARS2 expression in AD cases compared with controls. We believe that our findings provide important information regarding the rs2373115 variant and expression of nearby genes with respect to AD risk.

Deciphering the genetic regulation of peripheral blood transcriptome in pigs through expression genome-wide association study and allele-specific expression analysis

BackgroundEfforts to improve sustainability in livestock production systems have focused on two objectives: investigating the genetic control of immune function as it pertains to robustness and disease resistance, and finding predictive markers for use in breeding programs. In this context, the peripheral blood transcriptome represents an important source of biological information about an individual’s health and immunological status, and has been proposed for use as an intermediate phenotype to measure immune capacity. The objective of this work was to study the genetic architecture of variation in gene expression in the blood of healthy young pigs using two approaches: an expression genome-wide association study (eGWAS) and allele-specific expression (ASE) analysis.ResultsThe blood transcriptomes of 60-day-old Large White pigs were analyzed by expression microarrays for eGWAS (242 animals) and by RNA-Seq for ASE analysis (38 animals). Using eGWAS, the expression levels of 1901 genes were found to be associated with expression quantitative trait loci (eQTLs). We recovered 2839 local and 1752 distant associations (Single Nucleotide Polymorphism or SNP located less or more than 1 Mb from expression probe, respectively). ASE analyses confirmed the extensive cis-regulation of gene transcription in blood, and revealed allelic imbalance in 2286 SNPs, which affected 763 genes. eQTLs and ASE-genes were widely distributed on all chromosomes. By analyzing mutually overlapping eGWAS results, we were able to describe putative regulatory networks, which were further refined using ASE data. At the functional level, genes with genetically controlled expression that were detected by eGWAS and/or ASE analyses were significantly enriched in biological processes related to RNA processing and immune function. Indeed, numerous distant and local regulatory relationships were detected within the major histocompatibility complex region on chromosome 7, revealing ASE for most class I and II genes.ConclusionsThis study represents, to the best of our knowledge, the first genome-wide map of the genetic control of gene expression in porcine peripheral blood. These results represent an interesting resource for the identification of genetic markers and blood biomarkers associated with variations in immunity traits in pigs, as well as any other complex traits for which blood is an appropriate surrogate tissue.

Integration of liver gene co-expression networks and eGWAs analyses highlighted candidate regulators implicated in lipid metabolism in pigs

In the present study, liver co-expression networks and expression Genome Wide Association Study (eGWAS) were performed to identify DNA variants and molecular pathways implicated in the functional regulatory mechanisms of meat quality traits in pigs. With this purpose, the liver mRNA expression of 44 candidates genes related with lipid metabolism was analysed in 111 Iberian x Landrace backcross animals. The eGWAS identified 92 eSNPs located in seven chromosomal regions and associated with eight genes: CROT, CYP2U1, DGAT1, EGF, FABP1, FABP5, PLA2G12A, and PPARA. Remarkably, cis-eSNPs associated with FABP1 gene expression which may be determining the C18:2(n-6)/C18:3(n-3) ratio in backfat through the multiple interaction of DNA variants and genes were identified. Furthermore, a hotspot on SSC8 associated with the gene expression of eight genes was identified and the TBCK gene was pointed out as candidate gene regulating it. Our results also suggested that the PI3K-Akt-mTOR pathway plays an important role in the control of the analysed genes highlighting nuclear receptors as the NR3C1 or PPARA. Finally, sex-dimorphism associated with hepatic lipid metabolism was identified with over-representation of female-biased genes. These results increase our knowledge of the genetic architecture underlying fat composition traits.

Expression-based GWAS identifies variants, gene interactions and key regulators affecting intramuscular fatty acid content and composition in porcine meat.

The aim of this work is to better understand the genetic mechanisms determining two complex traits affecting porcine meat quality: intramuscular fat (IMF) content and its fatty acid (FA) composition. With this purpose, expression Genome-Wide Association Study (eGWAS) of 45 lipid-related genes associated with meat quality traits in swine muscle (Longissimus dorsi) of 114 Iberian × Landrace backcross animals was performed. The eGWAS identified 241 SNPs associated with 11 genes: ACSM5, CROT, FABP3, FOS, HIF1AN, IGF2, MGLL, NCOA1, PIK3R1, PLA2G12A and PPARA. Three expression Quantitative Trait Loci (eQTLs) for IGF2, ACSM5 and MGLL were identified, showing cis-acting effects, whereas 16 eQTLs had trans regulatory effects. A polymorphism in the ACSM5 promoter region associated with its expression was identified. In addition, strong candidate genes regulating ACSM5, FOS, PPARA, PIK3R1, PLA2G12A and HIF1AN gene expression were also seen. Notably, the analysis highlighted the NR3C1 transcription factor as a strong candidate gene involved in the regulation of the 45 genes analysed. Finally, the IGF2, MGLL, MC2R, ARHGAP6, and NR3C1 genes were identified as potential regulators co-localizing within QTLs for fatness and growth traits in the IBMAP population. The results obtained increase our knowledge in the functional regulatory mechanisms involved in these complex traits.

Combined meta-analysis of systemic effects of allogeneic stem cell transplantation and systemic sclerosis

BackgroundChronic graft-versus-host disease (cGVHD) is a major factor of morbidity and mortality for allogeneic stem cell transplantation (aSCT). The skin and internal organ involvement is the most common systemic complication of cGVHD and closely resembles systemic sclerosis (SSc). Circulating lymphocytes characterize the autoimmune nature of both conditions. Therefore we hypothesized that the common clinical manifestation (systemic organ and skin injury) and the common underlying players (lymphocytes) justify the combined meta-analysis of these diseases.ResultsThe aSCT and SSc datasets were uploaded from Gene Expression Omnibus (GEO), a public functional genomics data repository. The available microarray studies of peripheral blood mononuclear cells (PBMCs) and isolated lymphocytes were limited to well established microarray platforms (Affymetrix, Agilent, Canvac, and Illumina) and experimental settings with ≥10 patients per group. The resulting pools of data were merged by unique gene identifier and analyzed by the expression genome-wide association studies (eGWAS) coupled with the subtraction of the cGVHD+ and cGVHD− molecular signatures. The eGWAS was applied to 47 and 50 lymphocyte profiles from aSCT and SSc patients, respectively. The identified 35 candidates were represented by 8 known cGVHD genes (including CXCR4, LTBR and PML) and 28 new candidate genes (including SEPX1 and DNJGB1). The further mutual subtraction of cGVHD+ and cGVHD− candidates and pathway analysis identified a list of 25 genes. Seven of these genes belong to the fibroblast development and function pathway, consisting of the well known cGVHD genes CCND1, JUN, and FOS, and the new molecular targets MMP2, FOSB, TNFAIP8, and DUSP1. These genes become primary candidates for a potential link of systemic effects of cGVHD and SSc.ConclusionsWe designed a new approach for meta-analysis by combining data from different diseases using common clinical manifestation as a linker. This allowed us to power up the insufficient standalone meta-analysis of aSCT microarray studies, by adding SSc samples to the data pool. This new method has successfully identified novel molecular targets for systemic effects of both aSCT and SSc. We believe that this approach is generalizable and can be applied to an array of diseases with common clinical manifestations.

Brain expression genome-wide association study (eGWAS) identifies human disease-associated variants.

Genetic variants that modify brain gene expression may also influence risk for human diseases. We measured expression levels of 24,526 transcripts in brain samples from the cerebellum and temporal cortex of autopsied subjects with Alzheimer's disease (AD, cerebellar n = 197, temporal cortex n = 202) and with other brain pathologies (non–AD, cerebellar n = 177, temporal cortex n = 197). We conducted an expression genome-wide association study (eGWAS) using 213,528 cisSNPs within ±100 kb of the tested transcripts. We identified 2,980 cerebellar cisSNP/transcript level associations (2,596 unique cisSNPs) significant in both ADs and non–ADs (q<0.05, p = 7.70×10−5–1.67×10−82). Of these, 2,089 were also significant in the temporal cortex (p = 1.85×10−5–1.70×10−141). The top cerebellar cisSNPs had 2.4-fold enrichment for human disease-associated variants (p<10−6). We identified novel cisSNP/transcript associations for human disease-associated variants, including progressive supranuclear palsy SLCO1A2/rs11568563, Parkinson's disease (PD) MMRN1/rs6532197, Paget's disease OPTN/rs1561570; and we confirmed others, including PD MAPT/rs242557, systemic lupus erythematosus and ulcerative colitis IRF5/rs4728142, and type 1 diabetes mellitus RPS26/rs1701704. In our eGWAS, there was 2.9–3.3 fold enrichment (p<10−6) of significant cisSNPs with suggestive AD–risk association (p<10−3) in the Alzheimer's Disease Genetics Consortium GWAS. These results demonstrate the significant contributions of genetic factors to human brain gene expression, which are reliably detected across different brain regions and pathologies. The significant enrichment of brain cisSNPs among disease-associated variants advocates gene expression changes as a mechanism for many central nervous system (CNS) and non–CNS diseases. Combined assessment of expression and disease GWAS may provide complementary information in discovery of human disease variants with functional implications. Our findings have implications for the design and interpretation of eGWAS in general and the use of brain expression quantitative trait loci in the study of human disease genetics.