Integrating genetic and gene expression data: application to cardiovascular and metabolic traits in mice

The millions of common DNA variations that occur in the human population, or among inbred strains of mice and rats, perturb the expression (transcript levels) of a large fraction of the genes expressed in a particular tissue. The hundreds or thousands of common cis-acting variations that occur in the population may in turn affect the expression of thousands of other genes by affecting transcription factors, signaling molecules, RNA processing, and other processes that act in trans. The levels of transcripts are conveniently quantitated using expression arrays, and the cis- and trans-acting loci can be mapped using quantitative trait locus (QTL) analysis, in the same manner as loci for physiologic or clinical traits. Thousands of such expression QTL (eQTL) have been mapped in various crosses in mice, as well as other experimental organisms, and less detailed maps have been produced in studies of cells from human pedigrees. Such an integrative genetics approach (sometimes referred to as “genetical genomics”) is proving useful for identifying genes and pathways that contribute to complex clinical traits. The coincidence of clinical trait QTL and eQTL can help in the prioritization of positional candidate genes. More importantly, mathematical modeling of correlations between levels of transcripts and clinical traits in genetic crosses can allow prediction of causal interactions and the identification of “key driver” genes. An important objective of such studies will be to model biological networks in physiologic processes. When combined with high-density single nucleotide polymorphism (SNP) mapping, it should be feasible to identify genes that contribute to transcript levels using association analysis in outbred populations. In this review we discuss the basic concepts and applications of this integrative genomic approach to cardiovascular and metabolic diseases

Detection of Molecular Paths Associated with Insulitis and Type 1 Diabetes in Non-Obese Diabetic Mouse

Recent clinical evidence suggests important role of lipid and amino acid metabolism in early pre-autoimmune stages of type 1 diabetes pathogenesis. We study the molecular paths associated with the incidence of insulitis and type 1 diabetes in the Non-Obese Diabetic (NOD) mouse model using available gene expression data from the pancreatic tissue from young pre-diabetic mice. We apply a graph-theoretic approach by using a modified color coding algorithm to detect optimal molecular paths associated with specific phenotypes in an integrated biological network encompassing heterogeneous interaction data types. In agreement with our recent clinical findings, we identified a path downregulated in early insulitis involving dihydroxyacetone phosphate acyltransferase (DHAPAT), a key regulator of ether phospholipid synthesis. The pathway involving serine/threonine-protein phosphatase (PP2A), an upstream regulator of lipid metabolism and insulin secretion, was found upregulated in early insulitis. Our findings provide further evidence for an important role of lipid metabolism in early stages of type 1 diabetes pathogenesis, as well as suggest that such dysregulation of lipids and related increased oxidative stress can be tracked to beta cells

A Generalized Framework for Quantifying the Dynamics of EEG Event-Related Desynchronization

Brains were built by evolution to react swiftly to environmental challenges. Thus, sensory stimuli must be processed ad hoc, i.e., independent—to a large extent—from the momentary brain state incidentally prevailing during stimulus occurrence. Accordingly, computational neuroscience strives to model the robust processing of stimuli in the presence of dynamical cortical states. A pivotal feature of ongoing brain activity is the regional predominance of EEG eigenrhythms, such as the occipital alpha or the pericentral mu rhythm, both peaking spectrally at 10 Hz. Here, we establish a novel generalized concept to measure event-related desynchronization (ERD), which allows one to model neural oscillatory dynamics also in the presence of dynamical cortical states. Specifically, we demonstrate that a somatosensory stimulus causes a stereotypic sequence of first an ERD and then an ensuing amplitude overshoot (event-related synchronization), which at a dynamical cortical state becomes evident only if the natural relaxation dynamics of unperturbed EEG rhythms is utilized as reference dynamics. Moreover, this computational approach also encompasses the more general notion of a “conditional ERD,” through which candidate explanatory variables can be scrutinized with regard to their possible impact on a particular oscillatory dynamics under study. Thus, the generalized ERD represents a powerful novel analysis tool for extending our understanding of inter-trial variability of evoked responses and therefore the robust processing of environmental stimuli

Inhalation of β2 agonists impairs the clearance of nontypable Haemophilus influenzae from the murine respiratory tract

BACKGROUND: Nontypable Haemophilus influenzae (NTHi) is a common bacterial pathogen causing human respiratory tract infections under permissive conditions such as chronic obstructive pulmonary disease. Inhalation of β2-receptor agonists is a widely used treatment in patients with chronic obstructive pulmonary disease. The aim of this study was to determine the effect of inhalation of β2 agonists on the host immune response to respiratory tract infection with NTHi. METHODS: Mouse alveolar macrophages were stimulated in vitro with NTHi in the presence or absence of the β2 receptor agonists salmeterol or salbutamol. In addition, mice received salmeterol or salbutamol by inhalation and were intranasally infected with NTHi. End points were pulmonary inflammation and bacterial loads. RESULTS: Both salmeterol and salbutamol inhibited NTHi induced tumor necrosis factor-α (TNFα) release by mouse alveolar macrophages in vitro by a β receptor dependent mechanism. In line, inhalation of either salmeterol or salbutamol was associated with a reduced early TNFα production in lungs of mice infected intranasally with NTHi, an effect that was reversed by concurrent treatment with the β blocker propranolol. The clearance of NTHi from the lungs was impaired in mice treated with salmeterol or salbutamol, an adverse effect that was prevented by propranolol and independent of the reduction in TNFα. CONCLUSION: These data suggest that inhalation of salmeterol or salbutamol may negatively influence an effective clearance of NTHi from the airways

The malignant phenotype in breast cancer is driven by eIF4A1-mediated changes in the translational landscape

Human mRNA DeXD/H-box helicases are ubiquitous molecular motors that are required for the majority of cellular processes that involve RNA metabolism. One of the most abundant is eIF4A, which is required during the initiation phase of protein synthesis to unwind regions of highly structured mRNA that would otherwise impede the scanning ribosome. Dysregulation of protein synthesis is associated with tumorigenesis, but little is known about the detailed relationships between RNA helicase function and the malignant phenotype in solid malignancies. Therefore, immunohistochemical analysis was performed on over 3000 breast tumors to investigate the relationship among expression of eIF4A1, the helicase-modulating proteins eIF4B, eIF4E and PDCD4, and clinical outcome. We found eIF4A1, eIF4B and eIF4E to be independent predictors of poor outcome in ER-negative disease, while in contrast, the eIF4A1 inhibitor PDCD4 was related to improved outcome in ER-positive breast cancer. Consistent with these data, modulation of eIF4A1, eIF4B and PCDC4 expression in cultured MCF7 cells all restricted breast cancer cell growth and cycling. The eIF4A1-dependent translatome of MCF7 cells was defined by polysome profiling, and was shown to be highly enriched for several classes of oncogenic genes, including G-protein constituents, cyclins and protein kinases, and for mRNAs with G/C-rich 5′UTRs with potential to form G-quadruplexes and with 3′UTRs containing microRNA target sites. Overall, our data show that dysregulation of mRNA unwinding contributes to the malignant phenotype in breast cancer via preferential translation of a class of genes involved in pro-oncogenic signaling at numerous levels. Furthermore, immunohistochemical tests are promising biomarkers for tumors sensitive to anti-helicase therapies

Findings from the Peutz-Jeghers Syndrome Registry of Uruguay

Background: Peutz-Jeghers syndrome (PJS) is characterized by intestinal polyposis, mucocutaneous pigmentation and an increased cancer risk, usually caused by mutations of the STK11 gene. This study collected epidemiological, clinical and genetic data from all Uruguayan PJS patients. Methods: Clinical data were obtained from public and private medical centers and updated annually. Sequencing of the STK11 gene in one member of each family was performed. Results and discussion: 25 cases in 11 unrelated families were registered (15 males, 10 females). The average age of diagnosis and death was 18 and 41 years respectively. All patients had characteristic PJS pigmentation and gastrointestinal polyps. 72% required urgent surgery due to intestinal obstruction. 3 families had multiple cases of seizure disorder, representing 20% of cases. 28% developed cancer and two patients had more than one cancer. An STK11 mutation was found in 8 of the 9 families analyzed. A unique M136K missense mutation was noted in one family. Comparing annual live births and PJS birth records from 1970 to 2009 yielded an incidence of 1 in 155,000. Conclusion: The Uruguayan Registry for Peutz-Jeghers patients showed a high chance of emergent surgery, epilepsy, cancer and shortened life expectancy. The M136K missense mutation is a newly reported STK 11 mutation

Identification of QTLs controlling gene expression networks defined a priori

BACKGROUND: Gene expression microarrays allow the quantification of transcript accumulation for many or all genes in a genome. This technology has been utilized for a range of investigations, from assessments of gene regulation in response to genetic or environmental fluctuation to global expression QTL (eQTL) analyses of natural variation. Current analysis techniques facilitate the statistical querying of individual genes to evaluate the significance of a change in response, also known as differential expression. Since genes are also known to respond as groups due to their membership in networks, effective approaches are needed to investigate transcriptome variation as related to gene network responses. RESULTS: We describe a statistical approach that is capable of assessing higher-order a priori defined gene network response, as measured by microarrays. This analysis detected significant network variation between two Arabidopsis thaliana accessions, Bay-0 and Shahdara. By extending this approach, we were able to identify eQTLs controlling network responses for 18 out of 20 a priori-defined gene networks in a recombinant inbred line population derived from accessions Bay-0 and Shahdara. CONCLUSION: This approach has the potential to be expanded to facilitate direct tests of the relationship between phenotypic trait and transcript genetic architecture. The use of a priori definitions for network eQTL identification has enormous potential for providing direction toward future eQTL analyses

Moving toward a system genetics view of disease

Testing hundreds of thousands of DNA markers in human, mouse, and other species for association to complex traits like disease is now a reality. However, information on how variations in DNA impact complex physiologic processes flows through transcriptional and other molecular networks. In other words, DNA variations impact complex diseases through the perturbations they cause to transcriptional and other biological networks, and these molecular phenotypes are intermediate to clinically defined disease. Because it is also now possible to monitor transcript levels in a comprehensive fashion, integrating DNA variation, transcription, and phenotypic data has the potential to enhance identification of the associations between DNA variation and diseases like obesity and diabetes, as well as characterize those parts of the molecular networks that drive these diseases. Toward that end, we review methods for integrating expression quantitative trait loci (eQTLs), gene expression, and clinical data to infer causal relationships among gene expression traits and between expression and clinical traits. We further describe methods to integrate these data in a more comprehensive manner by constructing coexpression gene networks that leverage pairwise gene interaction data to represent more general relationships. To infer gene networks that capture causal information, we describe a Bayesian algorithm that further integrates eQTLs, expression, and clinical phenotype data to reconstruct whole-gene networks capable of representing causal relationships among genes and traits in the network. These emerging network approaches, aimed at processing high-dimensional biological data by integrating data from multiple sources, represent some of the first steps in statistical genetics to identify multiple genetic perturbations that alter the states of molecular networks and that in turn push systems into disease states. Evolving statistical procedures that operate on networks will be critical to extracting information related to complex phenotypes like disease, as research goes beyond a single-gene focus. The early successes achieved with the methods described herein suggest that these more integrative genomics approaches to dissecting disease traits will significantly enhance the identification of key drivers of disease beyond what could be achieved by genetic association studies alone