Weighted Interaction SNP Hub (WISH) network method for building genetic networks for complex diseases and traits using whole genome genotype data

Genome-Wide Epistatic Network Analyses of Semantic Fluency in Older Adults

Semantic fluency impairment has been attributed to a wide range of neurocognitive and psychiatric conditions, especially in the older population. Moderate heritability estimates on semantic fluency were obtained from both twin and family-based studies suggesting genetic contributions to the observed variation across individuals. Currently, effort in identifying the genetic variants underlying the heritability estimates for this complex trait remains scarce. Using the semantic fluency scale and genome-wide SNP genotype data from the Long Life Family Study (LLFS), we performed a genome-wide association study (GWAS) and epistasis network analysis on semantic fluency in 2289 individuals aged over 60 years from the American LLFS cohorts and replicated the findings in 1129 individuals aged over 50 years from the Danish LLFS cohort. In the GWAS, two SNPs with genome-wide significance (rs3749683, p = 2.52 × 10-8; rs880179, p = 4.83 × 10-8) mapped to the CMYAS gene on chromosome 5 were detected. The epistasis network analysis identified five modules as significant (4.16 × 10-5 < p < 7.35 × 10-3), of which two were replicated (p < 3.10 × 10-3). These two modules revealed significant enrichment of tissue-specific gene expression in brain tissues and high enrichment of GWAS catalog traits, e.g., obesity-related traits, blood pressure, chronotype, sleep duration, and brain structure, that have been reported to associate with verbal performance in epidemiological studies. Our results suggest high tissue specificity of genetic regulation of gene expression in brain tissues with epistatic SNP networks functioning jointly in modifying individual verbal ability and cognitive performance.

The effect of STAT1, miR-99b, and MAP2K1 in alcoholic liver disease (ALD) mouse model and hepatocyte

Alcoholic liver disease (ALD) serves as the leading cause of chronic liver diseases-related morbidity and mortality, which threatens the life of millions of patients in the world. However, the molecular mechanisms underlying ALD progression remain unclear. Here, we applied microarray analysis and experimental approaches to identify miRNAs and related regulatory signaling that associated with ALD. Microarray analysis identified that the expression of miR-99b was elevated in the ALD mouse model. The AML-12 cells were treated with EtOH and the expression of miR-99b was enhanced in the cells. The expression of miR-99b was positively correlated with ALT levels in the ALD mice. The microarray analysis identified the abnormally expressed mRNAs in ALD mice and the overlap analysis was performed with based on the differently expressed mRNAs and the transcriptional factors of miR-99b, in which STAT1 was identified. The elevated expression of STAT1 was validated in ALD mice. Meanwhile, the treatment of EtOH induced the expression of STAT1 in the AML-12 cells. The expression of STAT1 was positively correlated with ALT levels in the ALD mice. The positive correlation of STAT1 and miR-99b expression was identified in bioinformatics analysis and ALD mice. The expression of miR-99b and pri-miR-99b was promoted by the overexpression of STAT1 in AML-12 cells. ChIP analysis confirmed the enrichment of STAT1 on miR-99b promoter in AML-12 cells. Next, we found that the expression of mitogen-activated protein kinase kinase 1 (MAP2K1) was negatively associated with miR-99b. The expression of MAP2K1 was downregulated in ALD mice. Consistently, the expression of MAP2K1 was reduced by the treatment of EtOH in AML-12 cells. The expression of MAP2K1 was negative correlated with ALT levels in the ALD mice. We identified the binding site of MAP2K1 and miR-99b. Meanwhile, the treatment of miR-99b mimic repressed the luciferase activity of MAP2K1 in AML-12 cells. The expression of MAP2K1 was suppressed by miR-99b in the cells. We observed that the expression of MAP2K1 was inhibited by the overexpression of STAT1 in AML-12 cells. Meanwhile, the apoptosis of AML-12 cells was induced by the treatment of EtOH, while miR-99b mimic promoted but the overexpression of MAP2K1 attenuated the effect of EtOH in the cells. In conclusion, we identified the correlation and effect of STAT1, miR-99b, and MAP2K1 in ALD mouse model and hepatocyte. STAT1, miR-99b, and MAP2K1 may serve as potential therapeutic target of ALD.

Orchestrating information across tissues via a novel multitask GAT framework to improve quantitative gene regulation relation modeling for survival analysis

Survival analysis is critical to cancer prognosis estimation. High-throughput technologies facilitate the increase in the dimension of genic features, but the number of clinical samples in cohorts is relatively small due to various reasons, including difficulties in participant recruitment and high data-generation costs. Transcriptome is one of the most abundantly available OMIC (referring to the high-throughput data, including genomic, transcriptomic, proteomic and epigenomic) data types. This study introduced a multitask graph attention network (GAT) framework DQSurv for the survival analysis task. We first used a large dataset of healthy tissue samples to pretrain the GAT-based HealthModel for the quantitative measurement of the gene regulatory relations. The multitask survival analysis framework DQSurv used the idea of transfer learning to initiate the GAT model with the pretrained HealthModel and further fine-tuned this model using two tasks i.e. the main task of survival analysis and the auxiliary task of gene expression prediction. This refined GAT was denoted as DiseaseModel. We fused the original transcriptomic features with the difference vector between the latent features encoded by the HealthModel and DiseaseModel for the final task of survival analysis. The proposed DQSurv model stably outperformed the existing models for the survival analysis of 10 benchmark cancer types and an independent dataset. The ablation study also supported the necessity of the main modules. We released the codes and the pretrained HealthModel to facilitate the feature encodings and survival analysis of transcriptome-based future studies, especially on small datasets. The model and the code are available at http://www.healthinformaticslab.org/supp/.

Genome-Wide Association Studies in Sunflower: Towards Sclerotinia sclerotiorum and Diaporthe/Phomopsis Resistance Breeding

Diseases caused by necrotrophic fungi, such as the cosmopolitan Sclerotinia sclerotiorum and the Diaporthe/Phomopsis complex, are among the most destructive diseases of sunflower worldwide. The lack of complete resistance combined with the inefficiency of chemical control makes assisted breeding the best strategy for disease control. In this work, we present an integrated genome-wide association (GWA) study investigating the response of a diverse panel of sunflower inbred lines to both pathogens. Phenotypic data for Sclerotinia head rot (SHR) consisted of five disease descriptors (disease incidence, DI; disease severity, DS; area under the disease progress curve for DI, AUDPCI, and DS, AUDPCS; and incubation period, IP). Two disease descriptors (DI and DS) were evaluated for two manifestations of Diaporthe/Phomopsis: Phomopsis stem canker (PSC) and Phomopsis head rot (PHR). In addition, a principal component (PC) analysis was used to derive transformed phenotypes as inputs to a univariate GWA (PC-GWA). Genotypic data comprised a panel of 4269 single nucleotide polymorphisms (SNP), generated via genotyping-by-sequencing. The GWA analysis revealed 24 unique marker-trait associations for SHR, 19 unique marker-trait associations for Diaporthe/Phomopsis diseases, and 7 markers associated with PC1 and PC2. No common markers were found for the response to the two pathogens. Nevertheless, epistatic interactions were identified between markers significantly associated with the response to S. sclerotiorum and Diaporthe/Phomopsis. This suggests that, while the main determinants of resistance may differ for the two pathogens, there could be an underlying common genetic basis. The exploration of regions physically close to the associated markers yielded 364 genes, of which 19 were predicted as putative disease resistance genes. This work presents the first simultaneous evaluation of two manifestations of Diaporthe/Phomopsis in sunflower, and undertakes a comprehensive GWA study by integrating PSC, PHR, and SHR data. The multiple regions identified, and their exploration to identify candidate genes, contribute not only to the understanding of the genetic basis of resistance, but also to the development of tools for assisted breeding.

Genome-wide epistatic interactions of litter size at birth in Chinese indigenous pigs

Improving litter size at birth (TNB) and the number of piglets born alive (NBA) are the main breeding goals related to litter traits, which are economically important. A better understanding of genetic architecture underlying TNB and NBA traits could increase pig production efficiency. However, most previous studies on these traits focus on additive genetic effects, while epistatic interactions underlying TNB and NBA traits has not yet been well investigated, which are essential to understand how traits-related genes interact. Herein, we conducted genome scans of epistatic interactions underlying TNB and NBA traits in a total of 150 Chinese indigenous pigs (75 Jinhua and 75 Shengxian Spotted pigs) with high throughput genomic data. Based on SNPs with high interaction values and connectivity scores, we identified eight promising candidate genes (AKT2, TSC1, MTOR, PIK3R5, TIAM1, FGF14, RALB and ROR2) potentially associated with litter traits in pigs. Moreover, the underlying pathways, e.g., calcium ion transport, pointed out their roles in litter size-related traits. Our findings provide new insight into genetic architecture of litter traits in pigs and will benefit economic profits in pig production.

Specific Differentially Methylated and Expressed Genes in People with Longevity Family History

Background

We attempt to identify specific differentially methylated and expressed genes in people with longevity family history, it will contribute to discover significant features about human longevity.

Methods

A prevalence study was conducted during October 2017 to January 2019 in Bama County of Guangxi, China and individuals were recruited and grouped into longevity family (n=60) and non-longevity family (n=60) to identify differentially methylated genes (DMGs). The expression profile dataset GSE16717 was downloaded from the GEO database in which individuals were divided into 3 groups, namely longevity (n=50), longevity offspring (n=50) and control (n=50) for identifying differentially expressed genes (DEGs). It was considered significantly different when P or adjusted P≤0.05.

Results

In total, 117 longevity-related hypermethylated genes enriched in interleukin secretion/production regulation, chemokine signaling pathway and natural killer cell-mediated cytotoxicity. Another 296 significant key longevity-related DEGs primarily involved in protein binding, nucleus, cytoplasm, T cell receptor signaling pathway and Metabolic pathway, H19 and PFKFB4 were found to be both methylated and downregulated in people with longevity family history.

Conclusion

Human longevity-specific genes involve in many immunity regulations and cellular immunity pathways, H19 and PFKFB4 show hypermethylated and suppressed status in people with longevity family history and might serve as longevity candidate genes.

Extensive haplotypes are associated with population differentiation and environmental adaptability in Upland cotton (Gossypium hirsutum)

Three extensive eco-haplotypes associated with population differentiation and environmental adaptability in Upland cotton were identified, with A06_85658585, A08_43734499 and A06_113104285 considered the eco-loci for environmental adaptability. Population divergence is suggested to be the primary force driving the evolution of environmental adaptability in various species. Chromosome inversion increases reproductive isolation between subspecies and accelerates population divergence to adapt to new environments. Although modern cultivated Upland cotton (Gossypium hirsutum L.) has spread worldwide, the noticeable phenotypic differences still existed among cultivars grown in different areas. In recent years, the long-distance migration of cotton cultivation areas throughout China has demanded that breeders better understand the genetic basis of environmental adaptability in Upland cotton. Here, we integrated the genotypes of 419 diverse accessions, long-term environment-associated variables (EAVs) and environment-associated traits (EATs) to evaluate subgroup differentiation and identify adaptive loci in Upland cotton. Two highly divergent genomic regions were found on chromosomes A06 and A08, which likely caused by extensive chromosome inversions. The subgroups could be geographically classified based on distinct haplotypes in the divergent regions. A genome-wide association study (GWAS) also confirmed that loci located in these regions were significantly associated with environmental adaptability in Upland cotton. Our study first revealed the cause of population divergence in Upland cotton, as well as the consequences of variation in its environmental adaptability. These findings provide new insights into the genetic basis of environmental adaptability in Upland cotton, which could accelerate the development of molecular markers for adaptation to climate change in future cotton breeding.

Co-expression network analysis identifies a gene signature as a predictive biomarker for energy metabolism in osteosarcoma

Background

Osteosarcoma (OS) is a common malignant bone tumor originating in the interstitial tissues and occurring mostly in adolescents and young adults. Energy metabolism is a prerequisite for cancer cell growth, proliferation, invasion, and metastasis. However, the gene signatures associated with energy metabolism and their underlying molecular mechanisms that drive them are unknown.

Methods

Energy metabolism-related genes were obtained from the TARGET database. We applied the “NFM” algorithm to classify putative signature gene into subtypes based on energy metabolism. Key genes related to progression were identified by weighted co-expression network analysis (WGCNA). Based on least absolute shrinkage and selection operator (LASSO) Cox proportional regression hazards model analyses, a gene signature for the predication of OS progression and prognosis was established. Robustness and estimation evaluations and comparison against other models were used to evaluate the prognostic performance of our model.

Results

Two subtypes associated with energy metabolism was determined using the “NFM” algorithm, and significant modules related to energy metabolism were identified by WGCNA. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) suggested that the genes in the significant modules were enriched in kinase, immune metabolism processes, and metabolism-related pathways. We constructed a seven-gene signature consisting of SLC18B1, RBMXL1, DOK3, HS3ST2, ATP6V0D1, CCAR1, and C1QTNF1 to be used for OS progression and prognosis. Upregulation of CCAR1, and C1QTNF1 was associated with augmented OS risk, whereas, increases in the expression SCL18B1, RBMXL1, DOK3, HS3ST2, and ATP6VOD1 was correlated with a diminished risk of OS. We confirmed that the seven-gene signature was robust, and was superior to the earlier models evaluated; therefore, it may be used for timely OS diagnosis, treatment, and prognosis.

Conclusions

The seven-gene signature related to OS energy metabolism developed here could be used in the early diagnosis, treatment, and prognosis of OS.

Genome-Wide Epistatic Interaction Networks Affecting Feed Efficiency in Duroc and Landrace Pigs

Interactions among genomic loci have often been overlooked in genome-wide association studies, revealing the combinatorial effects of variants on phenotype or disease manifestation. Unexplained genetic variance, interactions among causal genes of small effects, and biological pathways could be identified using a network biology approach. The main objective of this study was to determine the genome-wide epistatic variants affecting feed efficiency traits [feed conversion ratio (FCR) and residual feed intake (RFI)] based on weighted interaction SNP hub (WISH-R) method. Herein, we detected highly interconnected epistatic SNP modules, pathways, and potential biomarkers for the FCR and RFI in Duroc and Landrace purebreds considering the whole population, and separately for low and high feed efficient groups. Highly interacting SNP modules in Duroc (1,247 SNPs) and Landrace (1,215 SNPs) across the population and for low feed efficient (Duroc-80 SNPs, Landrace-146 SNPs) and high feed efficient group (Duroc-198 SNPs, Landrace-232 SNPs) for FCR and RFI were identified. Gene and pathway analyses identified ABL1, MAP3K4, MAP3K5, SEMA6A, KITLG, and KAT2B from chromosomes 1, 2, 5, and 13 underlying ErbB, Ras, Rap1, thyroid hormone, axon guidance pathways in Duroc. GABBR2, GNA12, and PRKCG genes from chromosomes 1, 3, and 6 pointed towards thyroid hormone, cGMP-PKG and cAMP pathways in Landrace. From Duroc low feed efficient group, the TPK1 gene was found involved with thiamine metabolism, whereas PARD6G, DLG2, CRB1 were involved with the hippo signaling pathway in high feed efficient group. PLOD1 and SETD7 genes were involved with lysine degradation in low feed efficient group in Landrace, while high feed efficient group pointed to genes underpinning valine, leucine, isoleucine degradation, and fatty acid elongation. Some SNPs and genes identified are known for their association with feed efficiency, others are novel and potentially provide new avenues for further research. Further validation of epistatic SNPs and genes identified here in a larger cohort would help to establish a framework for modelling epistatic variance in future methods of genomic prediction, increasing the accuracy of estimated genetic merit for FE and helping the pig breeding industry.

Identification of the Biomarkers and Pathological Process of Heterotopic Ossification: Weighted Gene Co-Expression Network Analysis

Heterotopic ossification (HO) is the formation of abnormal mature lamellar bone in extra-skeletal sites, including soft tissues and joints, which result in high rates of disability. The understanding of the mechanism of HO is insufficient. The aim of this study was to explore biomarkers and pathological processes in HO+ samples. The gene expression profile GSE94683 was downloaded from the Gene Expression Omnibus database. Sixteen samples from nine HO- and seven HO+ subjects were analyzed. After data preprocessing, 3,529 genes were obtained for weighted gene co-expression network analysis. Highly correlated genes were divided into 13 modules. Finally, the cyan and purple modules were selected for further study. Gene ontology functional annotation and Kyoto Encyclopedia of Genes and Genomes pathway enrichment indicated that the cyan module was enriched in a variety of components, including protein binding, membrane, nucleoplasm, cytosol, poly(A) RNA binding, biosynthesis of antibiotics, carbon metabolism, endocytosis, citrate cycle, and metabolic pathways. In addition, the purple module was enriched in cytosol, mitochondrion, protein binding, structural constituent of ribosome, rRNA processing, oxidative phosphorylation, ribosome, and non-alcoholic fatty liver disease. Finally, 10 hub genes in the cyan module [actin related protein 3 (ACTR3), ADP ribosylation factor 4 (ARF4), progesterone receptor membrane component 1 (PGRMC1), ribosomal protein S23 (RPS23), mannose-6-phosphate receptor (M6PR), WD repeat domain 12 (WDR12), synaptosome associated protein 23 (SNAP23), actin related protein 2 (ACTR2), siah E3 ubiquitin protein ligase 1 (SIAH1), and glomulin (GLMN)] and 2 hub genes in the purple module [proteasome 20S subunit alpha 3 (PSMA3) and ribosomal protein S27 like (RPS27L)] were identified. Hub genes were validated through quantitative real-time polymerase chain reaction. In summary, 12 hub genes were identified in two modules that were associated with HO. These hub genes could provide new biomarkers, therapeutic ideas, and targets in HO.

Data Integration in Poplar: 'Omics Layers and Integration Strategies

Populus trichocarpa is an important biofuel feedstock that has been the target of extensive research and is emerging as a model organism for plants, especially woody perennials. This research has generated several large ‘omics datasets. However, only few studies in Populus have attempted to integrate various data types. This review will summarize various ‘omics data layers, focusing on their application in Populus species. Subsequently, network and signal processing techniques for the integration and analysis of these data types will be discussed, with particular reference to examples in Populus.

Prediction of VEGF-C as a Key Target of Pure Total Flavonoids From Citrus Against NAFLD in Mice via Network Pharmacology

Pure total flavonoids from Citrus (PTFC) effectively reduce the symptoms of non-alcoholic fatty liver disease (NAFLD). Our previous microarray analysis uncovered the alterations of important signaling pathways in the treatment of NAFLD with PTFC. However, the underlying core genes that might be targeted by PTFC, which play important roles in the progression of NALFD are yet to be identified. In this study, we predicted the vascular endothelial growth factor-C (VEGF-C) as potential key molecular target of PTFC against NAFLD via network pharmacology analysis. The network pharmacology approach presented here provided important clues for understanding the mechanisms of PTFC treatment in the development of NAFLD.

Identification of epistasis loci underlying rice flowering time by controlling population stratification and polygenic effect

Flowering time is an important agronomic trait, attributed by multiple genes, gene-gene interactions and environmental factors. Population stratification and polygenic effects might confound genetic effects of the causal loci underlying this complex trait. We proposed a two-step approach for detecting epistasis interactions underlying rice flowering time by accounting population structure and polygenic effects. Simulation studies showed that the approach used in this study performs better than classical and PC-linear approaches in terms of powers and false discovery rates in the case of population stratification and polygenic effects. Whole genome epistasis analyses identified 589 putative genetic interactions for flowering time. Eighteen of these interactions are located within 10 kilobases of regions of known protein-protein interactions. Thirty-seven SNPs near to twenty-five genes involve in rice or/and Arabidopsis (orthologue) flowering pathway. Bioinformatics analysis showed that 66.55% pairwise genes of the identified interactions (392 out of the 589 interactions) have similarity in various genomic features. Moreover, significant numbers of detected epistatic genes have high expression in different floral tissues. Our findings highlight the importance of epistasis analysis by controlling population stratification and polygenic effect and provided novel insights into the genetic architecture of rice flowering which could assist breeding programmes.

Age-related but not longevity-related genes are found by weighted gene co-expression network analysis in the peripheral blood cells of humans

Human lifespan is determined by genetic and environmental factors. Potential longevity genes are neither specific nor reproducible, and longevity-related genes are constantly confused with age-related genes. To distinguish specific age- and longevity-related genes, we analyzed a Gene Expression Omnibus (GEO) dataset established by the Leiden Longevity Study. The individuals were classified into longevity (mean age, 93.4 ± 3.0 years), longevity offspring (60.8 ± 6.1) and control (61.9 ± 6.9) groups. The series matrix files were downloaded, and average expression values were calculated. Differentially expressed genes (DEGs) between longevity and control groups and those between longevity and their offspring were identified by GEO2R online. A total of 507 longevity- and 755 age-related DEGs were visualized using a Venn diagram. Weighted gene co-expression network analysis (WGCNA) was performed on the longevity- and age-related DEGs. Age-related color modules and genes were identified. However, no longevity-related modules or genes were found. The green module, with 46 age-related DEGs, was the most biologically significant to age and aging. Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, and protein-protein interaction pathway analyses were conducted on these 46 DEGs, which are mainly enriched in B cell activation and receptor signaling pathways. CR2, VPREB3, MS4A1 and CCR6 were considered the most crucial candidate genes for aging.

WISH-R- a fast and efficient tool for construction of epistatic networks for complex traits and diseases

Background

Genetic epistasis is an often-overlooked area in the study of the genomics of complex traits. Genome-wide association studies are a useful tool for revealing potential causal genetic variants, but in this context, epistasis is generally ignored. Data complexity and interpretation issues make it difficult to process and interpret epistasis. As the number of interaction grows exponentially with the number of variants, computational limitation is a bottleneck. Gene Network based strategies have been successful in integrating biological data and identifying relevant hub genes and pathways related to complex traits. In this study, epistatic interactions and network-based analysis are combined in the Weighted Interaction SNP hub (WISH) method and implemented in an efficient and easy to use R package.

Results

The WISH R package (WISH-R) was developed to calculate epistatic interactions on a genome-wide level based on genomic data. It is easy to use and install, and works on regular genomic data. The package filters data based on linkage disequilibrium and calculates epistatic interaction coefficients between SNP pairs based on a parallelized efficient linear model and generalized linear model implementations. Normalized epistatic coefficients are analyzed in a network framework, alleviating multiple testing issues and integrating biological signal to identify modules and pathways related to complex traits. Functions for visualizing results and testing runtimes are also provided.

Conclusion

The WISH-R package is an efficient implementation for analyzing genome-wide epistasis for complex diseases and traits. It includes methods and strategies for analyzing epistasis from initial data filtering until final data interpretation. WISH offers a new way to analyze genomic data by combining epistasis and network based analysis in one method and provides options for visualizations. This alleviates many of the existing hurdles in the analysis of genomic interactions.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2291-2) contains supplementary material, which is available to authorized users.

Identification of genome-wide SNP-SNP interactions associated with important traits in chicken

BACKGROUND

In addition to additive genetic effects, epistatic interactions can play key roles in the control of phenotypic variation of traits of interest. In the current study, 475 male birds from lean and fat chicken lines were utilized as a resource population to detect significant epistatic effects associated with growth and carcass traits.

RESULTS

A total of 421 significant epistatic effects were associated with testis weight (TeW), from which 11 sub-networks (Sub-network1 to Sub-network11) were constructed. In Sub-network1, which was the biggest network, there was an interaction between GGA21 and GGAZ. Three genes on GGA21 (SDHB, PARK7 and VAMP3) and nine genes (AGTPBP1, CAMK4, CDC14B, FANCC, FBP1, GNAQ, PTCH1, ROR2 and STARD4) on GGAZ that might be potentially important candidate genes for testis growth and development were detected based on the annotated gene function. In Sub-network2, there was a SNP on GGA19 that interacted with 8 SNPs located on GGA10. The SNP (Gga_rs15834332) on GGA19 was located between C-C motif chemokine ligand 5 (CCL5) and MIR142. There were 32 Refgenes on GGA10, including TCF12 which is predicted to be a target gene of miR-142-5p. We hypothesize that miR-142-5p and TCF12 may interact with one another to regulate testis growth and development. Two genes (CDH12 and WNT8A) in the same cadherin signaling pathway were implicated as potentially important genes in the control of metatarsus circumference (MeC). There were no significant epistatic effects identified for the other carcass and growth traits, e.g. heart weight (HW), liver weight (LW), spleen weight (SW), muscular and glandular stomach weight (MGSW), carcass weight (CW), body weight (BW1, BW3, BW5, BW7), chest width (ChWi), metatarsus length (MeL).

CONCLUSIONS

The results of the current study are helpful to better understand the genetic basis of carcass and growth traits, especially for testis growth and development in broilers.

Systems genetics analysis of pharmacogenomics variation during antidepressant treatment

Selective serotonin reuptake inhibitors (SSRIs) are the most widely used antidepressants, but the efficacy of the treatment varies significantly among individuals. It is believed that complex genetic mechanisms play a part in this variation. We have used a network based approach to unravel the involved genetic components. Moreover, we investigated the potential difference in the genetic interaction networks underlying SSRI treatment response over time. We found four hub genes (ASCC3, PPARGC1B, SCHIP1 and TMTC2) with different connectivity in the initial SSRI treatment period (baseline to week 4) compared with the subsequent period (4-8 weeks after initiation), suggesting that different genetic networks are important at different times during SSRI treatment. The strongest interactions in the initial SSRI treatment period involved genes encoding transcriptional factors, and in the subsequent period genes involved in calcium homeostasis. In conclusion, we suggest a difference in genetic interaction networks between initial and subsequent SSRI response.

Integrative Analysis of Genomics and Transcriptome Data to Identify Potential Functional Genes of BMDs in Females

Osteoporosis is known to be highly heritable. However, to date, the findings from more than 20 genome-wide association studies (GWASs) have explained less than 6% of genetic risks. Studies suggest that the missing heritability data may be because of joint effects among genes. To identify novel heritability for osteoporosis, we performed a system-level study on bone mineral density (BMD) by weighted gene coexpression network analysis (WGCNA), using the largest GWAS data set for BMD in the field, Genetic Factors for Osteoporosis Consortium (GEFOS-2), and a transcriptomic gene expression data set generated from transiliac bone biopsies in women. A weighted gene coexpression network was generated for 1574 genes with GWAS nominal evidence of association (p ≤ 0.05) based on dissimilarity measurement on the expression data. Twelve distinct gene modules were identified, and four modules showed nominally significant associations with BMD (p ≤ 0.05), but only one module, the yellow module, demonstrated a good correlation between module membership (MM) and gene significance (GS), suggesting that the yellow module serves an important biological role in bone regulation. Interestingly, through characterization of module content and topology, the yellow module was found to be significantly enriched with contractile fiber part (GO:044449), which is widely recognized as having a close relationship between muscle and bone. Furthermore, detailed submodule analyses of important candidate genes (HOMER1, SPTBN1) by all edges within the yellow module implied significant enrichment of functional connections between bone and cytoskeletal protein binding. Our study yielded novel information from system genetics analyses of GWAS data jointly with transcriptomic data. The findings highlighted a module and several genes in the model as playing important roles in the regulation of bone mass in females, which may yield novel insights into the genetic basis of osteoporosis. © 2016 American Society for Bone and Mineral Research.

Multi-omic data integration and analysis using systems genomics approaches: methods and applications in animal production, health and welfare

In the past years, there has been a remarkable development of high-throughput omics (HTO) technologies such as genomics, epigenomics, transcriptomics, proteomics and metabolomics across all facets of biology. This has spearheaded the progress of the systems biology era, including applications on animal production and health traits. However, notwithstanding these new HTO technologies, there remains an emerging challenge in data analysis. On the one hand, different HTO technologies judged on their own merit are appropriate for the identification of disease-causing genes, biomarkers for prevention and drug targets for the treatment of diseases and for individualized genomic predictions of performance or disease risks. On the other hand, integration of multi-omic data and joint modelling and analyses are very powerful and accurate to understand the systems biology of healthy and sustainable production of animals. We present an overview of current and emerging HTO technologies each with a focus on their applications in animal and veterinary sciences before introducing an integrative systems genomics framework for analysing and integrating multi-omic data towards improved animal production, health and welfare. We conclude that there are big challenges in multi-omic data integration, modelling and systems-level analyses, particularly with the fast emerging HTO technologies. We highlight existing and emerging systems genomics approaches and discuss how they contribute to our understanding of the biology of complex traits or diseases and holistic improvement of production performance, disease resistance and welfare.

Review: High-performance computing to detect epistasis in genome scale data sets

It is becoming clear that most human diseases have a complex etiology that cannot be explained by single nucleotide polymorphisms (SNPs) or simple additive combinations; the general consensus is that they are caused by combinations of multiple genetic variations. The limited success of some genome-wide association studies is partly a result of this focus on single genetic markers. A more promising approach is to take into account epistasis, by considering the association of multiple SNP interactions with disease. However, as genomic data continues to grow in resolution, and genome and exome sequencing become more established, the number of combinations of variants to consider increases rapidly. Two potential solutions should be considered: the use of high-performance computing, which allows us to consider a larger number of variables, and heuristics to make the solution more tractable, essential in the case of genome sequencing. In this review, we look at different computational methods to analyse epistatic interactions within disease-related genetic data sets created by microarray technology. We also review efforts to use epistatic analysis results to produce biomarkers for diagnostic tests and give our views on future directions in this field in light of advances in sequencing technology and variants in non-coding regions.

Systems genetics of obesity in an F2 pig model by genome-wide association, genetic network, and pathway analyses

Obesity is a complex condition with world-wide exponentially rising prevalence rates, linked with severe diseases like Type 2 Diabetes. Economic and welfare consequences have led to a raised interest in a better understanding of the biological and genetic background. To date, whole genome investigations focusing on single genetic variants have achieved limited success, and the importance of including genetic interactions is becoming evident. Here, the aim was to perform an integrative genomic analysis in an F2 pig resource population that was constructed with an aim to maximize genetic variation of obesity-related phenotypes and genotyped using the 60K SNP chip. Firstly, Genome Wide Association (GWA) analysis was performed on the Obesity Index to locate candidate genomic regions that were further validated using combined Linkage Disequilibrium Linkage Analysis and investigated by evaluation of haplotype blocks. We built Weighted Interaction SNP Hub (WISH) and differentially wired (DW) networks using genotypic correlations amongst obesity-associated SNPs resulting from GWA analysis. GWA results and SNP modules detected by WISH and DW analyses were further investigated by functional enrichment analyses. The functional annotation of SNPs revealed several genes associated with obesity, e.g., NPC2 and OR4D10. Moreover, gene enrichment analyses identified several significantly associated pathways, over and above the GWA study results, that may influence obesity and obesity related diseases, e.g., metabolic processes. WISH networks based on genotypic correlations allowed further identification of various gene ontology terms and pathways related to obesity and related traits, which were not identified by the GWA study. In conclusion, this is the first study to develop a (genetic) obesity index and employ systems genetics in a porcine model to provide important insights into the complex genetic architecture associated with obesity and many biological pathways that underlie it.