Obesity is a complex metabolic condition in strong association with various diseases, like type 2 diabetes, resulting in major public health and economic implications. Obesity is the result of ...environmental and genetic factors and their interactions, including genome-wide genetic interactions. Identification of co-expressed and regulatory genes in RNA extracted from relevant tissues representing lean and obese individuals provides an entry point for the identification of genes and pathways of importance to the development of obesity. The pig, an omnivorous animal, is an excellent model for human obesity, offering the possibility to study in-depth organ-level transcriptomic regulations of obesity, unfeasible in humans. Our aim was to reveal adipose tissue co-expression networks, pathways and transcriptional regulations of obesity using RNA Sequencing based systems biology approaches in a porcine model.
We selected 36 animals for RNA Sequencing from a previously created F2 pig population representing three extreme groups based on their predicted genetic risks for obesity. We applied Weighted Gene Co-expression Network Analysis (WGCNA) to detect clusters of highly co-expressed genes (modules). Additionally, regulator genes were detected using Lemon-Tree algorithms.
WGCNA revealed five modules which were strongly correlated with at least one obesity-related phenotype (correlations ranging from -0.54 to 0.72, P < 0.001). Functional annotation identified pathways enlightening the association between obesity and other diseases, like osteoporosis (osteoclast differentiation, P = 1.4E-7), and immune-related complications (e.g. Natural killer cell mediated cytotoxity, P = 3.8E-5; B cell receptor signaling pathway, P = 7.2E-5). Lemon-Tree identified three potential regulator genes, using confident scores, for the WGCNA module which was associated with osteoclast differentiation: CCR1, MSR1 and SI1 (probability scores respectively 95.30, 62.28, and 34.58). Moreover, detection of differentially connected genes identified various genes previously identified to be associated with obesity in humans and rodents, e.g. CSF1R and MARC2.
To our knowledge, this is the first study to apply systems biology approaches using porcine adipose tissue RNA-Sequencing data in a genetically characterized porcine model for obesity. We revealed complex networks, pathways, candidate and regulatory genes related to obesity, confirming the complexity of obesity and its association with immune-related disorders and osteoporosis.
RNA-seq is a useful next-generation sequencing (NGS) technology that has been widely used to understand mammalian transcriptome architecture and function. In this study, a breed-specific RNA-seq ...experiment was utilized to detect putative single nucleotide polymorphisms (SNPs) in liver tissue of young bulls of the Polish Red, Polish Holstein-Friesian (HF) and Hereford breeds, and to understand the genomic variation in the three cattle breeds that may reflect differences in production traits.
The RNA-seq experiment on bovine liver produced 107,114,4072 raw paired-end reads, with an average of approximately 60 million paired-end reads per library. Breed-wise, a total of 345.06, 290.04 and 436.03 million paired-end reads were obtained from the Polish Red, Polish HF, and Hereford breeds, respectively. Burrows-Wheeler Aligner (BWA) read alignments showed that 81.35%, 82.81% and 84.21% of the mapped sequencing reads were properly paired to the Polish Red, Polish HF, and Hereford breeds, respectively. This study identified 5,641,401 SNPs and insertion and deletion (indel) positions expressed in the bovine liver with an average of 313,411 SNPs and indel per young bull. Following the removal of the indel mutations, a total of 195,3804, 152,7120 and 205,3184 raw SNPs expressed in bovine liver were identified for the Polish Red, Polish HF, and Hereford breeds, respectively. Breed-wise, three highly reliable breed-specific SNP-databases (SNP-dbs) with 31,562, 24,945 and 28,194 SNP records were constructed for the Polish Red, Polish HF, and Hereford breeds, respectively. Using a combination of stringent parameters of a minimum depth of ≥10 mapping reads that support the polymorphic nucleotide base and 100% SNP ratio, 4,368, 3,780 and 3,800 SNP records were detected in the Polish Red, Polish HF, and Hereford breeds, respectively. The SNP detections using RNA-seq data were successfully validated by kompetitive allele-specific PCR (KASPTM) SNP genotyping assay. The comprehensive QTL/CG analysis of 110 QTL/CG with RNA-seq data identified 20 monomorphic SNP hit loci (CARTPT, GAD1, GDF5, GHRH, GHRL, GRB10, IGFBPL1, IGFL1, LEP, LHX4, MC4R, MSTN, NKAIN1, PLAG1, POU1F1, SDR16C5, SH2B2, TOX, UCP3 and WNT10B) in all three cattle breeds. However, six SNP loci (CCSER1, GHR, KCNIP4, MTSS1, EGFR and NSMCE2) were identified as highly polymorphic among the cattle breeds.
This study identified breed-specific SNPs with greater SNP ratio and excellent mapping coverage, as well as monomorphic and highly polymorphic putative SNP loci within QTL/CGs of bovine liver tissue. A breed-specific SNP-db constructed for bovine liver yielded nearly six million SNPs. In addition, a KASPTM SNP genotyping assay, as a reliable cost-effective method, successfully validated the breed-specific putative SNPs originating from the RNA-seq experiments.
In cattle, the in vitro production (IVP) of embryos is becoming more relevant than embryos produced in vivo, i.e. after multiple ovulation and embryo transfer (MOET). However, the effects of IVP on ...the developmental programming of specific organs in the postnatal calves are yet unknown. Previously, we reported an epigenomic and transcriptomic profile of the hypothalamus–pituitary–testicular axis compatible with its earlier activation in IVP calves compared to MOET animals. Here, we studied the hepatic and muscular epigenome and transcriptome of those same male dairy calves (n = 4 per group). Tissue samples from liver and semitendinosus muscle were obtained at 3 months of age, and the extracted gDNA and RNA were sequenced through whole-genome bisulfite sequencing and RNA-sequencing, respectively. Next, bioinformatic analyses determined differentially methylated cytosines or differentially expressed genes false discovery rate (FDR) < 0.05 for each Omic dataset; and nonparametrically combined genes (NPCG) for both integrated omics (P < 0.05). KEGG pathways enrichment analysis showed that NPCG upregulated in the liver and the muscle of the IVP calves were involved in oxidative phosphorylation and the tricarboxylic acid cycle. In contrast, ribosome and translation were upregulated in the liver but downregulated in the muscle of the IVP calves compared to the MOET calves (FDR < 0.05). A model considering the effect of the methylation levels and the group on the expression of all the genes involved in these pathways confirmed these findings. In conclusion, the multiomics data integration approach indicated an altered hepatic and muscular energy regulation in phenotypically normal IVP calves compared to MOET calves. Summary Sentence Transcriptomic and epigenomic results suggest that aerobic respiration was upregulated in both liver and muscle, while protein synthesis was increased in the liver but downregulated in the muscle of in vitro produced calves compared to in vivo counterparts. Graphical Abstract
In cattle, several calves born after IVP (“in vitro” embryo production) present similar birthweight to those generated after MOET (multiple ovulation and embryo transfer). However, the underlying ...molecular patterns in organs involved in the developmental process are unknown and could indicate physiological programming. The objectives of this study were: (1) to compare epigenomic and transcriptomic modifications in the hypothalamus, pituitary, gonadal and adrenal organs between 3 months old ovum pick‐up‐IVP and MOET male calves (n = 4 per group) and (2) to use blood epigenomic data to proxy methylation of the inner organs. Extracted gDNA and RNA were sequenced through whole‐genome bisulfite sequencing and RNA sequencing, respectively. Next, bioinformatic analyses determined differentially methylated cytosines (DMC) and differentially expressed genes (DEG) (FDR < 0.05) in IVP versus MOET samples and the KEGG pathways that were overrepresented by genes associated with DMC or DEG (FDR < 0.1). Pathways related to hypothalamus, pituitary, gonadal (HPG) axis activation (GnRH secretion in the hypothalamus, GnRH signaling in the pituitary, and steroidogenesis in the testicle) were enriched in IVP calves. Modeling the effect of the methylation levels and the group on the expression of all the genes involved in these pathways confirmed their upregulation in HPG organs in IVP calves. The application of the DIABLO method allowed the identification of 15 epigenetic and five transcriptomic biomarkers, which were able to predict the embryo origin using the epigenomic data from the blood. In conclusion, the use of an integrated epigenomic‐transcriptomic approach suggested an early activation of the HPG axis in male IVP calves compared to MOET counterparts, and the identification of potential biomarkers allowed the use of blood samples to proxy methylation levels of the relevant internal organs.
Hypothyroidism is associated with obesity, and thyroid hormones are involved in the regulation of body composition, including fat mass. Genome-wide association studies (GWAS) in adults have ...identified 19 and 6 loci associated with plasma concentrations of thyroid stimulating hormone (TSH) and free thyroxine (fT4), respectively.
This study aimed to identify and characterize genetic variants associated with circulating TSH and fT4 in Danish children and adolescents and to examine whether these variants associate with obesity.
Genome-wide association analyses of imputed genotype data with fasting plasma concentrations of TSH and fT4 from a population-based sample of Danish children, adolescents, and young adults, and a group of children, adolescents, and young adults with overweight and obesity were performed (N = 1,764, mean age = 12.0 years range 2.5-24.7). Replication was performed in additional comparable samples (N = 2,097, mean age = 11.8 years 1.2-22.8). Meta-analyses, using linear additive fixed-effect models, were performed on the results of the discovery and replication analyses.
No novel loci associated with TSH or fT4 were identified. Four loci previously associated with TSH in adults were confirmed in this study population (PDE10A (rs2983511: β = 0.112SD, p = 4.8 ∙ 10-16), FOXE1 (rs7847663: β = 0.223SD, p = 1.5 ∙ 10-20), NR3C2 (rs9968300: β = 0.194SD), p = 2.4 ∙ 10-11), VEGFA (rs2396083: β = 0.088SD, p = 2.2 ∙ 10-10)). Effect sizes of variants known to associate with TSH or fT4 in adults showed a similar direction of effect in our cohort of children and adolescents, 11 of which were associated with TSH or fT4 in our study (p<0.0002). None of the TSH or fT4 associated SNPs were associated with obesity in our cohort, indicating no pleiotropic effects of these variants on obesity.
In a group of Danish children and adolescents, four loci previously associated with plasma TSH concentrations in adults, were associated with plasma TSH concentrations in children, suggesting comparable genetic determinants of thyroid function in adults and children.
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.
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.
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.
DNA methylation of different gene components, including different exons and introns, or different lengths of exons and introns is associated with differences in gene expression. To investigate the ...methylation of porcine gene components associated with the boar taint (BT) trait, this study used reduced representation bisulfite sequencing (RRBS) data from nine porcine testis samples in three BT groups (low, medium and high BT). The results showed that the methylation levels of the first exons and first introns were lower than those of the other exons and introns. The first exons/introns of CpG island regions had even lower levels of methylation. A total of 123 differentially methylated promoters (DMPs), 194 differentially methylated exons (DMEs) and 402 differentially methylated introns (DMIs) were identified, of which 80 DMPs (DMP-CpGis), 112 DMEs (DME-CpGis) and 166 DMIs (DMI-CpGis) were discovered in CpG islands. Importantly,
contained one each of DMP, DME, DMI, DMP-CpGi, DME-CpGi and DMI-CpGi. Gene-GO term relationships and pathways analysis showed DMP-CpGi-related genes are mainly involved in methylation-related biological functions. In addition, gene-gene interaction networks consisted of nodes that were hypo-methylated
, hypo-methylated
, hypo-methylated
, hyper-methylated
, hyper-methylated
and hyper-methylated
. They could be used as candidate biomarkers for reducing boar taint in pigs, after further validation in large cohorts.
Feed efficiency (FE) is an economically important trait. Thus, reliable predictors would help to reduce the production cost and provide sustainability to the pig industry. We carried out ...metabolome-transcriptome integration analysis on 40 purebred Duroc and Landrace uncastrated male pigs to identify potential gene-metabolite interactions and explore the molecular mechanisms underlying FE. To this end, we applied untargeted metabolomics and RNA-seq approaches to the same animals. After data quality control, we used a linear model approach to integrate the data and find significant differently correlated gene-metabolite pairs separately for the breeds (Duroc and Landrace) and FE groups (low and high FE) followed by a pathway over-representation analysis. We identified 21 and 12 significant gene-metabolite pairs for each group. The valine-leucine-isoleucine biosynthesis/degradation and arginine-proline metabolism pathways were associated with unique metabolites. The unique genes obtained from significant metabolite-gene pairs were associated with sphingolipid catabolism, multicellular organismal process, cGMP, and purine metabolic processes. While some of the genes and metabolites identified were known for their association with FE, others are novel and provide new avenues for further research. Further validation of genes, metabolites, and gene-metabolite interactions in larger cohorts will elucidate the regulatory mechanisms and pathways underlying FE.
Resorption of bones and cartilage coupled with structural changes in the inflamed joints are the major hallmark of rheumatoid arthritis (RA). Genetic polymorphisms in pro-inflammatory interleukins ...(ILs) appear to play an important role in the susceptibility towards progressive RA. We therefore aimed to investigate the association of single nucleotide polymorphisms (SNP), present in the hotspot coding/promoter regions of IL-6, −17 and −18, with RA susceptibility or severity in a larger study cohort from Pakistan together with finding clues as to how a functional SNP impacts the predisposition towards RA. TaqMan SNP genotyping approach was first used to assess IL-6 (rs1800795), IL-17 F (rs763780), IL-17A (rs2275913), and IL-18 (rs1946518) polymorphisms in 310 subjects (150 RA and 160 control). Molecular dynamic simulations (MDS) of wild- and mutant-type IL-17A with corresponding receptor were thereafter performed using AMBER-16; Chimera 1.13 was used for analyses. Our results showed the association of two SNPs, namely IL-6 − 174 G/C allelic (OR = 0.960, 95% CI = 0.929-0.992, p = .009) and IL-17 F 7488 T/C allelic (OR = 0.907, 95%CI = 0.861-0.954, p = .000) with increased RA risk in Pakistani subjects. When mapped, IL-17 F 7488 T/C was found involved in His
161
→Arg
161
change near the C-terminus of IL-17 F. Comparative MDS revealed enhanced stability of the mutant hence advocating a potential role of IL-17F functional SNP in RA susceptibility and/or severity. This study provides a novel structural insight for SNP-derived functional mutation and its overall impact on binding with heterotrimeric receptor complex of IL-17 receptor thereby opening new avenues for understanding the biochemical basis of the disease.