Selection on MHC-linked microsatellite loci in sheep populations

Polymorphism of Selected Regions of Ovar-MHC and the Health Status of the Ovine Mammary Gland

Udder diseases (mastitis) are a serious cause of economic losses in sheep breeding as they have a negative impact on lamb rearing and the quality of dairy products. Thus far, progress in treatment and prevention of these diseases has been insufficient-giving ground for searching possibilities of using natural immunity to combat mastitis. This study aims to assess the relationship between the microsatellite polymorphism of selected Ovar-MHC genes and the health status of the mammary gland of sheep. The research is carried out on sheep of the Polish Heath and Polish Lowland breeds. In ovine milk, the SCC and the percentage of the lymphocyte subpopulation are assessed. Based on genomic DNA, molecular analysis of the Ovar-MHC gene fragments (OLADRB1, OLADRB2, OMHC1) polymorphism is performed. Significant differences are found in SCC level and the percentage of lymphocytes (CD4, CD8, CD19) in the milk depending on the alleles of the Ovar-MHC genes. Alleles of 488 bp (DRB1) and 284 bp (DRB2) are found more frequently in sheep, which milk contained <200 × 10³/mL SCC, while in carriers of the 508 bp (DRB1) and 272 bp (DRB2) alleles, SCC level in milk is significantly higher (>200 × 10³/mL). The obtained results justify the need for further research to better understand the genetic basis of mastitis, and to search for effective molecular markers that can be used in breeding practice.

Genetic variability and ontogeny predict microbiome structure in a disease-challenged montane amphibian

Amphibian populations worldwide are at risk of extinction from infectious diseases, including chytridiomycosis caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd). Amphibian cutaneous microbiomes interact with Bd and can confer protective benefits to the host. The composition of the microbiome itself is influenced by many environment- and host-related factors. However, little is known about the interacting effects of host population structure, genetic variation and developmental stage on microbiome composition and Bd prevalence across multiple sites. Here we explore these questions in Amietia hymenopus, a disease-affected frog in southern Africa. We use microsatellite genotyping and 16S amplicon sequencing to show that the microbiome associated with tadpole mouthparts is structured spatially, and is influenced by host genotype and developmental stage. We observed strong genetic structure in host populations based on rivers and geographic distances, but this did not correspond to spatial patterns in microbiome composition. These results indicate that demographic and host genetic factors affect microbiome composition within sites, but different factors are responsible for host population structure and microbiome structure at the between-site level. Our results help to elucidate complex within- and among- population drivers of microbiome structure in amphibian populations. That there is a genetic basis to microbiome composition in amphibians could help to inform amphibian conservation efforts against infectious diseases.

Microsatellite markers for evaluating the diversity of the natural killer complex and major histocompatibility complex genomic regions in domestic horses

Genotyping microsatellite markers represents a standard, relatively easy, and inexpensive method of assessing genetic diversity of complex genomic regions in various animal species, such as the major histocompatibility complex (MHC) and/or natural killer cell receptor (NKR) genes. MHC-linked microsatellite markers have been identified and some of them were used for characterizing MHC polymorphism in various species, including horses. However, most of those were MHC class II markers, while MHC class I and III sub-regions were less well covered. No tools for studying genetic diversity of NKR complex genomic regions are available in horses. Therefore, the aims of this work were to establish a panel of markers suitable for analyzing genetic diversity of the natural killer complex (NKC), and to develop additional microsatellite markers of the MHC class I and class III genomic sub-regions in horses. Nine polymorphic microsatellite loci were newly identified in the equine NKC. Along with two previously reported microsatellites flanking this region, they constituted a panel of 11 loci allowing to characterize genetic variation in this functionally important part of the horse genome. Four newly described MHC class I/III-linked markers were added to 11 known microsatellites to establish a panel of 15 MHC markers with a better coverage of the class I and class III sub-regions. Major characteristics of the two panels produced on a group of 65 horses of 13 breeds and on five Przewalski's horses showed that they do reflect genetic variation within the horse species.

MHC haplotype diversity in Persian Arabian horses determined using polymorphic microsatellites

Previous research on the equine major histocompatibility complex (MHC) demonstrated strong correlations between haplotypes defined by polymorphic intra-MHC microsatellites and haplotypes defined using classical serology. Here, we estimated MHC diversity in a sample of 124 Arabian horses from an endangered strain native to Iran (Persian Asil Arabians), using a validated 10-marker microsatellite panel. In a group of 66 horses related as parent-offspring pairs or half-sibling groups, we defined 51 MHC haplotypes, 49 of which were new. In 47 of the remaining 58 unrelated horses, we could assign one previously identified MHC haplotype, and by default, we gave provisional haplotype status to the remaining constellation of microsatellite alleles. In these horses, we found 21 haplotypes that we had previously defined and 31 provisional haplotypes, two of which had been identified in an earlier study. This gave a total of 78 new MHC haplotypes. The final 11 horses were MHC heterozygotes that we could not phase using information from any of the previously validated or provisional haplotypes. However, we could determine that these horses carried a total of 22 different undefined haplotypes. In the overall population sample, we detected three homozygous horses and one maternally inherited recombinant from 21 informative segregations. Virtually all of the horses tested were MHC heterozygotes, and most unrelated horses (98%) were heterozygous for rare microsatellite-defined haplotypes found less than three times in the sampled horses. This is evidence for a very high level of MHC haplotype variation in the Persian Asil Arabian horse.

Estimation of genetic diversity in Gute sheep: pedigree and microsatellite analyses of an ancient Swedish breed

Background

Breeds with small population size are in danger of an increased inbreeding rate and loss of genetic diversity, which puts them at risk for extinction. In Sweden there are a number of local breeds, native breeds which have adapted to specific areas in Sweden, for which efforts are being made to keep them pure and healthy over time. One example of such a breed is the Swedish Gute sheep. The objective of this study was to estimate inbreeding and genetic diversity of Swedish Gute sheep.

Results

Three datasets were analysed: pedigree information of the whole population, pedigree information for 100 animals of the population, and microsatellite genotypes for 94 of the 100 animals. The average inbreeding coefficient for lambs born during a six year time period (2007–2012) did not increase during that time period. The inbreeding calculated from the entire pedigree (0.038) and for a sample of the population (0.018) was very low. Sheep were more heterozygous at the microsatellite markers than expected (average multilocus heterozygosity and Ritland inbreeding estimates 1.01845 and -0.03931) and five of seven microsatellite markers were not in Hardy Weinberg equilibrium due to heterozygosity excess. The total effective population size estimated from the pedigree information was 155.4 and the average harmonic mean effective population size estimated from microsatellites was 88.3. Pedigree and microsatellite genotype estimations of inbreeding were consistent with a breeding program with the purpose of reducing inbreeding.

Conclusion

Our results showed that current breeding programs of the Swedish Gute sheep are consistent with efforts of keeping this breed viable and these breeding programs are an example for other small local breeds in conserving breeds for the future.

Electronic supplementary material

The online version of this article (doi:10.1186/s41065-017-0026-4) contains supplementary material, which is available to authorized users.

Genotyping of major histocompatibility complex Class II DRB gene in Rohilkhandi goats by polymerase chain reaction-restriction fragment length polymorphism and DNA sequencing

Aim:

To study the major histocompatibility complex (MHC) Class II DRB1 gene polymorphism in Rohilkhandi goat using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) and nucleotide sequencing techniques.

Materials and Methods:

DNA was isolated from 127 Rohilkhandi goats maintained at sheep and goat farm, Indian Veterinary Research Institute, Izatnagar, Bareilly. A 284 bp fragment of exon 2 of DRB1 gene was amplified and digested using BsaI and TaqI restriction enzymes. Population genetic parameters were calculated using Popgene v 1.32 and SAS 9.0. The genotypes were then sequenced using Sanger dideoxy chain termination method and were compared with related breeds/species using MEGA 6.0 and Megalign (DNASTAR) software.

Results:

TaqI locus showed three and BsaI locus showed two genotypes. Both the loci were found to be in Hardy–Weinberg equilibrium (HWE), however, population genetic parameters suggest that heterozygosity is still maintained in the population at both loci. Percent diversity and divergence matrix, as well as phylogenetic analysis revealed that the MHC Class II DRB1 gene of Rohilkhandi goats was found to be in close cluster with Garole and Scottish blackface sheep breeds as compared to other goat breeds included in the sequence comparison.

Conclusion:

The PCR-RFLP patterns showed population to be in HWE and absence of one genotype at one locus (BsaI), both the loci showed excess of one or the other homozygote genotype, however, effective number of alleles showed that allelic diversity is present in the population. Sequence comparison of DRB1 gene of Rohilkhandi goat with other sheep and goat breed assigned Rohilkhandi goat in divergence with Jamanupari and Angora goats.

Admixture and local breed marginalization threaten Algerian sheep diversity

Due to its geo-climatic conditions, Algeria represents a biodiversity hotspot, with sheep breeds well adapted to a patchwork of extremely heterogeneous harsh habitats. The importance of this peculiar genetic reservoir increases as climate change drives the demand for new adaptations. However, the expansion of a single breed (Ouled-Djellal) which occurred in the last decades has generated a critical situation for the other breeds; some of them are being subjected to uncontrolled cross-breeding with the favored breed and/or to marginalization (effective size contraction). This study investigated genetic diversity within and among six of the nine Algerian breeds, by use of 30 microsatellite markers. Our results showed that, in spite of the census contraction experienced by most of the considered breeds, genetic diversity is still substantial (average gene diversity ranging 0.68 to 0.76) and inbreeding was not identified as a problem. However, two breeds (Rembi and Taâdmit) appeared to have lost most of their genetic originality because of intensive cross-breeding with Ouled-Djellal. Based on the above evidence, we suggest Hamra, Sidaoun, and D’man as breeds deserving the highest priority for conservation in Algeria.

Selection methods for resistance to and tolerance of helminths in livestock

Helminthiases are among the most important livestock diseases worldwide, in particular for small ruminants, which are the focus of this review. Resource Allocation Theory implies that high-productivity farm animals proportionate insufficient resources for adequate coping with stressful conditions. Significant differences between breeds and within breeds are seen, as well as genotype vs. environment interactions. With improvement of genetic host resistance to infection, transmission of infection will be impacted. On the other hand, genetic improvement of resilience can lead to a reduction in clinical signs of disease, but not necessarily reduce transmission of infection to other animals. Faecal egg count (FEC) is the main measurement used to evaluate helminthiasis load, despite the fact that the protocols and analytical methods can affect the results, and the FEC data frequently shows aggregative, negative skewed distribution, and a high coefficient of variation. Mass selection where heritability is generally medium to low generally produces slow results and low economic returns. Many studies have been published linking resistance to nematodes in livestock to Quantitative Trait Loci and most studies have concentrated on chromosomes where the major histocompatibility complex region is located. Nevertheless, these complex traits have been seen to be affected by thousands of variants that each has a small effect. More recent studies have shown that genome-wide selection strategies can be useful in selecting animals for improved production and resistance traits in this case.

Development of MHC-Linked Microsatellite Markers in the Domestic Cat and Their Use to Evaluate MHC Diversity in Domestic Cats, Cheetahs, and Gir Lions

Diversity within the major histocompatibility complex (MHC) reflects the immunological fitness of a population. MHC-linked microsatellite markers provide a simple and an inexpensive method for studying MHC diversity in large-scale studies. We have developed 6 MHC-linked microsatellite markers in the domestic cat and used these, in conjunction with 5 neutral microsatellites, to assess MHC diversity in domestic mixed breed (n = 129) and purebred Burmese (n = 61) cat populations in Australia. The MHC of outbred Australian cats is polymorphic (average allelic richness = 8.52), whereas the Burmese population has significantly lower MHC diversity (average allelic richness = 6.81; P < 0.01). The MHC-linked microsatellites along with MHC cloning and sequencing demonstrated moderate MHC diversity in cheetahs (n = 13) and extremely low diversity in Gir lions (n = 13). Our MHC-linked microsatellite markers have potential future use in diversity and disease studies in other populations and breeds of cats as well as in wild felid species.

Isolation and characterization of 10 MHC Class I-associated microsatellite loci in tammar wallaby (Macropus eugenii)

The major histocompatibility complex (MHC) contain genes which play a key role in immune response and mate choice, and are therefore of functional importance to molecular ecologists. Here we describe the design of 10 MHC Class I-associated microsatellite loci from the tammar wallaby. All 10 loci are highly polymorphic, with the expected heterozygosity ranging from 0.547 to 0.919. Six loci successfully cross-amplify in other macropodid species. These microsatellites will serve as useful tools for studying the level of MHC diversity, the impact of selection on genetic variation and the unique structure of the tammar wallaby MHC.

Contrasted patterns of selection on MHC-linked microsatellites in natural populations of the Malagasy plague reservoir

Plague (Yersinia pestis infection) is a highly virulent rodent disease that persists in many natural ecosystems. The black rat (Rattus rattus) is the main host involved in the plague focus of the central highlands of Madagascar. Black rat populations from this area are highly resistant to plague, whereas those from areas in which the disease is absent (low altitude zones of Madagascar) are susceptible. Various lines of evidence suggest a role for the Major Histocompatibility Complex (MHC) in plague resistance. We therefore used the MHC region as a candidate for detecting signatures of plague-mediated selection in Malagasy black rats, by comparing population genetic structures for five MHC-linked microsatellites and neutral markers in two sampling designs. We first compared four pairs of populations, each pair including one population from the plague focus and one from the disease-free zone. Plague-mediated selection was expected to result in greater genetic differentiation between the two zones than expected under neutrality and this was observed for one MHC-class I-linked locus (D20Img2). For this marker as well as for four other MHC-linked loci, a geographic pattern of genetic structure was found at local scale within the plague focus. This pattern would be expected if plague selection pressures were spatially variable. Finally, another MHC-class I-linked locus (D20Rat21) showed evidences of balancing selection, but it seems more likely that this selection would be related to unknown pathogens more widely distributed in Madagascar than plague.

A novel approach for mining polymorphic microsatellite markers in silico

An important emerging application of high-throughput 454 sequencing is the isolation of molecular markers such as microsatellites from genomic DNA. However, few studies have developed microsatellites from cDNA despite the added potential for targeting candidate genes. Moreover, to develop microsatellites usually requires the evaluation of numerous primer pairs for polymorphism in the focal species. This can be time-consuming and wasteful, particularly for taxa with low genetic diversity where the majority of primers often yield monomorphic polymerase chain reaction (PCR) products. Transcriptome assemblies provide a convenient solution, functional annotation of transcripts allowing markers to be targeted towards candidate genes, while high sequence coverage in principle permits the assessment of variability in silico. Consequently, we evaluated fifty primer pairs designed to amplify microsatellites, primarily residing within transcripts related to immunity and growth, identified from an Antarctic fur seal (Arctocephalus gazella) transcriptome assembly. In silico visualization was used to classify each microsatellite as being either polymorphic or monomorphic and to quantify the number of distinct length variants, each taken to represent a different allele. The majority of loci (n = 36, 76.0%) yielded interpretable PCR products, 23 of which were polymorphic in a sample of 24 fur seal individuals. Loci that appeared variable in silico were significantly more likely to yield polymorphic PCR products, even after controlling for microsatellite length measured in silico. We also found a significant positive relationship between inferred and observed allele number. This study not only demonstrates the feasibility of generating modest panels of microsatellites targeted towards specific classes of gene, but also suggests that in silico microsatellite variability may provide a useful proxy for PCR product polymorphism.

A microsatellite-based analysis for the detection of selection on BTA1 and BTA20 in northern Eurasian cattle (Bos taurus) populations

BACKGROUND

Microsatellites surrounding functionally important candidate genes or quantitative trait loci have received attention as proxy measures of polymorphism level at the candidate loci themselves. In cattle, selection for economically important traits is a long-term strategy and it has been reported that microsatellites are linked to these important loci.

METHODS

We have investigated the variation of seven microsatellites on BTA1 (Bos taurus autosome 1) and 16 on BTA20, using bovine populations of typical production types and horn status in northern Eurasia. Genetic variability of these loci and linkage disequilibrium among these loci were compared with those of 28 microsatellites on other bovine chromosomes. Four different tests were applied to detect molecular signatures of selection.

RESULTS

No marked difference in locus variability was found between microsatellites on BTA1, BTA20 and the other chromosomes in terms of different diversity indices. Average D' values of pairwise syntenic markers (0.32 and 0.28 across BTA 1 and BTA20 respectively) were significantly (P < 0.05) higher than for non-syntenic markers (0.15). The Ewens-Watterson test, the Beaumont and Nichol's modified frequentist test and the Bayesian FST-test indicated elevated or decreased genetic differentiation, at SOD1 and AGLA17 markers respectively, deviating significantly (P < 0.05) from neutral expectations. Furthermore, lnRV, lnRH and lnRtheta' statistics were used for the pairwise population comparison tests and were significantly less variable in one population relative to the other, providing additional evidence of selection signatures for two of the 51 loci. Moreover, the three Finnish native populations showed evidence of subpopulation divergence at SOD1 and AGLA17. Our data also indicate significant intergenic linkage disequilibrium around the candidate loci and suggest that hitchhiking selection has played a role in shaping the pattern of observed linkage disequilibrium.

CONCLUSION

Hitchhiking due to tight linkage with alleles at candidate genes, e.g. the POLL gene, is a possible explanation for this pattern. The potential impact of selective breeding by man on cattle populations is discussed in the context of selection effects. Our results also suggest that a practical approach to detect loci under selection is to simultaneously apply multiple neutrality tests based on different assumptions and estimations.

Rapid selection against inbreeding in a wild population of a rare frog

Populations that are small and isolated can be threatened through loss of fitness due to inbreeding. Nevertheless, an increased frequency of recessive homozygotes could increase the efficiency of selection against deleterious mutants, thus reducing inbreeding depression. In wild populations, observations of evolutionary changes determined by selection against inbreeding are few. We used microsatellite DNA markers to compare the genetic features of tadpoles immediately after hatch with those of metamorphosing froglets belonging to the same cohort in a small, isolated population of the threatened frog Rana latastei. Within a generation, the inbreeding coefficient (F_IS) decreased: at hatch, F_IS was significantly >0, whereas F_IS was <0 after metamorphosis. Furthermore, heterozygosity increased and allelic frequencies changed over time, resulting in the loss of genotypes at metamorphosis that were present in hatchlings. One microsatellite locus exhibited atypically large F_ST values, suggesting it might be linked to a locus under selection. These results support the hypothesis that strong selection against the most inbred genotypes occurred among early life-history stages in our population. Selective forces can promote changes that can affect population dynamics and should be considered in conservation planning.

Does genetic diversity predict health in humans?

Genetic diversity, especially at genes important for immune functioning within the Major Histocompatibility Complex (MHC), has been associated with fitness-related traits, including disease resistance, in many species. Recently, genetic diversity has been associated with mate preferences in humans. Here we asked whether these preferences are adaptive in terms of obtaining healthier mates. We investigated whether genetic diversity (heterozygosity and standardized mean d²) at MHC and nonMHC microsatellite loci, predicted health in 153 individuals. Individuals with greater allelic diversity (d²) at nonMHC loci and at one MHC locus, linked to HLA-DRB1, reported fewer symptoms over a four-month period than individuals with lower d². In contrast, there were no associations between MHC or nonMHC heterozygosity and health. NonMHC-d² has previously been found to predict male preferences for female faces. Thus, the current findings suggest that nonMHC diversity may play a role in both natural and sexual selection acting on human populations.

High levels of genetic variation at MHC class II DBB loci in the tammar wallaby (Macropus eugenii)

High levels of MHC diversity are crucial for immunological fitness of populations, with island populations particularly susceptible to loss of genetic diversity. In this study, the level of MHC class II DBB diversity was examined in tammar wallabies (Macropus eugenii) from Kangaroo Island by genotyping class II-linked microsatellite loci and sequencing of DBB genes. Here we show that the tammar wallaby has at least four expressed MHC class II DBB loci and extensive genetic variation in the peptide-binding region of the DBB genes. These results contradict early studies which suggested that wallabies lacked MHC class II diversity and demonstrate that, in spite of the long-term isolation on an offshore island, this population of wallabies has a high level of DBB diversity.

Genetic structure of European sheep breeds

Large-scale evaluations of genetic diversity in domestic livestock populations are necessary so that region-specific conservation measures can be implemented. We performed the first such survey in European sheep by analysing 820 individuals from 29 geographically and phenotypically diverse breeds and a closely related wild species at 23 microsatellite loci. In contrast to most other domestic species, we found evidence of widespread heterozygote deficit within breeds, even after removing loci with potentially high frequency of null alleles. This is most likely due to subdivision among flocks (Wahlund effect) and use of a small number of rams for breeding. Levels of heterozygosity were slightly higher in southern than in northern breeds, consistent with declining diversity with distance from the Near Eastern centre of domestication. Our results highlight the importance of isolation in terms of both geography and management in augmenting genetic differentiation through genetic drift, with isolated northern European breeds showing the greatest divergence and hence being obvious targets for conservation. Finally, using a Bayesian cluster analysis, we uncovered evidence of admixture between breeds, which has important implications for breed management.