Integrated systems approach identifies genetic nodes and networks in late-onset Alzheimer's disease.

The genetics of complex disease produce alterations in the molecular interactions of cellular pathways whose collective effect may become clear through the organized structure of molecular networks. To characterize molecular systems associated with late-onset Alzheimer's disease (LOAD), we constructed gene-regulatory networks in 1,647 postmortem brain tissues from LOAD patients and nondemented subjects, and we demonstrate that LOAD reconfigures specific portions of the molecular interaction structure. Through an integrative network-based approach, we rank-ordered these network structures for relevance to LOAD pathology, highlighting an immune- and microglia-specific module that is dominated by genes involved in pathogen phagocytosis, contains TYROBP as a key regulator, and is upregulated in LOAD. Mouse microglia cells overexpressing intact or truncated TYROBP revealed expression changes that significantly overlapped the human brain TYROBP network. Thus the causal network structure is a useful predictor of response to gene perturbations and presents a framework to test models of disease mechanisms underlying LOAD.

The Cytochrome P450 (CYP) 2C8*3 Variant (rs11572080) is Associated with Improved Asthma Symptom Control in Children and Altered Lipid Mediator Production and Inflammatory Response in Human Bronchial Epithelial Cells.

This study investigated an association between the cytochrome P450 (CYP) 2C8*3 polymorphism with asthma symptom control in children and changes in lipid metabolism and pro-inflammatory signaling by human bronchial epithelial cells (HBECs) treated with cigarette smoke condensate (CSC). CYP genes are inherently variable in sequence and while such variations are known to produce clinically relevant effects on drug pharmacokinetics and pharmacodynamics, the effects on endogenous substrate metabolism and associated physiological processes are less understood. In this study, CYP2C8*3 was associated with improved asthma symptom control among children: Mean asthma control scores were 3.68 [n=207] for patients with one or more copies of the CYP2C8*3 allele vs. 4.42 [n=965] for CYP2C8*1/*1 (p=0.0133). In vitro, CYP2C8*3 was associated with an increase in montelukast 36-hydroxylation and a decrease in linoleic acid (LA) metabolism despite lower mRNA and protein expression. Additionally, CYP2C8*3 was associated with reduced mRNA expression of interleukin-6 (IL-6) and C-X-C motif chemokine ligand 8 (CXCL-8) by HBECs in response to CSC, which was replicated using the soluble epoxide hydrolase inhibitor, AUDA. Interestingly, 9(10)- and 12(13)-DiHOME, the hydrolyzed metabolites of 9(10)- and 12(13)-EpOME, increased the expression of IL-6 and CXCL-8 mRNA by HBECs. This study reveals previously undocumented effects of the CYP2C8*3 variant on the response of HBECs to exogenous stimuli. Significance Statement These findings suggest a role for CYP2C8 in regulating the EpOME:DiHOME ratio leading to a change in cellular inflammatory responses elicited by environmental stimuli that exacerbate asthma.

Impact of CYP3A5 Polymorphisms on Pediatric Asthma Outcomes.

Genetic variation among inhaled corticosteroid (ICS)-metabolizing enzymes may affect asthma control, but evidence is limited. This study tested the hypothesis that single-nucleotide polymorphisms (SNPs) in Cytochrome P450 3A5 (CYP3A5) would affect asthma outcomes. Patients aged 2-18 years with persistent asthma were recruited to use the electronic AsthmaTracker (e-AT), a self-monitoring tool that records weekly asthma control, medication use, and asthma outcomes. A subset of patients provided saliva samples for SNP analysis and participated in a pharmacokinetic study. Multivariable regression analysis adjusted for age, sex, race, and ethnicity was used to evaluate the impact of CYP3A5 SNPs on asthma outcomes, including asthma control (measured using the asthma symptom tracker, a modified version of the asthma control test or ACT), exacerbations, and hospital admissions. Plasma corticosteroid and cortisol concentrations post-ICS dosing were also assayed using liquid chromatography-tandem mass spectrometry. Of the 751 patients using the e-AT, 166 (22.1%) provided saliva samples and 16 completed the PK study. The e-AT cohort was 65.1% male, and 89.6% White, 6.0% Native Hawaiian, 1.2% Black, 1.2% Native American, 1.8% of unknown race, and 15.7% Hispanic/Latino; the median age was 8.35 (IQR: 5.51-11.3) years. CYP3A5*3/*3 frequency was 75.8% in White subjects, 50% in Native Hawaiians and 76.9% in Hispanic/Latino subjects. Compared with CYP3A5*3/*3, the CYP3A5*1/*x genotype was associated with reduced weekly asthma control (OR: 0.98; 95% CI: 0.97-0.98; p < 0.001), increased exacerbations (OR: 6.43; 95% CI: 4.56-9.07; p < 0.001), and increased asthma hospitalizations (OR: 1.66; 95% CI: 1.43-1.93; p < 0.001); analysis of 3/*3, *1/*1 and *1/*3 separately showed an allelic copy effect. Finally, PK analysis post-ICS dosing suggested muted changes in cortisol concentrations for patients with the CYP3A5*3/*3 genotype, as opposed to an effect on ICS PK. Detection of CYP3A5*3/3, CYPA35*1/*3, and CYP3A5*1/*1 could impact inhaled steroid treatment strategies for asthma in the future.

Limits and potential of combined folding and docking.

MOTIVATION: In the last decade, de novo protein structure prediction accuracy for individual proteins has improved significantly by utilising deep learning (DL) methods for harvesting the co-evolution information from large multiple sequence alignments (MSAs). The same approach can, in principle, also be used to extract information about evolutionary-based contacts across protein-protein interfaces. However, most earlier studies have not used the latest DL methods for inter-chain contact distance prediction. This article introduces a fold-and-dock method based on predicted residue-residue distances with trRosetta. RESULTS: The method can simultaneously predict the tertiary and quaternary structure of a protein pair, even when the structures of the monomers are not known. The straightforward application of this method to a standard dataset for protein-protein docking yielded limited success. However, using alternative methods for generating MSAs allowed us to dock accurately significantly more proteins. We also introduced a novel scoring function, PconsDock, that accurately separates 98% of correctly and incorrectly folded and docked proteins. The average performance of the method is comparable to the use of traditional, template-based or ab initio shape-complementarity-only docking methods. Moreover, the results of conventional and fold-and-dock approaches are complementary, and thus a combined docking pipeline could increase overall docking success significantly. This methodology contributed to the best model for one of the CASP14 oligomeric targets, H1065. AVAILABILITY AND IMPLEMENTATION: All scripts for predictions and analysis are available from https://github.com/ElofssonLab/bioinfo-toolbox/ and https://gitlab.com/ElofssonLab/benchmark5/. All models joined alignments, and evaluation results are available from the following figshare repository https://doi.org/10.6084/m9.figshare.14654886.v2. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The evolutionary history of topological variations in the CPA/AT transporters.

CPA/AT transporters are made up of scaffold and a core domain. The core domain contains two non-canonical helices (broken or reentrant) that mediate the transport of ions, amino acids or other charged compounds. During evolution, these transporters have undergone substantial changes in structure, topology and function. To shed light on these structural transitions, we create models for all families using an integrated topology annotation method. We find that the CPA/AT transporters can be classified into four fold-types based on their structure; (1) the CPA-broken fold-type, (2) the CPA-reentrant fold-type, (3) the BART fold-type, and (4) a previously not described fold-type, the Reentrant-Helix-Reentrant fold-type. Several topological transitions are identified, including the transition between a broken and reentrant helix, one transition between a loop and a reentrant helix, complete changes of orientation, and changes in the number of scaffold helices. These transitions are mainly caused by gene duplication and shuffling events. Structural models, topology information and other details are presented in a searchable database, CPAfold (cpafold.bioinfo.se).

DisProt: intrinsic protein disorder annotation in 2020.

The Database of Protein Disorder (DisProt, URL: https://disprot.org) provides manually curated annotations of intrinsically disordered proteins from the literature. Here we report recent developments with DisProt (version 8), including the doubling of protein entries, a new disorder ontology, improvements of the annotation format and a completely new website. The website includes a redesigned graphical interface, a better search engine, a clearer API for programmatic access and a new annotation interface that integrates text mining technologies. The new entry format provides a greater flexibility, simplifies maintenance and allows the capture of more information from the literature. The new disorder ontology has been formalized and made interoperable by adopting the OWL format, as well as its structure and term definitions have been improved. The new annotation interface has made the curation process faster and more effective. We recently showed that new DisProt annotations can be effectively used to train and validate disorder predictors. We believe the growth of DisProt will accelerate, contributing to the improvement of function and disorder predictors and therefore to illuminate the 'dark' proteome.

Dynamic Expression of Transient Receptor Potential Vanilloid-3 and Integrated Signaling with Growth Factor Pathways during Lung Epithelial Wound Repair following Wood Smoke Particle and Other Forms of Lung Cell Injury.

Prior studies revealed increased expression of the transient receptor potential vanilloid-3 (TRPV3) ion channel after wood smoke particulate matter (WSPM) treatment of human bronchial epithelial cells (HBECs). TRPV3 attenuated pathologic endoplasmic reticulum stress and cytotoxicity mediated by transient receptor potential ankyrin-1. Here, the basis for how TRPV3 expression is regulated by cell injury and the effects this has on HBEC physiology and WSPM-induced airway remodeling in mice was investigated. TRPV3 mRNA was rapidly increased in HBECs treated with WSPM and after monolayer damage caused by tryptic disruption, scratch wounding, and cell passaging. TRPV3 mRNA abundance varied with time, and stimulated expression occurred independent of new protein synthesis. Overexpression of TRPV3 in HBECs reduced cell migration and wound repair while enhancing cell adhesion. This phenotype correlated with disrupted mRNA expression of ligands of the epidermal growth factor, tumor growth factor-ß, and frizzled receptors. Accordingly, delayed wound repair by TRPV3 overexpressing cells was reversed by growth factor supplementation. In normal HBECs, TRPV3 upregulation was triggered by exogenous growth factor supplementation and was attenuated by inhibitors of growth factor receptor signaling. In mice, subacute oropharyngeal instillation with WSPM also promoted TRPV3 mRNA expression and epithelial remodeling, which was attenuated by TRPV3 antagonist pre- and cotreatment. This latter effect may be the consequence of antagonist-induced TRPV3 expression. These findings provide insights into the roles of TRPV3 in lung epithelial cells under basal and dynamic states, as well as highlight potential roles for TRPV3 ligands in modulating epithelial damage/repair. SIGNIFICANCE STATEMENT: Coordinated epithelial repair is essential for the maintenance of the airways, with deficiencies and exaggerated repair associated with adverse consequences to respiratory health. This study shows that TRPV3, an ion channel, is involved in coordinating repair through integrated repair signaling pathways, wherein TRPV3 expression is upregulated immediately after injury and returns to basal levels as cells complete the repair process. TRPV3 may be a novel target for understanding and/or treating conditions in which airway/lung epithelial repair is not properly orchestrated.

Guide Swap enables genome-scale pooled CRISPR-Cas9 screening in human primary cells.

CRISPR-Cas9 screening allows genome-wide interrogation of gene function. Currently, to achieve the high and uniform Cas9 expression desirable for screening, one needs to engineer stable and clonal Cas9-expressing cells-an approach that is not applicable in human primary cells. Guide Swap permits genome-scale pooled CRISPR-Cas9 screening in human primary cells by exploiting the unexpected finding that editing by lentivirally delivered, targeted guide RNAs (gRNAs) occurs efficiently when Cas9 is introduced in complex with nontargeting gRNA. We validated Guide Swap in depletion and enrichment screens in CD4+ T cells. Next, we implemented Guide Swap in a model of ex vivo hematopoiesis, and identified known and previously unknown regulators of CD34+ hematopoietic stem and progenitor cell (HSPC) expansion. We anticipate that this platform will be broadly applicable to other challenging cell types, and thus will enable discovery in previously inaccessible but biologically relevant human primary cell systems.

Tau overexpression exacerbates neuropathology after repeated mild head impacts in male mice.

Repeated mild traumatic brain injury (rmTBI) can lead to development of chronic traumatic encephalopathy (CTE), which is characterized by progressive neurodegeneration with presence of white matter damage, gliosis and hyper-phosphorylated tau. While animal models of rmTBI have been documented, few characterize the molecular pathogenesis and expression profiles of relevant injured brain regions. Additionally, while the usage of transgenic tau mice in rmTBI is prevalent, the effects of tau on pathological outcomes has not been well studied. Here we characterized a 42-impact closed-head rmTBI paradigm on 3-4 month old male C57BL/6 (WT) and Tau-overexpressing mice (Tau58.4). This injury paradigm resulted in chronic gliosis, T-cell infiltration, and demyelination of the optic nerve and associated white matter tracts at 1-month post-injury. At 3-months post-injury, Tau58.4 mice showed progressive neuroinflammation and neurodegeneration in multiple brain regions compared to WT mice. Corresponding to histopathology, RNAseq of the optic nerve tract at 1-month post-injury showed significant upregulation of inflammatory pathways and downregulation of myelin synthetic pathways in both genotypes. However, Tau58.4 mice showed additional changes in neurite development, protein processing, and cell stress. Comparisons with published transcriptomes of human Alzheimer's Disease and CTE revealed common signatures including neuroinflammation and downregulation of protein phosphatases. We next investigated the demyelination and T-cell infiltration phenotypes to determine whether these offer potential avenues for therapeutic intervention. Tau58.4 mice were treated with the histamine H3 receptor antagonist GSK239512 for 1-month post-injury to promote remyelination of white matter lesions. This restored myelin gene expression to sham levels but failed to repair the histopathologic lesions. Likewise, injured T-cell-deficient Rag2/Il2rg (R2G2) mice also showed evidence for inflammation and loss of myelin. However, unlike immune-competent mice, R2G2 mice had altered myeloid cell gene expression and fewer demyelinated lesions. Together this data shows that rmTBI leads to chronic white matter inflammatory demyelination and axonal loss exacerbated by human tau overexpression but suggests that immune-suppression and remyelination alone are insufficient to reverse damage.

In vivo discovery of immunotherapy targets in the tumour microenvironment.

Recent clinical trials showed that targeting of inhibitory receptors on T cells induces durable responses in a subset of cancer patients, despite advanced disease. However, the regulatory switches controlling T-cell function in immunosuppressive tumours are not well understood. Here we show that such inhibitory mechanisms can be systematically discovered in the tumour microenvironment. We devised an in vivo pooled short hairpin RNA (shRNA) screen in which shRNAs targeting negative regulators became highly enriched in murine tumours by releasing a block on T-cell proliferation upon tumour antigen recognition. Such shRNAs were identified by deep sequencing of the shRNA cassette from T cells infiltrating tumour or control tissues. One of the target genes was Ppp2r2d, a regulatory subunit of the PP2A phosphatase family. In tumours, Ppp2r2d knockdown inhibited T-cell apoptosis and enhanced T-cell proliferation as well as cytokine production. Key regulators of immune function can therefore be discovered in relevant tissue microenvironments.

Cultivar and phosphorus effects on switchgrass yield and rhizosphere microbial diversity.

Switchgrass (Panicum virgatum L.) is a native perennial grass identified as a promising biofuel crop for production on marginal agricultural lands. As such, research into switchgrass fertility and the switchgrass rhizosphere microbiome has been ongoing in an effort to increase production sustainability. We examined the effects of cultivar and phosphorus (P) fertilization on biomass yield, P removal, and rhizosphere bacterial and fungal community structure in three switchgrass cultivars: Sunburst, Shawnee, and Liberty. The Liberty cv. is the first lowland-type bioenergy switchgrass adapted to USDA hardiness zones 4, 5, and 6. On a medium soil test P clay loam soil, biomass yield response to applied P was linear, increasing 135 kg ha-1 for every kilogram of P applied prior to establishment. Average post-frost biomass yield was 9.6 Mg ha-1 year-1 when unfertilized, and maximum biomass yield was 10.3 Mg ha-1 year-1 when fertilized at 58.6 kg ha-1 P, suggesting that P application on medium soil test P soils is beneficial for switchgrass establishment and early growth. Switchgrass cv. Shawnee was more productive than cvs. Liberty or Sunburst (11.3, 10.2, and 8.6 Mg ha-1 year-1, respectively). Both bacterial and fungal communities were significantly shaped by cultivar. These shifts, while inconsistent between year and cultivar, may reflect a selection of the microbial community from that present in soil to maximize total nutrient uptake, regardless of additional P amendments. Phosphorus fertilization did not affect microbial community structure. Results of this study suggest that the cultivar-associated selection of particular microbial taxa may have implications for increased productivity.

Activation of Human Transient Receptor Potential Melastatin-8 (TRPM8) by Calcium-Rich Particulate Materials and Effects on Human Lung Cells.

To better understand how adverse health effects are caused by exposure to particulate materials, and to develop preventative measures, it is important to identify the properties of particles and molecular targets that link exposure with specific biologic outcomes. Coal fly ash (CFA) is a by-product of coal combustion that can affect human health. We report that human transient receptor potential melastatin-8 (TRPM8) and an N-terminally truncated TRPM8 variant (TRPM8-Δ801) are activated by CFA and calcium-rich nanoparticles and/or soluble salts within CFA. TRPM8 activation by CFA was potentiated by cold temperature involving the phosphatidylinositol 4,5-bisphosphate binding residue (L1008), but was independent of the icilin and menthol binding site residue Y745 and, essentially, the N-terminal amino acids 1-800. CFA, calcium nanoparticles, and calcium salts also activated transient receptor potential vanilloid-1 (TRPV1) and transient receptor potential ankyrin-1 (TRPA1), but not TRPV4. CFA treatment induced CXCL1 and interleukin-8 mRNA in BEAS-2B and primary human bronchial epithelial cells through activation of both TRPM8 and TRPV1. However, neither mouse nor rat TRPM8 was activated by these materials, and Trpm8 knockout had no effect on cytokine induction in the lungs of CFA-instilled mice. Amino acids S921 and S927 in mouse Trpm8 were identified as important for the lack of response to CFA. These results imply that TRPM8, in conjunction with TRPV1 and TRPA1, might sense selected forms of inhaled particulate materials in human airways, shaping cellular responses to these materials, and improving our understanding of how and why certain particulate materials elicit different responses in biologic systems, affecting human health.

Identification of Cyclic Depsipeptides and Their Dedicated Synthetase from Hapsidospora irregularis.

Seven cyclic depsipeptides were isolated from Hapsidospora irregularis and structurally characterized as the calcium channel blocker leualacin and six new analogues based on the NMR and HRESIMS data. These new compounds were named leualacins B-G. The absolute configurations of the amino acids and 2-hydroxyisocaproic acids were determined by recording the optical rotation values. Biological studies showed that calcium influx elicited by leualacin F in primary human lobar bronchial epithelial cells involves the TRPA1 channel. Through genome sequencing and targeted gene disruption, a noniterative nonribosomal peptide synthetase was found to be involved in the biosynthesis of leualacin. A comparison of the structures of leualacin and its analogues indicated that the A2 and A4 domains of the leualacin synthetase are substrate specific, while A1, A3, and A5 can accept alternative precursors to yield new molecules.

A new approach for trace analysis of guanidine compounds in surface water with resorcinarene-based ion chromatography columns.

Trace levels of pharmaceuticals have been detected in surface water and may pose a health risk to humans and other organisms. New chromatographic materials will help identify and quantify these contaminants. We introduce a new ion chromatographic (IC) material designed to separate cationic pharmaceuticals and report its ability to separate a group of guanidine compounds. Guanidine moieties are strongly basic and protonated under acid conditions, and therefore can potentially be separated on the newly designed stationary phase and detected by ion exchange chromatography. The new column packing material is based on glutamic acids bonded to resorcinarene moieties that in turn are bound to divinylbenzene macroporous resin. Detection limits in the range of 5-30 µg L(-1) were achieved using integrated pulsed amperometric detection (IPAD) for guanidine (G), methylguanidine (MG), 1,1-dimethylbiguanide (DMG), agmatine (AGM), guanidinobenzoic acid (GBA) and cimetidine (CIM). Suppressed conductivity (CD) and UV-vis detection resulted in limits of detection similar to IPAD, in the range of 2-66 µg L(-1), but were not able to detect all of the analytes. Three water sources, river, lake, and marsh, were analyzed and despite matrix effects, sensitivity for guanidine compounds was in the 100 µg L(-1) range and apparent recoveries were 80-96%. The peak area precision was 0.01-2.89% for IPAD, CD and UV-vis detection.

Nitrogen Management for Corn and Groundwater Quality in Upper Midwest Irrigated Sands.

Groundwater contamination from NO-N leaching in corn ( L.) production with coarse-textured soils poses an environmental concern. Our objectives were to evaluate NO-N leaching in continuous corn (CC), corn after soybean ( L.) (CSb), and soybean after corn (SbC) in irrigated sandy soils in Minnesota related to (i) N rate using best management practices of split-N application, (ii) a split-N application and single preplant applications of enhanced-efficiency fertilizers (EEF), and (iii) residual N treatment in SbC. Urea (0-315 kg N ha in 45-kg increments) was broadcast as a split application (half at preplant and half at the V4 development stage) and polymer-coated urea (ESN), ESN/urea, and SuperU at preplant at a rate of 180 kg N ha on an Arvilla sandy loam soil. In May and June, 75% of the total drainage and 73% of the total NO-N leached occurred. At the economic optimum N rate (EONR), season-long NO-N leaching rates were 86 and 106 kg NO-N ha for CC and CSb, respectively. In CC, reducing the EONR by 20% reduced grain yield by 4% and NO-N leached by 9%, and a 25% reduction in EONR resulted in an additional 2% reduction for both, whereas no significant reductions occurred for CSb. Similar NO-N leaching occurred with EEFs and the split-N application. After 4 yr of no N application, we measured 9 to 20 mg NO-N L and leaching of 21 to 51 kg NO-N ha, highlighting the difficulty of meeting drinking water quality standards in corn cropping systems.

Common dysregulation network in the human prefrontal cortex underlies two neurodegenerative diseases.

Using expression profiles from postmortem prefrontal cortex samples of 624 dementia patients and non-demented controls, we investigated global disruptions in the co-regulation of genes in two neurodegenerative diseases, late-onset Alzheimer's disease (AD) and Huntington's disease (HD). We identified networks of differentially co-expressed (DC) gene pairs that either gained or lost correlation in disease cases relative to the control group, with the former dominant for both AD and HD and both patterns replicating in independent human cohorts of AD and aging. When aligning networks of DC patterns and physical interactions, we identified a 242-gene subnetwork enriched for independent AD/HD signatures. This subnetwork revealed a surprising dichotomy of gained/lost correlations among two inter-connected processes, chromatin organization and neural differentiation, and included DNA methyltransferases, DNMT1 and DNMT3A, of which we predicted the former but not latter as a key regulator. To validate the inter-connection of these two processes and our key regulator prediction, we generated two brain-specific knockout (KO) mice and show that Dnmt1 KO signature significantly overlaps with the subnetwork (P = 3.1 × 10(-12)), while Dnmt3a KO signature does not (P = 0.017).

Variations in DNA elucidate molecular networks that cause disease.

Identifying variations in DNA that increase susceptibility to disease is one of the primary aims of genetic studies using a forward genetics approach. However, identification of disease-susceptibility genes by means of such studies provides limited functional information on how genes lead to disease. In fact, in most cases there is an absence of functional information altogether, preventing a definitive identification of the susceptibility gene or genes. Here we develop an alternative to the classic forward genetics approach for dissecting complex disease traits where, instead of identifying susceptibility genes directly affected by variations in DNA, we identify gene networks that are perturbed by susceptibility loci and that in turn lead to disease. Application of this method to liver and adipose gene expression data generated from a segregating mouse population results in the identification of a macrophage-enriched network supported as having a causal relationship with disease traits associated with metabolic syndrome. Three genes in this network, lipoprotein lipase (Lpl), lactamase beta (Lactb) and protein phosphatase 1-like (Ppm1l), are validated as previously unknown obesity genes, strengthening the association between this network and metabolic disease traits. Our analysis provides direct experimental support that complex traits such as obesity are emergent properties of molecular networks that are modulated by complex genetic loci and environmental factors.

Genetics of gene expression and its effect on disease.

Common human diseases result from the interplay of many genes and environmental factors. Therefore, a more integrative biology approach is needed to unravel the complexity and causes of such diseases. To elucidate the complexity of common human diseases such as obesity, we have analysed the expression of 23,720 transcripts in large population-based blood and adipose tissue cohorts comprehensively assessed for various phenotypes, including traits related to clinical obesity. In contrast to the blood expression profiles, we observed a marked correlation between gene expression in adipose tissue and obesity-related traits. Genome-wide linkage and association mapping revealed a highly significant genetic component to gene expression traits, including a strong genetic effect of proximal (cis) signals, with 50% of the cis signals overlapping between the two tissues profiled. Here we demonstrate an extensive transcriptional network constructed from the human adipose data that exhibits significant overlap with similar network modules constructed from mouse adipose data. A core network module in humans and mice was identified that is enriched for genes involved in the inflammatory and immune response and has been found to be causally associated to obesity-related traits.

Integrating genotypic and expression data in a segregating mouse population to identify 5-lipoxygenase as a susceptibility gene for obesity and bone traits.

Forward genetic approaches to identify genes involved in complex traits such as common human diseases have met with limited success. Fine mapping of linkage regions and validation of positional candidates are time-consuming and not always successful. Here we detail a hybrid procedure to map loci involved in complex traits that leverages the strengths of forward and reverse genetic approaches. By integrating genotypic and expression data in a segregating mouse population, we show how clusters of expression quantitative trait loci linking to regions of the genome accurately reflect the underlying perturbation to the transcriptional network induced by DNA variations in genes that control the complex traits. By matching patterns of gene expression in a segregating population with expression responses induced by single-gene perturbation experiments, we show how genes controlling clusters of expression and clinical quantitative trait loci can be mapped directly. We demonstrate the utility of this approach by identifying 5-lipoxygenase as underlying previously identified quantitative trait loci in an F(2) cross between strains C57BL/6J and DBA/2J and showing that it has pleiotropic effects on body fat, lipid levels and bone density.

An integrative genomics approach to infer causal associations between gene expression and disease.

A key goal of biomedical research is to elucidate the complex network of gene interactions underlying complex traits such as common human diseases. Here we detail a multistep procedure for identifying potential key drivers of complex traits that integrates DNA-variation and gene-expression data with other complex trait data in segregating mouse populations. Ordering gene expression traits relative to one another and relative to other complex traits is achieved by systematically testing whether variations in DNA that lead to variations in relative transcript abundances statistically support an independent, causative or reactive function relative to the complex traits under consideration. We show that this approach can predict transcriptional responses to single gene-perturbation experiments using gene-expression data in the context of a segregating mouse population. We also demonstrate the utility of this approach by identifying and experimentally validating the involvement of three new genes in susceptibility to obesity.