Diet-Dependent Changes of the DNA Methylome Using a Göttingen Minipig Model for Obesity.

Objective: Environmental factors can influence obesity by epigenetic mechanisms. The aim of this study was to investigate obesity-related epigenetic changes and the potential for reversal of these changes in the liver of Göttingen minipigs subjected to diet interventions. Methods: High-throughput liquid hybridization capture-based bisulfite sequencing (LHC-BS) was used to quantify the methylation status of gene promotor regions in liver tissue in three groups of male castrated Göttingen minipigs: a standard chow group (SD, N = 7); a group fed high fat/fructose/cholesterol diet (FFC, N = 10) and a group fed high fat/fructose/cholesterol diet during 7 months and reversed to standard diet for 6 months (FFC/SD, N = 12). Expression profiling by qPCR of selected metabolically relevant genes was performed in liver tissue from all pigs. Results: The pigs in the FFC diet group became morbidly obese. The FFC/SD diet did not result in a complete reversal of the body weight to the same weight as in the SD group, but it resulted in reversal of all lipid related metabolic parameters. Here we identified widespread differences in the patterning of cytosine methylation of promoters between the different feeding groups. By combining detection of differentially methylated genes with a rank-based hypergeometric overlap algorithm, we identified 160 genes showing differential methylation in corresponding promoter regions in the FFC diet group when comparing with both the SD and FFC/SD groups. As expected, this differential methylation under FFC diet intervention induced de-regulation of several metabolically-related genes involved in lipid/cholesterol metabolism, inflammatory response and fibrosis generation. Moreover, five genes, of which one is a fibrosis-related gene (MMP9), were still perturbed after diet reversion. Conclusion: Our findings highlight the potential of exploring diet-epigenome interactions for treatment of obesity.

Unraveling molecular mechanisms involved in the development of leptin resistance using the pig as a model.

The increase in obesity worldwide underlines the need for research concerning its metabolic and genetic determinants. One of the most intriguing mechanisms regarding obesity involves leptin and its signaling cascade. Leptin is a key regulator contributing to the fine-tuned crosstalk between nutrient availability and appetite signaling in the central nervous system. Owing to ethical concerns, many human tissues are not readily available and pigs can serve as a good animal model owing to their comparable anatomy, metabolism and genetics. In the present study, we utilized the pig to investigate the possible impact of increased adiposity on the development of alterations within the leptin signaling pathway. Two divergent groups of pigs (High and Low) were defined based on a high and low amount of mesenteric fat. Cortex, cerebellum, hypothalamus, mesenteric, subcutaneous and retroperitoneal fat tissues were used to study changes in expression levels of 94 mRNA transcripts related to the leptin signaling pathway using the qPCR approach. No significant differences were found at the central nervous system, whereas the expression level of STAT1 was reduced in mesenteric fat and leptin (LEP) and interleukin 6 (IL6) were shown to be consistently increased in all analyzed fat compartments from pigs with a high amount of mesenteric fat. These results could imply the onset of leptin and pro-inflammatory cytokine overexpression at early stages of obesity in the analyzed pigs without affecting any key components in the central nervous system. Thus, these pigs showing a unique leptin deregulation in adipose tissues could be a useful translational resource for studies of obesity and leptin resistance phenotypes.

A targeted genotyping approach enhances identification of variants in taste receptor and appetite/reward genes of potential functional importance for obesity-related porcine traits.

Taste receptors (TASRs) and appetite and reward (AR) mechanisms influence eating behaviour, which in turn affects food intake and risk of obesity. In a previous study, we used next generation sequencing to identify potentially functional mutations in TASR and AR genes and found indications for genetic associations between identified variants and growth and fat deposition in a subgroup of animals (n = 38) from the UNIK resource pig population. This population was created for studying obesity and obesity-related diseases. In the present study we validated results from our previous study by investigating genetic associations between 24 selected single nucleotide variants in TASR and AR gene variants and 35 phenotypes describing obesity and metabolism in the entire UNIK population (n = 564). Fifteen variants showed significant association with specific obesity-related phenotypes after Bonferroni correction. Six of the 15 genes, namely SIM1, FOS, TAS2R4, TAS2R9, MCHR2 and LEPR, showed good correlation between known biological function and associated phenotype. We verified a genetic association between potentially functional variants in TASR/AR genes and growth/obesity and conclude that the combination of identification of potentially functional variants by next generation sequencing followed by targeted genotyping and association studies is a powerful and cost-effective approach for increasing the power of genetic association studies.

Circulating let-7g is down-regulated in Bernese Mountain dogs with disseminated histiocytic sarcoma and carcinomas - a prospective study.

Cancer is a prevalent cause of mortality in Bernese mountain dogs (BMDs). Circulating microRNAs (miRNAs) are found in blood and have been identified as promising biomarkers in various neoplastic diseases in humans. In the current study, the expression profile of different types of miRNAs was investigated in healthy BMDs and BMDs with cancer. Seven healthy and six non-treated BMDs with cancer [four with disseminated histiocytic sarcomas (DHS)] were enrolled in this study. Clinical evaluations including physical examination, blood analysis, urinalysis and diagnostic imaging were performed on all dogs. Twenty-four different miRNAs were profiled from RNA isolated from whole blood preserved in PAXgene® tubes using quantitative real-time PCR (qPCR). The miRNA let-7g was significantly down-regulated in dogs with cancer (P = 0.002) and dogs with DHS (P = 0.011) compared with healthy controls. This miRNA is a known tumour suppressor and further analyses are warranted to assess its value as a non-invasive biomarker for early detection of different types of cancer in BMDs.

Alpha-smooth muscle actin and serotonin receptors 2A and 2B in dogs with myxomatous mitral valve disease.

Canine Myxomatous mitral valve disease (MMVD) is an age-related disease. Serotonin (5-HT) is implicated in the pathogenesis as locally-produced or platelet-derived. Involvement of the 5-HT2A receptor (R) and 5-HT2BR in the induction of myxomatous-mediating valvular myofibroblasts (MF) has been suggested. In an age-matched population of dogs with non-clinical and clinical MMVD, the objectives were to investigate (1) gene expression of 5-HT2AR and 5-HT2BR, (2) protein expression and spatial relationship of 5-HT2AR, 5-HT2BR and MF in the mitral valve (MV) and the cardiac anterior papillary muscle (AP) and (3) serum 5-HT concentrations. Gene expression of 5-HT2BR was significantly higher in MV and AP among dogs with clinical MMVD. This was not found for 5-HT2BR protein expression, though association of 5-HT2BR with myxomatous pathology and co-localization of 5-HT2BR and MF in MV and AP support a functional relationship, perhaps perpetuation of clinical MMVD. 5-HT2AR-expression and serum 5-HT showed no differences between groups.

Serotonin markers show altered transcription levels in an experimental pig model of mitral regurgitation.

Serotonin (5-hydroxytryptamine, 5-HT) signalling is implicated in the pathogenesis of myxomatous mitral valve disease (MMVD) through 5-HT1B receptor (R), 5-HT2AR and 5-HT2BR-induced myxomatous pathology. Based on increased tryptophan hydroxylase-1 (TPH-1) and decreased serotonin re-uptake transporter (SERT) in MMVD-affected valves, increased valvular 5-HT synthesis and decreased clearance have been suggested. It remains unknown how haemodynamic changes associated with mitral regurgitation (MR) affect 5-HT markers in the mitral valve, myocardium and circulation. Twenty-eight pigs underwent surgically induced MR or sham-operation, resulting in three MR groups: control (CON, n = 12), mild MR (mMR, n = 10) and severe MR (sMR, n = 6). The gene expression levels of 5-HT1BR, 5-HT2AR, 5-HT2BR, SERT and TPH-1 were analysed using quantitative PCR (qPCR) in the mitral valve (MV), anterior papillary muscle (AP) and left ventricle (LV). MV 5-HT2BR was also analysed with immunohistochemistry (IHC) in relation to histological lesions and valvular myofibroblasts. All 5-HTR mRNAs were up-regulated in MV compared to AP and LV (P <0.01). In contrast, SERT and TPH-1 were up-regulated in AP and LV compared to MV (P <0.05). In MV, mRNA levels were increased for 5-HT2BR (P = 0.02) and decreased for SERT (P = 0.03) in sMR vs. CON. There were no group differences in 5-HT2BR staining (IHC) but co-localisation was found with α-SMA-positive cells in 91% of all valves and with 33% of histological lesions. In LV, 5-HT1BR mRNA levels were increased in sMR vs. CON (P = 0.01). In conclusion, these data suggest that MR may affect mRNA expression of valvular 5-HT2BR and SERT, and left ventricular 5-HT1BR in some pigs.

Extensive changes in innate immune gene expression in obese Göttingen minipigs do not lead to changes in concentrations of circulating cytokines and acute phase proteins.

The usefulness of Göttingen minipigs as models for obesity and obesity-related pathologies is well established. The low-grade inflammation associated with obesity involves a range of innate immune factors; however, to our knowledge, the impact of obesity on innate immune factor expression has not been studied in Göttingen minipigs. Therefore, we studied the expression of innate immune genes in liver and adipose tissues as well as serum concentrations of cytokines and acute phase proteins in obese vs. lean Göttingen minipigs. In the liver, of 35 investigated genes, the expression of nine was significantly different in obese pigs (three up-regulated, six down-regulated). Of 33 genes in adipose tissues, obesity was associated with changed expression of 12 genes in the visceral adipose tissue (VAT) (three up-regulated), 11 in the abdominal retroperitoneal adipose tissue (RPAT) (seven of these up-regulated) and eight in the subcutaneous adipose tissue (SAT) from the neck (five of which were up-regulated). Obesity-associated expression changes were observed for three genes in all adipose tissues, namely chemokine (C-C motif) ligand 3-like 1 (up-regulated), CD200 molecule (down-regulated) and interleukin 1 receptor antagonist (up-regulated) with interleukin 1 receptor antagonist being the most highly regulated gene in both VAT and RPAT. Looking at patterns of expression across the three types of adipose tissues, obesity was associated with an increased number of acute phase proteins differentially expressed between adipose tissues and a decreased tissue-specific expression of cytokines and chemokines. In contrast to obese humans, no changes in serum concentrations of haptoglobin, C-reactive protein, serum amyloid A, tumor necrosis factor-α and interleukin 6 were found in obese Göttingen minipigs.

Expression studies of six human obesity-related genes in seven tissues from divergent pig breeds.

Obesity has reached epidemic proportions globally and has become the cause of several major health risks worldwide. Presently, more than 100 loci have been related to obesity and metabolic traits in humans by genome-wide association studies. The complex genetic architecture behind obesity has triggered a need for the development of better animal models than rodents. The pig has emerged as a very promising biomedical model to study human obesity traits. In this study, we have characterized the expression patterns of six obesity-related genes, leptin (LEP), leptin receptor (LEPR), melanocortin 4 receptor (MC4R), fat mass and obesity associated (FTO), neuronal growth regulator 1 (NEGR)1 and adiponectin (ADIPOQ), in seven obesity-relevant tissues (liver; muscle; pancreas; hypothalamus; and retroperitoneal, subcutaneous and mesenteric adipose tissues) in two pig breeds (production pigs and Göttingen minipigs) that deviate phenotypically and genetically from each other with respect to obesity traits. We observe significant differential expression for LEP, LEPR and ADIPOQ in muscle and in all three adipose tissues. Interestingly, in pancreas, LEP expression is only detected in the fat minipigs. FTO shows significant differential expression in all tissues analyzed, and NEGR1 shows significant differential expression in muscle, pancreas, hypothalamus and subcutaneous adipose tissue. The MC4R transcript can be detected only in hypothalamus. In general, the expression profiles of the investigated genes are in accordance with those observed in human studies. Our study shows that both the differences between the investigated breeds and the phenotypic state with respect to obesity/leanness play a large role for differential expression of the obesity-related genes.

Potential role of the porcine superoxide dismutase 1 (SOD1) gene in pig reproduction.

In a recent study we confirmed that QTL regions on pig chromosomes 11, 13, and 15 are associated with reproduction traits in the pig. Within these regions the genetic variation was largest on chromosome 13. The QTL region on this chromosome was therefore studied further to identify genes known to contribute to litter size. The superoxide dismutase (SOD1) gene localized at around 200 Mb in the pig (Sscrofa10) was the most obvious candidate gene. In the present study, we have cloned and sequenced the porcine SOD1 gene. The SOD1 amino acid sequence is highly conserved between human, mouse, rat, and pig. Expression studies by quantitative PCR showed differential levels of the SOD1 transcript in all tissues investigated. Sequence comparison between sows with high and low estimated breeding value (EBV) for litter size, revealed a total of eight single nucleotide polymorphisms (SNPs) in the noncoding sequence and no SNPs in the coding region. One of the intronic SNPs was genotyped in 248 sows with high and low EBV for litter size. Allele frequency differed significantly between the two group of sows indicating that polymorphism in the chromosome 13 locus has an impact on litter size. The sows homozygous for the A/A genotype conceive three piglets more compared to the A/T genotype, making this SNP a possible marker for litter size. However, this genotype was negatively correlated with other important traits under selection in the Danish pig production.

Exonization of a LINE1 fragment implicated in X-linked hypohidrotic ectodermal dysplasia in cattle.

A case of X-linked hypohidrotic ectodermal dysplasia (XHED) was identified in a family of Danish Red Holstein cattle. The ectodysplasin-signalling protein (EDA) is known to be central in the normal development of ectodermal structures, and mutations in the ectodysplasin A (EDA) gene have been reported to cause XHED. In this study, we analysed different EDA transcript variants in affected and unaffected cattle and identified a new transcript variant including a LINE1-derived pseudoexon between EDA exons 1 and 2. The 161-bp-long pseudoexon introduces a shift in reading frame and a premature stop codon early in EDA exon 2 and is probably the cause of XHED in this Danish Red Holstein family.

Refined candidate region specified by haplotype sharing for Escherichia coli F4ab/F4ac susceptibility alleles in pigs.

Infection of the small intestine by enterotoxigenic Escherichia coli F4ab/ac is a major welfare problem and financial burden for the pig industry. Natural resistance to this infection is inherited as a Mendelian recessive trait, and a polymorphism in the MUC4 gene segregating for susceptibility/resistance is presently used in a selection programme by the Danish pig breeding industry. To elucidate the genetic background involved in E. coli F4ab/ac susceptibility in pigs, a detailed haplotype map of the porcine candidate region was established. This region covers approximately 3.7 Mb. The material used for the study is a three generation family, where the founders are two Wild boars and eight Large White sows. All pigs have been phenotyped for susceptibility to F4ab/ac using an adhesion assay. Their haplotypes are known from segregation analysis using flanking markers. By a targeted approach, the candidate region was subjected to screening for polymorphisms, mainly focusing on intronic sequences. A total of 18 genes were partially sequenced, and polymorphisms were identified in GP5, CENTB2, APOD, PCYT1A, OSTalpha, ZDHHC19, TFRC, ACK1, MUC4, MUC20, KIAA0226, LRCH3 and MUC13. Overall, 227 polymorphisms were discovered in the founder generation. The analysis revealed a large haplotype block, spanning at least 1.5 Mb around MUC4, to be associated with F4ab/ac susceptibility.

Refined localization of the Escherichia coli F4ab/F4ac receptor locus on pig chromosome 13.

Diarrhoea in newborn and weaned pigs caused by enterotoxigenic Escherichia coli (ETEC) expressing F4 fimbriae leads to considerable losses in pig production. In this study, we refined the mapping of the receptor locus for ETEC F4ab/F4ac adhesion (F4bcR) by joint analysis of Nordic and Swiss data. A total of 236 pigs from a Nordic experimental herd, 331 pigs from a Swiss experimental herd and 143 pigs from the Swiss performing station were used for linkage analysis. Genotyping data of six known microsatellite markers, two newly developed markers (MUC4gt and HSA125gt) and an intronic SNP in MUC4 (MUC4-8227) were used to create the linkage map. The region for F4bcR was refined to the interval SW207-S0075 on pig chromosome 13. The most probable position of F4bcR was in the SW207-MUC4 region. The order of six markers was supported by physical mapping on the BAC fingerprint contig from the Wellcome Trust Sanger Institute. Thus, the region for F4bcR could be reduced from 26 to 14 Mb.

Sequence assembly.

Despite the rapidly increasing number of sequenced and re-sequenced genomes, many issues regarding the computational assembly of large-scale sequencing data have remain unresolved. Computational assembly is crucial in large genome projects as well for the evolving high-throughput technologies and plays an important role in processing the information generated by these methods. Here, we provide a comprehensive overview of the current publicly available sequence assembly programs. We describe the basic principles of computational assembly along with the main concerns, such as repetitive sequences in genomic DNA, highly expressed genes and alternative transcripts in EST sequences. We summarize existing comparisons of different assemblers and provide a detailed descriptions and directions for download of assembly programs at: http://genome.ku.dk/resources/assembly/methods.html.

Porcine lung surfactant protein B gene (SFTPB): cDNA sequencing, chromosomal location, and expression studies.

The porcine surfactant protein B (SFTPB) is a single copy gene on chromosome 3. Three different cDNAs for the SFTPB have been isolated and sequenced. Nucleotide sequence comparison revealed six nonsynonymous single nucleotide polymorphisms (SNPs), four synonymous SNPs and an in-frame deletion of 69 bp in the region coding for the active protein. Northern analysis showed lung-specific expression of three different isoforms of the SFTPB transcript. The expression level for the SFTPB gene is low in 50 days-old fetus and it increases during lung development. Quantitative real-time polymerase chain reaction (qPCR) showed significant down regulation of the SFTPB in pigs with Actinobacillus pleuropneumoniae lung infection.

SNP-finding in pig mitochondrial ESTs.

The Sino-Danish pig genome project produced 685 851 ESTs (Gorodkin et al. 2007), of which 41 499 originated from the mitochondrial genome. In this study, the mitochondrial ESTs were assembled, and 374 putative SNPs were found. Chromatograms for the ESTs containing SNPs were manually inspected, and 112 total (52 non-synonymous) SNPs were found to be of high confidence (five of them are close to disease-causing SNPs in humans). Nine of the high-confidence SNPs were tested experimentally, and eight were confirmed. The SNPs can be accessed online at http://pigest.ku.dk/more/mito.

Characterization of the porcine carboxypeptidase E cDNA: comparative analysis and investigation of a non-synonymous SNP.

Carboxypeptidase E (CPE) is an important enzyme responsible for the proteolytic processing of prohormone intermediates. A naturally occurring point mutation, leading to an accumulation of many neuroendocrine peptides has been characterized within exon 5 of the CPE gene in mice. In the present study the sequence of the cDNA for the porcine CPE gene including all the coding region and the 3'-UTR region was generated. Comparisons with bovine, human, mouse, and rat CPE cDNA sequences showed that the coding regions of the gene are highly conserved both at the nucleotide and at the amino acid level. A very low nonsynonymous/synonymous substitution ratios between the proteins was found indicating that purifying selection is acting on the CPE gene. A nonsynonymous SNP identified at position 1272 in the transcript resulting in a codon change from TCA (Serine) to TTA (Leucine) was genotyped in the Danish pig populations. However, the mutated allele was only identified in one particular family. Exon 5 was sequenced in 45 unrelated pigs. No evidence of variation was found in this region. In accordance with previous results and in accordance with comparative mapping information the gene was accurately mapped to porcine chromosome 8 using the ImpRH panel.

Molecular characterization of the porcine surfactant, pulmonary-associated protein C gene.

The surfactant, pulmonary-associated protein C (SFTPC) is a peptide secreted by the alveolar type II pneumocytes of the lung. We have characterized the porcine SFTPC gene at genomic, transcriptional, and protein levels. The porcine SFTPC is a single-copy gene on pig chromosome 14. Two transcripts were found in a newborn pig lung cDNA library: a full-length clone and a clone missing exon 5. cDNA sequence comparison revealed four synonymous and two nonsynonymous substitutions and in-frame insertions at the beginning of exon 5. Comparison of the SFTPC coding region between several mammals showed high levels of conservation. Northern blot studies showed lung-specific expression of the full-length SFTPC transcript, appearing in 50-day-old fetus and increasing during lung development. Both SFTPC transcripts were detected mainly in lung by real-time RT-PCR and they were significantly down-regulated in necrotic lungs of pigs infected with Actinobacillus pleuropneumoniae. Additionally, the protein levels were also decreased or absent in the necrotic tissue.

Linkage mapping of gene-associated SNPs to pig chromosome 11.

Single nucleotide polymorphisms (SNPs) were discovered in porcine expressed sequence tags (ESTs) orthologous to genes from human chromosome 13 (HSA13) and predicted to be located on pig chromosome 11 (SSC11). The SNPs were identified as sequence variants in clusters of EST sequences from pig cDNA libraries constructed in the Sino-Danish pig genome project. In total, 312 human gene sequences from HSA13 were used for similarity searches in our pig EST database. Pig ESTs showing significant similarity with HSA13 genes were clustered and candidate SNPs were identified. Allele frequencies for 26 SNPs were estimated in a group of 80 unrelated pigs from Danish commercial pig breeds: Duroc, Hampshire, Landrace and Large White. Eighteen of the 26 SNPs genotyped in the PiGMaP Reference Families were mapped by linkage analysis to SSC11. The EST-based SNPs published here are new genetic markers useful for linkage and association studies in commercial and experimental pig populations. This study represents the first gene-associated SNP linkage map of pig chromosome 11 and adds new comparative mapping information between SSC11 and HSA13. Furthermore, our data facilitate future studies aimed at the identification of interesting regions on pig chromosome 11, positional cloning and fine mapping of quantitative trait loci in pig.