Genome-wide interaction-based association (GWIBA) analysis has the potential to identify novel susceptibility loci. These interaction effects could be missed with the prevailing approaches in genome-wide association studies (GWAS). However, no convincing loci have been discovered exclusively from GWIBA methods, and the intensive computation involved is a major barrier for application. Here, we developed a fast, multi-thread/parallel program named "pair-wise interaction-based association mapping" (PIAM) for exhaustive two-locus searches. With this program, we performed a complete GWIBA analysis on seven diseases with stringent control for false positives, and we validated the results for three of these diseases. We identified one pair-wise interaction between a previously identified locus, C1orf106, and one new locus, TEC, that was specific for Crohn's disease, with a Bonferroni corrected P < 0.05 (P = 0.039). This interaction was replicated with a pair of proxy linked loci (P = 0.013) on an independent dataset. Five other interactions had corrected P < 0.5. We identified the allelic effect of a locus close to SLC7A13 for coronary artery disease. This was replicated with a linked locus on an independent dataset (P = 1.09 × 10⁻⁷). Through a local validation analysis that evaluated association signals, rather than locus-based associations, we found that several other regions showed association/interaction signals with nominal P < 0.05. In conclusion, this study demonstrated that the GWIBA approach was successful for identifying novel loci, and the results provide new insights into the genetic architecture of common diseases. In addition, our PIAM program was capable of handling very large GWAS datasets that are likely to be produced in the future.

Aging involves a progressive physiological remodeling that is controlled by both genetic and environmental factors. Many of these factors impact also on white adipose tissue (WAT), which has been shown to be a determinant of lifespan. Interrogating a transcriptional network for predicted causal regulatory interactions in a collection of mouse WAT from F2 crosses with a seed set of 60 known longevity genes, we identified a novel transcriptional subnetwork of 742 genes which represent thus-far-unknown longevity genes. Within this subnetwork, one gene was Pparg (Nr1c3), an adipose-enriched nuclear receptor previously not associated with longevity. In silico, both the PPAR signaling pathway and the transcriptional signature of Ppargamma agonist rosiglitazone overlapped with the longevity subnetwork, while in vivo, lowered expression of Pparg reduced lifespan in both the lipodystrophic Pparg1/2-hypomorphic and the Pparg2-deficient mice. These results establish Ppargamma2 as one of the determinants of longevity and suggest that lifespan may be rather determined by a purposeful genetic program than a random process.

Chronic Obstructive Pulmonary Disease (COPD) is a complex disease. Genetic, epigenetic, and environmental factors are known to contribute to COPD risk and disease progression. Therefore we developed a systematic approach to identify key regulators of COPD that integrates genome-wide DNA methylation, gene expression, and phenotype data in lung tissue from COPD and control samples. Our integrative analysis identified 126 key regulators of COPD. We identified EPAS1 as the only key regulator whose downstream genes significantly overlapped with multiple genes sets associated with COPD disease severity. EPAS1 is distinct in comparison with other key regulators in terms of methylation profile and downstream target genes. Genes predicted to be regulated by EPAS1 were enriched for biological processes including signaling, cell communications, and system development. We confirmed that EPAS1 protein levels are lower in human COPD lung tissue compared to non-disease controls and that Epas1 gene expression is reduced in mice chronically exposed to cigarette smoke. As EPAS1 downstream genes were significantly enriched for hypoxia responsive genes in endothelial cells, we tested EPAS1 function in human endothelial cells. EPAS1 knockdown by siRNA in endothelial cells impacted genes that significantly overlapped with EPAS1 downstream genes in lung tissue including hypoxia responsive genes, and genes associated with emphysema severity. Our first integrative analysis of genome-wide DNA methylation and gene expression profiles illustrates that not only does DNA methylation play a 'causal' role in the molecular pathophysiology of COPD, but it can be leveraged to directly identify novel key mediators of this pathophysiology.

The ubiquitin-proteasome system plays important roles in various biological processes as it degrades the majority of cellular proteins. Adequate proteasomal degradation of crucial transcription regulators ensures the proper development of neutrophils. The ubiquitin E3 ligase of Growth factor independent 1 (GFI1), a key transcription repressor governing terminal granulopoiesis, remains obscure. Here we report that the deficiency of the ring finger protein Interferon regulatory factor 2 binding protein 2a (Irf2bp2a) leads to an impairment of neutrophils differentiation in zebrafish. Mechanistically, Irf2bp2a functions as a ubiquitin E3 ligase targeting Gfi1aa for proteasomal degradation. Moreover, irf2bp2a gene is repressed by Gfi1aa, thus forming a negative feedback loop between Irf2bp2a and Gfi1aa during neutrophils maturation. Different levels of GFI1 may turn it into a tumor suppressor or an oncogene in malignant myelopoiesis. Therefore, discovery of certain drug targets GFI1 for proteasomal degradation by IRF2BP2 might be an effective anti-cancer strategy.

Complex diseases result from molecular changes induced by multiple genetic factors and the environment. To derive a systems view of how genetic loci interact in the context of tissue-specific molecular networks, we constructed an F2 intercross comprised of >500 mice from diabetes-resistant (B6) and diabetes-susceptible (BTBR) mouse strains made genetically obese by the Leptin(ob/ob) mutation (Lep(ob)). High-density genotypes, diabetes-related clinical traits, and whole-transcriptome expression profiling in five tissues (white adipose, liver, pancreatic islets, hypothalamus, and gastrocnemius muscle) were determined for all mice. We performed an integrative analysis to investigate the inter-relationship among genetic factors, expression traits, and plasma insulin, a hallmark diabetes trait. Among five tissues under study, there are extensive protein-protein interactions between genes responding to different loci in adipose and pancreatic islets that potentially jointly participated in the regulation of plasma insulin. We developed a novel ranking scheme based on cross-loci protein-protein network topology and gene expression to assess each gene's potential to regulate plasma insulin. Unique candidate genes were identified in adipose tissue and islets. In islets, the Alzheimer's gene App was identified as a top candidate regulator. Islets from 17-week-old, but not 10-week-old, App knockout mice showed increased insulin secretion in response to glucose or a membrane-permeant cAMP analog, in agreement with the predictions of the network model. Our result provides a novel hypothesis on the mechanism for the connection between two aging-related diseases: Alzheimer's disease and type 2 diabetes.

The majority of the heritability of coronary artery disease (CAD) remains unexplained, despite recent successes of genome-wide association studies (GWAS) in identifying novel susceptibility loci. Integrating functional genomic data from a variety of sources with a large-scale meta-analysis of CAD GWAS may facilitate the identification of novel biological processes and genes involved in CAD, as well as clarify the causal relationships of established processes. Towards this end, we integrated 14 GWAS from the CARDIoGRAM Consortium and two additional GWAS from the Ottawa Heart Institute (25,491 cases and 66,819 controls) with 1) genetics of gene expression studies of CAD-relevant tissues in humans, 2) metabolic and signaling pathways from public databases, and 3) data-driven, tissue-specific gene networks from a multitude of human and mouse experiments. We not only detected CAD-associated gene networks of lipid metabolism, coagulation, immunity, and additional networks with no clear functional annotation, but also revealed key driver genes for each CAD network based on the topology of the gene regulatory networks. In particular, we found a gene network involved in antigen processing to be strongly associated with CAD. The key driver genes of this network included glyoxalase I (GLO1) and peptidylprolyl isomerase I (PPIL1), which we verified as regulatory by siRNA experiments in human aortic endothelial cells. Our results suggest genetic influences on a diverse set of both known and novel biological processes that contribute to CAD risk. The key driver genes for these networks highlight potential novel targets for further mechanistic studies and therapeutic interventions.

The application of deep sequencing to map 5' capped transcripts has confirmed the existence of at least two distinct promoter classes in metazoans: "focused" promoters with transcription start sites (TSSs) that occur in a narrowly defined genomic span and "dispersed" promoters with TSSs that are spread over a larger window. Previous studies have explored the presence of genomic features, such as CpG islands and sequence motifs, in these promoter classes, but virtually no studies have directly investigated the relationship with chromatin features. Here, we show that promoter classes are significantly differentiated by nucleosome organization and chromatin structure. Dispersed promoters display higher associations with well-positioned nucleosomes downstream of the TSS and a more clearly defined nucleosome free region upstream, while focused promoters have a less organized nucleosome structure, yet higher presence of RNA polymerase II. These differences extend to histone variants (H2A.Z) and marks (H3K4 methylation), as well as insulator binding (such as CTCF), independent of the expression levels of affected genes. Notably, differences are conserved across mammals and flies, and they provide for a clearer separation of promoter architectures than the presence and absence of CpG islands or the occurrence of stalled RNA polymerase. Computational models support the stronger contribution of chromatin features to the definition of dispersed promoters compared to focused start sites. Our results show that promoter classes defined from 5' capped transcripts not only reflect differences in the initiation process at the core promoter but also are indicative of divergent transcriptional programs established within gene-proximal nucleosome organization.

Pollen wall consists of several complex layers which form elaborate species-specific patterns. In Arabidopsis, the transcription factor ABORTED MICROSPORE (AMS) is a master regulator of exine formation, and another transcription factor, TRANSPOSABLE ELEMENT SILENCING VIA AT-HOOK (TEK), specifies formation of the nexine layer. However, knowledge regarding the temporal regulatory roles of TEK in pollen wall development is limited. Here, TEK-GFP driven by the AMS promoter was prematurely expressed in the tapetal nuclei, leading to complete male sterility in the pAMS:TEK-GFP (pat) transgenic lines with the wild-type background. Cytological observations in the pat anthers showed impaired callose synthesis and aberrant exine patterning. CALLOSE SYNTHASE5 (CalS5) is required for callose synthesis, and expression of CalS5 in pat plants was significantly reduced. We demonstrated that TEK negatively regulates CalS5 expression after the tetrad stage in wild-type anthers and further discovered that premature TEK-GFP in pat directly represses CalS5 expression through histone modification. Our findings show that TEK flexibly mediates its different functions via different temporal regulation, revealing that the temporal regulation of TEK is essential for exine patterning. Moreover, the result that the repression of CalS5 by TEK after the tetrad stage coincides with the timing of callose wall dissolution suggests that tapetum utilizes temporal regulation of genes to stop callose wall synthesis, which, together with the activation of callase activity, achieves microspore release and pollen wall patterning.