Towards identification of new genetic determinants for post-weaning diarrhea in piglets

Post-weaning diarrhea in pigs is a considerable challenge in the pig farming industry due to its effect on animal welfare and production costs, as well as the large volume of antibiotics, which are used to treat diarrhea in pigs after weaning. Previous studies have revealed loci on SSC6 and SSC13 associated with susceptibility to specific diarrhea causing pathogens. This study aimed to identify new genetic loci for resistance to diarrhea based on phenotypic data. In depth clinical characterization of diarrhea was performed in 257 pigs belonging to two herds during the first 14 days post weaning. The daily diarrhea assessments were used for the classification of pigs into case and control groups. Pigs were assigned to case and control groups based only on the incidence of diarrhea in the second week of the study in order to differentiate between differences in etiology. Genome-wide association studies and metabolomics association analysis were performed in order to identify new biological determinants for diarrhea susceptibility. With the present work, we revealed a new locus for diarrhea resistance on SSC16. Furthermore, studies of metabolomics in the same pigs revealed one metabolite associated with diarrhea.

Modeling microRNA-driven post-transcriptional regulation using exon-intron split analysis in pigs

The contribution of microRNAs (miRNAs) to mRNA post-transcriptional regulation has often been explored by the post hoc selection of downregulated genes and determining whether they harbor binding sites for miRNAs of interest. This approach, however, does not discriminate whether these mRNAs are also downregulated at the transcriptional level. Here, we have characterized the transcriptional and post-transcriptional changes in mRNA expression in two porcine tissues: gluteus medius muscle of fasted and fed Duroc gilts and adipose tissue of lean and obese Duroc-Göttingen minipigs. Exon-intron split analysis of RNA-seq data allowed us to identify downregulated mRNAs with high post-transcriptional signals in fed or obese states, and we assessed whether they harbor binding sites for upregulated miRNAs in any of these two physiological states. We found 26 downregulated mRNAs with high post-transcriptional signals in the muscle of fed gilts and 21 of these were predicted targets of miRNAs upregulated in fed pigs. For adipose tissue, 44 downregulated mRNAs in obese minipigs displayed high post-transcriptional signals, and 25 of these were predicted targets of miRNAs upregulated in the obese state. These results suggest that the contribution of miRNAs to mRNA repression is more prominent in the skeletal muscle system. Finally, we identified several genes that may play relevant roles in the energy homeostasis of the pig skeletal muscle (DKK2 and PDK4) and adipose (SESN3 and ESRRG) tissues. By differentiating transcriptional from post-transcriptional changes in mRNA expression, exon-intron split analysis provides a valuable view of the regulation of gene expression, complementary to canonical differential expression analyses.

Epigenetic and Transcriptomic Characterization of Pure Adipocyte Fractions From Obese Pigs Identifies Candidate Pathways Controlling Metabolism

Reprogramming of adipocyte function in obesity is implicated in metabolic disorders like type 2 diabetes. Here, we used the pig, an animal model sharing many physiological and pathophysiological similarities with humans, to perform in-depth epigenomic and transcriptomic characterization of pure adipocyte fractions. Using a combined DNA methylation capture sequencing and Reduced Representation bisulfite sequencing (RRBS) strategy in 11 lean and 12 obese pigs, we identified in 3529 differentially methylated regions (DMRs) located at close proximity to-, or within genes in the adipocytes. By sequencing of the transcriptome from the same fraction of isolated adipocytes, we identified 276 differentially expressed transcripts with at least one or more DMR. These transcripts were over-represented in gene pathways related to MAPK, metabolic and insulin signaling. Using a candidate gene approach, we further characterized 13 genes potentially regulated by DNA methylation and identified putative transcription factor binding sites that could be affected by the differential methylation in obesity. Our data constitute a valuable resource for further investigations aiming to delineate the epigenetic etiology of metabolic disorders.

Haplotypes on pig chromosome 3 distinguish metabolically healthy from unhealthy obese individuals

We have established a pig resource population specifically designed to elucidate the genetics involved in development of obesity and obesity related co-morbidities by crossing the obesity prone Göttingen Minipig breed with two lean production pig breeds. In this study we have performed genome wide association (GWA) to identify loci with effect on blood lipid levels. The most significantly associated single nucleotide polymorphisms (SNPs) were used for linkage disequilibrium (LD) and haplotype analyses. Three separate haploblocks which influence the ratio between high density lipoprotein cholesterol and total cholesterol (HDL-C/CT), triglycerides (TG) and low density lipoprotein cholesterol (LDL-C) levels respectively were identified on Sus Scrofa chromosome 3 (SSC3). Large additive genetic effects were found for the HDL-C/CT and LDL-C haplotypes. Haplotypes segregating from Göttingen Minipigs were shown to impose a positive effect on blood lipid levels. Thus, the genetic profile of the Göttingen Minipig breed seems to support a phenotype comparable to the metabolic healthy obese (MHO) phenotype in humans.

Transcriptional immune response in mesenteric lymph nodes in pigs with different levels of resistance to Ascaris suum

A single nucleotide polymorphism on chromosome 4 (SNP TXNIP) has been reported to be associated with roundworm (Ascaris suum) burden in pigs. The objective of the present study was to analyse the immune response to A. suum mounted by pigs with genotype AA (n = 24) and AB (n = 23) at the TXNIP locus. The pigs were repeatedly infected with A. suum from eight weeks of age until necropsy eight weeks later. An uninfected control group (AA; n = 5 and AB; n = 5) was also included. At post mortem, we collected mesenteric lymph nodes and measured the expression of 28 selected immune-related genes. Recordings of worm burdens confirmed our previous results that pigs of the AA genotype were more resistant to infection than AB pigs. We estimated the genotype difference in relative expression levels in infected and uninfected animals. No significant change in expression levels between the two genotypes due to infection was observed for any of the genes, although IL-13 approached significance (P = 0.08; Punadjusted = 0.003). Furthermore, statistical analysis testing for the effect of infection separately in each genotype showed significant up-regulation of IL-13 (P<0.05) and CCL17 (P<0.05) following A. suum infection in the 'resistant' AA genotype and not in the 'susceptible' AB genotype. Pigs of genotype AB had higher expression of the high-affinity IgG receptor (FCGR1A) than AA pigs in both infected and non-infected animals (P = 1.85*10-11).

Altered Methylation Profile of Lymphocytes Is Concordant with Perturbation of Lipids Metabolism and Inflammatory Response in Obesity

Obesity is associated with immunological perturbations that contribute to insulin resistance. Epigenetic mechanisms can control immune functions and have been linked to metabolic complications, although their contribution to insulin resistance still remains unclear. In this study, we investigated the link between metabolic dysfunction and immune alterations with the epigenetic signature in leukocytes in a porcine model of obesity. Global DNA methylation of circulating leukocytes, adipose tissue leukocyte trafficking, and macrophage polarisation were established by flow cytometry. Adipose tissue inflammation and metabolic function were further characterised by quantification of metabolites and expression levels of genes associated with obesity and inflammation. Here we show that obese pigs showed bigger visceral fat pads, higher levels of circulating LDL cholesterol, and impaired glucose tolerance. These changes coincided with impaired metabolism, sustained macrophages infiltration, and increased inflammation in the adipose tissue. Those immune alterations were linked to global DNA hypermethylation in both B-cells and T-cells. Our results provide novel insight into the possible contribution of immune cell epigenetics into the immunological disturbances observed in obesity. The dramatic changes in the transcriptomic and epigenetic signature of circulating lymphocytes reinforce the concept that epigenetic processes participate in the increased immune cell activation and impaired metabolic functions in obesity.

Genome-wide association study reveals a locus for nasal carriage of Staphylococcus aureus in Danish crossbred pigs

Background

Staphylococcus aureus is an important human opportunistic pathogen residing on skin and mucosae of healthy people. Pigs have been identified as a source of human colonization and infection with methicillin-resistant Staphylococcus aureus (MRSA) and novel measures are needed to control zoonotic transmission. A recent longitudinal study indicated that a minority of pigs characterized by high nasal load and stable carriage may be responsible for the maintenance of S. aureus within farms. The primary objective of the present study was to detect genetic loci associated with nasal carriage of S. aureus in Danish crossbred pigs (Danish Landrace/Yorkshire/Duroc).

Results

Fifty-six persistent carriers and 65 non-carriers selected from 15 farms surveyed in the previous longitudinal study were genotyped using Illumina’s Porcine SNP60 beadchip. In addition, spa typing was performed on 126 S. aureus isolates from 37 pigs to investigate possible relationships between host and S. aureus genotypes. A single SNP (MARC0099960) on chromosome 12 was found to be associated with nasal carriage of S. aureus at a genome-wide level after permutation testing (p = 0.0497) whereas the association of a neighboring SNP was found to be borderline (p = 0.114). Typing of S. aureus isolates led to detection of 11 spa types belonging to the three main S. aureus clonal complexes (CC) previously described in pigs (CC9, CC30 and CC398). Individual carriers often harbored multiple S. aureus genotypes and the host-pathogen interaction seems to be independent of S. aureus genotype.

Conclusion

Our results suggest it may be possible to select pigs genetically resistant to S. aureus nasal colonization as a tool to control transmission of livestock-associated MRSA to humans.

Electronic supplementary material

The online version of this article (doi:10.1186/s12917-015-0599-y) contains supplementary material, which is available to authorized users.

Comparative Analyses of QTLs Influencing Obesity and Metabolic Phenotypes in Pigs and Humans

The pig is a well-known animal model used to investigate genetic and mechanistic aspects of human disease biology. They are particularly useful in the context of obesity and metabolic diseases because other widely used models (e.g. mice) do not completely recapitulate key pathophysiological features associated with these diseases in humans. Therefore, we established a F2 pig resource population (n = 564) designed to elucidate the genetics underlying obesity and metabolic phenotypes. Segregation of obesity traits was ensured by using breeds highly divergent with respect to obesity traits in the parental generation. Several obesity and metabolic phenotypes were recorded (n = 35) from birth to slaughter (242 ± 48 days), including body composition determined at about two months of age (63 ± 10 days) via dual-energy x-ray absorptiometry (DXA) scanning. All pigs were genotyped using Illumina Porcine 60k SNP Beadchip and a combined linkage disequilibrium-linkage analysis was used to identify genome-wide significant associations for collected phenotypes. We identified 229 QTLs which associated with adiposity- and metabolic phenotypes at genome-wide significant levels. Subsequently comparative analyses were performed to identify the extent of overlap between previously identified QTLs in both humans and pigs. The combined analysis of a large number of obesity phenotypes has provided insight in the genetic architecture of the molecular mechanisms underlying these traits indicating that QTLs underlying similar phenotypes are clustered in the genome. Our analyses have further confirmed that genetic heterogeneity is an inherent characteristic of obesity traits most likely caused by segregation or fixation of different variants of the individual components belonging to cellular pathways in different populations. Several important genes previously associated to obesity in human studies, along with novel genes were identified. Altogether, this study provides novel insight that may further the current understanding of the molecular mechanisms underlying human obesity.

Gender and Obesity Specific MicroRNA Expression in Adipose Tissue from Lean and Obese Pigs

Obesity is a complex condition that increases the risk of life threatening diseases such as cardiovascular disease and diabetes. Studying the gene regulation of obesity is important for understanding the molecular mechanisms behind the obesity derived diseases and may lead to better intervention and treatment plans. MicroRNAs (miRNAs) are short non-coding RNAs regulating target mRNA by binding to their 3'UTR. They are involved in numerous biological processes and diseases, including obesity. In this study we use a mixed breed pig model designed for obesity studies to investigate differentially expressed miRNAs in subcutaneous adipose tissue by RNA sequencing (RNAseq). Both male and female pigs are included to explore gender differences. The RNAseq study shows that the most highly expressed miRNAs are in accordance with comparable studies in pigs and humans. A total of six miRNAs are differentially expressed in subcutaneous adipose tissue between the lean and obese group of pigs, and in addition gender specific significant differential expression is observed for a number of miRNAs. The differentially expressed miRNAs have been verified using qPCR. The results of these studies in general confirm the trends found by RNAseq. Mir-9 and mir-124a are significantly differentially expressed with large fold changes in subcutaneous adipose tissue between lean and obese pigs. Mir-9 is more highly expressed in the obese pigs with a fold change of 10 and a p-value < 0.001. Mir-124a is more highly expressed in the obese pigs with a fold change of 114 and a p-value < 0.001. In addition, mir-124a is significantly higher expressed in abdominal adipose tissue in male pigs with a fold change of 119 and a p-value < 0.05. Both miRNAs are also significantly higher expressed in the liver of obese male pigs where mir-124a has a fold change of 12 and mir-9 has a fold change of 1.6, both with p-values < 0.05.

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.

An f2 pig resource population as a model for genetic studies of obesity and obesity-related diseases in humans: design and genetic parameters

Obesity is a rising worldwide public health problem. Difficulties to precisely measure various obesity traits and the genetic heterogeneity in human have been major impediments to completely disentangle genetic factors causing obesity. The pig is a relevant model for studying human obesity and obesity-related (OOR) traits. Using founder breeds divergent with respect to obesity traits we have created an F2 pig resource population (454 pigs), which has been intensively phenotyped for 36 OOR traits. The main rationale for our study is to characterize the genetic architecture of OOR traits in the F2 pig design, by estimating heritabilities, genetic, and phenotypic correlations using mixed- and multi-trait BLUP animal models. Our analyses revealed high coefficients of variation (15–42%) and moderate to high heritabilities (0.22–0.81) in fatness traits, showing large phenotypic and genetic variation in the F2 population, respectively. This fulfills the purpose of creating a resource population divergent for OOR traits. Strong genetic correlations were found between weight and lean mass at dual-energy x-ray absorptiometry scanning (0.56–0.97). Weight and conformation also showed strong genetic correlations with slaughter traits (e.g., r_g between abdominal circumference and leaf fat at slaughtering: 0.66). Genetic correlations between fat-related traits and the glucose level vary between 0.35 and 0.74 and show a strong correlation between adipose tissue and impaired glucose metabolism. Our power calculations showed a minimum of 80% power for QTL detection for all phenotypes. We revealed genetic correlations at population level, for the first time, for several difficult to measure and novel OOR traits and diseases. The results underpin the potential of the established F2 pig resource population for further genomic, systems genetics, and functional investigations to unravel the genetic background of OOR traits.

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.

A molecular study of congenital erythropoietic porphyria in cattle

Previous studies have shown that congenital erythropoietic porphyria (CEP) in cattle is caused by an inherited deficiency of the enzyme uroporphyrinogen III synthase (UROS) encoded by the UROS gene. In this study, we have established the pedigree of an extended Holstein family in which the disease is segregating in a manner consistent with autosomal recessive inheritance. Biochemical analyses demonstrated accumulation of uroporphyrin, thus confirming that it is indeed insufficient activity of UROS which is the cause of the disease. We have therefore sequenced all nine exons of UROS in affected and non-affected individuals without detecting any potential causative mutations. However, a single nucleotide polymorphism (SNP) located within the spliceosome attachment region in intron 8 of UROS is shown to segregate with the disease allele. Our study supports the hypothesis that CEP in cattle is caused by a mutation affecting UROS; however, additional functional studies are needed to identify the causative mutation.

Characterisation of five candidate genes within the ETEC F4ab/ac candidate region in pigs

Background

Enterotoxigenic Escherichia coli (ETEC) that express the F4ab and F4ac fimbriae is a major contributor to diarrhoea outbreaks in the pig breeding industry, infecting both newborn and weaned piglets. Some pigs are resistant to this infection, and susceptibility is inherited as a simple dominant Mendelian trait. Indentifying the genetics behind this trait will greatly benefit pig welfare as well as the pig breeding industry by providing an opportunity to select against genetically susceptible animals, thereby reducing the number of diarrhoea outbreaks. The trait has recently been mapped by haplotype sharing to a 2.5 Mb region on pig chromosome 13, a region containing 18 annotated genes.

Findings

The coding regions of five candidate genes for susceptibility to ETEC F4ab/ac infection (TFRC, ACK1, MUC20, MUC4 and KIAA0226), all located in the 2.5 Mb region, were investigated for the presence of possible causative mutations. A total of 34 polymorphisms were identified in either coding regions or their flanking introns. The genotyping data for two of those were found to perfectly match the genotypes at the ETEC F4ab/ac locus, a G to C polymorphism in intron 11 of TFRC and a C to T silent polymorphism in exon 22 of KIAA0226. Transcriptional profiles of the five genes were investigated in a porcine tissue panel including various intestinal tissues. All five genes were expressed in intestinal tissues at different levels but none of the genes were found differentially expressed between ETEC F4ab/ac resistant and ETEC F4ab/ac susceptible animals in any of the tested tissues.

Conclusions

None of the identified polymorphisms are obvious causative mutations for ETEC F4ab/ac susceptibility, as they have no impact on the level of the overall mRNA expression nor predicted to influence the composition of the amino acids composition. However, we cannot exclude that the five tested genes are bona fide candidate genes for susceptibility to ETEC F4ab/ac infection since the identified polymorphism might affect the translational apparatus, alternative splice forms may exist and post translational mechanisms might contribute to disease susceptibility.

A study of alternative splicing in the pig

Background

Since at least half of the genes in mammalian genomes are subjected to alternative splicing, alternative pre-mRNA splicing plays an important contribution to the complexity of the mammalian proteome. Expressed sequence tags (ESTs) provide evidence of a great number of possible alternative isoforms. With the EST resource for the domestic pig now containing more than one million porcine ESTs, it is possible to identify alternative splice forms of the individual transcripts in this species from the EST data with some confidence.

Results

The pig EST data generated by the Sino-Danish Pig Genome project has been assembled with publicly available ESTs and made available in the PigEST database. Using the Distiller package 2,515 EST clusters with candidate alternative isoforms were identified in the EST data with high confidence. In agreement with general observations in human and mouse, we find putative splice variants in about 30% of the contigs with more than 50 ESTs. Based on the criteria that a minimum of two EST sequences confirmed each splice event, a list of 100 genes with the most distinct tissue-specific alternative splice events was generated from the list of candidates. To confirm the tissue specificity of the splice events, 10 genes with functional annotation were randomly selected from which 16 individual splice events were chosen for experimental verification by quantitative PCR (qPCR). Six genes were shown to have tissue specific alternatively spliced transcripts with expression patterns matching those of the EST data. The remaining four genes had tissue-restricted expression of alternative spliced transcripts. Five out of the 16 splice events that were experimentally verified were found to be putative pig specific.

Conclusions

In accordance with human and rodent studies we estimate that approximately 30% of the porcine genes undergo alternative splicing. We found a good correlation between EST predicted tissue-specificity and experimentally validated splice events in different porcine tissue. This study indicates that a cluster size of around 50 ESTs is optimal for in silico detection of alternative splicing. Although based on a limited number of splice events, the study supports the notion that alternative splicing could have an important impact on species differentiation since 31% of the splice events studied appears to be species specific.

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.

A study of the occurrence of monochorionic and monozygotic twinning in the pig

In humans as well as in most farm animals, monozygotic twins have been described. Nevertheless, only a few reports of twinning in the pig have been published. It has been suggested that monozygotic twins are formed during the first 14 days of pregnancy. This monozygotic twin study includes the investigation of porcine monochorionic embryos from 76 sows at days 26-29 post-insemination (p.i.), as well as an examination of 10 whole litters at days 21-22 p.i. In the former group, 29% of the sows carried monochorionic embryos. Based on DNA profiling using microsatellite markers, one monozygotic twin pair was found among these embryos. In the latter group, three monozygotic twin pairs were identified. Thus, it can be concluded that although the occurrence of monozygotic twins in pigs is a sporadic event, the fusion of extra-embryonic membranes is relatively common.

Phenotypic and genetic characterization of a novel phenotype in pigs characterized by juvenile hairlessness and age dependent emphysema

Background

A pig phenotype characterized by juvenile hairlessness, thin skin and age dependent lung emphysema has been discovered in a Danish pig herd. The trait shows autosomal co-dominant inheritance with all three genotypes distinguishable. Since the phenotype shows resemblance to the integrin β₆ -/- knockout phenotype seen in mice, the two genes encoding the two subunits of integrin α_vβ₆, i.e. ITGB6 and ITGAV, were considered candidate genes for this trait.

Results

The mutated pig phenotype is characterized by hairlessness until puberty, thin skin with few hair follicles and absence of musculi arrectores pili, and at puberty or later localized areas of emphysema are seen in the lungs. Comparative mapping predicted that the porcine ITGB6 andITGAV orthologs map to SSC15. In an experimental family (n = 113), showing segregation of the trait, the candidate region was confirmed by linkage analysis with four microsatellite markers. Mapping of the porcine ITGB6 and ITGAV in the IMpRH radiation hybrid panel confirmed the comparative mapping information. Sequencing of the ITGB6 and ITGAV coding sequences from affected and normal pigs revealed no evidence of a causative mutation, but alternative splicing of the ITGB6 pre-mRNA was detected. For both ITGB6 and ITGAV quantitative PCR revealed no significant difference in the expression levels in normal and affected animals. In a western blot, ITGB6 was detected in lung protein samples of all three genotypes. This result was supported by flow cytometric analyses which showed comparable reactions of kidney cells from affected and normal pigs with an integrin α_vβ₆ monoclonal antibody. Also, immunohistochemical staining of lung tissue with an integrin β₆ antibody showed immunoreaction in both normal and affected pigs.

Conclusion

A phenotype resembling the integrin β₆ -/- knockout phenotype seen in mice has been characterized in the pig. The candidate region on SSC15 has been confirmed by linkage analysis but molecular and functional analyses have excluded that the mutated phenotype is caused by structural mutations in or ablation of any of the two candidate genes.

Genome-wide identification of quantitative trait loci in a cross between Hampshire and Landrace II: meat quality traits

Background

Meat quality traits are important in pig breeding programs, but they are difficult to include in a traditional selection program. Marker assisted selection (MAS) of meat quality traits is therefore of interest in breeding programs and a Quantitative Trait Locus (QTL) analysis is the key to identifying markers that can be used in MAS. In this study, Landrace and Hampshire intercross and backcross families were used to investigate meat quality traits. Hampshire pigs are commonly used as the sire line in commercial pig breeding. This is the first time a pedigree including Hampshire pigs has been used for a QTL analysis of meat quality traits.

Results

In total, we analyzed 39 meat quality traits and identified eight genome-wide significant QTL peaks in four regions: one on chromosome 3, two on chromosome 6 and one on chromosome 16. At least two of the QTLs do not appear to have been detected in previous studies. On chromosome 6 we identified QTLs for water content in M. longissimus dorsi (LD), drip loss in LD and post mortem pH decline in LD. On chromosomes 3 and 16 we identified previously undetected QTLs for protein content in LD and for freezing and cooking loss respectively.

Conclusion

We identified at least two new meat quality trait QTLs at the genome-wide significance level. We detected two QTLs on chromosome 6 that possibly coincide with QTLs detected in other studies. We were also able to exclude the C1843T mutation in the ryanodine receptor (RYR1) as a causative mutation for one of the chromosome 6 QTLs in this cross.

Selection of reference genes for gene expression studies in pig tissues using SYBR green qPCR

BACKGROUND

Real-time quantitative PCR (qPCR) is a method for rapid and reliable quantification of mRNA transcription. Internal standards such as reference genes are used to normalise mRNA levels between different samples for an exact comparison of mRNA transcription level. Selection of high quality reference genes is of crucial importance for the interpretation of data generated by real-time qPCR.

RESULTS

In this study nine commonly used reference genes were investigated in 17 different pig tissues using real-time qPCR with SYBR green. The genes included beta-actin (ACTB), beta-2-microglobulin (B2M), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), hydroxymethylbilane synthase (HMBS), hypoxanthine phosphoribosyltransferase 1 (HPRT1), ribosomal protein L4 (RPL4), succinate dehydrogenase complex subunit A (SDHA), TATA box binding protein (TPB)and tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein zeta polypeptide (YWHAZ). The stability of these reference genes in different pig tissues was investigated using the geNorm application. The range of expression stability in the genes analysed was (from the most stable to the least stable): ACTB/RPL4, TBP, HPRT, HMBS, YWHAZ, SDHA, B2M and GAPDH.

CONCLUSION

Expression stability varies greatly between genes. ACTB, RPL4, TPB and HPRT1 were found to have the highest stability across tissues. Based on both expression stability and expression level, our data suggest that ACTB and RPL4 are good reference genes for high abundant transcripts while TPB and HPRT1 are good reference genes for low abundant transcripts in expression studies across different pig tissues.

SNP mining porcine ESTs with MAVIANT, a novel tool for SNP evaluation and annotation

MOTIVATION

Single nucleotide polymorphisms (SNPs) analysis is an important means to study genetic variation. A fast and cost-efficient approach to identify large numbers of novel candidates is the SNP mining of large scale sequencing projects. The increasing availability of sequence trace data in public repositories makes it feasible to evaluate SNP predictions on the DNA chromatogram level. MAVIANT, a platform-independent Multipurpose Alignment VIewing and Annotation Tool, provides DNA chromatogram and alignment views and facilitates evaluation of predictions. In addition, it supports direct manual annotation, which is immediately accessible and can be easily shared with external collaborators.

RESULTS

Large-scale SNP mining of polymorphisms bases on porcine EST sequences yielded more than 7900 candidate SNPs in coding regions (cSNPs), which were annotated relative to the human genome. Non-synonymous SNPs were analyzed for their potential effect on the protein structure/function using the PolyPhen and SIFT prediction programs. Predicted SNPs and annotations are stored in a web-based database. Using MAVIANT SNPs can visually be verified based on the DNA sequencing traces. A subset of candidate SNPs was selected for experimental validation by resequencing and genotyping. This study provides a web-based DNA chromatogram and contig browser that facilitates the evaluation and selection of candidate SNPs, which can be applied as genetic markers for genome wide genetic studies.

AVAILABILITY

The stand-alone version of MAVIANT program for local use is freely available under GPL license terms at http://snp.agrsci.dk/maviant.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

Evaluation of the porcine melanocortin 4 receptor (MC4R) gene as a positional candidate for a fatness QTL in a cross between Landrace and Hampshire

Melanocortin 4 receptor (MC4R) is expressed in the appetite-regulating areas of the brain where it is central in the regulation of feed intake and energy balance. A mutation in MC4R causing an Asp298Asn substitution has been associated with fatness, high daily gain and feed intake in the pig. In a previously performed genome scan based on a Hampshire x Landrace cross, we detected one quantitative trait loci (QTL) affecting carcass fat/meat ratio and one QTL affecting the biceps femoris muscle, both close to the position of MC4R on porcine chromosome 1. In this study, the two lines were found to be close to fixation for alternative alleles of the Asp298Asn polymorphism. Additional QTL analyses supported our hypothesis of MC4R as a positional candidate gene but only for the fat/meat QTL. The Asp298Asn polymorphism was also evaluated as a selection target for daily gain in a Danish pig breeding population that included four breeds (Hampshire, Duroc, Landrace and Yorkshire). Over a 12-year period (1990-2002), a significant increase in the allele frequency of 298Asn was found in Landrace and Duroc, whereas a non-significant decrease in the 298Asn allele frequency was observed in Yorkshire. The Hampshire breed was fixed for the 298Asn allele in 1990. The high 298Asn allele frequencies in Hampshire, Landrace and Duroc are most likely due to selection for daily gain, whereas selection for daily gain in the Yorkshire breed apparently focuses on other loci.

Genome-wide identification of quantitative trait loci in a cross between Hampshire and Landrace I: carcass traits

We report the identification of quantitative trait loci (QTL) affecting carcass composition, carcass length, fat deposition and lean meat content using a genome scan across 462 animals from a combined intercross and backcross between Hampshire and Landrace pigs. Data were analysed using multiple linear regression fitting additive and dominance effects. This model was compared with a model including a parent-of-origin effect to spot evidence of imprinting. Several precisely defined muscle phenotypes were measured in order to dissect body composition in more detail. Three significant QTL were detected in the study at the 1% genome-wide level, and twelve significant QTL were detected at the 5% genome-wide level. These QTL comprise loci affecting fat deposition and lean meat content on SSC1, 4, 9, 10, 13 and 16, a locus on SSC2 affecting the ratio between weight of meat and bone in back and weight of meat and bone in ham and two loci affecting carcass length on SSC12 and 17. The well-defined phenotypes in this study enabled us to detect QTL for sizes of individual muscles and to obtain information of relevance for the description of the complexity underlying other carcass traits.

Experimental infection with Escherichia coli O149:F4ac in weaned piglets

The outcome of experimental intestinal infections with enterotoxigenic Escherichia coli (ETEC) is dependent on several factors. An important factor is adhesion of the challenge strain to the intestinal mucosa. The test for susceptibility towards ETEC adhesion has so far been made by an intestinal adhesion test made after slaughter of piglets. However, in an experimental infection study with the purpose to obtain diarrhoeic piglets, it would be an advantage to test for susceptibility prior to experimentation. The Mucin 4 gene on porcine chromosome 13 has been proposed as a candidate gene for the production of the specific ETEC F4ab/ac receptor, and a DNA marker-based test has been developed to allow genotyping for ETEC F4ab/ac resistance/susceptibility [Jørgensen, C.B., Cirera, S., Archibald, A.L., Anderson, L., Fredholm, M., Edfors-Lilja, I., 2004. Porcine polymorphisms and methods for detecting them. International application published under the patent cooperation treaty (PCT). PCT/DK2003/000807 or WO2004/048606-A2]. The aim of this study was to test an experimental model for ETEC O149:F4ac-induced diarrhoea in piglets, selected for susceptibility towards ETEC O149:F4ac adhesion prior to experimentation using a DNA marker-based test. Sixty-two healthy 25-32 days old recently weaned Danish crossbred piglets were used. All piglets were tested prior to experimentation for susceptibility or resistance towards ETEC O149:F4ac adhesion. Thirty-nine piglets, both susceptible and resistant, were oro-gastric intubated with 10(9)CFU of ETEC O149:F4ac and 23 age-matched piglets, both susceptible and resistant, were used as non-infected controls. Of susceptible piglets, challenged with ETEC O149:F4ac, 74% had ETEC O149:F4ac-associated diarrhoea first day after first challenge, which were significantly higher relatively to the resistant and challenged piglets where 20% had diarrhoea (p=0.04). This study suggests a model for experimental ETEC induced diarrhoea.

A high-resolution comparative map between pig chromosome 17 and human chromosomes 4, 8, and 20: identification of synteny breakpoints

We report on the construction of a high-resolution comparative map of porcine chromosome 17 (SSC17) focusing on evolutionary breakpoints with human chromosomes. The comparative map shows high homology with human chromosome 20 but suggests more limited homologies with other human chromosomes. SSC17 is of particular interest in studies of chromosomal organization due to the presence of QTLs that affect meat quality and carcass composition. A total of 158 pig ESTs available in databases or developed by the Sino-Danish Pig Genome Sequencing Consortium were mapped using the INRA-University of Minnesota porcine radiation hybrid panel. The high-resolution map was further anchored by fluorescence in situ hybridization. This study confirmed the extensive conservation between SSC17 and HSA20 and enabled the gene order to be determined. The homology of the SSC17 pericentromeric region was extended to other human chromosomes (HSA4, HSA8) and the chromosomal breakpoint boundaries were accurately defined. In total 15 breakpoints were identified.

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.

Genetic diversity of eleven European pig breeds

A set of eleven pig breeds originating from six European countries, and including a small sample of wild pigs, was chosen for this study of genetic diversity. Diversity was evaluated on the basis of 18 microsatellite markers typed over a total of 483 DNA samples collected. Average breed heterozygosity varied from 0.35 to 0.60. Genotypic frequencies generally agreed with Hardy-Weinberg expectations, apart from the German Landrace and Schwäbisch-Hällisches breeds, which showed significantly reduced heterozygosity. Breed differentiation was significant as shown by the high among-breed fixation index (overall F_ST = 0.27), and confirmed by the clustering based on the genetic distances between individuals, which grouped essentially all individuals in 11 clusters corresponding to the 11 breeds. The genetic distances between breeds were first used to construct phylogenetic trees. The trees indicated that a genetic drift model might explain the divergence of the two German breeds, but no reliable phylogeny could be inferred among the remaining breeds. The same distances were also used to measure the global diversity of the set of breeds considered, and to evaluate the marginal loss of diversity attached to each breed. In that respect, the French Basque breed appeared to be the most "unique" in the set considered. This study, which remains to be extended to a larger set of European breeds, indicates that using genetic distances between breeds of farm animals in a classical taxonomic approach may not give clear resolution, but points to their usefulness in a prospective evaluation of diversity.

FISH mapping and sequence analysis of 87 porcine BAC clones

Ninety-one bacterial artificial chromosomes (BAC) clones, selected effectively at random from our library, were used as probes for fluorescence in situ hybridization. Of these, 87 clones gave a specific signal in one, two or three different pair(s) of swine metaphase chromosomes. The ends of 35 BAC clones were sequenced in order to obtain information for comparative mapping. Fifteen of them gave useful comparative mapping information.

Comparative mapping in the pig: localization of 214 expressed sequence tags

In total, 214 ESTs (Expressed Sequence Tags) were assigned to the porcine gene map by using somatic cell hybrid mapping, radiation hybrid mapping, and FISH. The ESTs were isolated from a porcine small intestine cDNA library on the basis of significant sequence identity with human annotated genes. In total, 390 primer pairs were designed primarily in the 3' UTR of the sequences. Overall, 58.6% of the ESTs were successfully mapped by this approach. In total, 191 of the localizations are in agreement with the human comparative map, strongly indicating that these represent true orthologous genes. The remaining 23 ESTs provide new comparative mapping data, which should be considered as preliminary until confirmed by other studies. Our mapping efforts provide a significant contribution to the porcine map as well as to the comparative map for human and pig.

The GENETPIG database: a tool for comparative mapping in pig (Sus scrofa)

The GENETPIG database has been established for storing and disseminating the results of the European project: 'GENETPIG: identification of genes controlling economic traits in pig'. The partners of this project have mapped about 630 porcine and human ESTs onto the pig genome. The database collects the mapping results and links them to other sources of mapping data; this includes pig maps as well as available comparative mapping information. Functional annotation of the mapped ESTs is also given when a significant similarity to cognate genes was established. The database is accessible for consultation via the Internet at http://www.infobiogen.fr/services/Genetpig/.

Amphibian respiration and olfaction and their relationships: from Robert Townson (1794) to the present

The present review examines the developments in the elucidation of the mechanisms of amphibian respiration and olfaction. Research in these two areas has largely proceeded along independent lines, despite the fact that ventilation of the nasobuccopharyngeal cavity is a basic element in both functions. The English naturalist Robert Townson demonstrated, in the 1790s, that amphibians, contrary to general belief, ventilated the lungs by a pressure-pump mechanism. Frogs and other amphibians respire by alternatively dilating and contracting the buccopharyngeal cavity. During dilatation, with the mouth and glottis closed, air is sucked in through the open nostrils to fill the cavity. During contraction of the throat, with nostrils closed and glottis open, the air in the buccopharyngeal cavity is pressed into the lungs. During expiration, the glottis and nostrils open and air is expelled from the lungs 'by their own contraction from a state of distention'. Herholdt (1801), a Danish army surgeon, independently described the buccal pressure-pump mechanism in frogs, his experiments being confirmed by the commissioners of the Société Philomatique in Paris. Haro (1842) reintroduced a suction mechanism for amphibian respiration, which Panizza (1845) refuted: excision of the tympanic membranes prevented lung inflation, the air in the buccopharyngeal cavity leaving through the tympanum holes. Closure of the holes with the fingers restored lung inflation. The importance of cutaneous respiration in frogs and other amphibians was discovered by Spallanzani (1803), who found that frogs might survive excision of the lungs and that the amounts of exhaled carbon dioxide were small compared with those eliminated through the skin. Edwards (1824) confirmed and extended Spallanzani's findings, and Regnault & Reiset (1849) attempted to establish the relative importance of skin and lungs as respiratory organs in frogs. The problem was solved by Krogh (1904a) who measured respiration through the skin and lungs separately and simultaneously. Krogh (1904a) confirmed that carbon dioxide was chiefly eliminated through the skin, correlated with its high diffusion rate in water and tissue, whereas the pattern of oxygen uptake varied seasonally, the pulmonary uptake being lower than the cutaneous during autumn and winter, but substantially higher during the breeding period. Dolk & Postma (1927) confirmed this respiratory pattern. More recently, Hutchison and coworkers have examined the relative role of pulmonary and cutaneous gas exchange in a large number of amphibians, equipped with head masks for the separate measurement of the lung respiration in normally ventilating animals (Vinegar & Hutchison, 1965; Guimond & Hutchison, 1968; Hutchison, Whitford & Kohl, 1968; Whitford & Hutchison, 1963, 1965, 1966). As early as 1758, Rösel von Rosenhof suggested that the lungs of frogs in water functioned as hydrostatic organs that permitted the animal to float at the surface or rest on the bottom of the pond. The suggestion was inspired by observations made in the second half of the seventeenth century by members of the Royal Academy of Sciences in Paris. The French anatomists demonstrated that a tortoise, presumably the European freshwater turtle Emys orbicularis, could regulate its buoyancy by changing the volume of the lungs, to descend passively or ascend in the water. The hydrostatic function of the lungs has been repeatedly rediscovered, by Emery (1869) in the frog, by Marcacci (1895) in frogs, toads and salamanders, by Whipple (1906b) in a newt, by Willem (1920, 1931) in frogs and Xenopus laevis, by Speer 1942) in several anurans and urodeles, and finally by de Jongh (1972) in Xenopus laevis. In the second half of the nineteenth century a number of important papers appeared which confirmed and extended Townson's (1794) and Panizza's (1845) analysis of the normal respiratory movements in frogs. (ABSTRACT TRUNCATED)

Role of urinary and cloacal bladders in chelonian water economy: historical and comparative perspectives

The Parisian comparative anatomist Claude Perrault, dissecting an Indian giant tortoise in 1676, was the first to observe that the urinary bladder is of an extraordinary size in terrestrial tortoises. In 1799, the English comparative physiologist Robert Townson suggested that the bladder functioned as a water reservoir, as he had shown previously for frogs and toads. However, these observations went unnoticed in subsequent reports on tortoise water economy that were made by travellers and naturalists visiting the Galapagos Archipelago and marvelling over the huge numbers of giant tortoises that inhabited these desert-like islands. The first such report was by an American naval officer, David Porter, who was a privateer in the 1812-15 war with England. In his journal he referred to the constant supply of water which the Galapagos tortoises carried with them. References to the location in the body, as well as the amounts and quality of the water stored, were, however, contradictory. The confusion concerning the anatomical identity of the water reservoir in the Galapagos tortoise, Geochelone elephantopus, persisted throughout the nineteenth century, and continued when studies of tortoise water economy and drinking behaviour in arid environments were taken up independently in the desert tortoise, Gopherus agassizii, which inhabits the desert regions in the south-western United States. In 1881 Cox found large sacs filled with clear water under the carapace, but it was half a century later that these sacs were identified as the large bilobed bladder; references to specific water sacs continued to appear in the literature until the 1960s. Since 1970, information on the water economy of desert tortoises has been obtained from extensive field studies. Rates of disappearance of tritiated water injected into the body have shown that during the drought periods of the summer, water turnover (intake) rates do not differ from the rates of metabolic water production. Under these conditions urine is not voided, but is stored in the large bladder. During a drought period the bladder urine increases from initially low osmolality finally to reach isosmolality with the blood plasma. Soluble K+ is the major cation of the urine, but large amounts of K+ are also present as precipitated urates. During a drought period the body is in negative water balance, but despite substantial losses of total body water, the plasma concentrations of Na+ and Cl- can remain constant for many months, indicating regulation of the extracellular fluid and water content of the body tissues by reabsorption of water from the urinary bladder. The bladder thus acts both as a store for nitrogenous waste and K+ and as a water reservoir during droughts. Following rain showers, there is a sharp decline in tritium activity correlated with copious drinking from temporary pools of rain water. The old bladder urine is voided and most of the water drunk is stored as a highly dilute urine. In 1676 Perrault observed that in a freshwater turtle, Emys orbicularis, but not in the giant tortoise, two other bladders opened into the cloaca. By the mid-twentieth century it had been established that these cloacal bladders typically were restricted to species of chelonians that led a semi-terrestrial or semi-aquatic life. The function of the bladders has been debated since Townson observed in 1799 that dehydrated freshwater turtles took up water by anal drinking, suggesting that anal drinking served in the water economy of semi-terrestrial turtles. Since then, the bladders have been ascribed hydrostatic and respiratory functions, but the recent literature mostly argues for a respiratory function. The possible role of the cloacal bladders as a water reservoir in amphibious turtles is still open. Terrestrial amphibians and tortoises are unique among vertebrates in possessing large urinary bladders that may function as water reservoirs in dry environments. (ABSTRACT TRUNCATED)

Characterisation of a bovine/murine hybrid cell panel informative for all bovine autosomes

A bovine/murine hybrid cell panel consisting of 57 cell lines was typed with 124 markers by PCR. Southern hybridisation and isozyme analysis in order to establish its utility as a resource for genome mapping. All bovine chromosomes, including the sex chromosomes were represented in the panel. Computerised analysis of syntenies indicated that there are no cell lines containing only a single bovine chromosome. The panel was used to map 10 new bovine microsatellite markers, and the MYL6 and CPE genes. This panel is informative for all bovine chromosomes other than the sex-specific region of the X chromosome and can be used in synteny mapping studies. At present, due to the relatively small number of markers typed, the resolution of the panel does not go beyond the chromosomal level.

Urea and amphibian water economy

Accumulation of urea in the body fluids enables some amphibians to tolerate high ambient salinities (Bufo viridis, Xenopus laevis, Rana cancrivora, Ambystoma tigrinum, Batrachoseps spp.) or to estivate in soil with low water potentials (Scaphiopus spp.). These species are assumed not only to accumulate urea produced in the normal metabolism, but to synthesize urea in response to water shortage. Re-examination of the data did not support the view of an osmoregulatory urea synthesis. Increased urea synthesis on exposure to high salinities in X. laevis, R. cancrivora and Batrachoseps spp. seemed to reflect reactions to an adverse environment. It is suggested that in amphibians, solute concentration in the plasma and rate of excretion of urea are coordinated so that at a certain plasma concentration, urea is excreted at the same rate at which it is produced. The higher the level of urea in the body fluids at balance between production and excretion, the higher the tolerance of the species of low external water potentials. The mechanisms that integrate the relationship between plasma solute concentration and handling of urea by the kidneys are not known.

200 years of amphibian water economy: from Robert Townson to the present

In the 1790s, Robert Townson established the main features of the water economy of terrestrial amphibians: rapid evaporative water loss in dry surroundings, 'drinking' by absorption of water through the abdominal skin pressed against moist substrates, and use of the urinary bladder as a reservoir from which water is reabsorbed on land. This knowledge was of little interest to the establishment in the first half of the nineteenth century of experimental physiology as a basic medical discipline, when frogs became models in the elucidation of general physiological processes. Townson's pioneer contributions to amphibian physiology were forgotten for 200 years (Jørgensen 1994b). During (1901) and particularly Overton (1904) restored knowledge about amphibian water economy to the level reached by Townson, but the papers had little impact on the young science of animal physiology because they primarily aimed at elucidating the transport of fluids across membranes. Frog skin remained a model preparation in such studies throughout the century. With the establishment of terrestrial ecology early in the century, the relations of animals, including amphibians, to water became a central theme. Concurrently with comparative studies of amphibian water economy in an ecological setting, the subject proceeded as an aspect of animal osmoregulation. Adolph (1920-1930) and Rey (1937a) established the highly dynamic nature of water balance in amphibians in water and on land. Their observations indicated functional links between environment, skin and kidneys, the nature of which remained to be explored. Thorson & Svihla (1943) reopened the ecological approach in a comparative study of the relations between amphibian habitat and tolerance of dehydration. By mid-century, the central themes of amphibian adaptations to terrestrial modes of life were re-established, except for the function of the bladder as a water-depot. During the following decades, a rich literature appeared, particularly focusing on adaptations of amphibians to arid environments. Thus, in the 1970s, it was found that 'waterproofing' of the highly permeable skins by means of skin secretions had evolved independently in several families of tropical arboreal frogs, and that a number of amphibians that aestivate whilst burrowed in dry soil could reduce evaporation by forming cocoons from shed strata cornea. In 1950-1970 the role of bladder urine as a water depot in terrestrial amphibians was recognized: this did not change the established view of water balance in terrestrial amphibians as alternating between dehydration on land and rehydration in response to the deficit in body water. Amphibians may, however, maintain normal water balance whether the ambient medium is water or air by means of little understood integrated mechanisms in control of cutaneous drinking behaviour, water permeability of the skin and bladder wall, and urine production.