Genetic diversity of three Chinese native sheep breeds

The sheep (Ovis aries L.) has been an important farm animal species since its domestication. A wide array of indigenous sheep breeds with abundant phenotypic diversity exists for domestication and selection. Therefore, assessing the genetic diversity of a local sheep resource using a multi-molecular system is helpful for maintaining and conserving those breeds. This study aimed to investigate the genetic diversity of three native Chinese sheep breeds (Tibetan sheep, Sishui Fur sheep, and Small-tailed Han sheep) using 15 microsatellite markers and the second exon of the DRA gene. In regards to the microsatellites, on average, 19 alleles per loci were observed among all individuals. Across loci, the HO within the population was 0.652 ± 0.022 in Tibetan sheep, 0.603 ± 0.023 in Small-tailed Han sheep and 0.635 ± 0.022 in SFS, and for most populations, the H _E and H _O were inconsistent. In addition, affluent private alleles within the breed indicated that the breeds have different domestication histories or sites. In regards to the 2 exon of the DRA gene, three haplotypes were constructed by seven single-nucleotide polymorphisms (SNPs), which were identified in the second DRA exon and inferred the potential for phenotypic variety in these Chinese native sheep. In summary, the current study reveals the importance of implementing effective conservation strategies for these three native Chinese sheep.     

Sequence-based genotyping of the sheep MHC class II DRB1 locus

The immunopolymorphism database (IPD) provides a single nomenclature for alleles at the major histocompatibility complex (MHC) loci for a range of different species. The minimum requirements for inclusion of a sheep class II DRB1 sequence is a submission that includes all polymorphic sites within the second exon from at least two independent polymerase chain reactions (PCR). In order to meet these requirements, we have developed a DNA-based genotyping method for the rapid analysis of allelic diversity at the DRB1 locus in domestic sheep, Ovis aries. Using a series of primers located within introns flanking exon 2 and genomic DNA from a cohort of 214 sheep representing 15 different breeds and crossbreeds, the complete exon 2 sequences of 38 Ovar-DRB1 alleles were obtained. This sequence resource allowed the development of a generic set of locus-specific primers which amplify a fragment that includes all polymorphic sites within the second exon. Bidirectional sequence analysis of the PCR product provides a composite sequence where each polymorphic site is represented by the corresponding International Union of Biochemistry nucleotide code. A Basic Local Alignment Search Tool search of alleles held within the IPD or National Center for Biotechnology Information databases allows individual allele sequences to be identified. Low levels of homozygosity (7.48%) within the cohort and verification of previously genotyped samples confirmed the broad allelic specificity of this method. It improves on currently available methods and is broadly applicable to the analysis of MHC diversity in studies investigating linkages with resistance or susceptibility to disease.     

Beta‐globin gene evolution in the ruminants: evidence for an ancient origin of sheep haplotype B

Domestic sheep (Ovis aries) can be divided into two groups with significantly different responses to hypoxic environments, determined by two allelic beta‐globin haplotypes. Haplotype A is very similar to the goat beta‐globin locus, whereas haplotype B has a deletion spanning four globin genes, including beta‐C globin, which encodes a globin with high oxygen affinity. We surveyed the beta‐globin locus using resequencing data from 70 domestic sheep from 42 worldwide breeds and three Ovis canadensis and two Ovis dalli individuals. Haplotype B has an allele frequency of 71.4% in O. aries and was homozygous (BB) in all five wild sheep. This shared ancestry indicates haplotype B is at least 2–3 million years old. Approximately 40 kb of the sequence flanking the ~37‐kb haplotype B deletion had unexpectedly low identity between haplotypes A and B. Phylogenetic analysis showed that the divergent region of sheep haplotype B is remarkably distinct from the beta‐globin loci in goat and cattle but still groups with the Ruminantia. We hypothesize that this divergent ~40‐kb region in haplotype B may be from an unknown ancestral ruminant and was maintained in the lineage to O. aries, but not other Bovidae, evolving independently of haplotype A. Alternatively, the ~40‐kb sequence in haplotype B was more recently acquired by an ancestor of sheep from an unknown non‐Bovidae ruminant, replacing part of haplotype A. Haplotype B has a lower nucleotide diversity than does haplotype A, suggesting a recent bottleneck, whereas the higher frequency of haplotype B suggests a subsequent spread through the global population of O. aries.     

Identification of the amino acids in the Major Histocompatibility Complex class II region of Scottish Blackface sheep that are associated with resistance to nematode infection

Lambs with the Major Histocompatibility Complex DRB1*1101 allele have been shown to produce fewer nematode eggs following natural and deliberate infection. These sheep also possess fewer adult Teladorsagia circumcincta than sheep with alternative alleles at the DRB1 locus. However, it is unclear if this allele is responsible for the reduced egg counts or merely acts as a marker for a linked gene. This study defined the MHC haplotypes in a population of naturally infected Scottish Blackface sheep by PCR amplification and sequencing, and examined the associations between MHC haplotypes and faecal egg counts by generalised linear mixed modelling. The DRB1*1101 allele occurred predominately on one haplotype and a comparison of haplotypes indicated that the causal mutation or mutations occurred in or around this locus. Additional comparisons with another resistant haplotype indicated that mutations in or around the DQB2*GU191460 allele were also responsible for resistance to nematode infections. Further analyses identified six amino acid substitutions in the antigen binding site of DRB1*1101 that were significantly associated with reductions in the numbers of adult T. circumcincta.     

Diversity of Extended HLA-DRB1 Haplotypes in the Finnish Population

The Major Histocompatibility Complex (MHC, 6p21) codes for traditional HLA and other host response related genes. The polymorphic HLA-DRB1 gene in MHC Class II has been associated with several complex diseases. In this study we focus on MHC haplotype structures in the Finnish population. We explore the variability of extended HLA-DRB1 haplotypes in relation to the other traditional HLA genes and a selected group of MHC class III genes. A total of 150 healthy Finnish individuals were included in the study. Subjects were genotyped for HLA alleles (HLA-A, -B, -DRB1, -DQB1, and -DPB1). The polymorphism of TNF, LTA, C4, BTNL2 and HLA-DRA genes was studied with 74 SNPs (single nucleotide polymorphism). The C4A and C4B gene copy numbers and a 2-bp silencing insertion at exon 29 in C4A gene were analysed with quantitative genomic realtime-PCR. The allele frequencies for each locus were calculated and haplotypes were constructed using both the traditional HLA alleles and SNP blocks. The most frequent Finnish A∼B∼DR -haplotype, uncommon in elsewhere in Europe, was A*03∼B*35∼DRB1*01∶01. The second most common haplotype was a common European ancestral haplotype AH 8.1 (A*01∼B*08∼DRB1*03∶01). Extended haplotypes containing HLA-B, TNF block, C4 and HLA-DPB1 strongly increased the number of HLA-DRB1 haplotypes showing variability in the extended HLA-DRB1 haplotype structures. On the contrary, BTNL2 block and HLA-DQB1 were more conserved showing linkage with the HLA-DRB1 alleles. We show that the use of HLA-DRB1 haplotypes rather than single HLA-DRB1 alleles is advantageous when studying the polymorphisms and LD patters of the MHC region. For disease association studies the HLA-DRB1 haplotypes with various MHC markers allows us to cluster haplotypes with functionally important gene variants such as C4 deficiency and cytokines TNF and LTA, and provides hypotheses for further assessment. Our study corroborates the importance of studying population-specific MHC haplotypes.     

Major histocompatibility complex class II genetic variation in forest musk deer (Moschus berezovskii) in China

The major histocompatibility complex (MHC) plays an important role in the immune system of vertebrates. We used the second exon of four MHC class II genes (DRA, DQA1, DQA2 and DRB3) to assess the overall MHC variation in forest musk deer (Moschus berezovskii). We also compared the MHC variation in captive and wild populations. We observed 22 alleles at four loci (four at DRA, four at DQA1, four at DQA2 and 10 at DRB3), 15 of which were newly identified alleles. Results suggest that forest musk deer maintain relatively high MHC variation, which may result from balancing selection. Moreover, considerable diversity was observed at the DRA locus. We found a high frequency of Mobe‐DRA*02, Mobe‐DQA1*01 and Mobe‐DQA2*05 alleles, which may be important for pathogen resistance. A Ewens–Watterson test showed that the DRB3 locus in the wild population had experienced recent balancing selection. We detected a small divergence at the DRA locus, suggesting the effect of weak positive selection on the DRA gene. Alternatively, this locus may be young and not yet adapted a wide spectrum of alleles for pathogen resistance. The significant heterozygosity deficit observed at the DQA1 and DRB3 loci in the captive population and at all four loci in the wild population may be the result of a population bottleneck. Additionally, MHC genetic diversity was higher in the wild population than in the captive, suggesting that the wild population may have the ability to respond to a wider range of pathogens.     

High-Density SNP Screening of the Major Histocompatibility Complex in Systemic Lupus Erythematosus Demonstrates Strong Evidence for Independent Susceptibility Regions

A substantial genetic contribution to systemic lupus erythematosus (SLE) risk is conferred by major histocompatibility complex (MHC) gene(s) on chromosome 6p21. Previous studies in SLE have lacked statistical power and genetic resolution to fully define MHC influences. We characterized 1,610 Caucasian SLE cases and 1,470 parents for 1,974 MHC SNPs, the highly polymorphic HLA-DRB1 locus, and a panel of ancestry informative markers. Single-marker analyses revealed strong signals for SNPs within several MHC regions, as well as with HLA-DRB1 (global p = 9.99×10⁻¹⁶). The most strongly associated DRB1 alleles were: *0301 (odds ratio, OR = 2.21, p = 2.53×10⁻¹²), *1401 (OR = 0.50, p = 0.0002), and *1501 (OR = 1.39, p = 0.0032). The MHC region SNP demonstrating the strongest evidence of association with SLE was rs3117103, with OR = 2.44 and p = 2.80×10⁻¹³. Conditional haplotype and stepwise logistic regression analyses identified strong evidence for association between SLE and the extended class I, class I, class III, class II, and the extended class II MHC regions. Sequential removal of SLE–associated DRB1 haplotypes revealed independent effects due to variation within OR2H2 (extended class I, rs362521, p = 0.006), CREBL1 (class III, rs8283, p = 0.01), and DQB2 (class II, rs7769979, p = 0.003, and rs10947345, p = 0.0004). Further, conditional haplotype analyses demonstrated that variation within MICB (class I, rs3828903, p = 0.006) also contributes to SLE risk independent of HLA-DRB1*0301. Our results for the first time delineate with high resolution several MHC regions with independent contributions to SLE risk. We provide a list of candidate variants based on biologic and functional considerations that may be causally related to SLE risk and warrant further investigation.     

Mitochondrial DNA Diversity of Modern,Ancient and Wild Sheep (Ovis gmelinii anatolica) from Turkey: New Insights on the Evolutionary History of Sheep

In the present study, to contribute to the understanding of the evolutionary history of sheep, the mitochondrial (mt) DNA polymorphisms occurring in modern Turkish native domestic (n = 628), modern wild (Ovis gmelinii anatolica) (n = 30) and ancient domestic sheep from Oylum Höyük in Kilis (n = 33) were examined comparatively with the accumulated data in the literature. The lengths (75 bp/76 bp) of the second and subsequent repeat units of the mtDNA control region (CR) sequences differentiated the five haplogroups (HPGs) observed in the domestic sheep into two genetic clusters as was already implied by other mtDNA markers: the first cluster being composed of HPGs A, B, D and the second cluster harboring HPGs C, E.To manifest genetic relatedness between wild Ovis gmelinii and domestic sheep haplogroups, their partial cytochrome B sequences were examined together on a median-joining network. The two parallel but wider aforementioned clusters were observed also on the network of Ovis gmelenii individuals, within which domestic haplogroups were embedded. The Ovis gmelinii wilds of the present day appeared to be distributed on two partially overlapping geographic areas parallel to the genetic clusters that they belong to (the first cluster being in the western part of the overall distribution). Thus, the analyses suggested that the domestic sheep may be the products of two maternally distinct ancestral Ovis gmelinii populations.Furthermore, Ovis gmelinii anatolica individuals exhibited a haplotype of HPG A (n = 22) and another haplotype (n = 8) from the second cluster which was not observed among the modern domestic sheep. HPG E, with the newly observed members (n = 11), showed signs of expansion. Studies of ancient and modern mtDNA suggest that HPG C frequency increased in the Southeast Anatolia from 6% to 22% some time after the beginning of the Hellenistic period, 500 years Before Common Era (BCE).     

Expression of the Multiple Sclerosis-Associated MHC Class II Allele HLA-DRB1*1501 Is Regulated by Vitamin D

Multiple sclerosis (MS) is a complex trait in which allelic variation in the MHC class II region exerts the single strongest effect on genetic risk. Epidemiological data in MS provide strong evidence that environmental factors act at a population level to influence the unusual geographical distribution of this disease. Growing evidence implicates sunlight or vitamin D as a key environmental factor in aetiology. We hypothesised that this environmental candidate might interact with inherited factors and sought responsive regulatory elements in the MHC class II region. Sequence analysis localised a single MHC vitamin D response element (VDRE) to the promoter region of HLA-DRB1. Sequencing of this promoter in greater than 1,000 chromosomes from HLA-DRB1 homozygotes showed absolute conservation of this putative VDRE on HLA-DRB1*15 haplotypes. In contrast, there was striking variation among non–MS-associated haplotypes. Electrophoretic mobility shift assays showed specific recruitment of vitamin D receptor to the VDRE in the HLA-DRB1*15 promoter, confirmed by chromatin immunoprecipitation experiments using lymphoblastoid cells homozygous for HLA-DRB1*15. Transient transfection using a luciferase reporter assay showed a functional role for this VDRE. B cells transiently transfected with the HLA-DRB1*15 gene promoter showed increased expression on stimulation with 1,25-dihydroxyvitamin D3 (P = 0.002) that was lost both on deletion of the VDRE or with the homologous “VDRE” sequence found in non–MS-associated HLA-DRB1 haplotypes. Flow cytometric analysis showed a specific increase in the cell surface expression of HLA-DRB1 upon addition of vitamin D only in HLA-DRB1*15 bearing lymphoblastoid cells. This study further implicates vitamin D as a strong environmental candidate in MS by demonstrating direct functional interaction with the major locus determining genetic susceptibility. These findings support a connection between the main epidemiological and genetic features of this disease with major practical implications for studies of disease mechanism and prevention.     

Applicability of SSCP analysis for MHC genotyping: fingerprinting of Ovar-DRB1 exon 2 alleles from Finnish and Russian breeds

The applicability of single strand conformation polymorphism (SSCP) analysis for major histocompatibility complex (MHC) genotyping in sheep was studied. A panel of Ovar-DRB1 exon 2 'allele fingerprints' was defined. The panel could accelerate DRB1 genotyping of new breeds when already existing sequences are used as references in SSCP analysis. In this study, seven new exon 2 sequences and 19 different alleles in total were detected from 31 animals of Finnish and Russian sheep breeds. Ovar-DRB1 * 0201 was detected in all the six grey Finnsheep animals included in this study, suggesting reduced MHC diversity within these animals.     

The chimpanzee major histocompatibility complex class II DR subregion contains an unexpectedly high number of beta-chain genes

The major histocompatibility complex (MHC) class II DR subregion of the chimpanzee was studied by restriction fragment length polymorphism (RFLP) analysis. Genomic DNA obtained from a panel of 94 chimpanzees was digested with the restriction enzyme Taq I and hybridized with an HLA-DR probe specific for the 3' untranslated (UT) region. Such a screening revealed the existence of 14 distinct DRB/Taq I gene-associated fragments allowing the definition of 11 haplotypes. Segregation studies proved that the number of chimpanzee class II DRB/Taq I fragments is not constant and varies from three to six depending on the haplotype. Comparison of these data with a human reference panel manifested that some MHC DRB/Taq I fragments are shared by man and chimpanzee. Moreover, the number of HLA-DRB/Taq I gene-associated fragments detected in a panel of homozygous typing cells varies from one to three and corresponds with the number of HLA-DRB genes present for most haplotypes. However, a discrepancy is observed for the HLA-DR4,-DR7, and -DR9 haplotypes since a fourth HLA-DRB pseudogene present within these haplotypes lacks its 3' UT region and thus is not detected with the probe used. These results suggest that chimpanzees have a higher maximum number of DRB genes per haplotype than man. As a consequence, some chimpanzee haplotypes must show a dissimilar organization of the MHC DR subregion compared to their human equivalents. The implications of these findings are discussed in the context of the trans-species theory of MHC polymorphism. Address correspondence and offprint requests to: R. E. Bontrop.     

Mitochondrial haplotypes reveal a strong genetic structure for three Indian sheep breeds

This survey represents the first characterization of mitochondrial DNA diversity within three breeds of Indian sheep (two strains of the Deccani breed, as well as the Bannur and Garole breeds) from different geographic regions and with divergent phenotypic characteristics. A 1061-bp fragment of the mitochondrial genome spanning the control region, a portion of the 12S rRNA gene and the complete phenyl tRNA gene, was sequenced from 73 animals and compared with the corresponding published sequence from European and Asian breeds and the European Mouflon (Ovis musimon). Analysis of all 156 sequences revealed 73 haplotypes, 52 of which belonged to the Indian breeds. The three Indian breeds had no haplotypes in common, but one Indian haplotype was shared with European and other Asian breeds. The highest nucleotide and haplotype diversity was observed in the Bannur breed (0.00355 and 0.981 respectively), while the minimum was in the Sangamneri strain of the Deccani breed (0.00167 and 0.882 respectively). All 52 Indian haplotypes belonged to mitochondrial lineage A. Therefore, these Indian sheep are distinct from other Asian and European breeds studied so far. The relationships among the haplotypes showed strong breed structure and almost no introgression among these Indian breeds, consistent with Indian sheep husbandry, which discourages genetic exchange between breeds. These results have implications for the conservation of India's ovine biodiversity and suggest a common origin for the breeds investigated.     

Sequence Polymorphism and Evolution of Three Cetacean MHC Genes

Sequence variability at three major histocompatibility complex (MHC) genes (DQB, DRA, and MHC-I) of cetaceans was investigated in order to get an overall understanding of cetacean MHC evolution. Little sequence variation was detected at the DRA locus, while extensive and considerable variability were found at the MHC-I and DQB loci. Phylogenetic reconstruction and sequence comparison revealed extensive sharing of identical MHC alleles among different species at the three MHC loci examined. Comparisons of phylogenetic trees for these MHC loci with the trees reconstructed only based on non-PBR sites revealed that allelic similarity/identity possibly reflected common ancestry and were not due to adaptive convergence. At the same time, trans-species evolution was also evidenced that the allelic diversity of the three MHC loci clearly pre-dated species divergence events according to the relaxed molecular clock. It may be the forces of balancing selection acting to maintain the high sequence variability and identical alleles in trans-specific manner at the MHC-I and DQB loci.     

Dominant Sequences of Human Major Histocompatibility Complex Conserved Extended Haplotypes from HLA-DQA2 to DAXX

We resequenced and phased 27 kb of DNA within 580 kb of the MHC class II region in 158 population chromosomes, most of which were conserved extended haplotypes (CEHs) of European descent or contained their centromeric fragments. We determined the single nucleotide polymorphism and deletion-insertion polymorphism alleles of the dominant sequences from HLA-DQA2 to DAXX for these CEHs. Nine of 13 CEHs remained sufficiently intact to possess a dominant sequence extending at least to DAXX, 230 kb centromeric to HLA-DPB1. We identified the regions centromeric to HLA-DQB1 within which single instances of eight “common” European MHC haplotypes previously sequenced by the MHC Haplotype Project (MHP) were representative of those dominant CEH sequences. Only two MHP haplotypes had a dominant CEH sequence throughout the centromeric and extended class II region and one MHP haplotype did not represent a known European CEH anywhere in the region. We identified the centromeric recombination transition points of other MHP sequences from CEH representation to non-representation. Several CEH pairs or groups shared sequence identity in small blocks but had significantly different (although still conserved for each separate CEH) sequences in surrounding regions. These patterns partly explain strong calculated linkage disequilibrium over only short (tens to hundreds of kilobases) distances in the context of a finite number of observed megabase-length CEHs comprising half a population''s haplotypes. Our results provide a clearer picture of European CEH class II allelic structure and population haplotype architecture, improved regional CEH markers, and raise questions concerning regional recombination hotspots.     

Mitochondrial DNA and Y‐chromosome diversity in East Adriatic sheep

Variation in mitochondrial DNA (mtDNA) and Y‐chromosome haplotypes was analysed in nine domestic sheep breeds (159 rams) and 21 mouflon ( Ovis musimon) sampled in the East Adriatic. Mitochondrial DNA analyses revealed a high frequency of type B haplotypes, predominantly in European breeds, and a very low frequency of type A haplotypes, which are more frequent in some Asian breeds. Mitochondrial haplotype Hmt‐3 was the most frequent (26.4%), and 37.1%, 20.8% and 7.6% of rams had haplotypes one, two and three mutations remote from Hmt‐3 respectively. In contrast, Y‐chromosome analyses revealed extraordinary paternal allelic richness: HY‐6, 89.3%; HY‐8, 5.0%; HY‐18, 3.1%; HY‐7, 1.3%; and HY‐5, 1.3%. In fact, the number of haplotypes observed is comparable to the number found in Turkish breeds and greater than the number found in European breeds so far. Haplotype HY‐18 (A‐oY1/135‐SRYM18), identified here for the first time, provides a link between the haplotype HY‐12 (A‐oY1/139‐SRYM18) found in a few rams in Turkey and haplotype HY‐9 (A‐oY1/131‐SRYM18) found in one ram in Ethiopia. All mouflons had type B mtDNA haplotypes, including the private haplotype (Hmt‐55), and all were paternally monomorphic for haplotype HY‐6. Our data support a quite homogeneous maternal origin of East Adriatic sheep, which is a characteristic of European breeds. At the same time, the high number of haplotypes found was surprising and intriguing, and it begs for further analysis. Simultaneous analysis of mtDNA and Y‐chromosome information allowed us to detect a large discrepancy between maternal and paternal lineages in some populations. This is most likely the result of breeder efforts to ‘upgrade’ local populations using rams with different paternal origins.     

Transcription and diversity of immunoglobulin lambda chain variable genes in the rat

The DRB region of the human major histocompatibility complex displays length polymorphism: Five major haplotypes differing in the number and type of genes they contain have been identified, each at appreciable frequency. In an attempt to determine whether this haplotype polymorphism, like the allelic polymorphism, predates the divergence of humansfrom great apes, we have worked out the organization of the DRB region of the chimpanzee Hugo using a combination of chromosome walking, pulsed-field gel electrophoresis, and sequencing. Hugo is a DRB homozygote whose single DRB haplotype is some 440 kilobases (kb) long and contains five genes. At least one and possibly two of these are pseudogenes, while three are presumably active genes. The genes are designated DRB ^* A0201, DRB2 ^* 0101, DRB3 ^* 0201, DRB6 ^* 0105, and DRB5 ^* 0301, and are arranged in this order on the chromosome. The DRB2 and DRB3 genes are separated by approximately 250 kb of sequence that does not seem to contain any additional DRB genes. The DRB ^* A0201 gene is related to the DRB1 gene of the human DR2 haplotype; the DRB2 ^* 0101 and DRB3 ^* 0201 genes are related to the DRB2 and DRB3 genes of the human DR3 haplotype, respectively; the DRB6 ^* 0105 and DRB5 ^* 0301 genes are related to the DRBVI and DRB5 genes of the human DR2 haplotype, respectively. Thus the Hugo haplotype appears to correspond to the entire human DR2 haplotype, into which a region representing a portion of the human DR3 haplotype has been inserted. Since other chimpanzees have their DRB regions organized in different ways, we conclude that, first, the chimpanzee DRB region, like the human DRB region, displays length polymorphism; second, some chimpanzee DRB haplotypes are longer than the longest known human DRB haplotypes; third, in some chimpanzee haplotypes at least, the DRB genes occur in combinations different from those of the human haplotypes; fourth, and most importantly, certain DRB gene combinations have been conserved in the evolution of chimpanzees and humans from their common ancestors. These data thus provide evidence that not only allelic but also haplotype polymorphism can be passed on from one species to another in a given evolutionary lineage.     

Application of Selection Mapping to Identify Genomic Regions Associated with Dairy Production in Sheep

In Europe, especially in Mediterranean areas, the sheep has been traditionally exploited as a dual purpose species, with income from both meat and milk. Modernization of husbandry methods and the establishment of breeding schemes focused on milk production have led to the development of “dairy breeds.” This study investigated selective sweeps specifically related to dairy production in sheep by searching for regions commonly identified in different European dairy breeds. With this aim, genotypes from 44,545 SNP markers covering the sheep autosomes were analysed in both European dairy and non-dairy sheep breeds using two approaches: (i) identification of genomic regions showing extreme genetic differentiation between each dairy breed and a closely related non-dairy breed, and (ii) identification of regions with reduced variation (heterozygosity) in the dairy breeds using two methods. Regions detected in at least two breeds (breed pairs) by the two approaches (genetic differentiation and at least one of the heterozygosity-based analyses) were labeled as core candidate convergence regions and further investigated for candidate genes. Following this approach six regions were detected. For some of them, strong candidate genes have been proposed (e.g. ABCG2, SPP1), whereas some other genes designated as candidates based on their association with sheep and cattle dairy traits (e.g. LALBA, DGAT1A) were not associated with a detectable sweep signal. Few of the identified regions were coincident with QTL previously reported in sheep, although many of them corresponded to orthologous regions in cattle where QTL for dairy traits have been identified. Due to the limited number of QTL studies reported in sheep compared with cattle, the results illustrate the potential value of selection mapping to identify genomic regions associated with dairy traits in sheep.     

MHC class II association study in eight breeds of dog with hypoadrenocorticism

Canine hypoadrenocorticism is an endocrine disorder characterised by inadequate secretion of steroid hormones from the adrenal glands. Pathology results from immune-mediated destruction of the adrenal cortex, which is similar to that seen in the human Addison’s disease. Both the canine and human diseases have similar clinical presentation, with the diagnosis based on performing a dynamic adrenocorticotropic hormone stimulation test. MHC class II has previously been associated with the human and canine diseases. In the current study, we conducted an MHC class II association study in eight breeds of dog with diagnoses of hypoadrenocorticism. We demonstrated significant differences in dog leukocyte antigen (DLA) haplotype frequencies in six of these breeds: Cocker spaniel, Springer spaniel, Labrador, West Highland white terrier (WHWT), Bearded collie, and Standard poodle. In the Springer spaniel, the DLA-DRB1*015:01--DQA1*006:01--DQB1*023:01 haplotype was significantly associated with disease risk (p?=?0.014, odds ratio (OR)?=?5.14) and showed a similar trend in the Cocker spaniel. This haplotype is related to one associated with hypoadrenocorticism in the Nova Scotia duck tolling retriever. Similar haplotypes shared between breeds were demonstrated, with DLA-DRB1*001:01--DQA1*001:01--DQB1*002:01 more prevalent in both affected Labrador (p?=?0.0002, OR?=?3.06) and WHWT (p?=?0.01, OR?=?2.11). Other haplotypes that have not previously been associated with the disease were identified. The inter-breed differences in DLA haplotypes associated with susceptibility to canine hypoadrenocorticism could represent divergent aetiologies. This could have implications for clinical diagnosis and future comparative studies. Alternatively, it may suggest that the gene of interest is closely linked to the MHC.     

Large herbivores facilitate an insect herbivore by modifying plant community composition in a temperate grassland

Large herbivores often co‐occur and share plant resources with herbivorous insects in grassland ecosystems; yet, how they interact with each other remains poorly understood. We conducted a series of field experiments to investigate whether and how large domestic herbivores (sheep; Ovis aries) may affect the abundance of a common herbivorous insect (aphid; Hyalopterus pruni) in a temperate grassland of northeast China. Our exclosure experiment showed that 3 years (2010–2012) of sheep grazing had led to 86% higher aphid abundance compared with ungrazed sites. Mechanistically, this facilitative effect was driven by grazing altering the plant community, rather than by changes in food availability and predator abundance for aphids. Sheep significantly altered plant community by reducing the abundance of unpalatable forbs for the aphids. Our small‐scale forb removal experiment revealed an “associational plant defense” by forbs which protect the grass Phragmites australis from being attacked by the aphids. However, selective grazing on forbs by sheep indirectly disrupted such associational plant defense, making P. australis more susceptible to aphids, consequentially increasing the density of aphids. These findings provide a novel mechanistic explanation for the effects of large herbivores on herbivorous insects by linking selective grazing to plant community composition and the responses of insect populations in grassland ecosystems.     

Extensive polymorphism of the major histocompatibility complex DRA gene in Balkan donkeys: perspectives on selection and genealogy

The major histocompatibility complex (MHC) contains genes important for immune response in mammals, and these genes exhibit high polymorphism and diversity. The DRA gene, a member of the MHC class II family, is highly conserved across a large number of mammalian species, but it displays exceptionally rich sequence variations in Equidae members. We analyzed allelic polymorphism of the DRA locus in 248 donkeys sampled across the Balkan Peninsula (Albania, Bulgaria, Croatia, Macedonia, Greece and Montenegro). Five known alleles and two new alleles were identified. The new allele Eqas‐DRA*0601 was found to carry a synonymous mutation, and new allele Eqas‐DRA*0701, a non‐synonymous mutation. We further analyzed the historical selection and allele genealogy at the DRA locus in equids. Signals of positive selection obtained by various tests were ambiguous. A conservative conclusion is that DRA polymorphism occurred relatively recently and that positive selection has been acting on the DRA locus for a relatively brief period.