Next Article in Journal
Evaluation of Water Scarcity Footprint for Taiwanese Dairy Farming
Next Article in Special Issue
Prophylactic and Therapeutic Efficacy of Prebiotic Supplementation against Intestinal Coccidiosis in Rabbits
Previous Article in Journal
Growth Performance and Carcass Traits Responses to Dried Distillers’ Grain with Solubles Feeding of Growing Awassi Ram Lambs
Previous Article in Special Issue
Effect of Silver Nanoparticle Administration on Productive Performance, Blood Parameters, Antioxidative Status, and Silver Residues in Growing Rabbits under Hot Climate
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Diversity of IGHM and IGHE in the Leporids Revealed Different Patterns of Diversity in the Two European Rabbit Subspecies (O. cuniculus algirus and O. c. cuniculus)

1
CIBIO Centro de Investigação em Biodiversidade e Recursos Genéticos, InBio Laboratório Associado, Universidade do Porto, Campus Agrário de Vairão, 4485-661 Vairão, Portugal
2
Departamento de Biologia, Faculdade de Ciências, Universidade do Porto, 4169-007 Porto, Portugal
3
CESPU, Instituto de Investigação e Formação Avançada em Ciências e Tecnologias da Saúde, 4585-116 Gandra PRD, Portugal
*
Author to whom correspondence should be addressed.
Animals 2019, 9(11), 955; https://doi.org/10.3390/ani9110955
Submission received: 7 October 2019 / Revised: 8 November 2019 / Accepted: 8 November 2019 / Published: 12 November 2019
(This article belongs to the Special Issue Disease and Immunology of Rabbits)

Abstract

:

Simple Summary

The study of European rabbit immunoglobulin genes has contributed decisively to the current knowledge on antibody structure and diversification. The European rabbit has also been increasingly used as an animal model for the study of many human diseases, such as syphilis, tuberculosis, and AIDS. As such, the study of its immune system genes is of crucial relevance, but the study of rabbit immunoglobulins has focused only on the IgG and IgA antibodies. In this study, we added to the knowledge of the rabbit immune system by investigating the genetic diversity of two antibodies, IgM and IgE, in wild and domestic rabbits as well as other rabbit close species. With the data obtained in this study, we showed a high similarity between the different rabbit close species studied and we pointed out important genetic differences in the wild and domestic rabbits. Our findings are a valuable tool for the management of rabbit wild populations and domestic breeds and may contribute to the identification of immunoglobulins genetic variants with greater efficiency against pathogens.

Abstract

The European rabbit (Oryctolagus cuniculus) has been an important model for immunological studies but the study of its immunoglobulins (Ig) has been restricted to its unique IgA and IgG. Here, we studied the genetic diversity of IgM and IgE in several species of leporids and performed population genetics studies on European rabbit wild populations and domestic breeds. The leporids sequencing showed that these Ig are well conserved (98% sequence similarity among leporids), For IgM the Cµ1 and Cµ4 were the most diverse and most conserved domains, respectively, while for IgE the Cε1 was the most diverse domain and Cε2 and Cε3 the most conserved domains. The differences in the pattern of most conserved and most diverse domain between the Ig isotypes are most likely related to each isotype function. The genetic population data showed contrasting results for IgM and IgE. For both Ig, as expected, a greater diversity was observed in the original species range, the Iberian Peninsula. However, unexpectedly the genetic diversity found for IgE in the domestic animals is higher than that for the French wild populations. These results will increase knowledge of the genetic diversity of leporids and wild and domestic rabbit populations and are important tools for the management of wild populations and rabbitries.

1. Introduction

Immunoglobulins (Igs) are one of the most emblematic components of the jawed vertebrate immune system. These heterodimeric glycoproteins act as a critical part of the immune system by recognizing and binding to antigens and subsequently triggering mechanisms that lead to its destruction [1]. In eutherian mammals, five Ig isotypes have been described, IgM, IgD, IgE, IgG and IgA. Differences between Ig isotypes reside in the constant domains of the heavy chain (CH) and determine that each isotype has specific effector functions [1].
The European rabbit (Oryctolagus cuniculus) has been used as an animal model for diseases and immunological studies (reviewed in [2]). One of the peculiarities of the European rabbit is the unique features of Ig isotypes. As some other mammals the European rabbit has no IgD and no functional IGHD gene has been found in the rabbit germline [3,4]. IgM and IgE are encoded by single genes, IGHM and IGHE, as well as IgG, a surprising fact among mammals that typically have multiple IGHG copies encoding multiple IgG subclasses [4,5]. For IgA, on the other hand, the European rabbit has at least 15 IGHA genes encoding 15 IgA subclasses [6,7] being the most complex IgA system of all studied mammals. These unique aspects may, perhaps, explain that the study of European rabbit Igs has focused on IgG and IgA (e.g., [6,7,8,9,10,11]) forgetting other Ig isotypes. In fact, the information available for European rabbit IgM and IgE is extremely scarce with only a few sequences available in public databases (Figure 1).
IgM, in its transmembrane form, is the first isotype to be expressed during B-cell development. As such the VH and VL regions associated with this Ig have not undergone much antigenic driven diversification being able to recognize a wider range of antigens than other isotypes, although showing low affinity [1]. In mammals, IgM is secreted in pentamers, in which the IgM monomers are linked by disulfide bonds in the CH4 domain and further bound a J chain. This pentameric form is very effective in opsonizing pathogens and IgM antibodies are also very efficient at fixing complement [12]. Immunoglobulin E (IgE) has the lowest serum concentration and shortest half-life of all Ig subclasses. Despite this, it is a very potent immunoglobulin associated with hypersensitivity, allergic reactions and response to parasitic worm infections [1].
The European rabbit originated on the Iberian Peninsula (IP) [13], where two subspecies coexist: O. cuniculus algirus in the south-west area of the IP and O. cuniculus cuniculus in the north east area. After the last glacial peak, O. c. cuniculus expanded its range north, across the Pyrennes, towards France [14], from where the species was domesticated in the last 1500 years [15,16]. It belongs to the family Leporidae which comprises 11 genera: Brachylagus, Bunolagus, Caprolagus, Lepus, Nesolagus, Oryctolagus, Pentalagus, Poelagus, Pronolagus, Romerolagus and Sylvilagus) [17,18]. Most are monotypic genera with restricted distributions, with the exception of Lepus, a polytypic and cosmopolitan genus, and Nesolagus, Pronolagus and Sylvilagus, which comprise two, four and more than 16 species, respectively [17,19,20,21,22,23,24,25]. Despite some discrepancies in leporid taxonomy are still unresolved, the widest consensus agrees that the Leporidae first basal group included Pronolagus, and that Romerolagus is at the second basal position [26,27,28].
In this study, we extend the knowledge on IgM and IgE isotypes in leporids by sequencing the complete IGHM and IGHE genes for six additional extant leporid genera: Bunolagus, Brachylagus, Lepus, Pentalagus, Romerolagus and Sylvilagus. Additionally, we compared the genetic diversity of these Igs between wild rabbit populations of O. c. algirus and wild and domestic populations of O. c. cuniculus.

2. Material and Methods

For the study of IGHE and IGHM evolution in leporids, we analysed specimens of Bunolagus, Brachylagus, Lepus, Pentalagus, Romerolagus and Sylvilagus genera. To assess the genetic diversity of European rabbit we analysed 10 individuals for each of three wild populations: O. c. algirus collected in Portugal and O. c. cuniculus collected in Spain and in France (Figure 2); as well as 15 domestic animals from 6 breeds: French Lop, Flemish Giant, Argent Champagne, Angora, New Zealand White and Netherland Dwarf. Total genomic DNA was extracted from frozen liver or ear tissue using an EasySpin Genomic DNA Minipreps Tissue Kit (Citomed). No animals were killed for this study. All samples used in this work belong to the Cibio/InBio sample collection and had been previously used in other works [11,29,30].
PCR amplifications of the four IGHE and IGHM exons was conducted using primers designed on the basis of European rabbit available sequences (GenBank accession number AY386696) [4]. A fragment containing the IGHE CH1 and CH2 domains was amplified using primers FE12 5’AGGTGGAAGCCAGGGTGAGG 3’ and RE21 5’ CGCCTCGCGGTTCGTAATCT 3’ under the following conditions: 15 min at 95 °C followed by 35 cycles at 95 °C (30 s), 60 °C (30 s) and 72 °C (45 s), with a final extension at 60 °C (20 min). Another fragment containing IGHE CH3 and CH4 exons was amplified using primers FE31 5’ TGACACCCAGAAAGACAGGG 3’ and RE41 5’ CCGAGTGACACTGCAGTGTT 3’ under the following conditions: 15 min at 95 °C followed by 35 cycles at 95 °C (30 s), 62 °C (30 s) and 72 °C (45 s), with a final extension at 60 °C (20 min). A fragment containing the IGHM CH1 and CH2 domains was amplified using primers FM11 5’ AGCTTTTCACACCTCCCCTT 3’ and RM21 5’ AAACCCATGAGGACGCCTGT 3’ under the following conditions: 15 min at 95 °C followed by 35 cycles at 95 °C (30 s), 62 °C (30 s) and 72 °C (45 s), with a final extension at 60 °C (20 min). Another fragment containing IGHM CH3 and CH4 exons was amplified using primers FM31 5’ TCTGGGTGAAACCACCCCTT 3’ and RM41 5’ CTGACAGGGTTAGTTTGCAT 3’ under the following conditions: 15 min at 95 °C followed by 35 cycles at 95 °C (30 s), 58 °C (30 s) and 72 °C (45 s), with a final extension at 60 °C (20 min). Sequences were determined by automated sequencing following the Big Dye Terminator Cycle Sequencing protocol (Perkin Elmer, Warrington, UK) using the referred primers.
Sequences obtained in this study were edited and aligned using CLUSTAL W [31] as implemented in BioEdit software [32] and the amino acid sequences were inferred using the software BioEdit [32]. The obtained sequences were also aligned and compared to European rabbit sequences available in GenBank. Codon numbering is according to the IMGT unique numbering for C-DOMAIN [33].
Sequence nucleotide diversity was estimated using DnaSP version 5.10 [34]. Haplotype reconstruction was obtained for each individual using PHASE version 2.1 [35] as implemented in DnaSP version 5.10 [34]. DnaSP was used to estimate population genetic parameters (S, π), Tajima’s D, number of haplotypes, haplotype diversity, variance of haplotype diversity and standard deviation of haplotype diversity.
Sequences obtained in this study were deposited in GenBank (Accession Numbers: MN650919-MN651079).

3. Results

3.1. Leporids IGHE and IGHM

We newly sequenced the IGHM and IGHE genes for six leporid genera: Bunolagus, Brachylagus, Lepus, Pentalagus, Romerolagus and Sylvilagus. The gene structure is similar to the European rabbit IGHM and IGHE and we thus inferred the intron-exon organization from published rabbit IGHM and IGHE sequences. Splicing signals are present at the intron boundaries and therefore we assumed that all studied leporids share the same IGHM and IGHE exon organization as the European rabbit.

3.1.1. IGHM

Cµ1 Domain

The Cµ1 domain is the most varied of leporid IGHM domains with 24 amino acid variable positions, the majority of which involve one substitution and 11 are conservative regarding the amino acid properties. For the Cµ1 domain some genera specific residues are observed: Arg 117 for the European rabbit, Asn 3 and Ala 122 for Pentalagus and Asn 115, Gln 117 and Asn 118 for Romerolagus (Figure 3).

Cµ2 Domain

The leporid Cµ2 domain is fairly conserved with only 12 amino acid variable positions. The great majority of these polymorphic positions involve one substitution but seven encode for changes in amino acid properties. The European rabbit has specific Ser residues at positions 100, while Lepus have specific Lys 45 residues and Sylvilagus have Ala 90 residues (Figure 3).

Cµ3 Domain

The variability observed for the leporid Cµ3 domain is similar to that of the Cµ2 domain. Of the 10 variable positions present in the CH3 domain all but one involve one substitution but seven encode for changes in amino acid properties. Interestingly, the European rabbit has two Pro residues at positions 1.3 and 1.1 and Met 86 residues whereas other leporids have 1.3 GluSer 1.1 and Thr 86. Lepus IgM specifically has Val 22 and Ser 83 (Figure 3).

Cµ4 domain

This is the most conserved of the leporid IGHM domains with 8 variable positions, half of which contain changes in amino acid properties. In this domain only one diagnostic position was observed for two genera: Lepus at Ile105 and Brachylagus at Leu 115 (Figure 3).

3.1.2. IGHE

Cε1 Domain

Similarly, to what was observed for IGHM, the Cε1 domain is the most diverse of leporid IGHE domains. In this domain 27 amino acid variable positions were observed. Of these, only four involve more than one substitution and six encode changes in amino acid properties. Specific residues were observed for Bunolagus, Ala 13, Lepus, Ser 43, and Romerolagus, Pro 1.3 and Phe 1.2 (Figure 4).

Cε2 Domain

In the leporid Cε2 domain 16 amino acid variable positions were observed. These involve one substitution but the majority encodes for changes in amino acid properties. Two diagnostic positions were observed for Romerolagus Gln15 and Phe21 (Figure 4).

Cε3 Domain

The leporid Cε3 domain has 14 amino acid variable positions, two of which involve more than one substitution and five encode for changes in amino acid properties. Two genera specific positions were observed for each of three genera: Lepus has specific Glu84 and Tyr102, Brachylagus has characteristic Asn45.4 and Asn84 and Pentalagus has as distinctive residues Thr85.1 and Thr116 (Figure 4).

Cε4 Domain

Contrarily to that observed for leporid IGHM, the Cε4 domain has higher variability than the Cε2 and Cε3 domains. The Cε4 domain has 19 amino acid variable positions, 13 of which involve changes in amino acid properties. Lepus has specific Glu 113 residues, Pentalagus has specific Cys 43 residues and Romerolagus has specific Glu13, Gln 45.2, Lys 45.5, and Gln 101 residues (Figure 4).

3.2. Wild and Domestic European Rabbits IgM and IgE Diversity

3.2.1. IGHM

European rabbits IGHM diversity seems to be higher for the O. c. cuniculus subspecies: of the 31 SNP’s we found, 20 were observed for the Spanish O. c. cuniculus population and 15 for the French O. c. cuniculus population, while 13 and 11 were found for the O. c. algirus subspecies and domestic breeds, respectively (Table 1). The nucleotide diversity agrees with this pattern (Table 1). Forty haplotypes were recovered (Table 1) most of which are rare, occurring once in heterozygosity, or uncommon (0.01–0.10 frequencies). The most common haplotype (0.178 frequency) is shared between O. c. algirus and the domestic breeds. Other haplotypes are shared between O. c. algirus and French wild O. c. cuniculus, O. c. cuniculus and domestic breeds and between French wild O. c. cuniculus and domestic breeds. The haplotype diversity values are in accordance with the nucleotide diversity, being highest for the French wild O. c. cuniculus and lowest for the domestic breeds (Table 1).
Less than one quarter of the recovered SNPs involves nonsynonymous substitutions resulting in 7 polymorphic amino acid positions (Table 2). Unique amino acids were found for O. c. algirus, G at CH1 residue 95, and Spanish wild O. c. cuniculus, I at CH3 residue 103, T and L at CH4 residues 82 and 114, respectively (Table 2).

3.2.2. IGHE

Our results show a greater diversity of IGHE in the IP populations: of the 26 SNPs observed, 16 were found for the O. c. algirus and 16 other for the Spanish wild O. c. cuniculus (Table 1). Surprisingly, for the Domestic breeds 9 SNPs were observed but only three SNPs were observed for the French O. c. cuniculus (Table 1). The nucleotide diversity values are in agreement with the observed SNPs (Table 1). We recovered 24 haplotypes, two of which are common: one haplotype is shared among all populations (0.322 frequency) and another haplotype is shared between O. c. cuniculus French populations and domestic breed’s (0.267 frequency); the remaining haplotypes are rare or uncommon (0.01–0.10 frequencies).
Approximately half of the recovered SNP’s involve nonsynonymous substitutions resulting in 12 polymorphic amino acid positions (Table 3). Unique amino acids were found for O. c. algirus: M and F at CH2 residues 81 and 116, respectively; Spanish wild O. c. cuniculus: I and M at CH1 residues 4 and 84, respectively, and K at CH4 residue 82; and for the domestic breeds: K, A and T at CH1 residues 12, 44 and 45, respectively, G at CH2 residue 1.7 and V at CH3 residue 85 (Table 3)

4. Discussion

The study of European rabbit immunoglobulin genes has contributed to elucidating numerous fundamental aspects of antibody structure and diversification mechanisms (reviewed in [36]). The rabbit has also been increasingly used as an animal model to many human diseases [2]. Despite this, the study of rabbit immunoglobulins has focused on the IgG and IgA isotypes. The knowledge on other leporids immune system genes is also very scarce.

4.1. IGHM and IGHE Evolution

The results obtained in this study show that throughout Leporid evolution IGHM and IGHE have been fairly conserved (each Ig has 98% sequence similarity among leporids), similarly to what has been found for IGHG [11]. A gradient in IGHM constant domains diversity was observed, with the Cµ1 domain being the most diverse and the Cµ4 the most conserved. This gradient in domain diversity was not observed for IGHE, the Cε1 is the most diverse domain but Cε4 shows greater diversity than Cε2 or Cε3 domains, the most conserved domains. For leporid IGHG a gradient in the constant domains diversity was described from the most conserved Cγ1 domain to the most diverse Cγ3 domain [11] and for the European rabbit IgA subclasses the Cα1 domain has twice more variable amino acid sites than Cα2 and Cα3 domains [6,7]. The differences in the pattern of most conserved and most diverse domain between the different Ig isotypes are most likely related to each isotype function. To trigger the effector functions that lead to pathogen elimination Igs must bind to other effector molecules such as Fc receptors and the complement. IgM interacts with three receptors, FcµR, Fcα/µR and pIgR, the binding sites of which have been mapped to the Cµ4 domain [37,38,39] and with the complement which binds to sites in the Cµ3 domain [40,41,42], the domains we found to be the most conserved during leporid IGHM evolution. IgE two receptors, FcεRI and FcεRII, binding sites locate to the Cε3 [43,44], again the domain we found to be the most conserved. IgA ligands, both Fc receptors, FcαRI and pIgR, and some pathogen evasion proteins interact with sites on the Cα2 and Cα3 domains (reviewed in [45], the most conserved IgA domains in rabbit IgA subclasses [6,7] and during mammalian evolution [45]. As for rabbit IgG, it’s receptor, FcγR, and the complement C1q protein have been shown to bind residues in the Cγ2 domain [46,47,48] a very conserved domain in leporid IgG evolution [11]. Thus, Ig domains which are likely essential for interaction with host ligands have been conserved during leporid evolution. Curiously, for the IgE Cε3 domain all but two of the 14 amino acid variable positions we found, including the five positions that encode for changes in amino acid properties, lie in or in direct vicinity of the binding sites to FcεRI and FcεRII. Variation at these codons may prove adaptive by improving the binding of IgE to its Fc receptors. It would be interesting to study the rabbit IgE Fc receptors which may also show specific differences reflecting the host adaptation to different selective pressures imposed by the different environments that each leporid species inhabit.

4.2. Wild and Domestic European Rabbits IgM and IgE Diversity

The IP is the region were the European rabbit originated [13]. Later, during the 7th century the rabbit was domesticated from French wild European rabbit populations [15,16]. Consistent with this, immunological markers genetic diversity has been shown to decrease from the IP to the French wild European rabbit populations and that the domestic animals are, to a great extent, a subset of the diversity found in the French wild O. c. cuniculus [29,30,36,49,50,51,52]. The results we obtained for IgE show a greater nucleotide and haplotype diversity in the IP O.c. algirus populations, in agreement with the previously described pattern of highest genetic diversity in IP populations. Interestingly, the genetic diversity we found for the domestic animals IgE is higher than that for the French wild O. c. cuniculus (domestic animals nucleotide diversity is twice as high as for French O. c. cuniculus and haplotype diversity is 1,5 times higher for domestic animals, Table 1), similar to what was found for European rabbit IgG [29]. This may reflect changes in allele frequencies caused by the founder effect on the origin of domestic breeds [15]. Striking for IgE is also the fact that despite the highest number of SNPs was found for the wild Iberian populations, the higher number of unique amino acids was found for the domestic animals, again similarly to what was found for IgG [29]. This indicates that the domestic breeds are under a higher selective pressure, possibly due to higher densities leading to exposure to different parasitic worms, leading to the emergence and maintenance of mutations of functional significance. In fact, despite the prophylaxis adopted in rabbitries, the prevalence of parasitic infections in domestic rabbits is high and mixed infestations predominate (approx. 80%; [53,54]) as an excessive stocking density of animals favors invasion mechanisms [54]. On the other hand, the more ancient Iberian European rabbit populations have accumulated a higher number of random mutations. We cannot discard that the higher variability in the domestic rabbits could be associated with a compensatory mechanism already proposed by van der Loo [55]: to compensate the lack of genetic diversity shown by the rabbit domestic breeds for most immunologic genes, some particular genes show a higher diversity when compared to the wild populations.
The results we obtained for IgM, on the other hand, revealed a higher nucleotide and haplotype diversity for the O. c. cuniculus subspecies, particularly for the French populations, and that the domestic animals have the lowest genetic diversity (Table 1). Considering the results previously obtained for most immunological markers [29,30,36,49,50,51,52] the higher genetic diversity observed for the French O. c. cuniculus was unexpected. Nevertheless, this higher nucleotide and haplotype diversity does not translate into higher number of polymorphic amino acid sites for the French O. c. cuniculus. Noteworthy for IgM is also the fact that only a small proportion of the observed SNP’s (7 out of 31 SNPs) involve non-synonymous substitutions, in spite of the observed 31 SNPs only seven polymorphic amino acid positions were detected. IgM is the most ancient Ig isotype being present with the same function in all gnathostomes [56]. As such it is a very conserved gene that most likely suffers stochastic nucleotide substitutions and is subject to strong purifying selection.

5. Conclusions

Our findings show that throughout Leporid evolution IGHM and IGHE have been fairly conserved. Despite this, variable amino acid residues were identified for IgE which may improve the binding of IgE to its Fc receptors. The observed pattern of genetic diversity of IGHM and IGHE along the wild rabbit populations from Iberia and France is contradictory, being, as expected, higher for the Iberian populations for IGHE but surprisingly higher for French populations for IGHM. This is the first time that these two immunoglobulin isotypes were fully characterized for Leporids and for European rabbit wild populations and domestic breeds. These findings will be useful for the management of rabbit populations and for the understanding of the rabbit adaptive immunity against pathogens.

Author Contributions

conceptualization, A.P. and P.J.E.; methodology, A.P. and T.A.; formal analysis, A.P. and T.A.; writing—original draft preparation, A.P.; writing—review and editing, T.A. and P.J.E.

Funding

The Portuguese Foundation for Science and Technology (FCT) supported the postdoctoral fellowship of Ana Pinheiro (SFRH/BPD/117451/2016), doctoral fellowship of Tereza Almeida (PD/BD/114542/2016) and the FCT Investigator grant of Pedro J. Esteves (IF/00376/2015). This article is a result of the project AGRIGEN – NORTE-01-0145-FEDER-000007, supported by Norte Portugal Regional Operational Programme (NORTE2020), under the PORTUGAL 2020 Partnership Agreement, through the European Regional Development Fund (ERDF).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Schroeder, H.W., Jr.; Cavacini, L. Structure and function of immunoglobulins. J. Allergy Clin. Immunol. 2010, 125, S41–S52. [Google Scholar] [CrossRef] [PubMed]
  2. Esteves, P.J.; Abrantes, J.; Baldauf, H.M.; BenMohamed, L.; Chen, Y.; Christensen, N.; Gonzalez-Gallego, J.; Giacani, L.; Hu, J.; Kaplan, G.; et al. The wide utility of rabbits as models of human diseases. Exp. Mol. Med. 2018, 50, 66. [Google Scholar] [CrossRef] [PubMed]
  3. Lanning, D.K.; Zhai, S.K.; Knight, K.L. Analysis of the 3′ Cmu region of the rabbit Ig heavy chain locus. Gene 2003, 309, 135–144. [Google Scholar] [CrossRef]
  4. Ros, F.; Puels, J.; Reichenberger, N.; van Schooten, W.; Buelow, R.; Platzer, J. Sequence analysis of 0.5 Mb of the rabbit germline immunoglobulin heavy chain locus. Gene 2004, 330, 49–59. [Google Scholar] [CrossRef] [PubMed]
  5. Gertz, E.M.; Schaffer, A.A.; Agarwala, R.; Bonnet-Garnier, A.; Rogel-Gaillard, C.; Hayes, H.; Mage, R.G. Accuracy and coverage assessment of Oryctolagus cuniculus (rabbit) genes encoding immunoglobulins in the whole genome sequence assembly (OryCun2.0) and localization of the IGH locus to chromosome 20. Immunogenetics 2013, 65, 749–762. [Google Scholar] [CrossRef] [PubMed]
  6. Burnett, R.C.; Hanly, W.C.; Zhai, S.K.; Knight, K.L. The IgA heavy-chain gene family in rabbit: Cloning and sequence analysis of 13 C alpha genes. EMBO J. 1989, 8, 4041–4047. [Google Scholar] [CrossRef] [PubMed]
  7. Pinheiro, A.; de Sousa-Pereira, P.; Strive, T.; Knight, K.L.; Woof, J.M.; Esteves, P.J.; Abrantes, J. Identification of a new European rabbit IgA with a serine-rich hinge region. PLoS ONE 2018, 13, e0201567. [Google Scholar] [CrossRef] [PubMed]
  8. Esteves, P.J.; Alves, P.C.; Ferrand, N.; van der Loo, W. Hotspot variation at the CH2-CH3 interface of leporid IgG antibodies (Oryctolagus, Sylvilagus and Lepus). Eur. J. Immunogenet. 2002, 29, 529–535. [Google Scholar] [CrossRef] [PubMed]
  9. Girardi, E.; Holdom, M.D.; Davies, A.M.; Sutton, B.J.; Beavil, A.J. The crystal structure of rabbit IgG-Fc. Biochem. J. 2009, 417, 77–83. [Google Scholar] [CrossRef] [PubMed]
  10. Hamers, R.; Hamers-Casterman, C. Molecular localization of a chain allotypic specificities in rabbit IgG (7S gamma-globulin). J. Mol. Biol. 1965, 14, 288–289. [Google Scholar] [CrossRef]
  11. Pinheiro, A.; Woof, J.M.; Almeida, T.; Abrantes, J.; Alves, P.C.; Gortazar, C.; Esteves, P.J. Leporid immunoglobulin G shows evidence of strong selective pressure on the hinge and CH3 domains. Open Biol. 2014, 4, 140088. [Google Scholar] [CrossRef] [PubMed]
  12. Owen, J.A.; Punt, J.; Stranford, S.A.; Jones, P.P.; Kuby, J. Kuby Immunology, 7th ed.; W.H. Freeman: New York, NY, USA, 2013; p. 692. [Google Scholar]
  13. Lopez-Martinez, N. Revisión Sistemática y Biostratigráfica de los Lagomorpha (Mammalia) del Terciario y Cuatemario de España; Diputación General de Aragón: Zaragoza, Spain, 1989; Volume 3. [Google Scholar]
  14. Queney, G.; Ferrand, N.; Weiss, S.; Mougel, F.; Monnerot, M. Stationary Distributions of Microsatellite Loci between Divergent Population Groups of the European Rabbit (Oryctolagus cuniculus). Mol. Biol. Evol. 2001, 18, 2169–2178. [Google Scholar] [CrossRef] [PubMed]
  15. Carneiro, M.; Afonso, S.; Geraldes, A.; Garreau, H.; Bolet, G.; Boucher, S.; Tircazes, A.; Queney, G.; Nachman, M.W.; Ferrand, N. The genetic structure of domestic rabbits. Mol. Biol. Evol. 2011, 28, 1801–1816. [Google Scholar] [CrossRef] [PubMed]
  16. Carneiro, M.; Albert, F.W.; Afonso, S.; Pereira, R.J.; Burbano, H.; Campos, R.; Melo-Ferreira, J.; Blanco-Aguiar, J.A.; Villafuerte, R.; Nachman, M.W.; et al. The genomic architecture of population divergence between subspecies of the European rabbit. PLoS Genet 2014, 10, e1003519. [Google Scholar] [CrossRef] [PubMed]
  17. Alves, P.C.; Hackländer, K. Lagomorph species: Geographical distribution and conservation status. In Lagomorph Biology; Alves, P.C., Ferrand, N., Hackländer, K., Eds.; Springer: Berlin/Heidelberg, Germany, 2008; pp. 395–405. [Google Scholar]
  18. Chapman, J.; Flux, J. Introduction and overview of the lagomorphs. In Rabbits, Hares and Pikas, Status, Survey and Conservation Action Plan; Chapman, J., Flux, J., Eds.; UICN/SSC Lagomorph Specialist Group: Gland, Switzerland, 1990; pp. 1–6. [Google Scholar]
  19. Corbet, G.B.; Hill, J.E. A World List of Mammalian Species; Oxford University Press: Oxford, UK, 1980. [Google Scholar]
  20. Can, D.N.; Abramov, A.V.; Tikhonov, A.N.; Averianov, A.O. Annamite striped rabbit Nesolagus timminsi in Vietnam. Acta Theriol. 2001, 46, 437–440. [Google Scholar] [CrossRef]
  21. Chapman, J.A.; Cramer, K.L.; Dippenaar, N.J.; Robinson, T.J. Systematics and Biogeography of the New-England Cottontail, Sylvilagus transitionalis (Bangs, 1895), with the Description of a New Species from the Appalachian Mountains. Proc. Biol. Soc. Wash. 1992, 105, 841–866. [Google Scholar]
  22. Flux, J.E.C.; Angermann, R. The hares and jackrabbits. In Rabbits, Hares and Pikas: Status Conservation Action Plan; Chapman, J.A., Flux, J.E.C., Eds.; International Union for Conservation of Nature and Natural Resources: Gland, Switzerland, 1990; pp. 61–94. [Google Scholar]
  23. Frey, J.K.; Fisher, R.D.; Ruedas, L.A. Identification and restriction of the type locality of the Manzano Mountains cottontail, Sylvilagus cognatus Nelson, 1907 (Mammalia: Lagomorpha: Leporidae). Proc. Biol. Soc. Wash. 1997, 110, 329–331. [Google Scholar]
  24. Matthee, C.A.; Robinson, T.J. Mitochondrial DNA differentiation among geographical populations of Pronolagus rupestris, Smith’s red rock rabbit (Mammalia: Lagomorpha). Heredity 1996, 76 Pt 5, 514–523. [Google Scholar] [CrossRef]
  25. Surridge, A.K.; Timmins, R.J.; Hewitt, G.M.; Bell, D.J. Striped rabbits in Southeast Asia. Nature 1999, 400, 726. [Google Scholar] [CrossRef]
  26. Matthee, C.A.; van Vuuren, B.J.; Bell, D.; Robinson, T.J. A molecular supermatrix of the rabbits and hares (Leporidae) allows for the identification of five intercontinental exchanges during the Miocene. Syst. Biol. 2004, 53, 433–447. [Google Scholar] [CrossRef] [PubMed]
  27. Sparwel, M.; Doronina, L.; Churakov, G.; Stegemann, A.; Brosius, J.; Robinson, T.J.; Schmitz, J. The Volcano Rabbit in the Phylogenetic Network of Lagomorphs. Genome Biol. Evol. 2019, 11, 11–16. [Google Scholar] [CrossRef] [PubMed]
  28. Robinson, T.J.; Matthee, C.A. Phylogeny and evolutionary origins of the Leporidae: A review of cytogenetics, molecular analyses and a supermatrix analysis. Mammal Rev. 2005, 35, 231–247. [Google Scholar] [CrossRef]
  29. Pinheiro, A.; Almeida, T.; Esteves, P.J. Survey of genetic diversity of IgG in wild and domestic rabbits. Int. J. Immunogenet. 2015, 42, 364–367. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Abrantes, J.; Areal, H.; Esteves, P.J. Insights into the European rabbit (Oryctolagus cuniculus) innate immune system: Genetic diversity of the toll-like receptor 3 (TLR3) in wild populations and domestic breeds. BMC Genet. 2013, 14, 73. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Thompson, J.D.; Higgins, D.G.; Gibson, T.J. 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, 4673–4680. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Hall, T.A. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp. Ser. 1999, 41, 3. [Google Scholar]
  33. Lefranc, M.P.; Pommie, C.; Kaas, Q.; Duprat, E.; Bosc, N.; Guiraudou, D.; Jean, C.; Ruiz, M.; Da Piedade, I.; Rouard, M.; et al. IMGT unique numbering for immunoglobulin and T cell receptor constant domains and Ig superfamily C-like domains. Dev. Comp. Immunol. 2005, 29, 185–203. [Google Scholar] [CrossRef] [PubMed]
  34. Librado, P.; Rozas, J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics 2009, 25, 1451–1452. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Stephens, M.; Donnelly, P. A comparison of bayesian methods for haplotype reconstruction from population genotype data. Am. J. Hum. Genet. 2003, 73, 1162–1169. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Pinheiro, A.; Neves, F.; Lemos de Matos, A.; Abrantes, J.; van der Loo, W.; Mage, R.; Esteves, P.J. An overview of the lagomorph immune system and its genetic diversity. Immunogenetics 2016, 68, 83–107. [Google Scholar] [CrossRef] [PubMed]
  37. Ghumra, A.; Shi, J.; McIntosh, R.S.; Rasmussen, I.B.; Braathen, R.; Johansen, F.E.; Sandlie, I.; Mongini, P.K.; Areschoug, T.; Lindahl, G.; et al. Structural requirements for the interaction of human IgM and IgA with the human Fcalpha/mu receptor. Eur. J. Immunol. 2009, 39, 1147–1156. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Braathen, R.; Sorensen, V.; Brandtzaeg, P.; Sandlie, I.; Johansen, F.E. The carboxyl-terminal domains of IgA and IgM direct isotype-specific polymerization and interaction with the polymeric immunoglobulin receptor. J. Biol. Chem. 2002, 277, 42755–42762. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Lloyd, K.A.; Wang, J.; Urban, B.C.; Czajkowsky, D.M.; Pleass, R.J. Glycan-independent binding and internalization of human IgM to FCMR, its cognate cellular receptor. Sci. Rep. 2017, 7, 42989. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Sharp, T.H.; Boyle, A.L.; Diebolder, C.A.; Kros, A.; Koster, A.J.; Gros, P. Insights into IgM-mediated complement activation based on in situ structures of IgM-C1-C4b. Proc. Natl. Acad. Sci. USA 2019, 116, 11900–11905. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Arya, S.; Chen, F.; Spycher, S.; Isenman, D.E.; Shulman, M.J.; Painter, R.H. Mapping of amino acid residues in the C mu 3 domain of mouse IgM important in macromolecular assembly and complement-dependent cytolysis. J. Immunol. 1994, 152, 1206–1212. [Google Scholar] [PubMed]
  42. Taylor, B.; Wright, J.F.; Arya, S.; Isenman, D.E.; Shulman, M.J.; Painter, R.H. C1q binding properties of monomer and polymer forms of mouse IgM mu-chain variants. Pro544Gly and Pro434Ala. J. Immunol. 1994, 153, 5303–5313. [Google Scholar] [PubMed]
  43. Holdom, M.D.; Davies, A.M.; Nettleship, J.E.; Bagby, S.C.; Dhaliwal, B.; Girardi, E.; Hunt, J.; Gould, H.J.; Beavil, A.J.; McDonnell, J.M.; et al. Conformational changes in IgE contribute to its uniquely slow dissociation rate from receptor FcepsilonRI. Nat. Struct. Mol. Biol. 2011, 18, 571–576. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Garman, S.C.; Wurzburg, B.A.; Tarchevskaya, S.S.; Kinet, J.P.; Jardetzky, T.S. Structure of the Fc fragment of human IgE bound to its high-affinity receptor FcepsilonRI alpha. Nature 2000, 406, 259–266. [Google Scholar] [CrossRef] [PubMed]
  45. Pinheiro, A.; Woof, J.M.; Abi-Rached, L.; Parham, P.; Esteves, P.J. Computational analyses of an evolutionary arms race between mammalian immunity mediated by immunoglobulin A and its subversion by bacterial pathogens. PLoS ONE 2013, 8, e73934. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Rayner, L.E.; Kadkhodayi-Kholghi, N.; Heenan, R.K.; Gor, J.; Dalby, P.A.; Perkins, S.J. The solution structure of rabbit IgG accounts for its interactions with the Fc receptor and complement C1q and its conformational stability. J. Mol. Biol. 2013, 425, 506–523. [Google Scholar] [CrossRef] [PubMed]
  47. Sondermann, P.; Huber, R.; Oosthuizen, V.; Jacob, U. The 3.2-A crystal structure of the human IgG1 Fc fragment-Fc gammaRIII complex. Nature 2000, 406, 267–273. [Google Scholar] [CrossRef] [PubMed]
  48. Schneider, S.; Zacharias, M. Atomic resolution model of the antibody Fc interaction with the complement C1q component. Mol. Immunol. 2012, 51, 66–72. [Google Scholar] [CrossRef] [PubMed]
  49. van der Loo, W.; Mougel, F.; Sanchez, M.S.; Bouton, C.; Castien, E.; Fonseca, A.; Ferrand, N.; Soriguer, R.; Monnerot, M. Cytonuclear disequilibria in wild populations of rabbit (Oryctolagus cuniculus L.) suggest unequal allele turnover rates at the b locus (IGKC1). Immunogenetics 1999, 49, 629–643. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Esteves, P.J.; Lanning, D.; Ferrand, N.; Knight, K.L.; Zhai, S.K.; van der Loo, W. Allelic variation at the VHa locus in natural populations of rabbit (Oryctolagus cuniculus, L.). J. Immunol. 2004, 172, 1044–1053. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Surridge, A.K.; van der Loo, W.; Abrantes, J.; Carneiro, M.; Hewitt, G.M.; Esteves, P.J. Diversity and evolutionary history of the MHC DQA gene in leporids. Immunogenetics 2008, 60, 515–525. [Google Scholar] [CrossRef] [PubMed]
  52. Pinheiro, A.; Lanning, D.; Alves, P.C.; Mage, R.G.; Knight, K.L.; van der Loo, W.; Esteves, P.J. Molecular bases of genetic diversity and evolution of the immunoglobulin heavy chain variable region (IGHV) gene locus in leporids. Immunogenetics 2011, 63, 397–408. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Ilic, T.; Stepanovic, P.; Nenadovic, K.; Dimitrijevic, S. Improving agricultural production of domestic rabbits in Serbia by follow-up study of their parasitic infections. Iran. J. Vet. Res. 2018, 19, 290–297. [Google Scholar] [PubMed]
  54. Szkucik, K.; Pyz-Lukasik, R.; Szczepaniak, K.O.; Paszkiewicz, W. Occurrence of gastrointestinal parasites in slaughter rabbits. Parasitol. Res. 2014, 113, 59–64. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. van der Loo, W. Variance analysis of immunoglobulin alleles in natural populations of rabbit (Oryctolagus cuniculus): The extensive interallelic divergence at the b locus could be the outcome of overdominance-type selection. Genetics 1993, 135, 171–187. [Google Scholar] [PubMed]
  56. Flajnik, M.F.; Kasahara, M. Origin and evolution of the adaptive immune system: Genetic events and selective pressures. Nat. Rev. Genet. 2010, 11, 47–59. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. European rabbit (Oryctolagus cuniculus) IgM and IgE DNA sequences of constant regions (Genbank accession numbers for IGHM: J00666, and IGHE: AY386696.1). The amino acid translation is shown. The IMGT unique numbering for the constant domain (13) is shown on top.
Figure 1. European rabbit (Oryctolagus cuniculus) IgM and IgE DNA sequences of constant regions (Genbank accession numbers for IGHM: J00666, and IGHE: AY386696.1). The amino acid translation is shown. The IMGT unique numbering for the constant domain (13) is shown on top.
Animals 09 00955 g001
Figure 2. Map with the geographical locations of the populations used in this study. Marked with (●) are O. c. algirus populations from Portugal, with (■) are O. c. cuniculus from Spain and with (▴) O. c. cuniculus from France.
Figure 2. Map with the geographical locations of the populations used in this study. Marked with (●) are O. c. algirus populations from Portugal, with (■) are O. c. cuniculus from Spain and with (▴) O. c. cuniculus from France.
Animals 09 00955 g002
Figure 3. Variable codons and amino acid physicochemical properties for leporid IGHM constant domains. The codon numbering is according to the IMGT unique numbering for the constant domain (13). The colours represent amino acid properties: polar neutral (green), polar positive (yellow), polar negative (orange), non-polar neutral (purple), non-polar aliphatic (blue) and non-polar aromatic (light pink). Occ—Oryctolagus cuniculus cuniculus (accession number J00666); Lcalif—Lepus californicus; Lcal—Lepus callotis; Lcas—Lepus castroviejoi; Lcap—Lepus capensis; Leur—Lepus europaeus; Lgra—Lepus granatensis; Lsax—Lepus saxatilis; Ltwn—Lepus townsendi; Lam—Lepus americanus; Lcap—Lepus capensis; Ltim—Lepus timidus; Pfur—Pentalagus furnessii; Sflo—Sylvilagus floridanus; Sbac—Sylvilagus bachmanii; Rdia—Romerolagus diazii; Bida—Brachylagus idahoensis; Bmon—Bunolagus monticularis.
Figure 3. Variable codons and amino acid physicochemical properties for leporid IGHM constant domains. The codon numbering is according to the IMGT unique numbering for the constant domain (13). The colours represent amino acid properties: polar neutral (green), polar positive (yellow), polar negative (orange), non-polar neutral (purple), non-polar aliphatic (blue) and non-polar aromatic (light pink). Occ—Oryctolagus cuniculus cuniculus (accession number J00666); Lcalif—Lepus californicus; Lcal—Lepus callotis; Lcas—Lepus castroviejoi; Lcap—Lepus capensis; Leur—Lepus europaeus; Lgra—Lepus granatensis; Lsax—Lepus saxatilis; Ltwn—Lepus townsendi; Lam—Lepus americanus; Lcap—Lepus capensis; Ltim—Lepus timidus; Pfur—Pentalagus furnessii; Sflo—Sylvilagus floridanus; Sbac—Sylvilagus bachmanii; Rdia—Romerolagus diazii; Bida—Brachylagus idahoensis; Bmon—Bunolagus monticularis.
Animals 09 00955 g003
Figure 4. Variable codons and amino acid physicochemical properties for leporid IGHE constant domains. The codon numbering is according to the IMGT unique numbering for the constant domain (13). The colours represent amino acid properties: polar neutral (green), polar positive (yellow), polar negative (orange), non-polar neutral (purple), non-polar aliphatic (blue) and non-polar aromatic (light pink). Occ—Oryctolagus cuniculus cuniculus (accession number AY386696.1); Lcalif—Lepus californicus; Lcall—Lepus callotis; Lcas—Lepus castroviejoi; Lcap—Lepus capensis; Leur—Lepus europaeus; Lgra—Lepus granatensis; Lsax—Lepus saxatilis; Ltwn—Lepus townsendi; Lam—Lepus americanus; Lcap—Lepus capensis; Ltim—Lepus timidus; Pfur—Pentalagus furnessii; Sflo—Sylvilagus floridanus; Sbac—Sylvilagus bachmanii; Rdia—Romerolagus diazii; Bida—Brachylagus idahoensis; Bmon—Bunolagus monticularis.
Figure 4. Variable codons and amino acid physicochemical properties for leporid IGHE constant domains. The codon numbering is according to the IMGT unique numbering for the constant domain (13). The colours represent amino acid properties: polar neutral (green), polar positive (yellow), polar negative (orange), non-polar neutral (purple), non-polar aliphatic (blue) and non-polar aromatic (light pink). Occ—Oryctolagus cuniculus cuniculus (accession number AY386696.1); Lcalif—Lepus californicus; Lcall—Lepus callotis; Lcas—Lepus castroviejoi; Lcap—Lepus capensis; Leur—Lepus europaeus; Lgra—Lepus granatensis; Lsax—Lepus saxatilis; Ltwn—Lepus townsendi; Lam—Lepus americanus; Lcap—Lepus capensis; Ltim—Lepus timidus; Pfur—Pentalagus furnessii; Sflo—Sylvilagus floridanus; Sbac—Sylvilagus bachmanii; Rdia—Romerolagus diazii; Bida—Brachylagus idahoensis; Bmon—Bunolagus monticularis.
Animals 09 00955 g004
Table 1. IGHM and IGHE nucleotide sequence polymorphism and estimates of haplotype diversity, haplotype diversity variance and haplotype standard deviation in the European rabbit populations/domestic animals studied.
Table 1. IGHM and IGHE nucleotide sequence polymorphism and estimates of haplotype diversity, haplotype diversity variance and haplotype standard deviation in the European rabbit populations/domestic animals studied.
LocusSubspeciesPopulationNSπDHHdHd VarianceHd Std Deviation
IGHM All45310.00317−0.98878400.9490.000170.013
O. c. algirusPortugal10130.002830.10630110.9110.001760.042
O. c. cuniculusSpain10200.00340−0.84715130.9110.002890.054
France10150.003450.39814130.9580.000650.025
Domestic15110.002430.5946580.7470.003630.060
IGHE All45260.00226−1.38819240.8170.000780.028
O. c. algirusPortugal10160.00348−0.02986110.8680.00410.064
O. c. cuniculusSpain10160.00222−1.3744980.7740.00510.071
France1030.000680.0886830.5420.010940.105
Domestic1590.00136−0.7176980.7430.0030.055
N, sample size; S, segregating sites; π, nucleotide diversity; D, Tajima’s D statistical test; H, number of haplotypes; Hd, haplotype diversity; Hd variance, variance of haplotype diversity; Hd std deviation, standard deviation of haplotype diversity.
Table 2. IgM amino acid polymorphic positions for European rabbit populations/domestic animals. The nucleotide variation is shown, dots represent identity. The codon numbering is according to the IMGT unique numbering for the constant domain [33].
Table 2. IgM amino acid polymorphic positions for European rabbit populations/domestic animals. The nucleotide variation is shown, dots represent identity. The codon numbering is according to the IMGT unique numbering for the constant domain [33].
SubspeciesPopulationDomain, Amino Acid Position and Most Common Triplet and Amino Acid
Cµ1Cµ3Cµ4
26377810910382114
cggRgtcVgtcVagcSacgTgccAcacH
O. c. algirusPortugalRVVga.D…/..aTAH
.a.Qa..I S
g..G
O. c. cuniculusSpainRVVga.D…/..aTAH
.a.Qa..I S.taIa..T.t.L
FranceRVVga.D…/..aTAH
a..Ia..IS
DomesticRVVga.D…/..aTA
a..Ia..IS
Table 3. IgE amino acid polymorphic positions for European rabbit populations/domestic animals. The nucleotide variation is shown, dots represent identity. The codon numbering is according to the IMGT unique numbering for the constant domain [33].
Table 3. IgE amino acid polymorphic positions for European rabbit populations/domestic animals. The nucleotide variation is shown, dots represent identity. The codon numbering is according to the IMGT unique numbering for the constant domain [33].
SubspeciesPopulationDomain, Amino Acid Position and Most Common Triplet and Amino Acid
Cε1Cε2Cε3Cε4
4124445.284.41.7811168485.4182
gtcVagaRggcGaacNacgTgcgAacgTctcLtatYgccAaccTgagE
O. c. algirusPortugalVRGNTAT LYATE
.t.Mt..Fg..D g..A
O. c. cuniculusSpainVRGNTATLg..DA….TE
a..I .t.M Y g..Aa..K
FranceVRGNTATLg..DATE
g..A
DomesticVRGNTATLg..D ATE
.a.K.c.A.c.T .g.G .t.Vg..A

Share and Cite

MDPI and ACS Style

Pinheiro, A.; Almeida, T.; Esteves, P.J. Genetic Diversity of IGHM and IGHE in the Leporids Revealed Different Patterns of Diversity in the Two European Rabbit Subspecies (O. cuniculus algirus and O. c. cuniculus). Animals 2019, 9, 955. https://doi.org/10.3390/ani9110955

AMA Style

Pinheiro A, Almeida T, Esteves PJ. Genetic Diversity of IGHM and IGHE in the Leporids Revealed Different Patterns of Diversity in the Two European Rabbit Subspecies (O. cuniculus algirus and O. c. cuniculus). Animals. 2019; 9(11):955. https://doi.org/10.3390/ani9110955

Chicago/Turabian Style

Pinheiro, Ana, Tereza Almeida, and Pedro J. Esteves. 2019. "Genetic Diversity of IGHM and IGHE in the Leporids Revealed Different Patterns of Diversity in the Two European Rabbit Subspecies (O. cuniculus algirus and O. c. cuniculus)" Animals 9, no. 11: 955. https://doi.org/10.3390/ani9110955

APA Style

Pinheiro, A., Almeida, T., & Esteves, P. J. (2019). Genetic Diversity of IGHM and IGHE in the Leporids Revealed Different Patterns of Diversity in the Two European Rabbit Subspecies (O. cuniculus algirus and O. c. cuniculus). Animals, 9(11), 955. https://doi.org/10.3390/ani9110955

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop