Skip to main content

Transcription variants of SLA-7, a swine non classical MHC class I gene

Abstract

In pig, very little information is available on the non classical class I (Ib) genes of the Major Histocompatibility Complex (MHC) i.e. SLA-6, -7 and -8. Our aim was to focus on the transcription pattern of the SLA-7 gene. RT-PCR experiments were carried out with SLA-7 specific primers targeting either the full coding sequence (CDS) from exon 1 to the 3 prime untranslated region (3UTR) or a partial CDS from exon 4 to the 3UTR. We show that the SLA-7 gene expresses a full length transcript not yet identified that refines annotation of the gene with eight exons instead of seven as initially described from the existing RefSeq RNA. These two RNAs encode molecules that differ in cytoplasmic tail length. In this study, another SLA-7 transcript variant was characterized, which encodes a protein with a shorter alpha 3 domain, as a consequence of a splicing site within exon 4. Surprisingly, a cryptic non canonical GA-AG splicing site is used to generate this transcript variant. An additional SLA-7 variant was also identified in the 3UTR with a splicing site occurring 31 nucleotides downstream to the stop codon. In conclusion, the pig SLA-7 MHC class Ib gene presents a complex transcription pattern with two transcripts encoding various molecules and transcripts that do not alter the CDS and may be subject to post-transcriptional regulation.

Background

The Major Histocompatibility Complex (MHC) class I gene family comprises classical (Ia) and non classical (Ib) genes. The highly polymorphic class Ia genes are widely expressed and encode membrane-bound glycoproteins that present self and viral peptides to cytotoxic T cells [1] and modulate the activity of natural killer cells [2]. In contrast, the class Ib genes display limited polymorphism, and are predominantly expressed in immunotolerant organ sites in human, notably at the feto-maternal interface [3]. In man, three MHC class Ib genes have been characterized, namely HLA-E, -F and -G[3] and HLA-G has been shown to express alternatively spliced variants encoding various membrane-bound as well as soluble proteins [4]. In mouse, the H2-QaI gene is orthologous to HLA-E[5] and functional homologies have been established between H2-Qa2 and HLA-G[6]. One to four MHC class Ib genes have been identified in rat according to haplotypes [7] and four MHC class Ib genes have been characterized in cattle [8]. There is a growing interest in addressing the role of the MHC class Ib genes in the species where they are characterized. Indeed, MHC class Ia genes seem to share similar functions across species but the MHC class Ib genes are good candidates to address questions on both shared and species-specific immunity-related roles.

In pig, very limited information is available on the MHC class Ib genes SLA-6, -7 and -8. The three genes have been fully sequenced from the homozygous Hp1a.0 haplotype [9, 10]. Nine allelic variants have been reported for SLA-6 and only two for SLA-7 or SLA-8[11]. It has been shown that SLA-Ib genes are expressed in a less restricted manner than the HLA-Ib genes [12, 13] despite a predominant transcription in the lymphoid organs, the lung and the digestive tract [13]. In addition, conversely to the SLA-Ia genes, transfection experiments have revealed that the promoters of SLA-7 and SLA-6 do not respond to interferon, suggesting distinct regulatory systems for pig MHC class Ia and Ib genes, as in human [14]. Our aim was to focus on the transcription of the SLA-7 gene known to have a unique reference transcript [12]. In this report, we show that the SLA-7 gene expresses a full-length transcript not yet identified as well as at least two additional alternative spliced variants that lead to either exon alteration in the resulting protein or modification of the 3’end of the transcript.

Methods

Animals, tissues and RNA extraction

Tissues from Melanoma-bearing Libechov Minipigs (MeLiM) [15] and French Large White pigs were used. The tissues from MeLiM pigs have been sampled on 13 months old animals. At the time of tissue sampling, all MeLiM animals had regressed, meaning that they were not bearing melanomas anymore [15]. Tissues included brain, thymus, tonsil, spleen and liver. Total RNA was extracted using QIAGEN RNeasy Mini Kits (Qiagen, France). All RNA samples were purified by on-column digestion of DNA with DNase I as recommended by the manufacturer (Qiagen, France).

Primer design

Three primers were designed from the SLA-7 reference cDNA [12] and genomic [9, 10] sequences, using the Primer3 online program [16]. The primer combinations were suitable to amplify the full coding sequence from exon 1 to the three prime untranslated region (3UTR) or a partial coding sequence from exon 4 to the 3UTR (Table 1 and Figure 1). Primers were also designed to amplify cDNAs of the RPL32 gene that has been used as control gene for expression levels (Table 1 and Figure 1).

Table 1 RT-PCR primers
Figure 1
figure 1

Schematic representation of the SLA-7 transcripts. Exons are numbered E1 to E8 and the three prime non coding region is indicated as 3UTR. Sizes of exons and 3UTR are given in nucleotides within the boxes. The number of aminoacids for each exon is indicated above the exon number. Exons represented by a dark blue box (E1) correspond to the leader sequences. Exons represented by orange boxes (E2 to E4) stand for the alpha1, 2 and 3 domains of the molecules. Exons represented in bright blue (E5) correspond to the transmembrane domain. Exons represented by green boxes encode the cytoplasmic tail of the molecule. The 3UTR is represented by a grey box. Positions of the primers used for RT-PCRs are indicated by arrows on top of the figure on E1 (SLA-7-e1-F), E4 (SLA-7-e4-F) and 3UTR (SLA-7-3UTR-R). The donor and acceptor splice sequences are positioned by arrows on E4 and 3UTR boxes.

RT-PCR and sequencing

Two micrograms of DNaseI-treated total RNA were reverse-transcribed (Superscript II enzyme, Invitrogen, USA) with Oligo (dT) primers in a final volume of 20 µL to which 30 µL of water were further added to prepare the stock solution of RT samples. PCRs were carried out in a final volume of 15 µL using 100 nM of each primer, 1 µL of the 1:10 RT sample and the GoTaq™ DNA polymerase (Promega, USA). Thermocycling conditions were as follows: 94°C for 3 min, followed by 35 amplification cycles at 94°C for 30 sec, 60°C for 30 sec and 72°C for 90 sec, followed by a final extension at 72°C for 5 min. The PCR products were purified using the JETQUICK Gel Extraction Spin Kit (Genomed, Germany) for further cloning into pCR2.1 vector (TA Cloning Kit, Invitrogen, USA) and sequencing (Eurofins MWG Operon, France).

Sequence analysis

Sequence similarities were searched with the BLAST tools [17]. Multiple alignments were carried out with CLUSTALW [18]. cDNA sequences were translated to protein by online DNA to Protein translation tool (http://bio.lundberg.gu.se/edu/translat.html).

Results and discussion

SLA-7 full coding sequences

Full length SLA-7 transcripts were characterized by RT-PCR from the thymus of MeLiM pigs using the primers SLA-7-e1-F and SLA-7-3UTR-R (Table 1 and Figure 1). A 1465 nucleotides long transcript was obtained and further named SLA-7-1465 (Accession number: GU322918). Annotation was carried out by aligning the cDNA sequence to the genomic reference sequence (GenBank accession number AJ251914) and eight exons were detected in this new transcript, in contrast to the reference full-length transcript (Accession number NM_213768) that harbours only seven exons [12] (Figure 1) and is referred to as SLA-7-001 (OTTSUST00000000782) in the Vertebrate Genome Annotation database [19]. The two RNAs encode proteins that differ in the cytoplasmic tail (Figure 2). The SLA-7-001 encoded protein contains a cytoplasmic tail that is defined by exons 6 and 7 and is 68 aminoacids long. The SLA-7-1465 encoded protein is characterized by a cytoplasmic tail that is defined by exons 6 to 8 and is 55 aminoacids long. It has been demonstrated that the cytoplasmic tail of MHC class I molecules contributes to their expression on the cell surface [20] and that mutations of cysteine residues in the cytoplasmic tail of MHC class Ia molecules modify extracellular recognition by Leukocyte Ig-Like receptor 1 [21]. Moreover, it has been reported that HLA-F molecules are entirely dependent on the cytoplasmic tail for export from the endoplasmic reticulum to the Golgi apparatus [22]. Altogether, these reports strongly support a major role for the cytoplasmic tail of MHC class I molecules in transport and function. Further experiments are required to study whether the SLA-7 molecules encoded by SLA-7-001 or SLA-7-1465 transcripts have distinct properties due to their different cytoplasmic tails.

Figure 2
figure 2

Multi-alignment of peptides encoded by the transcripts SLA-7-1465, SLA-7-1366 and SLA-7-001. The successive eight (SLA-7-1366 and SLA-7-1465) or seven (SLA-7-001) exons are alternatively indicated by black and blue font. Aminoacids at the junction between two exons are in grey boxes. Aminoacid similarities between two or three sequences are indicated below the sequence alignments by dots or stars, respectively.

An SLA-7 spliced variant encoding a protein with a shorter alpha 3 domain

A 1366 nucleotides long transcript was retrieved from brain RNA and further referred to as SLA-7-1366 (accession number: HQ224544). Surprisingly, annotation of the cDNA revealed the presence of nine exons due to a splicing site within exon 4 (Figures 1 and 2). The two exons matching to exon 4 were named exons 4a and 4b (Figure 1). Alignment of the SLA-7-1366 cDNA to the reference genomic sequence showed that between exons 4a and 4b, the donor and acceptor splice sites were GA and AG, respectively. This finding indicates that a cryptic non canonical splicing code is used to express this SLA-7 transcript variant. The general rule is the use of GT and AG for donor and acceptor splicing sites, respectively [23], but alternative codes may be functional [24]. It has been shown that the GA-AG splicing site is rarely used and a few cases have been reported among which splicing in the human parafibromin gene [25]. Our results suggest that the SLA-7 gene may be subject to subtle regulation resulting in the use of rarely used non canonical splicing sites. Additional studies are required to analyze whether this regulation is tissue-specific.

The SLA-7-1366 and SLA-7-1465 encoded molecules with different alpha 3 domain lengths (figures 1 and 2) i.e. 59 (39 from exon 4a and 20 from exon 4b) and 92 aminoacids long, respectively. The alpha 3 domain corresponds to the Immunoglobulin-like region and interacts with the cell surface CD8 glycoproteins that are expressed on cytotoxic T lymphocytes and function as a co-receptor with the T cell receptor [26]. Expression of SLA-7 molecules on the cell surface has not been demonstrated. However, the alpha 3 domain encoded by the SLA-7-1366 transcript is shortened by comparison to the full-length molecule, suggesting that such a modification may alter interactions with cell receptors.

A spliced variant in the 3’UTR with no alteration of the encoded protein

By using primers targeting the three prime end of the gene from exon 4 (SLA-7-e4-F) to the 3UTR (SLA-7-3UTR-R) (Table 1 and figure 1), two partial transcripts were recovered that differ in non coding sequence length (Figure 1). The 650 nucleotide long sequence (SLA-7-650) corresponds to the expected full-length sequence. The 464 nucleotide long sequence SLA-7-464 (accession number: GU322919) is the result of an alternative splicing within the 3UTR region, 31 nucleotides downstream to the stop codon. Modification of the 3UTR does not affect the encoded molecules. As indicated in figure 1, the canonical GT-AG rule was used for this splicing. A new category of transcripts has been recently characterized that can be subject to non sense mediated decay (NMD) [27, 28]. Variants targeted by NMD can present alternative splicing in the 3UTR and the distance between the stop codon and the splice site has been shown to be more than 50 nucleotides long [29]. The SLA-7-464 variant cannot fall into this category of transcripts subject to NMD because the distance between the stop codon and the splicing site in only 31 nucleotides. However, it is tempting to hypothesize that SLA-7-464 variants are subject to a post-transcriptional regulation that has to be explored.

The expression patterns of these two 3UTR variants were studied in four different tissues including spleen, thymus, tonsil, and liver from MeLiM and Large White pigs. Surprisingly, the SLA-7-650 band was detected in all tissues of both breeds but the SLA-7-464 band was detected only in MeLiM pigs (Figure 3). We cannot rule out a very weak expression of the short variant in Large White pigs but our results strongly suggest a co-expression of both variants in MeLiM pigs and a predominant expression of the SLA-7-650 variant in the Large White pigs included in our study.

Figure 3
figure 3

Tissue expression patterns in MeLiM and Large White pigs. Detection of the SLA-7 transcript variants in adult tissues from MeLiM (A) and Large White (B) pigs by RT-PCR using the primers SLA-7-e4-F and SLA-7-3UTR-R. The RPL32 gene was used as a control for expression levels as shown for four tissues of MeLiM pig (C).

Conclusion and perspectives

We have identified an SLA-7 full length transcript that had not been characterized before and that differs from the reference sequence by the length of the encoded cytoplasmic tail. In addition, we show that the SLA-7 gene is subject to alternative splicing transcription that leads to either a transcript encoding a molecule with a shortened alpha 3 domain or a transcript that is spliced in the 3UTR after the stop codon. In conclusion, the non classical MHC class Ib gene SLA-7 gene presents a complex transcription pattern, the regulation of which needs to be further investigated. The functions of the putative encoded molecules also need to be studied.

Abbreviations

MHC:

Major Histocompatibility Complex

HLA:

Human Leukocyte Antigen

SLA:

Swine Leukocyte Antigen

MeLiM:

Melanoblastoma-bearing Libechov Minipigs

3UTR:

3 prime untranslated region

Class Ib:

non classical class I

Class Ia:

classical class I

CDS:

coding sequence

NMD:

Non sense Mediated Decay

PCR:

Polymerase Chain Reaction

RT:

reverse transcription.

References

  1. Townsend A, Bodmer H: Antigen recognition by class I-restricted T lymphocytes. Annu Rev Immunol. 1989, 7: 601-624. 10.1146/annurev.iy.07.040189.003125.

    Article  CAS  PubMed  Google Scholar 

  2. Trowsdale J, Barten R, Haude A, Stewart CA, Beck S, Wilson MJ: The genomic context of natural killer receptor extended gene families. Immunol Rev. 2001, 181: 20-38. 10.1034/j.1600-065X.2001.1810102.x.

    Article  CAS  PubMed  Google Scholar 

  3. Moscoso J, Serrano-Vela JI, Pacheco R, Arnaiz-Villena A: HLA-G, -E and -F: allelism, function and evolution. Transpl Immunol. 2006, 17 (1): 61-64. 10.1016/j.trim.2006.09.010.

    Article  CAS  PubMed  Google Scholar 

  4. Carosella ED, Moreau P, Le Maoult J, Le Discorde M, Dausset J, Rouas-Freiss N: HLA-G molecules: from maternal-fetal tolerance to tissue acceptance. Adv Immunol. 2003, 81: 199-252. full_text.

    Article  CAS  PubMed  Google Scholar 

  5. Joly E, Rouillon V: The orthology of HLA-E and H2-Qa1 is hidden by their concerted evolution with other MHC class I molecules. Biol Direct. 2006, 1: 2-10.1186/1745-6150-1-2.

    Article  PubMed Central  PubMed  Google Scholar 

  6. Comiskey M, Goldstein CY, De-Fazio SR, Mammolenti M, Newmark JA, Warner CM: Evidence that HLA-G is the functional homolog of mouse Qa-2, the Ped gene product. Hum Immunol. 2003, 64: 999-1004. 10.1016/j.humimm.2003.08.352.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Lau P, Amadou C, Brun H, Rouillon V, McLaren F, Le Rolle AF, Graham M, Butcher GW, Joly E: Characterisation of RT1-E2, a multigenic family of highly conserved rat non-classical MHC class I molecules initially identified in cells from immunoprivileged sites. BMC Immunology. 2003, 4: 7-10.1186/1471-2172-4-7.

    Article  PubMed Central  PubMed  Google Scholar 

  8. Birch J, Codner G, Guzman E, Ellis SA: Genomic location and characterisation of nonclassical MHC class I genes in cattle. Immunogenetics. 2008, 60 (5): 267-273. 10.1007/s00251-008-0294-2.

    Article  CAS  PubMed  Google Scholar 

  9. Chardon P, Rogel-Gaillard C, Cattolico L, Duprat S, Vaiman M, Renard C: Sequence of the swine major histocompatibility complex region containing all non-classical class I genes. Tissue Antigens. 2001, 57 (1): 55-65. 10.1034/j.1399-0039.2001.057001055.x.

    Article  CAS  PubMed  Google Scholar 

  10. Renard C, Hart E, Sehra H, Beasley H, Coggill P, Howe K, Harrow J, Gilbert J, Sims S, Rogers J, Ando A, Shigenari A, Shiina T, Inoko H, Chardon P, Beck S: The genomic sequence and analysis of the swine major histocompatibility complex. Genomics. 2006, 88: 96-110. 10.1016/j.ygeno.2006.01.004.

    Article  CAS  PubMed  Google Scholar 

  11. Ho CS, Lunney JK, Ando A, Rogel-Gaillard C, Lee JH, Schook LB, Smith DM: Nomenclature for factors of the SLA system, update 2008. Tissue Antigens. 2009, 73: 307-15. 10.1111/j.1399-0039.2009.01213.x.

    Article  CAS  PubMed  Google Scholar 

  12. Crew MD, Phanavanh B, Garcia-Borges CN: Sequence and mRNA expression of nonclassical SLA class I genes SLA-7 and SLA-8. Immunogenetics. 2004, 56 (2): 111-114. 10.1007/s00251-004-0676-z.

    Article  CAS  PubMed  Google Scholar 

  13. Kusza S, Flori L, Gao Y, Teillaud A, Hu R, Lemonnier G, Bõsze Z, Bourneuf E, Vincent-Naulleau S, Rogel-Gaillard C: Transcription specificity of the class Ib genes SLA-6, SLA-7 and SLA-8 of the swine major histocompatibility complex and comparison with class Ia genes. Anim Genet. 2010.

    Google Scholar 

  14. Tennant LM, et al: Regulation of porcine classical and nonclassical MHC class I expression. Immunogenetics. 2007, 59 (5): 377-389. 10.1007/s00251-007-0206-x.

    Article  CAS  PubMed  Google Scholar 

  15. Vincent-Naulleau S, Le-Chalony C, Leplat JJ, Bouet S, Bailly C, Spatz A, Vielh P, Avril MF, Tricaud Y, Gruand J, et al: Clinical and histopathological characterization of cutaneous melanomas in the melanoblastoma-bearing Libechov minipig model. Pigment Cell Res. 2004, 17 (1): 24-35. 10.1046/j.1600-0749.2003.00101.x.

    Article  PubMed  Google Scholar 

  16. Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000, 132: 365-386.

    CAS  PubMed  Google Scholar 

  17. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.

    Article  CAS  PubMed  Google Scholar 

  18. Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22 (22): 4673-4680. 10.1093/nar/22.22.4673.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. VEGA: [http://vega.sanger.ac.uk/index.html]

  20. Davis DM, Reyburn HT, Pazmany L, Chiu I, Mandelboim O, Strominger JL: Impaired spontaneous endocytosis of HLA-G. Eur J Immunol. 1997, 27 (10): 2714-2719. 10.1002/eji.1830271035.

    Article  CAS  PubMed  Google Scholar 

  21. Gruda R, Achdout H, Stern-Ginossar N, Gazit R, Betser-Cohen G, Manaster I, Katz G, Gonen-Gross T, Tirosh B, Mandelboim O: Intracellular cysteine residues in the tail of MHC class I proteins are crucial for extracellular recognition by leukocyte Ig-like receptor 1. J Immunol. 2007, 179 (6): 3655-3661.

    Article  CAS  PubMed  Google Scholar 

  22. Boyle LH, Gillingham AK, Munro S, Trowsdale J: Selective export of HLA-F by its cytoplasmic tail. J Immunol. 2006, 176 (11): 6464-6472.

    Article  CAS  PubMed  Google Scholar 

  23. Mount SM: A catalogue of splice junction sequences. Nucleic Acids Res. 1982, 10 (2): 459-72. 10.1093/nar/10.2.459.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Burset M, Seledtsov IA, Solovyev VV: Analysis of canonical and non-canonical splice sites in mammalian genomes. Nucleic Acids Res. 2000, 28 (21): 4364-4375. 10.1093/nar/28.21.4364.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Bradley KJ, Cavaco BM, Bowl MR, Harding B, Young A, Thakker RV: Utilisation of a cryptic non-canonical donor splice site of the gene encoding PARAFIBROMIN is associated with familial isolated primary hyperparathyroidism. J Med Genet. 2005, 42 (8): e51-10.1136/jmg.2005.032201.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Sun J, Kavathas PB: Comparison of the roles of CD8 alpha alpha and CD8 alpha beta in interaction with MHC class I. J Immunol. 1997, 159 (12): 6077-6082.

    CAS  PubMed  Google Scholar 

  27. Barash Y, Calarco JA, Gao W, Pan Q, Wang X, Shai O, Blencowe BJ, Frey BJ: Deciphering the splicing code. Nature. 2010, 465 (7294): 53-59. 10.1038/nature09000.

    Article  CAS  PubMed  Google Scholar 

  28. Gardner LB: Nonsense-mediated RNA decay regulation by cellular stress: implications for tumorigenesis. Mol Cancer Res. 2010, 8 (3): 295-308. 10.1158/1541-7786.MCR-09-0502.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  29. Mendell JT, Sharifi NA, Meyers JL, Martinez-Murillo F, Dietz HC: Nonsense surveillance regulates expression of diverse classes of mammalian transcripts and mutes genomic noise. Nat Genet. 2004, 36 (10): 1073-1078. 10.1038/ng1429.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors would like to thank Dr. Patrick Chardon (GABI-GIS, INRA, France) for useful discussion and Hélène Hayes for English correction of the manuscript. Rui Hu’s PhD was supported by the Animal Genetics Department of INRA and Dadi Agriculture Comprehensive Exploitation Co. Ltd, Liaoning, China.

This article has been published as part of BMC Proceedings Volume 5 Supplement 4, 2011: Proceedings of the International Symposium on Animal Genomics for Animal Health (AGAH 2010). The full contents of the supplement are available online at http://www.biomedcentral.com/1753-6561/5?issue=S4.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Claire Rogel-Gaillard.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

RH participated in working out the experimental design, performed all experiments and sequence analysis, and drafted the manuscript. GL carried out RT-PCR experiments. EB and SVN contributed to tissue sampling. CRG coordinated the whole study, contributed to the experimental design, the analysis and interpretation of the results and corrected the manuscript. All authors have read and approved the final manuscript.

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Hu, R., Lemonnier, G., Bourneuf, E. et al. Transcription variants of SLA-7, a swine non classical MHC class I gene. BMC Proc 5 (Suppl 4), S10 (2011). https://doi.org/10.1186/1753-6561-5-S4-S10

Download citation

  • Published:

  • DOI: https://doi.org/10.1186/1753-6561-5-S4-S10

Keywords