Uncovering novel MHC alleles from RNA-Seq data: expanding the spectrum of MHC class I alleles in sheep.

BACKGROUND: Major histocompatibility complex (MHC) class I glycoproteins present selected peptides or antigens to CD8 + T cells that control the cytotoxic immune response. The MHC class I genes are among the most polymorphic loci in the vertebrate genome, with more than twenty thousand alleles known in humans. In sheep, only a very small number of alleles have been described to date, making the development of genotyping systems or functional studies difficult. A cost-effective way to identify new alleles could be to use already available RNA-Seq data from sheep. Current strategies for aligning RNA-Seq reads against annotated genome sequences or transcriptomes fail to detect the majority of class I alleles. Here, I combine the alignment of RNA-Seq reads against a specific reference database with de novo assembly to identify alleles. The method allows the comprehensive discovery of novel MHC class I alleles from RNA-Seq data (DinoMfRS). RESULTS: Using DinoMfRS, virtually all expressed MHC class I alleles could be determined. From 18 animals 75 MHC class I alleles were identified, of which 69 were novel. In addition, it was shown that DinoMfRS can be used to improve the annotation of MHC genes in the sheep genome sequence. CONCLUSIONS: DinoMfRS allows for the first time the annotation of unknown, more divergent MHC alleles from RNA-Seq data. Successful application to RNA-Seq data from 16 animals has approximately doubled the number of known alleles in sheep. By using existing data, alleles can now be determined very inexpensively for populations that have not been well studied. In addition, MHC expression studies or evolutionary studies, for example, can be greatly improved in this way, and the method should be applicable to a broader spectrum of other multigene families or highly polymorphic genes.

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.

Association of MHC class II haplotypes with reduced faecal nematode egg count and IgA activity in British Texel sheep.

Nematode infection is one of the principal diseases suffered by sheep and the class II region of the MHC has been repeatedly associated with differences in susceptibility and resistance to infection. The aim of this study was to examine the association of MHC class II haplotypes in a flock of Texel sheep with faecal egg counts and antibody responsiveness. Two haplotypes carried the DRB1*11:01 allele which has previously been associated with reduced egg counts in Scottish Blackface and Suffolk sheep. One of the two haplotypes was associated with reduced egg counts in the Texel breed, and both haplotypes were associated with reduced IgA activity against an extract from fourth-stage larvae. The reduced IgA activity is probably a consequence of reduced numbers of fourth-stage larvae in sheep carrying the resistance allele. The association of specific MHC alleles with reduced egg counts, reduced worm numbers and decreased IgA activity provides a mechanism for the density-dependent regulation of parasite growth and fecundity.

The genetic architecture of the MHC class II region in British Texel sheep.

Understanding the structure of the major histocompatibility complex, especially the number and frequency of alleles, loci and haplotypes, is crucial for efficient investigation of the way in which the MHC influences susceptibility to disease. Nematode infection is one of the most important diseases suffered by sheep, and the class II region has been repeatedly associated with differences in susceptibility and resistance to infection. Texel sheep are widely used in many different countries and are relatively resistant to infection. This study determined the number and frequency of MHC class II genes in a small flock of Texel sheep. There were 18 alleles at DRB1, 9 alleles at DQA1, 13 alleles at DQB1, 8 alleles at DQA2 and 16 alleles at DQB2. Several haplotypes had no detectable gene products at DQA1, DQB1 or DQB2, and these were defined as null alleles. Despite the large numbers of alleles, there were only 21 distinct haplotypes in the population. The relatively small number of observed haplotypes will simplify finding disease associations because common haplotypes provide more statistical power but complicate the discrimination of causative mutations from linked marker loci.

A novel RNAseq-assisted method for MHC class I genotyping in a non-model species applied to a lethal vaccination-induced alloimmune disease.

BACKGROUND: MHC class I genotyping is essential for a wide range of biomedical, immunological and biodiversity applications. Whereas in human a comprehensive MHC class I allele catalogue is available, respective data in non-model species is scarce in spite of decades of research. RESULTS: Taking advantage of the new high-throughput RNA sequencing technology (RNAseq), we developed a novel RNAseq-assisted method (RAMHCIT) for MHC class I typing at nucleotide level. RAMHCIT is performed on white blood cells, which highly express MHC class I molecules enabling reliable discovery of new alleles and discrimination of closely related alleles due to the high coverage of alleles with reads. RAMHCIT is more comprehensive than previous methods, because no targeted PCR pre-amplification of MHC loci is necessary, which avoids preselection of alleles as usually encountered, when amplification with MHC class I primers is performed prior to sequencing. In addition to allele identification, RAMHCIT also enables quantification of MHC class I expression at allele level, which was remarkably consistent across individuals. CONCLUSIONS: Successful application of RAMHCIT is demonstrated on a data set from cattle with different phenotype regarding a lethal, vaccination-induced alloimmune disease (bovine neonatal pancytopenia), for which MHC class I alleles had been postulated as causal agents.

Major Histocompatibility Complex class IIB polymorphism in an ancient Spanish breed.

Genes from the Major Histocompatibility Complex class II region are involved in the presentation of antigens. Therefore, they have the key role in regulating the immune response and in the resistance to infections. We investigated the Major Histocompatibility Complex class IIB genes, DRB and DQB, in Churra sheep, one of the most important indigenous breeds of Spain. These genes are among the most polymorphic in the mammalian genome. Furthermore, often different numbers of class IIB genes per haplotype exist, complicating the genotyping and sequencing of these genes. Especially the DQB region is only partially characterized in sheep and the repertoire of DRB and DQB alleles in Churra sheep, an ancient breed, is unknown. Here, we sequenced the class IIB genes for 15 rams that are the pedigree heads of a selection Nucleus herd. In total, we found 12 DRB and 25 DQB alleles. From these, 3 and 15 were new, respectively. Fourteen haplotypes carrying one or two DQB alleles could be deduced and the evolutionary relationship of these was investigated by phylogenetic trees. Based on the sequences of these most common class II alleles, a more efficient genotyping system for larger numbers of Churra sheep will be developed.

Arachnomelia syndrome in Simmental cattle is caused by a homozygous 2-bp deletion in the molybdenum cofactor synthesis step 1 gene (MOCS1).

BACKGROUND: Arachnomelia syndrome is an autosomal recessive inherited disease in cattle. Affected calves die around birth and show malformations of the skeleton mainly affecting the legs, the spinal column and the skull. A number of arachnomelia syndrome affected Simmental calves were recently detected by a surveillance system of anomalies with a peak of more than 120 recorded cases in the year 2006. The causative mutation was previously mapped to a 9 cM-region on bovine chromosome 23. We herein report the fine-mapping and identification of the gene causing arachnomelia syndrome in Simmental cattle. RESULTS: By using a dense set of markers, the arachnomelia syndrome linked region could be refined to 1.5 cM harbouring three protein coding genes. Comparative sequencing of these genes revealed a two-bp-deletion in the bovine MOCS1 gene resulting in a frame-shift and a premature termination codon. We genotyped affected calves and their ancestors and found that all affected were homozygous for the deletion whereas all carriers were heterozygous. Furthermore, cattle from the same population, but not directly related to known carriers mostly showed the wild type genotype. CONCLUSIONS: MOCS1 encodes two proteins that are involved in the first synthesis step of molybdenum cofactor. A non functional sulfite-oxydase, one of the enzymes requiring molybdenum cofactor, leads to a similar pathology in Brown Swiss cattle. In combination the perfect association of the mutation with the phenotype and the obvious disruption of protein translation provide strong evidence for the causality of the MOCS1 mutation. Our results are the first example for an oligogenic lethal inherited disease in cattle. Furthermore, they show the potential involvement of sulfite metabolism in aberrant bone development.

Bone marrow depletion with haemorrhagic diathesis in calves in Germany: characterization of the disease and preliminary investigations on its aetiology.

Since 2007 a new fatal haemorrhagic diathesis in calves has been observed in all areas of Germany. Analysis of 56 cases submitted for necropsy allowed its characterization. Calves fell ill within the first month of life independent of breed and sex. Only single or a few animals per herd were affected. Petechial and ecchymotic haemorrhages in many organs and tissues, particularly in skin, subcutis and gastrointestinal tract, were major findings in all animals. Microscopically a severe depletion of bone marrow cells was always observed. Lymphocytic depletion (43%) and inflammatory lesions (46%) were less frequently observed. Blood analysis of five animals indicated an aplastic pancytopenia. The resulting thrombocytopenia is regarded as major pathomechanism of this Haemorrhagic Disease Syndrome (HDS). Pedigree analysis gave no indication of hereditary disease. Tests for specific toxins such as S-(1,2-Dichlorovinyl)-L-cysteine (DCVC), furazolidone, or mycotoxins resulting in bone marrow depletion were negative. Bacterial infections, Bovine Viral Diarrhoea Virus, and Bluetongue Virus were ruled out as cause of the disease. HDS shares similarities with a circoviral infection in chickens (chicken infectious anaemia). A broad-spectrum PCR allowed detection of circoviral DNA in 5 of 25 HDS cases and in 1 of 8 non-HDS cases submitted for necropsy. Sequencing of the whole viral genome revealed a high similarity (up to 99%) with Porcine Circovirus type 2b. Single bone marrow cells stained weakly positive for PCV2 antigen by immunohistochemistry in 1 of 8 tested HDS animals. This is the first report of circovirus detection in cattle in Germany. The exact cause of HDS still remains unknown. A multifactorial aetiology involving infection, poisoning, immunopathy, or a genetic predisposition is conceivable. Additional research is necessary to clarify the pathogenesis and the potential role of PCV2 in HDS.

Syndrome of arachnomelia in Simmental cattle.

BACKGROUND: The syndrome of arachnomelia is an inherited malformation mainly of limbs, back and head in cattle. At present the arachnomelia syndrome has been well known mainly in Brown Swiss cattle. Nevertheless, the arachnomelia syndrome had been observed in the Hessian Simmental population during the decade 1964-1974. Recently, stillborn Simmental calves were observed having a morphology similar to the arachnomelia syndrome. The goal of this work was the characterization of the morphology and genealogy of the syndrome in Simmental to establish the basis for an effective management of the disease. RESULTS: The first pathologically confirmed arachnomelia syndrome-cases in the current Simmental population appeared in the year 2005. By 2007, an additional 140 calves with the arachnomelia syndrome were identified. The major pathological findings were malformed bones affecting the head, long bones of the legs and the vertebral column. It could be shown that, with the exception of two cases that were considered as phenocopies, all of the paternal and about two-third of the maternal pedigrees of the affected calves could be traced back to one common founder. Together with the data from experimental matings, the pedigree data support an autosomal recessive mutation being the etiology of the arachnomelia syndrome. The frequency of the mutation in the current population was estimated to be 3.32%. CONCLUSION: We describe the repeated occurrence of the arachnomelia syndrome in Simmental calves. It resembles completely the same defect occurring in the Brown Swiss breed. The mutation became relatively widespread amongst the current population. Therefore, a control system has to be established and it is highly desirable to map the disease and develop a genetic test system.

Comparative gene mapping in cattle, Indian muntjac, and Chinese muntjac by fluorescence in situ hybridization.

The Indian muntjac (Muntiacus muntjak vaginalis) has a karyotype of 2n = 6 in the female and 2n = 7 in the male. The karyotypic evolution of Indian muntjac via extensive tandem fusions and several centric fusions are well documented by molecular cytogenetic studies mainly utilizing chromosome paints. To achieve higher resolution mapping, a set of 42 different genomic clones coding for 37 genes and the nucleolar organizer region were used to examine homologies between the cattle (2n = 60), human (2n = 46), Indian muntjac (2n = 6/7) and Chinese muntjac (2n = 46) karyotypes. These genomic clones were mapped by fluorescence in situ hybridization (FISH). Localization of genes on all three pairs of M. m. vaginalis chromosomes and on the acrocentric chromosomes of M. reevesi allowed not only the analysis of the evolution of syntenic regions within the muntjac genus but also allowed a broader comparison of synteny with more distantly related species, such as cattle and human, to shed more light onto the evolving genome organization.

A robust, low- to medium-throughput prnp genotyping system in sheep.

BACKGROUND: In many countries breeding programs for resistance to scrapie in sheep are established. Therefore, the demand on genotyping capacities of the polymorphisms of the prion protein gene (prnp) relevant to presently known disease associations and EU regulations is steadily increasing. Most published typing methods are not well suited for routine typing of large sample numbers in smaller service laboratories for different reasons: they require partly manual data processing, sophisticated and sensitive protocols, high efforts regarding time and manpower, multiple step reactions or substantial hardware investments. To overcome these drawbacks, we developed a prnp typing method that is based on a 'multiplex amplification refractory mutation system' (ARMS) reaction. METHODS: In this study we combined the amplification refractory mutation system (ARMS) with standard fluorescent based fragment length analyses method to develop a prnp genotyping method (PRNP ARMS). RESULTS: By optimised primer design it was possible to type the 4 relevant single nucleotide polymorphisms (SNPs) in the prnp simultaneously in one multiplex reaction. Automated fragment length analysis enabled automated allele designation. Suitability of the PRNP ARMS for routine application was proven by typing samples with known genotypes and larger sample numbers from half-sib families. CONCLUSION: The ARMS PRNP typing method established in this study is universally suited for a broad range of typing projects with different requirements. It provides an efficient and inexpensive diagnostic mutation analysis that will improve the quality of prnp genotyping compared with other low-cost methods. It can be implemented by most molecular genetic laboratories using standard equipment.

Association of a lysine-232/alanine polymorphism in a bovine gene encoding acyl-CoA:diacylglycerol acyltransferase (DGAT1) with variation at a quantitative trait locus for milk fat content.

DGAT1 encodes diacylglycerol O-acyltransferase (EC ), a microsomal enzyme that catalyzes the final step of triglyceride synthesis. It became a functional candidate gene for lactation traits after studies indicated that mice lacking both copies of DGAT1 are completely devoid of milk secretion, most likely because of deficient triglyceride synthesis in the mammary gland. Our mapping studies placed DGAT1 close to the region of a quantitative trait locus (QTL) on bovine chromosome 14 for variation in fat content of milk. Sequencing of DGAT1 from pooled DNA revealed significant frequency shifts at several variable positions between groups of animals with high and low breeding values for milk fat content in different breeds (Holstein-Friesian, Fleckvieh, and Braunvieh). Among the variants was a nonconservative substitution of lysine by alanine (K232A), with the lysine-encoding allele being associated with higher milk fat content. Haplotype analysis indicated the lysine variant to be ancestral. Two animals that were typed heterozygous (Qq) at the QTL based on marker-assisted QTL-genotyping were heterozygous for the K232A substitution, whereas 14 animals that are most likely qq at the QTL were homozygous for the alanine-encoding allele. An independent association study in Fleckvieh animals confirmed the positive effect of the lysine variant on milk fat content. We consider the nonconservative K232A substitution to be directly responsible for the QTL variation, although our genetic studies cannot provide formal proof.