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.

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.

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.

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.

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.

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.

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.

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.

Evolutionary history of the sex-peptide (Acp70A) gene region in Drosophila melanogaster.

In Drosophila the products of the seminal fluid stimulate oviposition and suppress remating in the female. Of all the accessory gland peptides (Acp's) involved in these two responses, the sex-peptide (coded by the Acp70A gene) is among the best characterized at the functional level. A 1.2-kb fragment encompassing the Acp70A gene of nine lines from a natural population of D. melanogaster and one allele of D. sechellia was sequenced to study the forces shaping nucleotide variation within and between species. The coding region of D. simulans and D. mauritiana was also sequenced. A Ser to Ala replacement polymorphism at the last position of the signal peptide was detected in D. melanogaster. The Ser and Ala alleles are at intermediate frequencies. The level of nucleotide variation is lower for the derived Ala allele, which is compatible with a recent origin and an increase in frequency due to positive selection. Variation at the 5' flanking region is structured in two major highly differentiated haplotypes, whose distribution does not conform to neutral expectations. Selective and/or historical factors could contribute to the observed overall patterning of nucleotide variation at the Acp70A region.

Molecular characterization of the breakpoints of an inversion fixed between Drosophila melanogaster and D. subobscura.

The two breakpoints of a chromosomal inversion fixed since the split of Drosophila melanogaster and D. subobscura lineages have been isolated and sequenced in both species. The regions spanning the breakpoints initially were identified by the presence of two signals after interspecific in situ hybridization on polytene chromosomes. Interspecific comparison of the sequenced regions allowed us to delineate the location of the breakpoints. Close to one of these breakpoints a new transcription unit (bcn92) has been identified in both species. The inversion fixed between D. melanogaster and D. subobscura does not seem to have broken any transcription unit. Neither complete nor defective transposable elements were found in the regions encompassing the breakpoints. Short thymine-rich sequences (30-50 bp long) have been found bordering the breakpoint regions. Although alternating Pur-Pyr sequences were detected, these putative target sites for topoisomerase II were not differentially clustered in the breakpoints.

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.

The sex-peptide gene (Acp70A) is duplicated in Drosophila subobscura.

A 3.1-kb region of Drosphila subobscura homologous to the Acp70A region of D. melanogaster, which contains the sex-peptide gene, was cloned and sequenced. This region contains an approximately 600-bp duplication that includes the sex-peptide and its 5' and 3' flanking regions. The preproteins are 54 and 56 amino acids long, respectively (as compared to 55 amino acids in D. melanogaster), and each includes a 19-amino-acid-long signal peptide. The C-terminal part of the mature peptide is highly conserved between D. melanogaster and the two copies of D. subobscura. In this species, both copies of the gene are transcribed and, like in D. melangaster, only expressed in males. The duplicated region includes 300 bp upstream of the gene that would therefore seem sufficient for their expression in males. This region presents at its 5' end a stretch 93-bp that has a high similarity with the corresponding region of D. melanogaster and could be part of a still unidentified regulatory element of these genes.

Linkage and comparative mapping of the locus controlling susceptibility towards E. COLI F4ab/ac diarrhoea in pigs.

In 1995, Edfors-Lilja and coworkers mapped the locus for the E. COLI K88ab (F4ab) and K88ac (F4ac) intestinal receptor to pig chromosome 13 (SSC13). Using the same family material we have refined the map position to a region between the microsatellite markers Sw207 and Sw225. Primers from these markers were used to screen a pig BAC library and the positive clones were used for fluorescent in situ hybridization (FISH) analysis. The results of the FISH analysis helped to propose a candidate gene region in the SSC13q41-->q44 interval. Shotgun sequencing of the FISH-mapped BAC clones revealed that the candidate region contains an evolutionary breakpoint between human and pig. In order to further characterise the rearrangements between SSC13 and human chromosome 3 (HSA3), detailed gene mapping of SSC13 was carried out. Based on this mapping data we have constructed a detailed comparative map between SSC13 and HSA3. Two candidate regions on human chromosome 3 have been identified that are likely to harbour the human homologue of the gene responsible for susceptibility towards E. COLI F4ab/ac diarrhoea in pigs.