Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transcriptomics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

Transcriptomic Analysis of ISS Crewmembers’ Peripheral Blood Mononuclear Cells Reveals Homeostatic Regulations in Space

The impact of spaceflight on the immune system has been investigated for decades. Studies conducted in cell models, animals and humans suggest that the spaceflight environment affects the innate and acquired immune systems, as the ability to recognize antigens, defend against foreign invaders, and orchestrate repair is significantly hindered. However, the molecular mechanisms behind spaceflight-induced immune dysregulations are still unclear. In this study, blood from eleven (11) International Space Station (ISS) crewmembers was collected before, during and after long duration space missions, as well as from 11 matched ground control subjects. Transcriptomic analysis was performed in isolated peripheral blood mononuclear cells (PBMCs) using the RNA-sequencing technique. In comparison to the blood samples collected from the crewmembers pre-flight, a total of ~1000 genes were found to be upregulated and ~1000 genes downregulated in PBMC collected between 4 and 6 months after they were in space. The most significantly DEGs (differentially expressed genes) include activation of RUBCNL which is an autophagy enhancer and inhibition of GRASP which regulates cell trafficking. Genes involved in cell adhesion, cell cycle progression and other functions were also dysregulated. Pathway analysis of the DEGs indicates mitochondria dysfunction, particularly reduced ATP production in the electron transport chain. Other pathways impacted by spaceflight include glycolysis, autophagy and inflammatory response. Our results suggest that, in space, blood cells may have also experienced energy depletion and reduced metabolism. Consequently, the cells may become autophagic, which is a known homeostatic mechanism for blood cells to become quiescent, but to stay alive. Further analysis of the data shows recovery of the crewmembers after mission and potential differential responses between genders to the space environment. Our data potentially explains some of the physiological changes that have been observed in space such as mitochondria dysfunction, inhibition of T cell activation and telomere lengthening. Comparison of our results with other transcriptomics studies of ISS crewmembers’ blood cells will also be presented.

Maria Moreno Villanueva↗

Single Cell Genomics and Transcriptomics for Unicellular Eukaryotes

Despite their small size, unicellular eukaryotes have complex genomes with a high degree of plasticity that allow them to adapt quickly to environmental changes. Unicellular eukaryotes live with prokaryotes and higher eukaryotes, frequently in symbiotic or parasitic niches. To this day their contribution to the dynamics of the environmental communities remains to be understood. Unfortunately, the vast majority of eukaryotic microorganisms are either uncultured or unculturable, making genome sequencing impossible using traditional approaches. We have developed an approach to isolate unicellular eukaryotes of interest from environmental samples, and to sequence and analyze their genomes and transcriptomes. We have tested our methods with six species: an uncharacterized protist from cellulose-enriched compost identified as Platyophrya, a close relative of P. vorax; the fungus Metschnikowia bicuspidate, a parasite of water flea Daphnia; the mycoparasitic fungi Piptocephalis cylindrospora, a parasite of Cokeromyces and Mucor; Caulochytrium protosteloides, a parasite of Sordaria; Rozella allomycis, a parasite of the water mold Allomyces; and the microalgae Chlamydomonas reinhardtii. Here, we present the four components of our approach: pre-sequencing methods, sequence analysis for single cell genome assembly, sequence analysis of single cell transcriptomes, and genome annotation. This technology has the potential to uncover the complexity of single cell eukaryotes and their role in the environmental samples.

59 BASIC BIOLOGICAL SCIENCES↗

A Study on Contrastive Graph Neural Network Pretraining for Predicting Transcriptome Profiles

We study graph neural network learning for transcriptomics with limited amount of labeled data. Our study reveals that simple GNN architectures perform well and do not suffer from over-fitting as the more sophisticated ones. Our study shows that although contrastive learning as a pretraining strategy has been successful in predicting properties such as formation and binding energy, it is not effective for transcriptomics.

Ma, Jiaji [University of Virginia]↗

Quantitative trait locus (QTL) mapping and transcriptome profiling identify QTLs and candidate genes associated with heat stress response during reproductive development in Camelina sativa

Camelina sativa (L.) Crantz is a low-input oilseed crop that has great potential in providing sustainable feedstock for biofuels and bioproducts. Climate change is threatening production of camelina with rising global temperatures. Elucidating the genetic response to high temperatures is essential for successful breeding of heat-tolerant camelina varieties. Here, we report a combinatorial approach to identifying candidate genes associated with heat stress by quantitative trait locus (QTL) mapping and comparative transcriptome profiling. A population of recombinant inbred lines (RILs) was grown in a controlled growth chamber under the high-temperature regimes for 14 days beginning at the onset of the reproductive stage. Several traits related to seed production were evaluated at maturity. The QTL analysis identified several regions with co-located traits on chromosomes 8, 10, and 12. Two RILs with contrasting phenotypic responses to heat stress were chosen for gene expression profiling via RNA sequencing. Multiple pathways and genes were found to be strongly affected by heat stress, and many genes expressed differently between the two RILs. Several genes identified within the QTL regions were considered strong candidates that may control heat tolerance during reproduction in camelina. These studies provide resources for future studies that may assist in improving the heat tolerance of camelina.

60 APPLIED LIFE SCIENCES↗

Integrating transcriptomic and metabolomic analysis of the oleaginous yeast Rhodosporidium toruloides IFO0880 during growth under different carbon sources

Rhodosporidium toruloides is an oleaginous yeast capable of producing a variety of biofuels and bioproducts from diverse carbon sources. Despite numerous studies showing its promise as a platform microorganism, little is known about its metabolism and physiology. In this work, we investigated the central carbon metabolism in R. toruloides IFO0880 using transcriptomics and metabolomics during growth on glucose, xylose, acetate, or soybean oil. These substrates were chosen because they can be derived from plants. Significant changes in gene expression and metabolite concentrations were observed during growth on these four substrates. We mapped these changes onto the governing metabolic pathways to better understand how R. toruloides reprograms its metabolism to enable growth on these substrates. One notable finding concerns xylose metabolism, where poor expression of xylulokinase induces a bypass leading to arabitol production. Collectively, these results further our understanding of central carbon metabolism in R. toruloides during growth on different substrates. Furthermore, they may also help guide the metabolic engineering and development of better models of metabolism for R. toruloides.

59 BASIC BIOLOGICAL SCIENCES↗

Transcriptome-wide association analysis identifies candidate susceptibility genes for prostate-specific antigen levels in men without prostate cancer

Deciphering the genetic basis of prostate-specific antigen (PSA) levels may improve their utility for prostate cancer (PCa) screening. Using genome-wide association study (GWAS) summary statistics from 95,768 PCa-free men, we conducted a transcriptome-wide association study (TWAS) to examine impacts of genetically predicted gene expression on PSA. Analyses identified 41 statistically significant (p < 0.05/12,192 = 4.10 × 10 –6 ) associations in whole blood and 39 statistically significant (p < 0.05/13,844 = 3.61 × 10 –6 ) associations in prostate tissue, with 18 genes associated in both tissues. Cross-tissue analyses identified 155 statistically significantly (p < 0.05/22,249 = 2.25 × 10 –6 ) genes. Out of 173 unique PSA-associated genes across analyses, we replicated 151 (87.3%) in a TWAS of 209,318 PCa-free individuals from the Million Veteran Program. Based on conditional analyses, we found 20 genes (11 single tissue, nine cross-tissue) that were associated with PSA levels in the discovery TWAS that were not attributable to a lead variant from a GWAS. Ten of these 20 genes replicated, and two of the replicated genes had colocalization probability of >0.5: CCNA2 and HIST1H2BN. Six of the 20 identified genes are not known to impact PCa risk. Fine-mapping based on whole blood and prostate tissue revealed five protein-coding genes with evidence of causal relationships with PSA levels. Of these five genes, four exhibited evidence of colocalization and one was conditionally independent of previous GWAS findings. These results yield hypotheses that should be further explored to improve understanding of genetic factors underlying PSA levels.

60 APPLIED LIFE SCIENCES↗

Inference of Interactions in Cyanobacterial-Heterotrophic Co-Cultures via Transcriptome Sequencing

We employed deep sequencing technology to identify transcriptional adaptation of the euryhaline unicellular cyanobacterium Synechococcus sp. PCC 7002 and the marine facultative aerobe Shewanella putrefaciens W3-18-1 to growth in a co-culture and infer the effect of carbon flux distributions on photoautotroph-heterotroph interactions. The overall transcriptome response of both organisms to co-cultivation was shaped by their respective physiologies and growth constraints. Carbon limitation resulted in the expansion of metabolic capacities which was manifested through the transcriptional upregulation of transport and catabolic pathways. While growth coupling occurred via lactate oxidation or secretion of photosynthetically fixed carbon, there was evidence of specific metabolic interactions between the two organisms. On one hand, the production and excretion of specific amino acids (methionine and alanine) by the cyanobacterium correlated with the putative downregulation of the corresponding biosynthetic machinery of Shewanella W3-18-1. On the other hand, the broad and consistent decrease of mRNA levels for many Fe-regulated Synechococcus 7002 genes during co-cultivation suggested increased Fe availability as well as more facile and energy-efficient mechanisms for Fe acquisition by the cyanobacterium. Furthermore, evidence pointed at potentially novel interactions between oxygenic photoautotrophs and heterotrophs related to the oxidative stress response as transcriptional patterns suggested that Synechococcus 7002 rather than Shewanella W3-18-1 provided scavenging functions for reactive oxygen species under co-culture conditions. This study provides an initial insight into the complexity of photoautotrophic-heterotrophic interactions and brings new perspectives of their role in the robustness and stability of the association.

cyanobacterium↗

High-throughput single-cell transcriptomics of bacteria using combinatorial barcoding

Microbial split-pool ligation transcriptomics (microSPLiT) is a high-throughput single-cell RNA sequencing method for bacteria. With four combinatorial barcoding rounds, microSPLiT can profile transcriptional states in hundreds of thousands of Gram-negative and Gram-positive bacteria in a single experiment without specialized equipment. As bacterial samples are fixed and permeabilized before barcoding, they can be collected and stored ahead of time. During the first barcoding round, the fixed and permeabilized bacteria are distributed into a 96-well plate, where their transcripts are reverse transcribed into cDNA and labeled with the first well-specific barcode inside the cells. The cells are mixed and redistributed two more times into new 96-well plates, where the second and third barcodes are appended to the cDNA via in-cell ligation reactions. Finally, the cells are mixed and divided into aliquot sub-libraries, which can be stored until future use or prepared for sequencing with the addition of a fourth barcode. It takes 4 days to generate sequencing-ready libraries, including 1 day for collection and overnight fixation of samples. Here, the standard plate setup enables single-cell transcriptional profiling of up to 1 million bacterial cells and up to 96 samples in a single barcoding experiment, with the possibility of expansion by adding barcoding rounds. The protocol requires experience in basic molecular biology techniques, handling of bacterial samples and preparation of DNA libraries for next-generation sequencing. It can be performed by experienced undergraduate or graduate students. Data analysis requires access to computing resources, familiarity with Unix command line and basic experience with Python or R.

59 BASIC BIOLOGICAL SCIENCES↗

Dynamic Transcriptomic and Phosphoproteomic Analysis During Cell Wall Stress in Aspergillus nidulans

The fungal cell-wall integrity signaling (CWIS) pathway regulates cellular response to environmental stress to enable wall repair and resumption of normal growth. This complex, interconnected, pathway has been only partially characterized in filamentous fungi. To better understand the dynamic cellular response to wall perturbation, a β-glucan synthase inhibitor (micafungin) was added to a growing A. nidulans shake-flask culture. From this flask, transcriptomic and phosphoproteomic data were acquired over 10 and 120 min, respectively. To differentiate statistically-significant dynamic behavior from noise, a multivariate adaptive regression splines (MARS) model was applied to both data sets. Over 1800 genes were dynamically expressed and over 700 phosphorylation sites had changing phosphorylation levels upon micafungin exposure. Twelve kinases had altered phosphorylation and phenotypic profiling of all non-essential kinase deletion mutants revealed putative connections between PrkA, Hk-8–4, and Stk19 and the CWIS pathway. Our collective data implicate actin regulation, endocytosis, and septum formation as critical cellular processes responding to activation of the CWIS pathway, and connections between CWIS and calcium, HOG, and SIN signaling pathways.

59 BASIC BIOLOGICAL SCIENCES↗

Comparative Transcriptomics Analyses across Species, Organs, and Developmental Stages Reveal Functionally Constrained lncRNAs

The functionality of long noncoding RNAs (lncRNAs) is disputed. In general, lncRNAs are under weak selective pressures, suggesting that the majority of lncRNAs may be nonfunctional. However, although some surveys showed negligible phenotypic effects upon lncRNA perturbation, key biological roles were demonstrated for individual lncRNAs. Most lncRNAs with proven functions were implicated in gene expression regulation, in pathways related to cellular pluripotency, differentiation, and organ morphogenesis, suggesting that functional lncRNAs may be more abundant in embryonic development, rather than in adult organs. To test this hypothesis, we perform a multidimensional comparative transcriptomics analysis, across five developmental time points (two embryonic stages, newborn, adult, and aged individuals), four organs (brain, kidney, liver, and testes), and three species (mouse, rat, and chicken). We find that, overwhelmingly, lncRNAs are preferentially expressed in adult and aged testes, consistent with the presence of permissive transcription during spermatogenesis. LncRNAs are often differentially expressed among developmental stages and are less abundant in embryos and newborns compared with adult individuals, in agreement with a requirement for tighter expression control and less tolerance for noisy transcription early in development. For differentially expressed lncRNAs, we find that the patterns of expression variation among developmental stages are generally conserved between mouse and rat. Moreover, lncRNAs expressed above noise levels in somatic organs and during development show higher evolutionary conservation, in particular, at their promoter regions. Thus, we show that functionally constrained lncRNA loci are enriched in developing organs, and we suggest that many of these loci may function in an RNA-independent manner.

59 BASIC BIOLOGICAL SCIENCES↗

Transcriptomics reveal a mechanism of niche defense: two beneficial root endophytes deploy an antimicrobial GH18-CBM5 chitinase to protect their hosts

Effector secretion is crucial for root endophytes to establish and protect their ecological niche. We used time-resolved transcriptomics to monitor effector gene expression dynamics in two closely related Sebacinales, Serendipita indica and Serendipita vermifera, during symbiosis with three plant species, competition with the phytopathogenic fungus Bipolaris sorokiniana, and cooperation with root-associated bacteria. We observed increased effector gene expression in response to biotic interactions, particularly with plants, indicating their importance in host colonization. Some effectors responded to both plants and microbes, suggesting dual roles in intermicrobial competition and plant–microbe interactions. A subset of putative antimicrobial effectors, including a GH18-CBM5 chitinase, was induced exclusively by microbes. Functional analyses of this chitinase revealed its antimicrobial and plant-protective properties. We conclude that dynamic effector gene expression underpins the ability of Sebacinales to thrive in diverse ecological niches with a single fungal chitinase contributing substantially to niche defense.

59 BASIC BIOLOGICAL SCIENCES↗

Metabolic and transcriptomic study of pennycress natural variation identifies targets for oil improvement

Pennycress (Thlaspi arvense L.), a member of the Brassicaceae family, produces seed oil high in erucic acid, suitable for biodiesel and aviation fuel. Although pennycress, a winter annual, could be grown as a dedicated bioenergy crop, an increase in its seed oil content is required to improve its economic competitiveness. The success of crop improvement relies upon finding the right combination of biomarkers and targets, and the best genetic engineering and/or breeding strategies. In this work, we combined biomass composition with metabolomic and transcriptomic studies of developing embryos from 22 pennycress natural variants to identify targets for oil improvement. The selected accession collection presented diverse levels of fatty acids at maturity ranging from 29% to 41%. Pearson correlation analyses, weighted gene co‐expression network analysis and biomarker identifications were used as complementary approaches to detect associations between metabolite level or gene expression and oil content at maturity. The results indicated that improving seed oil content can lead to a concomitant increase in the proportion of erucic acid without affecting the weight of embryos. Processes, such as carbon partitioning towards the chloroplast, lipid metabolism, photosynthesis, and a tight control of nitrogen availability, were found to be key for oil improvement in pennycress. Besides identifying specific targets, our results also provide guidance regarding the best timing for their modification, early or middle maturation. Thus, this work lays out promising strategies, specific for pennycress, to accelerate the successful development of lines with increased seed oil content for biofuel applications.

59 BASIC BIOLOGICAL SCIENCES↗

Transcriptomic Analyses Elucidate Adaptive Differences of Closely Related Strains of Pseudomonas aeruginosa in Fuel

Pseudomonas aeruginosa can utilize hydrocarbons, but different strains have various degrees of adaptation despite their highly conserved genome. P. aeruginosa ATCC 33988 is highly adapted to hydrocarbons, while P. aeruginosa strain PAO1, a human pathogen, is less adapted and degrades jet fuel at a lower rate than does ATCC 33988. We investigated fuel-specific transcriptomic differences between these strains in order to ascertain the underlying mechanisms utilized by the adapted strain to proliferate in fuel. During growth in fuel, the genes related to alkane degradation, heat shock response, membrane proteins, efflux pumps, and several novel genes were upregulated in ATCC 33988. Overexpression of alk genes in PAO1 provided some improvement in growth, but it was not as robust as that of ATCC 33988, suggesting the role of other genes in adaptation. Expression of the function unknown gene PA5359 from ATCC 33988 in PAO1 increased the growth in fuel. Bioinformatic analysis revealed that PA5359 is a predicted lipoprotein with a conserved Yx(FWY)xxD motif, which is shared among bacterial adhesins. Overexpression of the putative resistance-nodulation-division (RND) efflux pump PA3521 to PA3523 increased the growth of the ATCC 33988 strain, suggesting a possible role in fuel tolerance. Interestingly, the PAO1 strain cannot utilize n -C 8 and n -C 10 . The expression of green fluorescent protein (GFP) under the control of alkB promoters confirmed that alk gene promoter polymorphism affects the expression of alk genes. Promoter fusion assays further confirmed that the regulation of alk genes was different in the two strains. Protein sequence analysis showed low amino acid differences for many of the upregulated genes, further supporting transcriptional control as the main mechanism for enhanced adaptation. Lastly, these results support that specific signal transduction, gene regulation, and coordination of multiple biological responses are required to improve the survival, growth, and metabolism of fuel in adapted strains. This study provides new insight into the mechanistic differences between strains and helpful information that may be applied in the improvement of bacterial strains for resistance to biotic and abiotic factors encountered during bioremediation and industrial biotechnological processes.

54 ENVIRONMENTAL SCIENCES↗

Transcriptomic data sets for Novosphingobium aromaticivorans DSM12444 and a ΔSARO_RS14285 mutant grown in the presence of glucose and either protocatechuic, vanillic, syringic, or 4-coumaric acid

The SARO_RS14285 gene, encoding a transcription factor, was deleted in Novosphingobium aromaticivorans DSM12444. The transcriptomes of the parent and ΔSARO_RS14285 strains were determined when grown in medium containing glucose with or without protocatechuic, vanillic, syringic, or 4-coumaric acid. We present the raw RNA sequencing data obtained from these cultures.

Novosphingobium aromaticivorans↗

De novo transcriptome in roots of switchgrass ( Panicum virgatum L. ) reveals gene expression dynamic and act network under alkaline salt stress

Background: Soil salinization is a major limiting factor for crop cultivation. Switchgrass is a perennial rhizomatous bunchgrass that is considered an ideal plant for marginal lands, including sites with saline soil. Here we investigated the physiological responses and transcriptome changes in the roots of Alamo (alkaline-tolerant genotype) and AM314/MS-155 (alkaline-sensitive genotype) under alkaline salt stress. Results: Alkaline salt stress significantly affected the membrane, osmotic adjustment and antioxidant systems in switchgrass roots, and the ASTTI values between Alamo and AM-314/MS-155 were divergent at different time points. A total of 108,319 unigenes were obtained after reassembly, including 73,636 unigenes in AM-314/MS-155 and 65,492 unigenes in Alamo. A total of 10,219 DEGs were identified, and the number of upregulated genes in Alamo was much greater than that in AM-314/MS-155 in both the early and late stages of alkaline salt stress. The DEGs in AM-314/MS-155 were mainly concentrated in the early stage, while Alamo showed greater advantages in the late stage. These DEGs were mainly enriched in plant-pathogen interactions, ubiquitin-mediated proteolysis and glycolysis/gluconeogenesis pathways. We characterized 1480 TF genes into 64 TF families, and the most abundant TF family was the C2H2 family, followed by the bZIP and bHLH families. A total of 1718 PKs were predicted, including CaMK, CDPK, MAPK and RLK. WGCNA revealed that the DEGs in the blue, brown, dark magenta and light steel blue 1 modules were associated with the physiological changes in roots of switchgrass under alkaline salt stress. The consistency between the qRT-PCR and RNA-Seq results confirmed the reliability of the RNA-seq sequencing data. A molecular regulatory network of the switchgrass response to alkaline salt stress was preliminarily constructed on the basis of transcriptional regulation and functional genes. Conclusions: Alkaline salt tolerance of switchgrass may be achieved by the regulation of ion homeostasis, transport proteins, detoxification, heat shock proteins, dehydration and sugar metabolism. These findings provide a comprehensive analysis of gene expression dynamic and act network induced by alkaline salt stress in two switchgrass genotypes and contribute to the understanding of the alkaline salt tolerance mechanism of switchgrass and the improvement of switchgrass germplasm.

59 BASIC BIOLOGICAL SCIENCES↗

Data for Integrating Transcriptomic and Metabolomic Analysis of the Oleaginous Yeast Rhodosporidium toruloides IFO0880 During Growth under Different Carbon Sources

Rhodosporidium toruloides is an oleaginous yeast capable of producing a variety of biofuels and bioproducts from diverse carbon sources. Despite numerous studies showing its promise as a platform microorganism, little is known about its metabolism and physiology. In this work, we investigated the central carbon metabolism in R. toruloides IFO0880 using transcriptomics and metabolomics during growth on glucose, xylose, acetate, or soybean oil. These substrates were chosen because they can be derived from plants. Significant changes in gene expression and metabolite concentrations were observed during growth on these four substrates. We mapped these changes onto the governing metabolic pathways to better understand how R. toruloides reprograms its metabolism to enable growth on these substrates. One notable finding concerns xylose metabolism, where poor expression of xylulokinase induces a bypass leading to arabitol production. Collectively, these results further our understanding of central carbon metabolism in R. toruloides during growth on different substrates. They may also help guide the metabolic engineering and development of better models of metabolism for R. toruloides .

Conversion↗

Genome assembly of the fungus Cochliobolus miyabeanus , and transcriptome analysis during early stages of infection on American wildrice ( Zizania palustris L.)

Here, the fungus Cochliobolus miyabeanus causes severe leaf spot disease on rice ( Oryza sativa ) and two North American specialty crops, American wildrice ( Zizania palustris ) and switchgrass ( Panicum virgatu ). Despite the importance of C. miyabeanus as a disease-causing agent in wildrice, little is known about either the mechanisms of pathogenicity or host defense responses. To start bridging these gaps, the genome of C. miyabeanus strain TG12bL2 was shotgun sequenced using Illumina technology. The genome assembly consists of 31.79 Mbp in 2,378 scaffolds with an N 50 = 74,921. It contains 11,000 predicted genes of which 94.5% were annotated. Approximately 10% of total gene number is expected to be secreted. The C. miyabeanus genome is rich in carbohydrate active enzymes, and harbors 187 small secreted peptides (SSPs) and some fungal effector homologs. Detoxification systems were represented by a variety of enzymes that could offer protection against plant defense compounds. The non-ribosomal peptide synthetases and polyketide synthases (PKS) present were common to other Cochliobolus species. Additionally, the fungal transcriptome was analyzed at 48 hours after inoculation in planta . A total of 10,674 genes were found to be expressed, some of which are known to be involved in pathogenicity or response to host defenses including hydrophobins, cutinase, cell wall degrading enzymes, enzymes related to reactive oxygen species scavenging, PKS, detoxification systems, SSPs, and a known fungal effector. This work will facilitate future research on C. miyabeanus pathogen-associated molecular patterns and effectors, and in the identification of their corresponding wildrice defense mechanisms.

59 BASIC BIOLOGICAL SCIENCES↗

­­­­Transcriptome-Wide Association Supplements Genome-Wide Association in Zea mays

Modern improvement of complex traits in agricultural species relies on successful associations of heritable molecular variation with observable phenotypes. Historically, this pursuit has primarily been based on easily measurable genetic markers. The recent advent of new technologies allows assaying and quantifying biological intermediates (hereafter endophenotypes) which are now readily measurable at a large scale across diverse individuals. The usefulness of endophenotypes for delineating the regulatory landscape of the genome and genetic dissection of complex trait variation remains underexplored in plants. The work presented here illustrated the utility of a large-scale (299-genotype and seven-tissue) gene expression resource to dissect traits across multiple levels of biological organization. Using single-tissue- and multi-tissue-based transcriptome-wide association studies (TWAS), we revealed that about half of the functional variation acts through altered transcript abundance for maize kernel traits, including 30 grain carotenoid abundance traits, 20 grain tocochromanol abundance traits, and 22 field-measured agronomic traits. Comparing the efficacy of TWAS with genomewide association studies (GWAS) and an ensemble approach that combines both GWAS and TWAS, we demonstrated that results of TWAS in combination with GWAS increase the power to detect known genes and aid in prioritizing likely causal genes. Using a variance partitioning approach in the largely independent maize Nested Association Mapping (NAM) population, we also showed that the most strongly associated genes identified by combining GWAS and TWAS explain more heritable variance for a majority of traits than the heritability captured by the random genes and the genes identified by GWAS or TWAS alone. As a result, this not only improves the ability to link genes to phenotypes, but also highlights the phenotypic consequences of regulatory variation in plants.

59 BASIC BIOLOGICAL SCIENCES↗