A physical map of large segments of pig chromosome 7q11-q14: comparative analysis with human chromosome 6p21

Wild Boar (Sus scrofa)—Fascioloides magna Interaction from the Perspective of the MHC Genes

Fascioloidosis is a parasitic disease caused by a trematode Fascioloides magna. Since major histocompatibility complex (MHC) genes play an important role in the immune response, the aim of this study was to compare the potential differences in MHC class II SLA-DRB1 exon 2 genes between wild boar populations from infected (cases) and non-infected areas (controls). During the winter of 2021, a total of 136 wild boar tissue samples were collected, 39 cases and 97 controls. DNA was extracted and sequenced using the Illumina platform. Differences in distributions of allele combinations were calculated using the Chi-Square test for homogeneity and between proportions using the large-sample test and Fisher–Irwin test. Analysis revealed 19 previously described swine leucocyte antigen (SLA) alleles. The number of polymorphic sites was 79 (29.6%), with 99 mutations in total. Nucleotide diversity π was estimated at 0.11. Proportions of the alleles SLA-DRB1*12:05 (p = 0.0008379) and SLA-DRB1*0101 (p = 0.0002825) were statistically significantly higher in controls, and proportions of the SLA-DRB1*0602 (p = 0.006059) and SLA-DRB1*0901 (p = 0.0006601) in cases. Alleles SLA-DRB1*04:09, SLA-DRB1*0501, SLA-DRB1*11:09, and SLA-DRB1*1301 were detected only in cases, while SLA-DRB1*0404, SLA-DRB1*0701, SLA-DRB1*02:10, and SLA-DRB1*04:08 were present only in controls. We did not confirm the existence of specific alleles that could be linked to F. magna infection. Detected high variability of the MHC class II SLA-DRB1 exon 2 genes indicate high resistance potential against various pathogens.

Evolution of MHC class II SLA-DRB1 locus in the Croatian wild boar (Sus scrofa) implies duplication and weak signals of positive selection

The wild boar is an ancestor of the domestic pig and an important game species with the widest geographical range of all ungulates. Although a large amount of data are available on major histocompatibility complex (MHC) variability in domestic pigs, only a few studies have been performed on wild boars. Due to their crucial role in appropriate immune responses and extreme polymorphism, MHC genes represent some of the best candidates for studying the processes of adaptive evolution. Here, we present the results on the variability and evolution of the entire MHC class II SLA-DRB1 locus exon 2 in 133 wild boars from Croatia. Using direct sequencing and cloning methods, we identified 20 SLA-DRB1 alleles, including eight new variants, with notable divergence. In some individuals, we documented functional locus duplication, and SLA-DRB1*04:10 was identified as the allele involved in the duplication. The expression of a duplicated locus was confirmed by cloning and sequencing cDNA-derived amplicons. Based on individual genotypes, we were able to assume that alleles SLA-DRB1*04:10 and SLA-DRB1*06:07 are linked as an allelic combination that co-evolves as a two-locus haplotype. Our investigation of evolutionary processes at the SLA-DRB1 locus confirmed the role of intralocus recombination in generating allelic variability, whereas tests of positive selection based on the dN/dS (non-synonymous/synonymous substitution rate ratio) test revealed atypically weak and ambiguous signals.

Localization, expression change in PRRSV infection and association analysis of the porcine TAP1 gene

The transporter associated with antigen processing (TAP) translocates antigenic peptides from the cytosol into the lumen of the endoplasmic reticular and plays a critical role in the major histocompatibility complex (MHC) class I molecule-mediated antigenic presentation pathway. In this study, the porcine TAP1 gene was mapped to the pig chromosome 7 (SSC7) and was closely linked to the marker SSC2B02 (retention fraction=43%, LOD=15.18). Subcellular localization of TAP1 by transient transfection of PK15 cells indicated that the TAP1 protein might be located in the endoplasmic reticulum (ER) in pig kidney epithelial cells (PK-15). Gene expression analysis by semi-quantitative RT-PCR revealed that TAP1 was selectively expressed in some immune and immune-related tissues. Quantitative real-time PCR (qRT-PCR) analysis revealed that this gene was up-regulated after treatments that mimic viral and bacterial infection (polyriboinosinic-polyribocytidylic acid (poly(I:C)) and lipopolysaccharide (LPS), respectively). In addition, elevated TAP1 expression was detected after porcine reproductive and respiratory syndrome virus (PRRSV) infection in porcine white blood cells (WBCs). One single nucleotide polymorphism (SNP) in exon 3 of TAP1 was detected in a Landrace pig population by Bsp143I restriction enzyme digestion. Different genotypes of this SNP had significant associations (P<0.05) with the red blood cell distribution width (RDW) of 1-day-old (1 d) pigs (P=0.0168), the PRRSV antibody level (PRRSV Ab) (P=0.0445) and the absolute lymphocyte count (LYM#) (P=0.024) of 17 d pigs. Our results showed that the TAP1 gene might have important roles in swine immune responses, and these results provide useful information for further functional studies.

Mapping quantitative trait loci for metabolic and cytological fatness traits of connected F2 crosses in pigs

In the present study 3 connected F(2) crosses were used to map QTL for classical fat traits as well as fat-related metabolic and cytological traits in pigs. The founder breeds were Chinese Meishan, European Wild Boar, and Pietrain with to some extent the same founder animals in the different crosses. The different selection history of the breeds for fatness traits as well as the connectedness of the crosses led to a high statistical power. The total number of F(2) animals varied between 694 and 966, depending on the trait. The animals were genotyped for around 250 genetic markers, mostly microsatellites. The statistical model was a multi-allele, multi-QTL model that accounted for imprinting. The model was previously introduced from plant breeding experiments. The traits investigated were backfat depth and fat area as well as relative number of fat cells with different sizes and 2 metabolic traits (i.e., soluble protein content as an indicator for the level of metabolic turnover and NADP-malate dehydrogenase as an indicator for enzyme activity). The results revealed in total 37 significant QTL on chromosomes 1, 2, 4, 5, 6, 7, 8, 9, 14, 17, and 18, with often an overlap of confidence intervals of several traits. These confidence intervals were in some cases remarkably small, which is due to the high statistical power of the design. In total, 18 QTL showed significant imprinting effects. The small and overlapping confidence intervals for the classical fatness traits as well as for the cytological and metabolic traits enabled positional and functional candidate gene identification for several mapped QTL.

Identification of the expressed MHC class II DQB gene of the Asiatic black bear, Ursus thibetanus, in Japan

Genetic diversity estimation of the major histocompatibility complex (MHC) gene may be an important tool in the assessment of immune response ability against infectious disease. We were able to identify a near full-length expressed DQB sequence by RACE-PCR method from the Asiatic black bear, Ursus thibetanus in Japan. This is the first such full length expression in the Ursidae. The bear had at least one functional DQB locus. In phylogenetic tree analysis its DQB amino acid sequence formed a monophyletic group with DQB sequences from members of the order Carnivora and had a 90% nucleotide sequence similarity with the DQB allele of the California sea lion, Zalophus californianus. We compared the DQB amino acid composition of U. thibetanus with those of several other mammalian species including Homo sapiens. Amino acid residues known to be functionally important for human MHC genes, tended to be also conserved among other mammalian species while PBRs in the beta1 domain were heterogeneous among mammalian species. The DQB sequence obtained from the bear had not only no putative frameshifts or deletions but also no abnormal amino acid mutations such as had been observed in human DQB molecules. This suggests that the bear DQB sequence was an apparently functional DQB allele. As a preliminary study, we sequenced the exon 2 region of DQB alleles from genomic DNA, and succeeded to amplify the exon 2 of DQB loci. Our study will provide useful information for conservation genetics of the U. thibetanus as well as more generally regarding the mammalian MHC region.

Haplotypes of the porcine peroxisome proliferator-activated receptor delta gene are associated with backfat thickness

BACKGROUND

Peroxisome proliferator-activated receptor delta belongs to the nuclear receptor superfamily of ligand-inducible transcription factors. It is a key regulator of lipid metabolism. The peroxisome proliferator-activated receptor delta gene (PPARD) has been assigned to a region on porcine chromosome 7, which harbours a quantitative trait locus for backfat. Thus, PPARD is considered a functional and positional candidate gene for backfat thickness. The purpose of this study was to test this candidate gene hypothesis in a cross of breeds that were highly divergent in lipid deposition characteristics.

RESULTS

Screening for genetic variation in porcine PPARD revealed only silent mutations. Nevertheless, significant associations between PPARD haplotypes and backfat thickness were observed in the F2 generation of the Mangalitsa x Piétrain cross as well as a commercial German Landrace population. Haplotype 5 is associated with increased backfat in F2 Mangalitsa x Piétrain pigs, whereas haplotype 4 is associated with lower backfat thickness in the German Landrace population. Haplotype 4 and 5 carry the same alleles at all but one SNP. Interestingly, the opposite effects of PPARD haplotypes 4 and 5 on backfat thickness are reflected by opposite effects of these two haplotypes on PPAR-delta mRNA levels. Haplotype 4 significantly increases PPAR-delta mRNA levels, whereas haplotype 5 decreases mRNA levels of PPAR-delta.

CONCLUSION

This study provides evidence for an association between PPARD and backfat thickness. The association is substantiated by mRNA quantification. Further studies are required to clarify, whether the observed associations are caused by PPARD or are the result of linkage disequilibrium with a causal variant in a neighbouring gene.

A robust linkage map of the porcine autosomes based on gene-associated SNPs

BACKGROUND

Genetic linkage maps are necessary for mapping of mendelian traits and quantitative trait loci (QTLs). To identify the actual genes, which control these traits, a map based on gene-associated single nucleotide polymorphism (SNP) markers is highly valuable. In this study, the SNPs were genotyped in a large family material comprising more than 5,000 piglets derived from 12 Duroc boars crossed with 236 Danish Landrace/Danish Large White sows. The SNPs were identified in sequence alignments of 4,600 different amplicons obtained from the 12 boars and containing coding regions of genes derived from expressed sequence tags (ESTs) and genomic shotgun sequences.

RESULTS

Linkage maps of all 18 porcine autosomes were constructed based on 456 gene-associated and six porcine EST-based SNPs. The total length of the averaged-sex whole porcine autosome was estimated to 1,711.8 cM resulting in an average SNP spacing of 3.94 cM. The female and male maps were estimated to 2,336.1 and 1,441.5 cM, respectively. The gene order was validated through comparisons to the cytogenetic and/or physical location of 203 genes, linkage to evenly spaced microsatellite markers as well as previously reported conserved synteny. A total of 330 previously unmapped genes and ESTs were mapped to the porcine autosome while ten genes were mapped to unexpected locations.

CONCLUSION

The linkage map presented here shows high accuracy in gene order. The pedigree family network as well as the large amount of meiotic events provide good reliability and make this map suitable for QTL and association studies. In addition, the linkage to the RH-map of microsatellites makes it suitable for comparison to other QTL studies.

Molecular characterization, expression profile and association analysis with fat deposition traits of the porcine APOM gene

Apolipoprotein M (APOM), a novel apolipoprotein presented mostly in high-density lipoprotein (HDL) in plasma, is involved in lipid and lipoprotein metabolism. Through comparative mapping, we have mapped this gene to SSC7 p1.1 in which many QTLs affecting fat deposition traits have been reported. As a candidate gene for fat deposition traits, in this study, we obtained the 742-bp mRNA sequence of porcine APOM including the full coding region and encoding a protein of 188 amino acids. The sequence was deposited into the GenBank under the accession no. DQ329240. Semi-quantitative RT-PCR results showed that the porcine APOM gene is expressed predominantly in liver and kidney tissue. The genomic sequence of this gene which contains six exons and five introns, is 3,621 bp in length (DQ272488). Bioinformatic analysis of the 5' regulatory region has revealed that classical TATA-box element and species conserved Hepatocyte nuclear factor-1a (HNF-1alpha) biding site were represented in this region. A G2289C single nucleotide polymorphism (SNP) in the intron 2 of porcine APOM gene detected as an Eco130I PCR-restriction fragment length polymorphism (PCR-RFLP) showed allele frequency differences among three purebreds. Association of the genotypes with fat deposition traits showed that different genotypes of porcine APOM gene were significantly associated with leaf fat weight (P < 0.05), backfat thickness at shoulder (P < 0.05), backfat thickness at thorax-waist (P < 0.05), backfat thickness at buttock (P < 0.01) and average backfat thickness over shoulder, thorax-waist and buttock (P < 0.01).

A BAC contig map over the proximal approximately 3.3 Mb region of horse chromosome 21

A total of 207 BAC clones containing 155 loci were isolated and arranged into a map of linearly ordered overlapping clones over the proximal part of horse chromosome 21 (ECA21), which corresponds to the proximal half of the short arm of human chromosome 19 (HSA19p) and part of HSA5. The clones form two contigs - each corresponding to the respective human chromosomes - that are estimated to be separated by a gap of approximately 200 kb. Of the 155 markers present in the two contigs, 141 (33 genes and 108 STS) were generated and mapped in this study. The BACs provide a 4-5x coverage of the region and span an estimated length of approximately 3.3 Mb. The region presently contains one mapped marker per 22 kb on average, which represents a major improvement over the previous resolution of one marker per 380 kb obtained through the generation of a dense RH map for this segment. Dual color fluorescence in situ hybridization on metaphase and interphase chromosomes verified the relative order of some of the BACs and helped to orient them accurately in the contigs. Despite having similar gene order and content, the equine region covered by the contigs appears to be distinctly smaller than the corresponding region in human (3.3 Mb vs. 5.5-6 Mb) because the latter harbors a host of repetitive elements and gene families unique to humans/primates. Considering limited representation of the region in the latest version of the horse whole genome sequence EquCab2, the dense map developed in this study will prove useful for the assembly and annotation of the sequence data on ECA21 and will be instrumental in rapid search and isolation of candidate genes for traits mapped to this region.

Genomic structure and gene order of swine chromosome 7q1.1q1.2

To clarify the structure of the porcine genomic region that contains quantitative trait loci (QTL) related to fat, we constructed a bacterial artificial chromosome (BAC) contig of the region from DST to SRPK1 on porcine chromosome 7 and performed low-redundancy 'skim' shotgun sequencing of the clones that composed a minimum tiling path of the contig. This analysis revealed that the gene order from VPS52 to SRPK1 is conserved between human and swine and that comparison with the human sequence identified a rearrangement in the swine genome at the proximal end of VPS52. Analysis of the nucleotide sequences of three BAC clones that included the rearrangement point demonstrated that COL21A1 and DST, which were not present in the corresponding human region, were located adjacent to the rearrangement point. These results provide useful information about the genomic region containing QTL for fat in pigs and help to clarify the structure of the so-called 'extended-class II' region distal to the porcine major histocompatibility complex class II region.

High resolution physical map of porcine chromosome 7 QTL region and comparative mapping of this region among vertebrate genomes

Background

On porcine chromosome 7, the region surrounding the Major Histocompatibility Complex (MHC) contains several Quantitative Trait Loci (QTL) influencing many traits including growth, back fat thickness and carcass composition. Previous studies highlighted that a fragment of ~3.7 Mb is located within the Swine Leucocyte Antigen (SLA) complex. Internal rearrangements of this fragment were suggested, and partial contigs had been built, but further characterization of this region and identification of all human chromosomal fragments orthologous to this porcine fragment had to be carried out.

Results

A whole physical map of the region was constructed by integrating Radiation Hybrid (RH) mapping, BAC fingerprinting data of the INRA BAC library and anchoring BAC end sequences on the human genome. 17 genes and 2 reference microsatellites were ordered on the high resolution IMNpRH2_12000rad Radiation Hybrid panel. A 1000:1 framework map covering 550 cR₁₂₀₀₀ was established and a complete contig of the region was developed. New micro rearrangements were highlighted between the porcine and human genomes. A bovine RH map was also developed in this region by mapping 16 genes. Comparison of the organization of this region in pig, cattle, human, mouse, dog and chicken genomes revealed that 1) the translocation of the fragment described previously is observed only on the bovine and porcine genomes and 2) the new internal micro rearrangements are specific of the porcine genome.

Conclusion

We estimate that the region contains several rearrangements and covers 5.2 Mb of the porcine genome. The study of this complete BAC contig showed that human chromosomal fragments homologs of this heavily rearranged QTL region are all located in the region of HSA6 that surrounds the centromere. This work allows us to define a list of all candidate genes that could explain these QTL effects.

Nomenclature for factors of the swine leukocyte antigen class II system, 2005

A systematic nomenclature for the genes and alleles of the swine major histocompatibility complex (MHC) is essential to the development and communication of research in swine immunology. The Swine Leukocyte Antigen (SLA) Nomenclature Committee of the International Society for Animal Genetics (ISAG) has reviewed all of the DNA-sequence information for MHC class II genes, available in GenBank/EMBL/DDBJ databases, and the associated published reports to develop such a systematic nomenclature. This article summarizes the proposed nomenclature, which parallels the World Health Organization's nomenclature for factors of the human MHC. The SLA class II genes expressed on the cell membrane will be noted as SLA-DRA, SLA-DRB1, SLA-DQA, and SLA-DQB1. Nomenclature assignments for all SLA class II GenBank sequences are now noted. The committee will add new SLA class II allele designations, as they are discovered, and will maintain a publicly available list of all recognized genes and alleles using the Immuno Polymorphism Database (IPD). The sequences will be available from the IPD-MHC section of the database which contains non-human MHC sequences (http://www.ebi.ac.uk/ipd/mhc/sla/).

Comparative analysis of a BAC contig of porcine chromosome 13q31-q32 and human chromosome 3q21-q22

Background

The gene(s) encoding the ETEC F4ab/ac receptors, involved in neonatal diarrhoea in pigs (a disease not yet described in humans), is located close to the TF locus on Sscr13. In order to reveal and characterize possible candidate genes encoding these receptors, a porcine physical map of the TF region is indispensable.

Results

A contig of 33 BAC clones, covering approximately 1.35 Mb surrounding the TF locus on Sscr13q31-q32, was built by chromosome walking. A total of 22,552 bp from the BAC contig were sequenced and compared with database sequences to identify genes, ESTs and repeat sequences, and to anchor the contig to the syntenic region of the human genome sequence (Hsap3q21-q22). The contig was further annotated based on this human/porcine comparative map, and was also anchored to the Sanger porcine framework map and the integrated map of Sscr13 by RH mapping.

Conclusion

The annotated contig, containing 10 genes and 2 ESTs, showed a complete conservation of linkage (gene order and orientation) with the human genome sequence, based on 46 anchor points. This underlines the importance of the human/porcine comparative map for the identification of porcine genes associated with genetic defects and economically important traits, and for assembly of the porcine genome sequence.