Expression Profiles of the Insulin-like Growth Factor System Components in Liver Tissue during Embryonic and Postnatal Growth of Erhualian and Yorkshire Reciprocal Cross F1 Pigs

DNA methylation analysis to differentiate reference, breed, and parent-of-origin effects in the bovine pangenome era

BACKGROUND

Most DNA methylation studies have used a single reference genome with little attention paid to the bias introduced due to the reference chosen. Reference genome artifacts and genetic variation, including single nucleotide polymorphisms (SNPs) and structural variants (SVs), can lead to differences in methylation sites (CpGs) between individuals of the same species. We analyzed whole-genome bisulfite sequencing data from the fetal liver of Angus (Bos taurus taurus), Brahman (Bos taurus indicus), and reciprocally crossed samples. Using reference genomes for each breed from the Bovine Pangenome Consortium, we investigated the influence of reference genome choice on the breed and parent-of-origin effects in methylome analyses.

RESULTS

Our findings revealed that ∼75% of CpG sites were shared between Angus and Brahman, ∼5% were breed specific, and ∼20% were unresolved. We demonstrated up to ∼2% quantification bias in global methylation when an incorrect reference genome was used. Furthermore, we found that SNPs impacted CpGs 13 times more than other autosomal sites (P < $5 \times {10}^{ - 324}$) and SVs contained 1.18 times (P < $5 \times {10}^{ - 324}$) more CpGs than non-SVs. We found a poor overlap between differentially methylated regions (DMRs) and differentially expressed genes (DEGs) and suggest that DMRs may be impacting enhancers that target these DEGs. DMRs overlapped with imprinted genes, of which 1, DGAT1, which is important for fat metabolism and weight gain, was found in the breed-specific and sire-of-origin comparisons.

CONCLUSIONS

This work demonstrates the need to consider reference genome effects to explore genetic and epigenetic differences accurately and identify DMRs involved in controlling certain genes.

A genomic estimated breeding value-assisted reduction method of single nucleotide polymorphism sets: a novel approach for determining the cutoff thresholds in genome-wide association studies and best linear unbiased prediction

Traditionally, the p-value is the criterion for the cutoff threshold to determine significant markers in genome-wide association studies (GWASs). Choosing the best subset of markers for the best linear unbiased prediction (BLUP) for improved prediction ability (PA) has become an interesting issue. However, when dealing with many traits having the same marker information, the p-values' themselves cannot be used as an obvious solution for having a confidence in GWAS and BLUP. We thus suggest a genomic estimated breeding value-assisted reduction method of the single nucleotide polymorphism (SNP) set (GARS) to address these difficulties. GARS is a BLUP-based SNP set decision presentation. The samples were Landrace pigs and the traits used were back fat thickness (BF) and daily weight gain (DWG). The prediction abilities (PAs) for BF and DWG for the entire SNP set were 0.8 and 0.8, respectively. By using the correlation between genomic estimated breeding values (GEBVs) and phenotypic values, selecting the cutoff threshold in GWAS and the best SNP subsets in BLUP was plausible as defined by GARS method. 6,000 SNPs in BF and 4,000 SNPs in DWG were considered as adequate thresholds. Gene Ontology (GO) analysis using the GARS results of the BF indicated neuron projection development as the notable GO term, whereas for the DWG, the main GO terms were nervous system development and cell adhesion.

The super-enhancer repertoire in porcine liver

The transcriptional initiation of genes is inextricably bound with the functions of cis-regulatory sequences. The pig is one of the most important livestock species and an ideal animal model for biomedical studies. At the same time, the liver is a critical organ with diverse and complex metabolic functions. Here, we performed Cleavage Under Targets and Tagmentation (CUT&Tag) coupled with high-throughput sequencing to profile the chromatin landscape of histone H3 lysine 27 acetylation (H3K27ac), histone H3 lysine 4 monomethylation (H3K4me1), and CCAAT enhancer-binding protein β (C-EBPβ) in the 70-d-old porcine liver, compared the different profiles among the three markers and their associated stitched-enhancers by stitching and sorting the peaks within 12.5 kb (Pott and Lieb, 2015) and generated the porcine liver-specific super-enhancers (SEs) by the combination of three markers. Compared to typical enhancers (TEs) and other stitched-enhancers, liver-specific SEs showed a higher density of cis-motifs and SNPs, which may recruit more tissue-specific vital TFs. The expression profiles in fetal and 70-d-old pigs proved that a large proportion of SE-associated genes were up-regulated and were more related to hepatic metabolisms and detoxification pathways. Our results illustrated the difference and connection among promoter and enhancer markers, identified the features of liver SEs and their associated genes, and provided novel insight into cis-element identification, function, and liver transcriptional regulation.

Transcriptomic profiling of Gh/Igf system reveals a prompted tissue-specific differentiation and novel hypoxia responsive genes in gilthead sea bream

A customized PCR-array was used for the simultaneous gene expression of the Gh/Igf system and related markers of muscle growth, and lipid and energy metabolism during early life stages of gilthead sea bream (60–127 days posthatching). Also, transcriptional reprogramming by mild hypoxia was assessed in fingerling fish with different history trajectories on O₂ availability during the same time window. In normoxic fish, the expression of almost all the genes in the array varied over time with a prompted liver and muscle tissue-specific differentiation, which also revealed temporal changes in the relative expression of markers of the full gilthead sea bream repertoire of Gh receptors, Igfs and Igf-binding proteins. Results supported a different contribution through development of ghr and igf subtypes on the type of action of GH via systemic or direct effects at the local tissue level. This was extensive to Igfbp1/2/4 and Igfbp3/5/6 clades that clearly evolved through development as hepatic and muscle Igfbp subtypes, respectively. This trade-off is however very plastic to cope changes in the environment, and ghr1 and igfbp1/3/4/5 emerged as hypoxic imprinting genes during critical early developmental windows leading to recognize individuals with different history trajectories of oxygen availability and metabolic capabilities later in life.

Long-distance chromatin interaction of IGF1 during embryonic and postnatal development in the liver of Sus scrofa

The dynamics of chromatin have been the focus of studies aimed at characterizing gene regulation. Among various chromosome conformation capture methods, 4C-seq is a powerful technique to identify genome-wide interactions with a single locus of interest. Insulin-like growth factor 1 (IGF1) is a member of the somatotropin axis that plays a significant role in cell proliferation and growth. Determining the IGF1-involved genome-wide chromatin interaction profile at different growth stages not only is important for understanding IGF1 transcriptional regulation but also provides a representation of genome-wide chromatin transformation during development. Using the IGF1 promoter as a "bait", we identified genome-wide interactomes of embryonic (E70) and postnatal (P1 and P70) pig liver cells by 4C-seq. The IGF1 promoter interactomes varied significantly among the three developmental stages. The most active chromatin interaction was observed in the P1 stage, while the highest interaction variability was observed in the P70 stage. The identified 4C sites were enriched around transcription start sites, CpG sites and functional pig QTLs. In addition, the genes located in the interacting regions and the involved pathways were also analysed. Overall, our work reveals a distinct long-distance regulatory pattern in pig liver during development for the first time, and the identified interacting sites and genes may serve as candidate targets in further transcriptional mechanism studies and effective molecular markers for functional traits.

Differential Expression of IGF1, IGFBP5, MSTN, and MYH1 Across Different Age Classes in American Quarter Horses

The objective of this study was to determine the influence of age on expression of insulin-like growth factor 1 (IGF1), insulin-like growth factor binding protein (IGFBP5), myostatin (MSTN), and myosin (MYH1) genes which are related to growth and muscle development in the American Quarter Horse. Thus, horses (n = 10) from weanling, yearling, 2-, 3-, and 10-year-old age classes were sampled and gene expression was assessed by RT-qPCR. ΔC_T was calculated using the hypoxanthine-guanine phosphoribosyltransferase gene as an internal normalizer. The generalized linear model was used to determine differentially expressed genes, by pairwise comparison between ages. Among technical replicates, the coefficient of variation ranged from 1.0 to 5.2% and was lower than the variation observed between biological replicates (2.1-12.9%). IGF1 demonstrated significantly lower expression in the 3-year-old age class than in weanlings and yearlings, but the 10-year-old age class displayed a significantly higher level than 2- and 3-year-old age classes. Expression of IGFBP5 was highest in weanlings compared with all other age classes. Expression of MSTN was significantly higher in weanlings than in other age classes, whereas 10-year-old horses had an intermediate level of expression, but significantly different from yearlings, 2- and 3-year-old fillies. Finally, expression of MYH1 was lower in 2- and 10-year-old horses than in weanlings and yearlings, whereas 3-year-old fillies demonstrated an intermediate level of expression. Differential expression patterns observed in this preliminary study provide insight into the physiological changes occurring throughout the life span of horses. These patterns could also help explain the variation in performance and endurance between individuals at different developmental stages.

Transcriptome comparison between prenatal and postnatal Large White livers identifies differences in the expression level of genes related to metabolism and postnatal growth

The current study examined the liver transcriptomic profiles of the Large White different in developmental periods. It was performed on pigs of two developmental stages: 70-day fetus (P70) and 70-day piglets (D70). The objective of the study was to identify genes associated with Large White liver lipid metabolism, growth and development. We sequenced eight sRNA libraries of liver samples from four Large White at P70 and D70 respectively. We totally obtained 19,202 genes. 4916 of them were found to be differentially expressed (DEGs) (p < 0.05, fold change ≥ 1), of which 2502 were up-regulated and 2414 were down-regulated. GO enrichment and KEGG pathway analysis indicated that ACACA, ACADM, ACAA2 and HADH were simultaneously enriched in diverse pathways related to lipid metabolism, and so they were considered to be the promising candidate genes which could affect the porcine liver lipid metabolism. Notably, the gene insulin-like growth factor 1 (IGF1) which participated in somatotropic axis signaling was found to be up-regulated in D70 compared with P70. miRWalk and TargetScan softwares were used to screen the miRNAs which bound to the 3' untranslated region (3'UTR) of IGF1. After integration analysis with miRNAs sequencing data, miR-18b and miR-130b-3p were selected for further study. MiR-18b and miR-130b-3p were down-regulated in D70 compared with P70. Dual luciferase assays indicated that miR-18b and miR-130b-3p could obviously decrease (p < 0.05) the fluorescence activity of the group transfected with the wild-type vector of IGF1 3'UTR, while the relative luciferase activity of the group transfected with the mutant vector of IGF1 3'UTR did not change significantly. Taken together, it indicated that miR-18b and miR-130b-3p could target IGF1 directly.

Somatotropic Axis Regulation Unravels the Differential Effects of Nutritional and Environmental Factors in Growth Performance of Marine Farmed Fishes

The Gh/Prl/Sl family has evolved differentially through evolution, resulting in varying relationships between the somatotropic axis and growth rates within and across fish species. This is due to a wide range of endogenous and exogenous factors that make this association variable throughout season and life cycle, and the present minireview aims to better define the nutritional and environmental regulation of the endocrine growth cascade over precisely defined groups of fishes, focusing on Mediterranean farmed fishes. As a result, circulating Gh and Igf-i are revitalized as reliable growth markers, with a close association with growth rates of gilthead sea bream juveniles with deficiency signs in both macro- or micro-nutrients. This, together with other regulated responses, promotes the use of Gh and Igf-i as key performance indicators of growth, aerobic scope, and nutritional condition in gilthead sea bream. Moreover, the sirtuin-energy sensors might modulate the growth-promoting action of somatotropic axis. In this scenario, transcripts of igf-i and gh receptors mirror changes in plasma Gh and Igf-i levels, with the ghr-i/ghr-ii expression ratio mostly unaltered over season. However, this ratio is nutritionally regulated, and enriched plant-based diets or diets with specific nutrient deficiencies downregulate hepatic ghr-i, decreasing the ghr-i/ghr-ii ratio. The same trend, due to a ghr-ii increase, is found in skeletal muscle, whereas impaired growth during overwintering is related to increase in the ghr-i/ghr-ii and igf-ii/igf-i ratios in liver and skeletal muscle, respectively. Overall, expression of insulin receptors and igf receptors is less regulated, though the expression quotient is especially high in the liver and muscle of sea bream. Nutritional and environmental regulation of the full Igf binding protein 1-6 repertoire remains to be understood. However, tissue-specific expression profiling highlights an enhanced and nutritionally regulated expression of the igfbp-1/-2/-4 clade in liver, whereas the igfbp-3/-5/-6 clade is overexpressed and regulated in skeletal muscle. The somatotropic axis is, therefore, highly informative of a wide-range of growth-disturbing and stressful stimuli, and multivariate analysis supports its use as a reliable toolset for the assessment of growth potentiality and nutrient deficiencies and requirements, especially in combination with selected panels of other nutritionally regulated metabolic biomarkers.