Expression quantitative trait loci are highly sensitive to cellular differentiation state

Blood cell development from multipotent hematopoietic stem cells to specialized blood cells is accompanied by drastic changes in gene expression for which the triggers remain mostly unknown. Genetical genomics is an approach linking natural genetic variation to gene expression variation, thereby allowing the identification of genomic loci containing gene expression modulators (eQTLs). In this paper, we used a genetical genomics approach to analyze gene expression across four developmentally close blood cell types collected from a large number of genetically different but related mouse strains. We found that, while a significant number of eQTLs (365) had a consistent “static” regulatory effect on gene expression, an even larger number were found to be very sensitive to cell stage. As many as 1,283 eQTLs exhibited a “dynamic” behavior across cell types. By looking more closely at these dynamic eQTLs, we show that the sensitivity of eQTLs to cell stage is largely associated with gene expression changes in target genes. These results stress the importance of studying gene expression variation in well-defined cell populations. Only such studies will be able to reveal the important differences in gene regulation between different ce

Genetical genomics: spotlight on QTL hotspots

No abstract available

A verification protocol for the probe sequences of Affymetrix genome arrays reveals high probe accuracy for studies in mouse, human and rat

BACKGROUND: The Affymetrix GeneChip technology uses multiple probes per gene to measure its expression level. Individual probe signals can vary widely, which hampers proper interpretation. This variation can be caused by probes that do not properly match their target gene or that match multiple genes. To determine the accuracy of Affymetrix arrays, we developed an extensive verification protocol, for mouse arrays incorporating the NCBI RefSeq, NCBI UniGene Unique, NIA Mouse Gene Index, and UCSC mouse genome databases. RESULTS: Applying this protocol to Affymetrix Mouse Genome arrays (the earlier U74Av2 and the newer 430 2.0 array), the number of sequence-verified probes with perfect matches was no less than 85% and 95%, respectively; and for 74% and 85% of the probe sets all probes were sequence verified. The latter percentages increased to 80% and 94% after discarding one or two unverifiable probes per probe set, and even further to 84% and 97% when, in addition, allowing for one or two mismatches between probe and target gene. Similar results were obtained for other mouse arrays, as well as for human and rat arrays. Based on these data, refined chip definition files for all arrays are provided online. Researchers can choose the version appropriate for their study to (re)analyze expression data. CONCLUSION: The accuracy of Affymetrix probe sequences is higher than previously reported, particularly on newer arrays. Yet, refined probe set definitions have clear effects on the detection of differentially expressed genes. We demonstrate that the interpretation of the results of Affymetrix arrays is improved when the new chip definition files are used

Generalized DNA Barcode Design Based on Hamming Codes

The diversity and scope of multiplex parallel sequencing applications is steadily increasing. Critically, multiplex parallel sequencing applications methods rely on the use of barcoded primers for sample identification, and the quality of the barcodes directly impacts the quality of the resulting sequence data. Inspection of the recent publications reveals a surprisingly variable quality of the barcodes employed. Some barcodes are made in a semi empirical fashion, without quantitative consideration of error correction or minimal distance properties. After systematic comparison of published barcode sets, including commercially distributed barcoded primers from Illumina and Epicentre, methods for improved, Hamming code-based sequences are suggested and illustrated. Hamming barcodes can be employed for DNA tag designs in many different ways while preserving minimal distance and error-correcting properties. In addition, Hamming barcodes remain flexible with regard to essential biological parameters such as sequence redundancy and GC content. Wider adoption of improved Hamming barcodes is encouraged in multiplex parallel sequencing applications

Intra- and inter-individual genetic differences in gene expression

Genetic variation is known to influence the amount of mRNA produced by a gene. Given that the molecular machines control mRNA levels of multiple genes, we expect genetic variation in the components of these machines would influence multiple genes in a similar fashion. In this study we show that this assumption is correct by using correlation of mRNA levels measured independently in the brain, kidney or liver of multiple, genetically typed, mice strains to detect shared genetic influences. These correlating groups of genes (CGG) have collective properties that account for 40-90% of the variability of their constituent genes and in some cases, but not all, contain genes encoding functionally related proteins. Critically, we show that the genetic influences are essentially tissue specific and consequently the same genetic variations in the one animal may up-regulate a CGG in one tissue but down-regulate the same CGG in a second tissue. We further show similarly paradoxical behaviour of CGGs within the same tissues of different individuals. The implication of this study is that this class of genetic variation can result in complex inter- and intra-individual and tissue differences and that this will create substantial challenges to the investigation of phenotypic outcomes, particularly in humans where multiple tissues are not readily available.

&#xa

Gene Set Enrichment in eQTL Data Identifies Novel Annotations and Pathway Regulators

Genome-wide gene expression profiling has been extensively used to generate biological hypotheses based on differential expression. Recently, many studies have used microarrays to measure gene expression levels across genetic mapping populations. These gene expression phenotypes have been used for genome-wide association analyses, an analysis referred to as expression QTL (eQTL) mapping. Here, eQTL analysis was performed in adipose tissue from 28 inbred strains of mice. We focused our analysis on “trans-eQTL bands”, defined as instances in which the expression patterns of many genes were all associated to a common genetic locus. Genes comprising trans-eQTL bands were screened for enrichments in functional gene sets representing known biological pathways, and genes located at associated trans-eQTL band loci were considered candidate transcriptional modulators. We demonstrate that these patterns were enriched for previously characterized relationships between known upstream transcriptional regulators and their downstream target genes. Moreover, we used this strategy to identify both novel regulators and novel members of known pathways. Finally, based on a putative regulatory relationship identified in our analysis, we identified and validated a previously uncharacterized role for cyclin H in the regulation of oxidative phosphorylation. We believe that the specific molecular hypotheses generated in this study will reveal many additional pathway members and regulators, and that the analysis approaches described herein will be broadly applicable to other eQTL data sets

Identification of expression QTL (eQTL) of genes expressed in porcine M. longissimus dorsi and associated with meat quality traits

<p>Abstract</p> <p>Background</p> <p>Genetic analysis of transcriptional profiles is a promising approach for identifying and dissecting the genetics of complex traits like meat performance. Accordingly, expression levels obtained by microarray analysis were taken as phenotypes in a linkage analysis to map eQTL. Moreover, expression levels were correlated with traits related to meat quality and principle components with high loadings of these traits. By using an up-to-date annotation and localization of the respective probe-sets, the integration of eQTL mapping data and information of trait correlated expression finally served to point to candidate genes for meat quality traits.</p> <p>Results</p> <p>Genome-wide transcriptional profiles of <it>M. longissimus dorsi </it>RNAs samples of 74 F2 animals of a pig resource population revealed 11,457 probe-sets representing genes expressed in the muscle. Linkage analysis of expression levels of these probe-sets provided 9,180 eQTL at the suggestive significance threshold of LOD > 2. We mapped 653 eQTL on the same chromosome as the corresponding gene and these were designated as 'putative <it>cis-</it>eQTL'. In order to link eQTL to the traits of interest, probe-sets were addressed with relative transcript abundances that showed correlation with meat quality traits at p ≤ 0.05. Out of the 653 'putative <it>cis-</it>eQTL', 262 transcripts were correlated with at least one meat quality trait. Furthermore, association of expression levels with composite traits with high loadings for meat quality traits generated by principle component analysis were taken into account leading to a list of 85 genes exhibiting <it>cis-</it>eQTL and trait dependent expression.</p> <p>Conclusion</p> <p>Holistic expression profiling was integrated with QTL analysis for meat quality traits. Correlations between transcript abundance and meat quality traits, combined with genetic positional information of eQTL allowed us to prioritise candidate genes for further study.</p

Combining transcriptional profiling and genetic linkage analysis to uncover gene networks operating in hematopoietic stem cells and their progeny

Stem cells are unique in that they possess both the capacity to self-renew and thereby maintain their original pool as well as the capacity to differentiate into mature cells. In the past number of years, transcriptional profiling of enriched stem cell populations has been extensively performed in an attempt to identify a universal stem cell gene expression signature. While stem-cell-specific transcripts were identified in each case, this approach has thus far been insufficient to identify a universal group of core “stemness” genes ultimately responsible for self-renewal and multipotency. Similarly, in the hematopoietic system, comparisons of transcriptional profiles between different hematopoietic cell stages have had limited success in revealing core genes ultimately responsible for the initiation of differentiation and lineage specification. Here, we propose that the combined use of transcriptional profiling and genetic linkage analysis, an approach called “genetical genomics”, can be a valuable tool to assist in the identification of genes and gene networks that specify “stemness” and cell fate decisions. We review past studies of hematopoietic cells that utilized transcriptional profiling and/or genetic linkage analysis, and discuss several potential future applications of genetical genomics

Expression profiling identifies novel candidate genes for ethanol sensitivity QTLs

The Inbred Long Sleep (ILS) and Inbred Short Sleep (ISS) mouse strains have a 16-fold difference in duration of loss of the righting response (LORR) following administration of a sedative dose of ethanol. Four quantitative trait loci (QTLs) have been mapped in these strains for this trait. Underlying each of these QTLs must be one or more genetic differences (polymorphisms in either gene coding or regulatory regions) influencing ethanol sensitivity. Because prior studies have tended to focus on differences in coding regions, genome-wide expression profiling in cerebellum was used here to identify candidate genes for regulatory region differences in these two strains. Fifteen differentially expressed genes were found that map to the QTL regions and polymorphisms were identified in the promoter regions of four of these genes by direct sequencing of ILS and ISS genomic DNA. Polymorphisms in the promoters of three of these genes, Slc22a4, Rassf2, and Tax1bp3, disrupt putative transcription factor binding sites. Slc22a4 and another candidate, Xrcc5, have human orthologs that map to genomic regions associated with human ethanol sensitivity in genetic linkage studies. These genes represent novel candidates for the LORR phenotype and provide new targets for future studies into the neuronal processes underlying ethanol sensitivity

Gemcitabine and Arabinosylcytosin Pharmacogenomics: Genome-Wide Association and Drug Response Biomarkers

Cancer patients show large individual variation in their response to chemotherapeutic agents. Gemcitabine (dFdC) and AraC, two cytidine analogues, have shown significant activity against a variety of tumors. We previously used expression data from a lymphoblastoid cell line-based model system to identify genes that might be important for the two drug cytotoxicity. In the present study, we used that same model system to perform a genome-wide association (GWA) study to test the hypothesis that common genetic variation might influence both gene expression and response to the two drugs. Specifically, genome-wide single nucleotide polymorphisms (SNPs) and mRNA expression data were obtained using the Illumina 550K® HumanHap550 SNP Chip and Affymetrix U133 Plus 2.0 GeneChip, respectively, for 174 ethnically-defined “Human Variation Panel” lymphoblastoid cell lines. Gemcitabine and AraC cytotoxicity assays were performed to obtain IC50 values for the cell lines. We then performed GWA studies with SNPs, gene expression and IC50 of these two drugs. This approach identified SNPs that were associated with gemcitabine or AraC IC50 values and with the expression regulation for 29 genes or 30 genes, respectively. One SNP in IQGAP2 (rs3797418) was significantly associated with variation in both the expression of multiple genes and gemcitabine and AraC IC50. A second SNP in TGM3 (rs6082527) was also significantly associated with multiple gene expression and gemcitabine IC50. To confirm the association results, we performed siRNA knock down of selected genes with expression that was associated with rs3797418 and rs6082527 in tumor cell and the knock down altered gemcitabine or AraC sensitivity, confirming our association study results. These results suggest that the application of GWA approaches using cell-based model systems, when combined with complementary functional validation, can provide insights into mechanisms responsible for variation in cytidine analogue response