Major Histocompatibility Complex as Marker Assisted Selection for Breeding Immunocompetent Animal

Introduction

During the past half-century, commercial breeding programs based on productive traits have improved animal growth rate and production. Still, high selection intensity for quantitative traits caused many physiological disorders as well as a reduction in overall immunocompetence. For breeding resistance animals, gathering knowledge regarding the population's genetic makeup is a fundamental tool that helps identify variations in immunological traits. On the other hand, knowledge of the genetic diversity at immune response genes is necessary for understanding the host-pathogen interactions, disease resistance, and susceptibility. In this regard, the combination of traditional genetic selection and molecular marker-assisted selection (MAS) based on the immunogenetics criteria is expected to be a more effective breeding approach (White & Knowles, 2013).

The major histocompatibility complex (MHC) is a cluster of genes that plays a central role in the immune system to distinguish self from non-self. It was first discovered in mice as the genetic locus responsible for rapid tissue allograft rejection and encoding some highly polymorphic alloantigens on the surface of cells. Graft rejection is an accidental consequence of the basic function of the polymorphic alloantigens, the MHC classical class I and II molecules, whose role is to bind and present peptide fragments to T lymphocytes of the immune system (Rock et al., 2016; Wieczorek et al., 2017; Wosen et al., 2018; Nakamura et al., 2019). It is now known that MHC also has important roles in autoimmunity and reproductive success. The extreme polymorphism and somatic variations in MHC genes enable the host to recognize many foreign peptides and direct immune reactions toward humoral and/or cell-mediated response (Zakharova et al., 2019). T cell epitope mapping could be applied in vaccine design and formulations for all species where functional MHC molecules have been analyzed. In veterinary science, the important role of MHC in disease resistance and production traits has made it a precious marker in selection programs (Behl et al., 2012; Q. H. Han et al., 2019; Manlik et al., 2019; Jaworska et al., 2020).

Several important characteristics of MHC genes that make them exceptional candidates for investigation are mentioned in Table 1.

Table 1. Substantial properties of MHC genes that make them precious targets of research

The MHC region, also known as leukocyte antigen (LA), comprises three clusters of genes, including classes I, II, and III. Class I and II genes express cell surface glycoproteins and are functionally divided into two types of classical and nonclassical MHC. Classical MHC molecules have polymorphic nature and present peptide ligands to T cell receptors (TcR). In contrast, nonclassical MHC molecules mediate non-restricted T cell functions such as inhibitory or activating stimuli in natural killer cells, peptide loading of class II molecules in the endosomal/lysosomal system, and displaying lipids. MHC class III encodes for other immune components like complement components (C2, C4, factor B) and cytokines (e.g., TNF-α) (Al Naqbi et al., 2021; Matzaraki et al., 2017; Radwan et al., 2020; Tumer et al., 2019).

MHC class I region includes three different loci, A/B/C, that encodes transmembrane α glycoproteins chains containing α1, α2, and mostly α3 domain, which is dimerized with soluble β-2-microglobulin. Classical MHC I molecules present endogenous antigens (e.g., viral and tumor peptides) to CD8+ T cells and mediate cellular immunity. MHC class II region includes four different loci, DP/DQ/DR/DW, that encode alpha (α) and beta (β) chains on the surface of antigen-presenting cells (APCs). Classical MHC II molecules present exogenous antigens (e.g., bacterial peptides) to CD4+ T cells and induce effective cytokines for humoral/cell-mediated immune responses (Figure 1) (Cruz-Tapias et al., 2013).

Figure 1. The three-dimensional (3D) structures and gene loci of human major histocompatibility complex (MHC). A) 3D structure of MHC class I, B) 3D structure of MHC class II, C) gene loci of MHC (Y. Li et al., 2010; McMahon et al., 2011; Sehnal et al., 2021)

The existence of the MHC system in different species was initially found as leukocyte antigens. Lymphocyte immunizations and the generation of monospecific antibodies against lymphocytes were commonly used for identifying the MHC haplotypes. Serological typing has several limitations, including the presence of highly similar MHC epitopes, the existence of novel haplotypes in a different population, variability of non-MHC antigens, and subjectivity in the interpretation of serological reactions (Janet E. Fulton et al., 2006).

Several molecular methods are now available to analyze MHC alleles and genotypes in domestic animals. Some of the PCR-based methods that are available for specified loci are sequence-specific primers (SSP), sequence-specific oligonucleotide probe analysis (SSOP), and single-strand conformational polymorphism (SSCP), restriction fragment length polymorphism (RFLP) analysis with identification by Southern blot analysis (Nikbakht et al., 2009). Direct DNA sequencing or cloning and sequencing the PCR products are the gold standard for most phylogenetic studies. However, the cost required for such analysis is relatively high and time-consuming. The MHC haplotypes will be determined using a microsatellite marker located within the MHC region. This approach has successfully been used for cattle and chicken (Abdurakhmonov, 2016; Janet E. Fulton et al., 2016; Iglesias et al., 2019; Manjula et al., 2021; Gholamreza Nikbakht et al., 2013).

The polymorphic nature of the classical MHC molecules is reflected by the high number of alleles detected for each locus. In most vertebrate species, a high degree of genetic polymorphism has been reported in the MHC II region. In this region, the DR locus is highly polymorphic and probably plays an essential role in developing MHC-restricted immune responses. The MHC class II alleles have also been associated with many immune and non-immune responses in humans and animals. They are considered a genetic risk factor for many autoimmune diseases. This review summarizes our findings over two decades of intense research on MHC polymorphism and its associations with immune responses in domestic animals.

Ruminants

Cattle

In cattle, the MHC system, called bovine leukocyte antigens (BoLA), was initially found by serological methods via lymphocyte immunizations and by the generation of monospecific antibodies. BoLA is located on chromosome 23 and contains about 154 predicted functional genes (Behl et al., 2012). Class I genes are composed of two BoLA‐A and BoLA‐B regions. Only the BoLA‐A locus seems to be functional, with nearly 100 known classical alleles listed on the Immuno Polymorphism Database (IPD) (https://www.ebi.ac.uk/ipd/mhc/group/BoLA/) (Figure 2).

Figure 2. Alleles of different regions of bovine leukocyte antigen (BoLA)

Class II genes are located in two separated subregions, IIa and IIb. The IIa comprises two gene clusters, DR and DQ (not DP gene), that express one DR molecule and one or two DQ molecules. Cattle have one monomorphic coding gene for α chain of the DR (DRA) and three genes that encode for the β chain of the DR (DRB 1-3) molecule. DRB1 is a pseudogene, DRB2 is poorly expressed, but DRB3 is strongly expressed. DRB3 is the most polymorphic locus for which, to date, 384 alleles have been identified (Figure 2). The DQ genes are composed of DQA (DQA 1–5) and DQB (DQB 1–5) genes, with at least 76 BoLA DQA and 91 BoLA DQB alleles. DR and DQ molecules are mainly involved in presenting a broad spectrum of antigens in cattle. The BoLA class IIb locus contains genes involved in antigen processing and transport (TAP) (Giovambattista et al., 2020; Kumari et al., 2021).

There is limited knowledge on variability, structure, interspecies phylogeny, and phylogenetic characteristics of cattle MHC. We have performed several projects on exon2 of BoLA-DRB3 (BoLA-DRB3.2) allelic polymorphism among Holstein cattle and river buffaloes. BoLA-DRB3.2 allelic frequencies detected in Iranian cattle are summarized in Table 2 (Alkafajy et al., 2020; Nikbakht Brujeni et al., 2016; Ranjbar et al., 2017).

The genetic diversity at BoLA-DRB3.2 was associated with the resistance/susceptibility of cattle to some infectious diseases (Table 3). In cattle, enzootic bovine leukosis (EBV) is one of the most important examples of the relationship between the MHC and infectious diseases. Nucleotide polymorphisms in the BoLA-DRB3.2 are associated with bovine leukemia virus (BLV) proviral load and different patterns of EBV, such as persistent lymphocytosis (PL) and leukemia or lymphosarcoma (LS). Specific amino acid sequences within the highly polymorphic peptide binding region of the BoLA-DRB3.2 influence the immune system's binding, orientation, and presentation. Due to the widespread of BLV, eradication programs based on the detection of BLV seropositive cattle and culling the infected animals are not feasible. In this case, breeding for disease resistance can be practiced while offering improved animal welfare and increased returns for breeders (Nikbakht Brujeni et al., 2016).

Even before immune molecules (especially MHC) had been structurally characterized by experimental and in-silico methods, it was hypothesized a long time ago that areas with greater amino acid sequence variability could be considered surrogate indicators of antigen interaction points in the receptor. This general correlation has been verified (Wu & Kabat, 2008). We further analyzed the associations between MHC alleles and calf diarrhea. Calf diarrhea was significantly associated with the BoLA-DRB3.2 genotype *1104/*0101*. A potential relationship between amino acid sequences in the antigen-binding groove of BoLA-DRB3 alleles and susceptibility or resistance to calf diarrhea is demonstrated. Among twelve different DRB3 alleles, 26 polymorphic positions were detected. A significant association was identified between diarrhea and the presence of Glutamic acid and Tyrosine in pocket 4 and Valine, Glutamine, and Leucine in pocket 9 of the peptide binding region. Thus, it can be concluded that pockets 4 and 9 of the BoLA-DRB3 molecule are involved in conferring susceptibility to calf diarrhea (Nikbakht et al., 2011).

Table 2. BoLA-DRB3.2 allelic frequencies detected in Iranian cattle

Table 3. BoLA-DRB3.2 alleles associated with resistance (R)/susceptibility (S) of cattle to infectious disease

Buffalo

The buffalo (Bubalus Bubalis) MHC genomic regions are named lymphocyte antigen (BuLA/Bubo), which have a similar overall organization to cattle (Tizard, 2018). Similarly, MHC class II genes of buffalo are located in two distinct locations but on chromosome 2p (C. Li et al., 2020). MHC IIa contains DR and DQ, and IIb contains DI, DY, and DO (Medina et al., 2019; Othman et al., 2018). The second exons of DRB, DQA, and DQB, coding peptide binding sites (PBS) of class II molecules, are highly polymorphic. The MHC class II DRB gene in buffalo is also orthologues to the cattle DRB3 (Behl et al., 2012). In buffalo, the polymorphism of DRB3 has been defined by using several PCR-based methods, including single-strand conformational polymorphism (SSCP) in Indian river buffalo, heteroduplex analysis (HA), polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP), and sequencing (De et al., 2002).

Polymorphism in DRB3 exon 2 of Iranian river buffalo (Bubalus bubalis) has been analyzed. The second exon of Bubu-DRB3 was amplified by polymerase chain reaction (PCR), and polymorphisms were detected by restriction fragment length polymorphism (RFLP) and heteroduplex analysis. A series of 8 RsaI and 11 HaeIII types of the second exon based on PCR-RFLP analysis was defined, while heteroduplex analysis of the same animals showed 5 types (Ranjbar et al., 2016). In different studies on BuLA-DRB3 polymorphism, some alleles were more frequent in the Holstein breed, and their frequency was between 70 to 89% (Behl et al., 2012).

A quantitative measure for conservation against the variability of aligned sequences is useful for identifying sites under constraints. Various methods to quantitatively evaluate the conservation or variability of alignment sites have been developed (Johansson & Toh, 2010). Therefore, the challenge of which regions in the MHC interact with antigen in BoLA/BuLA-DRB3.2 is still an open one, and clear identification of areas of greater amino acid sequence variability is still of paramount importance.

Sheep

Although considerable effort has been devoted to identifying genetic variations at the MHC in sheep, compared to other domestic animals, Ovar is still poorly characterized (Ilhan et al., 2016; Siva Subramaniam et al., 2015). Several studies have demonstrated that the Ovar-DRB1 and Ovar-DQB loci can be highly polymorphic and are of great importance for understanding the host-pathogen interactions, facilitating the selection of disease-resistant animals, and accelerating the development of effective vaccines (Scott et al., 1992; Sayers et al., 2005; Ballingall et al., 2008; Nikbakht et al., 2009; Jamshidi et al., 2011; Gholamreza Nikbakht et al., 2012; Shen et al., 2014). In the following, the results of the previous studies in this field are summarized.

The MHC gene of sheep is located on chromosome 20 and is called Ovar (Hediger et al., 1991). At least six MHC class II loci in sheep are transcribed (DRB1, DRA1, DQA1, DQB1, DQA2, and DQB2). The DRB1 locus encodes highly polymorphic glycoproteins composed of α and β-subunits. The first domain of DRB1 β-subunit (β1) is similar to HLA-DR1 β1, which contains the peptide-binding site. In particular, the polymorphism is concentrated in the pockets of the antigen-binding site (Herrmann-Hoesing et al., 2008).

In sheep, the polymorphism of Ovar-DRB1 has been defined using SSOP, SSCP, RFLP, DNA sequencing, and a Double Amplification Refractory Mutation System (ARMS) method, which recently has been successfully used for haplotyping (Konnai et al., 2003). There are over 100 different DRB1 DNA sequences reported in GenBank. We have previously reported genetic diversity of the exon 2 Ovar-DRB1 locus by either RFLP or direct sequencing in three different Iranian fat-tailed sheep breeds: Lori-Bakhtiari (LB), Shaul (Sh), and Zandi (Za) (Table 4) (Gholamreza Nikbakht et al., 2012).

Table 4. The genetic diversity of the exon 2 Ovar-DRB1 locus in Iranian fat-tailed sheep breeds

Recent studies have focused on the relationship between different Ovar alleles and resistance/susceptibility to parasitic diseases. In this regard, a study demonstrated that the DRB1*1101 allele was associated with reduced egg counts and reduced IgA activity against nematode parasites in British Texel sheep (Ali et al., 2019). Another study performed by Shen et al. revealed that Chinese Merino sheep with the genic haplotype DRB1-SacIab/DRB1-Mvalbb/DQB1- TaqIaa/DQB1-HaeIIInn were relatively resistant to cystic echinococcosis, while individuals with the genic haplotypes DRB1-Mvalbc/DQB1- Mvalyy/DQB1-TaqIab/DQB1-HaeIIImn and DRB1-Mvalbb/DQB1-Mvalcc/ DQB1-TaqIab/DQB1-HaeIIImn were more susceptible (Shen et al., 2014). Another factor that can discriminate resistant/susceptible individuals is the expression level of certain MHC genes. For instance, Keane et al. showed that the MHC class II gene Ovar-DQA1 was expressed 8.4-fold more highly in resistant sheep to the gastrointestinal nematodes than susceptible ones (Keane et al., 2007).

Previous studies have not been limited to evaluating the association between Ovar alleles and parasitic diseases. For instance, in a recently published study, the association of MHC class II alleles with inflammatory cytokines and acute-phase proteins was evaluated in sheep. According to the results, allele DRB1*2101 showed a negative influence on the IL-6 response and was associated with a lower serum level of IL-6. Also, compared to the homozygotes, haptoglobin concentration was higher in DRB1.2 heterozygous individuals (Atefeh Esmailnejad et al., 2021). Another research demonstrated that the DRB1*0404 haplotype was associated with increased mature weights and average daily gain in sheep, while the DRB1*0141 haplotype was associated with decreased growth (Cinar et al., 2016).

Bird MHC (B Complex)

In the late 1940s, Briles described two highly polymorphic blood group systems in chickens. One of the systems, designated B, was later found to determine erythroid alloantigens and control acute allograft rejection and genetic control of the immune response to natural and synthetic Antigens. The B system of the chicken was the second major histocompatibility complex (MHC) recognized in animals (W. E. Briles et al., 1982; W. Elwood Briles & Briles, 1982).

The chicken MHC is divided between two regions of micro chromosome 16; the typical or "classical" MHC genes are located in the B complex and "nonclassical genes" in the Rfp-Y region (W. E. Briles et al., 1982; W. Elwood Briles & Briles, 1982). The core B complex is simple and compact, with only 19 genes in approximately 92kb, and is composed of tightly linked polymorphic regions: BF (class I) and BL (class IIβ), and BG (class IV) (Miller & Taylor, 2016; Yuan et al., 2021).

There has been growing interest in genetic markers suitable for drawing population differences in qualitative and quantitative traits affecting immune responses. Besides the importance of identifying alleles in each locus, MAS for chicken has reached a precious finding on microsatellite markers located in the MHC region. Microsatellites are one of the DNA-based genetic markers that have been widely used as a standard technique for molecular genetic evaluation and mapping of a chicken's genome (Abebe et al., 2015). Microsatellites (or short tandem repeats) are tandem repeat loci with a core motif of 1 to 6 base pairs repeated several times. The variations in length define the alleles at a particular microsatellite locus within a population. Because of the relative ease of identification of microsatellites, exhibiting a high level of polymorphism, random distribution throughout the genome, and a high level of heterozygosity make its application a quite significant genetic tool. Microsatellite markers are highly informative and could be used for genotyping, individual DNA identification, parentage assignment, phylogeny, population and conservation genetics, molecular epidemiology, pathology, and Quantitative trait loci (QTL) mapping. Those microsatellites that are physically located within the MHC are important for MHC typing. Though microsatellite functions in the genome have not been identified, it has been asserted that they have a linkage with some functional genes (Abdurakhmonov, 2016; Vieira et al., 2016).

One of the most applicable immune microsatellite markers for chicken is LEI0258. It is a tetra-nucleotide tandem repeat that is located in the chromosome 16a region where there are genes of chicken MHC (Atefeh Esmailnejad et al., 2017; B. Han et al., 2013). An additional marker, MCW0371, is also located 10,560 bp downstream of LEI0258 in this region (Manjula et al., 2021). Interestingly, LEI0258 has many alleles and an extensive range in allele sizes. These alleles are easily distinguishable using relatively inexpensive electrophoretic size separation methods. Due to its close physical location to genes of the chicken MHC, LEI0258 was investigated as a genetic indicator for B haplotypes (Manjula et al., 2020; Mwambene et al., 2019). MHC genes play an important role in the immune system, autoimmunity, reproduction, economic traits, and life-history strategies. Therefore, MHC typing in the poultry industry could improve breeding strategies and help speculation about its future applications. Linkage of LEI0258 microsatellite with B-F haplotype and association of these microsatellite alleles with antibody response against Newcastle disease virus vaccine and body weight as well as ectoparasite were evaluated (Mpenda et al., 2020; Gholamreza Nikbakht et al., 2013).

The chicken MHC is known to have a very strong association with disease resistance/susceptibility to numerous pathogens (Marek's disease virus, Rous sarcoma virus, avian leukosis virus, fowl cholera, coccidiosis, and salmonella) and autoimmune diseases, obese strain spontaneous autoimmune thyroiditis (SAT) (Silva & Gallardo, 2020; J. E. Fulton, 2020; Tizard, 2018). Our studies have shown that polymorphisms in the chicken MHC region are associated with the immune response to vaccination and productive traits (Alkaragoly et al., 2018; Emam et al., 2014; Gholamreza Nikbakht & Esmailnejad, 2015). Moreover, economically important traits, such as juvenile and adult mortality, body weight, fertilization rate, embryonic mortality, hatchability, sexual maturity, and egg production, are influenced by MHC genotypes (Izadi et al., 2011; Lwelamira et al., 2008). These characteristics of the chicken B system have made it an important subject in commercial chicken breeding programs.

Merely selection based on the productive traits of broiler chicken has caused many physiological disorders as well as a reduction in overall immunocompetence. To simultaneously improve production and fitness traits, molecular markers associated with one or both sets of traits may be useful (Manjula et al., 2021). However, in the first decade after identifying chicken MHC, several studies were performed to examine its association with production traits. Most of these studies have been done on congenic or commercial inbred lines with low polymorphism. In addition, in most of these studies, MHC typing was carried out by serology that is now replaced with molecular methods (Livant et al., 2004; Martin et al., 1990).

MHC association with production traits can result from various causes, including:

• Differences between B genotypes in resistance to various sub-clinical infections affecting growth and production traits.
• Interaction between MHC molecules and peptide hormone receptors active in metabolism or affecting fertility.
• Linkage disequilibrium with causative genes controlling production traits (Ewald et al., 2007).

Production traits are multifactorial, controlled by a relatively large number of loci (QTLs), have low heritability, and are influenced by environmental factors. A candidate gene approach would be a powerful method for understanding the direct genetic basis involved in the expression of quantitative differences between individuals. So, we looked for a marker that doesn't have these limitations and is associated with production traits. Insulin-Like Growth Factors (IGFs) consist of a family of polypeptide hormones structurally associated with insulin with multiple metabolic and anabolic functions. IGF1 is an important regulator in stimulating growth, protein synthesis, and cell proliferation in various cell types. IGF1 polymorphism is associated with several productions and reproductive traits in poultry. They proposed that the IGF1 gene could link disequilibrium with the actual causative mutations that affect growth and carcass traits (A. Esmailnejad & Nikbakht Brujeni, 2016). Our results for the Khorasan indigenous chicken indicated significant linkage disequilibrium between LEI0258 and IGF1 loci. High MHC genetic diversity in indigenous breeds is a valuable genetic resource for associations and linkage studies (A. Esmailnejad & Nikbakht Brujeni, 2016).

Dog

The genes of canine MHC, also known as dog leukocyte antigen (DLA), are located close together on canine chromosome 12 (Niskanen et al., 2016). DLA comprises three gene clusters, including class I, II, and III genes, which have a high genetic similarity to human MHC. DLA class I and II molecules are involved in presenting self and non-self-antigens to the immune cells and the regulation of immune response. In contrast, the role of DLA class III molecules in the antigen presentation process is still unclear (Kennedy, Barnes, et al., 2006). Several studies have identified a large number of alleles in different loci of DLA. For instance, at DRB1, DQA1, and DQB1 loci of DLA-class II region 90, 22, and 54 alleles have been identified, respectively (Kennedy et al., 2007). Also, the association of these alleles with a wide range of immune and non-immune responses has been investigated (Kennedy et al., 2007; Vahedi et al., 2019). The results of some of these studies are summarized below.

Canine mammary gland tumors (CMTs) are considered the most common neoplasm in female dogs. Vahedi et al. evaluated the DLA-DRB1 genetic polymorphism in a dog population. Subsequently, they analyzed the association of obtained MHC genotypes with mammary gland tumor profiles. According to the results, there was a significant correlation between DLA-DRB1.2 genotypes and different CMT profiles. Some genotypes were significantly associated with the increased risk of carcinoma arising in a benign tumor. Some others represented a positive correlation with complex carcinoma. Also, heterozygote genotypes were associated with a lower risk of CMT (Vahedi et al., 2019). Another study evaluated the association between anal sac gland carcinoma and MHC class II loci in English Cocker Spaniels. The results revealed that the DLA-DQB1*00701 allele was associated with susceptibility to developing the carcinoma, while the allele DLA-DQB1*02001 was resistant (Aguirre-Hernández et al., 2010).

Recent research in this field has not been limited to tumors. For instance, Bozorgpanah et al. have evaluated the association between canine atopic dermatitis and DLA-DRB1 alleles. The results showed that the presence of the type D allele in the exon II of the DLA-DRB1 gene increases the risk of atopic dermatitis (Bozorgpanah et al., 2020). In another study, Hughes et al. evaluated the association of a specific type of DLA class II haplotype with hypoadrenocorticism in Nova Scotia Duck Tolling Retrievers (NSDTRs). The results demonstrated that the development of hypoadrenocorticism in NSDTRs is potentially greater if the dogs are homozygous for the DLA-DRB1*01502/DQA* 00601/DQB1*02301 haplotype (Hughes et al., 2010). Another study on NSDTR dogs revealed an elevated risk for immune-mediated rheumatic disease (IMRD) in NSDTR dogs that were homozygous for the DLA-DRB1*00601/DQA1*005011/ DQB1*02001 haplotype (Wilbe et al., 2009). Moreover, Kennedy et al. have examined the association of canine hypothyroidism and DLA alleles in several breeds of dogs. This study indicated a significant association between the occurrence of canine hypothyroidism and the DLA class II haplotype that contains the DLADQA1*00101 allele (Kennedy, Quarmby, et al., 2006). According to the results of another research, two potential DLA risk haplotypes for the canine primary immune-mediated hemolytic anemia were DLA-DRB1*00601/DQA1* 005011/ DQB1*00701 and DLA-DRB1*01501/ DQA1* 00601/DQB1*00301, while DLA-DRB1*00101/ DQA1* 00101/DQB1*00201 was a protective one (Kennedy, Barnes, et al., 2006). Several similar studies have evaluated the association between DLA haplotypes and alleles and a variety of pathological conditions in dogs, including chronic hepatitis (Bexfield et al., 2012), hypoadrenocorticism (Massey et al., 2013), chronic superficial keratitis (Barrientos et al., 2013), pancreatic acinar atrophy (Tsai et al., 2013), necrotizing meningoencephalitis (Greer et al., 2010), and anal furunculosis (Barnes et al., 2009).

Concluding Remarks

Due to antibiotic resistance and vaccination problems, more appropriate approaches are desired for breeding disease resistance animals, which is now practiced in cattle, sheep, chicken, and fish farming. Traditional high selection intensity for growth and production rate has caused many metabolic disorders and reduced overall immunocompetence and genetic diversity. Molecular markers associated with both resistance and fitness traits would be worthwhile in "Marker Assisted Selection" programs. The important role of MHC in disease resistance and production traits makes it a precious marker in predicting and controlling programs. It is worth noting that conserving genetic resources and diversity across populations are indispensable requirements for sustainable development. Conserving genetic resources and diversity across populations are essential requirements for sustainable development.

Acknowledgments

We would like to thank the veterinary research council of Tehran University.

Conflict of Interest

The authors declared no conflict of interest.