Combining microarrays and genetic analysis

An integrated genomics approach to identify genetic regions associated with neonatal growth trait in mice

Neonatal growth during the early post-partum period is closely associated with lactation performance. Neonatal growth reflects milk output and is a complex variable trait among inbred mouse strains, but few studies have compared this trait systematically across more than a few strains. In the present study, 11 inbred strains of mice were measured for a neonatal growth phenotype during the first eight days of lactation. Significant differences in neonatal growth trait were observed with QSi5 (3.71±0.05 g) and DBA/1J (2.67±0.06 g) strains defining the two extremes of the phenotype. In silico association analysis was performed for trait variability using the high density SNP information on inbred strains of mice. We found strong evidence to refine a previously identified large neonatal growth QTL on mouse chromosome 9, Neogq1. When an integrated strategy that combined fine mapping and analysis of mammary transcriptome expression profiles of lactating mice with divergent phenotypes was applied, we identified neogenin (Neo1), a gene important for mammary gland morphogenesis, as a likely quantitative trait gene (QTG) underlying the Neogq1 QTL in mice.

Expression QTL mapping in regulatory and helper T cells from the BXD family of strains reveals novel cell-specific genes, gene-gene interactions and candidate genes for auto-immune disease

BACKGROUND

Regulatory T cells (Tregs) play an essential role in the control of the immune response. Treg cells represent important targets for therapeutic interventions of the immune system. Therefore, it will be very important to understand in more detail which genes are specifically activated in Treg cells versus T helper (Th) cells, and which gene regulatory circuits may be involved in specifying and maintaining Treg cell homeostasis.

RESULTS

We isolated Treg and Th cells from a genetically diverse family of 31 BXD type recombinant inbred strains and the fully inbred parental strains of this family--C57BL/6J and DBA/2J. Subsequently genome-wide gene expression studies were performed from the isolated Treg and Th cells. A comparative analysis of the transcriptomes of these cell populations allowed us to identify many novel differentially expressed genes. Analysis of cis- and trans-expression Quantitative Trait Loci (eQTLs) highlighted common and unique regulatory mechanisms that are active in the two cell types. Trans-eQTL regions were found for the Treg functional genes Nrp1, Stat3 and Ikzf4. Analyses of the respective QTL intervals suggested several candidate genes that may be involved in regulating these genes in Treg cells. Similarly, possible candidate genes were found which may regulate the expression of F2rl1, Ctla4, Klrb1f. In addition, we identified a focused group of candidate genes that may be important for the maintenance of self-tolerance and the prevention of allergy.

CONCLUSIONS

Variation of expression across the strains allowed us to find many novel gene-interaction networks in both T cell subsets. In addition, these two data sets enabled us to identify many differentially expressed genes and to nominate candidate genes that may have important functions for the maintenance of self-tolerance and the prevention of allergy.

Genome-wide analysis of the mouse lung transcriptome reveals novel molecular gene interaction networks and cell-specific expression signatures

BACKGROUND

The lung is critical in surveillance and initial defense against pathogens. In humans, as in mice, individual genetic differences strongly modulate pulmonary responses to infectious agents, severity of lung disease, and potential allergic reactions. In a first step towards understanding genetic predisposition and pulmonary molecular networks that underlie individual differences in disease vulnerability, we performed a global analysis of normative lung gene expression levels in inbred mouse strains and a large family of BXD strains that are widely used for systems genetics. Our goal is to provide a key community resource on the genetics of the normative lung transcriptome that can serve as a foundation for experimental analysis and allow predicting genetic predisposition and response to pathogens, allergens, and xenobiotics.

METHODS

Steady-state polyA+ mRNA levels were assayed across a diverse and fully genotyped panel of 57 isogenic strains using the Affymetrix M430 2.0 array. Correlations of expression levels between genes were determined. Global expression QTL (eQTL) analysis and network covariance analysis was performed using tools and resources in GeneNetwork http://www.genenetwork.org.

RESULTS

Expression values were highly variable across strains and in many cases exhibited a high heritability factor. Several genes which showed a restricted expression to lung tissue were identified. Using correlations between gene expression values across all strains, we defined and extended memberships of several important molecular networks in the lung. Furthermore, we were able to extract signatures of immune cell subpopulations and characterize co-variation and shared genetic modulation. Known QTL regions for respiratory infection susceptibility were investigated and several cis-eQTL genes were identified. Numerous cis- and trans-regulated transcripts and chromosomal intervals with strong regulatory activity were mapped. The Cyp1a1 P450 transcript had a strong trans-acting eQTL (LOD 11.8) on Chr 12 at 36 ± 1 Mb. This interval contains the transcription factor Ahr that has a critical mis-sense allele in the DBA/2J haplotype and evidently modulates transcriptional activation by AhR.

CONCLUSIONS

Large-scale gene expression analyses in genetic reference populations revealed lung-specific and immune-cell gene expression profiles and suggested specific gene regulatory interactions.

Proteomic analysis of time-dependent difference of protein expression profile changes during neuronal differentiation of mouse embryonic stem cells

The study of ES cell-mediated neuronal differentiation allows elucidating the mechanism of neuronal development in spite of the complexity and the difficult accessibility. During the differentiation of embryonic stem cells into neuronal cell, the expression profiles in the level of protein were extensively investigated by proteomic analysis. These cells were analyzed for charges in proteome during the differentiation of ES cells by 2-dimensional electrophoresis (2-DE) and MALDI-TOF MS. Seven unique proteins were identified, some of which were differentially expressed at each stage. A complex system of neuronal differentiation can be activated in cultured embryonic stem cells and our two dimensional electrophoresis data should be useful for investigating some of the mechanism that regulates neuronal differentiation.

Complex nature of SNP genotype effects on gene expression in primary human leucocytes

Background

Genome wide association studies have been hugely successful in identifying disease risk variants, yet most variants do not lead to coding changes and how variants influence biological function is usually unknown.

Methods

We correlated gene expression and genetic variation in untouched primary leucocytes (n = 110) from individuals with celiac disease – a common condition with multiple risk variants identified. We compared our observations with an EBV-transformed HapMap B cell line dataset (n = 90), and performed a meta-analysis to increase power to detect non-tissue specific effects.

Results

In celiac peripheral blood, 2,315 SNP variants influenced gene expression at 765 different transcripts (< 250 kb from SNP, at FDR = 0.05, cis expression quantitative trait loci, eQTLs). 135 of the detected SNP-probe effects (reflecting 51 unique probes) were also detected in a HapMap B cell line published dataset, all with effects in the same allelic direction. Overall gene expression differences within the two datasets predominantly explain the limited overlap in observed cis-eQTLs. Celiac associated risk variants from two regions, containing genes IL18RAP and CCR3, showed significant cis genotype-expression correlations in the peripheral blood but not in the B cell line datasets. We identified 14 genes where a SNP affected the expression of different probes within the same gene, but in opposite allelic directions. By incorporating genetic variation in co-expression analyses, functional relationships between genes can be more significantly detected.

Conclusion

In conclusion, the complex nature of genotypic effects in human populations makes the use of a relevant tissue, large datasets, and analysis of different exons essential to enable the identification of the function for many genetic risk variants in common diseases.

Beyond standardization: dynamic software infrastructures for systems biology

Progress in systems biology is seriously hindered by slow production of suitable software infrastructures. Biologists need infrastructure that easily connects to work that is done in other laboratories, for which standardization is helpful. However, the infrastructure must also accommodate the specifics of their biological system, but appropriate mechanisms to support variation are currently lacking. We argue that a minimal computer language, and a software tool called a generator, can be used to quickly produce customized software infrastructures that 'systems biologists really want to have'.

Review of microarray experimental design strategies for genetical genomics studies

Genetical genomics approaches provide a powerful tool for studying the genetic mechanisms governing variation in complex traits. By combining information on phenotypic traits, pedigree structure, molecular markers, and gene expression, such studies can be used for estimating heritability of mRNA transcript abundances, for mapping expression quantitative trait loci (eQTL), and for inferring regulatory gene networks. Microarray experiments, however, can be extremely costly and time consuming, which may limit sample sizes and statistical power. Thus it is crucial to optimize experimental designs by carefully choosing the subjects to be assayed, within a selective profiling approach, and by cautiously controlling systematic factors affecting the system. Also, a rigorous strategy should be used for allocating mRNA samples across assay batches, slides, and dye labeling, so that effects of interest are not confounded with nuisance factors. In this presentation, we review some selective profiling strategies for genetical genomics studies, including the selection of individuals for increased genetic dissimilarity and for a higher number of recombination events. Efficient designs for studying epistasis are also discussed, as well as experiments for inferring heritability of transcriptional levels. It is shown that solving an optimal design problem generally requires a numerical implementation and that the optimality criteria should be intimately related to the goals of the experiment, such as the estimation of additive, dominance, and interacting effects, localizing putative eQTL, or inferring genetic and environmental variance components associated with transcriptional abundances.

GEDA: new knowledge base of gene expression in drug addiction

Abuse of drugs can elicit compulsive drug seeking behaviors upon repeated administration, and ultimately leads to the phenomenon of addiction. We developed a procedure for the standardization of microarray gene expression data of rat brain in drug addiction and stored them in a single integrated database system, focusing on more effective data processing and interpretation. Another characteristic of the present database is that it has a systematic flexibility for statistical analysis and linking with other databases. Basically, we adopt an intelligent SQL querying system, as the foundation of our DB, in order to set up an interactive module which can automatically read the raw gene expression data in the standardized format. We maximize the usability of this DB, helping users study significant gene expression and identify biological function of the genes through integrated up-to-date gene information such as GO annotation and metabolic pathway. For collecting the latest information of selected gene from the database, we also set up the local BLAST search engine and nonredundant sequence database updated by NCBI server on a daily basis. We find that the present database is a useful query interface and data-mining tool, specifically for finding out the genes related to drug addiction. We apply this system to the identification and characterization of methamphetamine-induced genes' behavior in rat brain.