Validating loci associated with bovine respiratory disease complex in pre-weaned Holstein calves.

Bovine respiratory disease (BRD) is considered one of the most economically important diseases in the cattle industry. Ultimately, the selection of cattle that are less susceptible to disease will allow producers to reduce the prevalence of BRD and lessen its economic impact. The objective of this study was to validate previously identified loci associated with susceptibility to BRD in an independent population of 140 pre-weaned Holstein calves from Wisconsin (WI). Using the McGuirk health scoring system, calves were classified as either clinically affected with BRD (n = 35) or healthy (n = 105). Additive genotypic tests were performed for genomic regions previously associated with susceptibility to BRD in calves from California (CA) and New Mexico (NM). Using this method, 4 loci (P < 0.01) consisting of 10 SNP were validated in the WI population, including 2 loci from CA, 1 locus from NM, and 1 locus from a combined CA + NM population. Most of the positional candidate genes and transcription factor binding site motifs associated with these loci have functions related to innate and adaptive immune responses. The validation of loci associated with susceptibility to BRD in independent populations allows producers to more reliably select cattle that are less susceptible to BRD, improving animal welfare, decreasing the annual revenue losses, and lowering the prevalence of the disease.

Genetic variation and gene conversions within the bovine NK-lysin gene family.

In contrast to a single copy of the NK-lysin gene in humans and many other mammals, we previously identified a family of four expressed NK-lysin genes arising by tandem duplications on cattle chromosome 11. Here, we report two genetic variants in the bovine NK-lysin complex with potential importance in the bovine innate immune system. The first one is a 9-bp deletion causing a three-amino-acid deletion in the pro-region of the NK1 gene product. The second is a deletion of NK2B in some Holstein cattle, resulting in copy number variation that is in disequilibrium with a SNP from the bovine 770K HD SNP array. We also show evidence for gene conversions within the three new NK2 genes, which at least partially accounts for their high degree of sequence identity.

Haplotype structure and copy number polymorphism of the beta-defensin 7 genes in diverse chicken breeds.

Beta-defensins is a family of avian peptides related to the innate immune system. Copy number variation was recently reported for the avian beta-defensin 7 gene (AvBD7) between the highly inbred Leghorn and Fayoumi lines. Here, we examined copy number variants in 35 different chicken breeds and found that 31 of them have at least the same representation of the duplicated AvBD7 allele. We also found haplotypes upstream of the AvBD6 regions that are strongly linked to the AvBD7 duplication. We observed a strong linkage disequilibrium spanning of the upstream region of the AvBD6 gene, with two SNPs being flanking markers to detect duplication of the AvBD7.

Description of bovine major histocompatibility complex class IIa haplotypes using parthenogenetic embryo-derived cells.

In this study, we report an approach to characterize individual BoLA haplotypes using cells from parthenogenetic bovine embryos derived from slaughterhouse ovaries. Eight of the 15 parthenogenetic embryos so obtained had not undergone meiotic recombination on the BoLA region and were suitable to describe BoLA haplotypes. Detailed analysis of the BoLA class IIa region identified seven different class IIa haplotypes, including six not previously described and two new alleles of BoLA-DQA and one BoLA-DQB. Our method provided reliable sources of homozygous DNA to describe BoLA haplotypes.

A high-resolution radiation hybrid map of the river buffalo major histocompatibility complex and comparison with BoLA.

The major histocompatibility complex (MHC) in mammals codes for antigen-presenting proteins. For this reason, the MHC is of great importance for immune function and animal health. Previous studies revealed this gene-dense and polymorphic region in river buffalo to be on the short arm of chromosome 2, which is homologous to cattle chromosome 23. Using cattle-derived STS markers and a river buffalo radiation hybrid (RH) panel (BBURH5000 ), we generated a high-resolution RH map of the river buffalo MHC region. The buffalo MHC RH map (cR5000 ) was aligned with the cattle MHC RH map (cR12000 ) to compare gene order. The buffalo MHC had similar organization to the cattle MHC, with class II genes distributed in two segments, class IIa and class IIb. Class IIa was closely associated with the class I and class III regions, and class IIb was a separate cluster. A total of 53 markers were distributed into two linkage groups based on a two-point LOD score threshold of ≥8. The first linkage group included 32 markers from class IIa, class I and class III. The second linkage group included 21 markers from class IIb. Bacterial artificial chromosome clones for seven loci were mapped by fluorescence in situ hybridization on metaphase chromosomes using single- and double-color hybridizations. The order of cytogenetically mapped markers in the region corroborated the physical order of markers obtained from the RH map and served as anchor points to align and orient the linkage groups.

Biallelic expression of the L-arginine:glycine amidinotransferase gene with different methylation status between male and female primordial germ cells in chickens.

The basic functions of DNA methylation include in gene silencing by methylation of specific gene promoters, defense of the host genome from retrovirus, and transcriptional suppression of transgenes. In addition, genomic imprinting, by which certain genes are expressed in a parent-of-origin-specific manner, has been observed in a wide range of plants and animals and has been associated with differential methylation. However, imprinting phenomena of DNA methylation effects have not been revealed in chickens. To analyze whether genomic imprinting occurs in chickens, methyl-DNA immunoprecipitation array analysis was applied across the entire genome of germ cells in early chick embryos. A differentially methylated region (DMR) was detected in the eighth intron of the l-arginine:glycine amidinotransferase (GATM) gene. When the DMR in GATM was analyzed by bisulfite sequencing, the methylation in male primordial germ cells (PGC) of 6-d-old embryos was higher than that in female PGC (57.5 vs. 35.0%). At 8 d, the DMR methylation of GATM in male PGC was 3.7-fold higher than that in female PGC (65.0 vs. 17.5%). Subsequently, to investigate mono- or biallelic expression of the GATM gene during embryo development, we found 2 indel sequences (GTTTAATGC and CAAAAA) within the GATM 3'-untranslated region in Korean Oge (KO) and White Leghorn (WL) chickens. When individual WL and KO chickens were genotyped for indel sequences, 3 allele combinations (homozygous insertion, homozygous deletion, and heterozygotes) were detected in both breeds using a gel shift assay and high-resolution melt assay. The deletion allele was predominant in KO, whereas the insertion allele was predominant in WL. Heterozygous animals were evenly distributed in both breeds (P < 0.01). Despite the different methylation status between male and female PGC, the GATM gene conclusively displayed biallelic expression in PGC as well as somatic embryonic, extraembryonic, and adult chicken tissues.

Comparative anchor tagged sequences (CATS) for integrative mapping of mammalian genomes.

Precise comparisons of mammalian gene maps require common anchor loci as landmarks for conserved chromosomal segments. Using a computer script that automates DNA sequence database alignments, we designed 410 evolutionarily conserved primer pair sequences which are specific for anchor locus gene amplification from any mammalian species' DNA. Primer pairs were designed to span introns for polymorphism ascertainment, and to include sufficient exonic sequence (25-400 bp) to allow for gene identification. A total of 318 primer pairs were optimized for domestic cats, and 86% of the sequenced feline PCR products showed homology to the gene of primer origin. A screen of 20 mammals from 11 orders revealed that 35-52% of the 318 primers yielded a single PCR product without further optimization suggesting that nearly 75% can be optimized for any eutherian mammal.

Anchored reference loci for comparative genome mapping in mammals.

Recent advances in gene mapping technologies have led to increased emphasis in developing representative genetic maps for several species, particularly domestic plants and animals. These maps are being compiled with two distinct goals: to provide a resource for genetic analysis, and to help dissect the evolution of genome organization by comparing linkage relationships of homologous genes. We propose here a list of 321 reference anchor loci suitable for comparative gene mapping in mammals and other vertebrate classes. We selected cloned mouse and human functional genes spaced an average of 5-10 centiMorgans throughout their respective genomes. We also attempted to include loci that are evolutionarily conserved and represented in comparative gene maps in other mammalian orders, particularly cattle and the domestic cat. We believe that the map may provide the basis for a unified approach to comparative analysis of mammalian species genomes.

Polymorphism and haplotype structure in River Buffalo (Bubalus bubalis) toll-like receptor 5 (TLR5) coding sequence.

Most of the 160 million river buffalo in the world are in Asia where they are used extensively, both as a food source and for draught power. Only recently have investigations begun exploring the buffalo genome for variation that might influence health and productivity of these economically important animals. This paper describes the sequence variability of the toll-like receptor 5 (TLR5) gene, which recognizes bacterial flagellin and is a key player in the immune system. TLR5 is comprised of a single exon that is 2577 bp and codes 858 amino acids. We examined single-nucleotide polymorphisms (SNPs) located within the coding region. Overall, 17 SNPs were discovered, seven of which are non-synonymous. Our study population yielded four different haplotypes. We examined predicted protein domain structure and found that river buffalo, swamp buffalo, and African Forest buffalo shared the same protein domain structure and are more similar to each other than they are to cattle and American bison, which are similar to each other. PolyPhen 2 analysis revealed one amino acid substitution in the river buffalo population with potential functional significance.

Construction and preliminary characterization of a river buffalo bacterial artificial chromosome library.

River buffalo genome analyses have advanced significantly in the last decade, and the genome sequence of Bubalus bubalis will be available shortly. Nonetheless, large-insert DNA library resources such as bacterial artificial chromosomes (BAC) are still required for validation and accurate assembly of the genome sequence. We constructed a river buffalo BAC library containing 52,224 clones with an average insert size of 97 kb, representing 1.7 × coverage of the genome. This genomic resource for river buffalo will facilitate further studies in this economically important species allowing for instance, whole genome physical mapping and isolation of genes and gene clusters, contributing to the elucidation of gene organization and identification of regulatory elements.

Polymorphism and gene organization of water buffalo MHC-DQB genes show homology to the BoLA DQB region.

In cattle (Bos taurus), there is evidence of more than 50 alleles of BoLA-DQB (bovine lymphocyte antigen DQB) that are distributed across at least five DQB loci, making this region one of the most complex in the BoLA gene family. In this study, DQB alleles were analysed for the water buffalo (Bubalus bubalis), another economically important bovine species. Twelve alleles for Bubu-DQB (Bubalis bubalis DQB) were determined by nucleotide sequence analysis. A phylogenetic analysis revealed numerous trans-species polymorphisms, with alleles from water buffalo assigned to at least three different loci (BoLA-DQB1, BoLA-DQB3 and BoLA-DQB4) that are also found in cattle. These presumptive loci were analysed for patterns of synonymous (d(S)) and non-synonymous (d(N)) substitution. Like BoLA-DQB1, Bubu-DQB1 was observed to be under strong positive selection for polymorphism. We conclude that water buffalo and cattle share the current arrangement of their DQB region because of their common ancestry.

Evidence for the horizontal acquisition of murine AKR virogenes by recent horizontal infection of the germ line.

Several recent reports (8, 10, 11, 13) have established the biological and molecular genetic similarity between the endogenous AKV virus of strain AKR, and an N-ecotropic endogenous virus found in the genome of feral Japanese mice, Mus musculus molossinus. The similarities are so striking as to suggest a common origin of these viruses, which are present in some, but not all, inbred mouse strains. The virogenes of AKR mice may have been acquired by either: (a) common descent of AKR (and other AKV(+) strains) from a common ancestor of AKR and molossinus animals, or (b) horizontal germ line infection of the AKR strains by molossinus virus at 1;he strain's inception followed by fixation through inbreeding. The sexual descent model carries with it a prediction of relative consanguinity of the AKR strain and molossinus, whereas the horizontal infection model does not. We have examined the polymorphic allozyme (allelic isozyme) genotype of 51 nonvirus-related loci in 17 strains of mice including AKR, C58, BALB/c, Swiss, and molossinus. By comparing the composite allozyme genotype of different inbred and outbred mouse strains, the "genetic distance" statistic was derived. Genetic distance measures the degree of allelic substitution between populations and increases proportionately with the amount of time the populations have been reproductively isolated. The genetic distance computed between molossinus and AKR is large, nearly 5-10 times the distance between known related populations and strains (e.g., C57L vs. C57BL/6). Molossinus had a similarly large distance from AKV negative strains (Swiss, C57L) as it did from AKV- positive strains. Cellular DNA sequences that flank the integrated AKV provirus were analyzed by restriction enzyme digestion of liver DNA from molossinus, AKR, and additional inbred strains that express ecotropic murine leukemia virus. The integration flanks of three AKR provirus sequences, Akv-1, Akv-2, and a third uncharacterized sequence, were not evident in molossinus cell DNA, which contained at least six different proviral integration fragments. These data effectively exclude the interpretation of consanguinity of AKR and molossinus and support the notion of acquisition of the endogenous virus in AKR by horizontal infection of the molossinus virus.

A high-resolution radiation hybrid map of sheep chromosome X and comparison with human and cattle.

A radiation hybrid (RH) map of sheep X chromosome (Ovisaries; OARX) containing 146 physically anchored loci was generated in this study, providing information for comparative X chromosome analysis between the maps of sheep, human, and cattle. Primers typed on the USUoRH5000 ovine whole-genome radiation hybrid panel were designed from sequences predicted to be on the ovine X chromosome, based on comparative mapping within the virtual sheep genome browser (v1.2). The resulting RH map for the ovine X chromosome consists of 4 linkage groups composed of 76 BAC end sequences (BES), 28 gene loci that were confirmed within ovine BAC clones in the CHORI-243 ovine BAC library, 28 additional gene loci from the ovine comparative map and 14 polymorphic sequence tagged sites (STS) from the OARX linkage map. This first-generation RH map of OARX contributes to the expansion of a comprehensive ovine genome map for sheep and provides evidence of rearrangements in loci order compared to the human and cattle orders.

Comparative RH maps of the river buffalo and bovine Y chromosomes.

Radiation hybrid maps were constructed for river buffalo and cattle Y chromosomes. A total of 41 cattle-derived Y-chromosome molecular markers were selected and tested with 2 previously described 5,000-rad whole-genome radiation hybrid (RH) panels (river buffalo - BBURH(5000) and cattle - BTARH(5000)) for generation of maps. Among the initial 41 selected markers, a subset of 26 markers generated PCR products suitable for scoring with the BBURH(5000) panel. Of these, 19 markers (73%) were distributed in 1 linkage group spanning 341.3 cR. Retention frequencies (RF) for individual markers ranged from 17.8% for SMCY to 56.7% for BTY1, with an average RF of 37.6%. From the selected markers, 37 generated reliable scores using the BTARH(5000) panel. The newly constructed BTAY RH map contains 28 markers distributed within 1 linkage group. Twenty-four of these markers had been previously mapped on BTAY using a 7,000-rad cattle-hamster WG-RH panel and 4 markers were mapped for the first time (ZFY, SeqRep, RepSeqS4 and BTY1). The length of the BTAY RH map was estimated to be 602.4 cR. Retention frequencies for individual mapped markers ranged from 10% (INRA126) to 63.3% (SeqRep), with an average RF of 35.3%. RH marker positions along the Y chromosome were compared between BBUY and BTAY, which revealed differences in the order of some of the markers. The BBUY pseudoautosomal region (PAR) is delineated by 3 BTAY PAR markers (MAF45, TGLA325 and UMN2008). These markers are telomeric in both species but are not found in the same order. Here we have demonstrated the effective use of bovine Y chromosome markers for the development of the first BBUY RH map. Likewise, these set of markers can be used for comparative assessment of Y chromosomes in other members of the Bovidae family.

Gene for neuraminidase activity on mouse chromosome 17 near h-2: pleiotropic effects on multiple hydrolases.

The low activity of liver neuraminidase that is characteristic of mouse strain SM/J is inherited as a single gene on chromosome 17, near the major histocompatibility complex. This gene, neuraminidase-1 (Neu-1), is represented by the low activity allele Neu-1s in SM/J and the high activity allele Neu-1b in C57BL/6J and most other strains. Previously described variations in the posttranslational processing of acid phosphatase, alpha-mannosidase, arylsulfatase-B, and alpha-glucosidase are attributed to pleiotropic effects of this gene.

The promise of comparative genomics in mammals.

Dense genetic maps of human, mouse, and rat genomes that are based on coding genes and on microsatellite and single-nucleotide polymorphism markers have been complemented by precise gene homolog alignment with moderate-resolution maps of livestock, companion animals, and additional mammal species. Comparative genetic assessment expands the utility of these maps in gene discovery, in functional genomics, and in tracking the evolutionary forces that sculpted the genome organization of modern mammalian species.

Analysis of sequence variability and protein domain architectures for bovine peptidoglycan recognition protein 1 and Toll-like receptors 2 and 6.

The mammalian Toll-like receptors (TLRs) recognize invading pathogens, thereafter provoking innate immune responses, whereas peptidoglycan recognition protein 1 (PGLYRP1) is directly microbicidal. The primary objective of this study was to characterize single-nucleotide polymorphisms (SNPs) and insertion-deletion polymorphisms (indels) within bovine TLR2, TLR6, and PGLYRP1, thereby facilitating future TLR signaling, association, and PGLYRP1 microbicidal assays relevant to bovine innate immunity. Comparative sequence analysis for 10 bovine breeds revealed 83 polymorphisms (82 SNPs, 1 indel), with 15 nonsynonymous SNPs located within predicted functional domains. Of the 83 polymorphisms detected, 72 (87%) are reported here for the first time. Several predicted amino acid replacements encoded by bovine TLR2 and TLR6, but not PGLYRP1, resulted in the confident prediction of protein domain alterations. Prediction and comparison of protein domain architectures for TLR2 and TLR6 revealed six regions of leucine-rich-repeat patterning that was conserved among multiple species. Collectively, differences in the patterns and frequencies of polymorphism were noted between bovine TLRs that predominantly recognize viral ligands (TLRs 3, 7, 8, 9) and those that recognize microbial and/or unknown ligands (TLRs 1, 2, 5, 6, 10).

Bovine genome mapping: evolutionary inference and the power of comparative genomics.

All bovine chromosomes are now represented by a syntenic group, a linkage map and at least one in situ hybridization. Almost 1,000 loci are mapped, about 300 of which are coding sequences useful for comparative mapping and evolutionary inference. Economically important loci are beginning to appear on bovine linkage maps and enhanced comparative maps are likely to be necessary to identify many of these genes by a comparative positional candidate gene approach.

First radiation hybrid map of the river buffalo X chromosome (BBUX) and comparison with BTAX.

We report the first radiation hybrid map of the river buffalo X chromosome generated from a recently constructed river buffalo (Bubalus bubalis) whole-genome radiation hybrid panel (BBURH(5000)). This map contains a total of 33 cattle-derived markers, including 10 genes, four ESTs and 19 microsatellites. The markers are distributed in two linkage groups: LG1 contains eight markers spanning 125.6 cR, and LG2 contains 25 markers spanning 366.3 cR. LG1 contains six markers in common with bovine sequence assembly build 3.1. With the exception of BMS2152, the order of these markers on our BBUX map is shuffled when compared to the cow X chromosome (Bos taurus; BTAX). From LG2, two markers (AMELX and BL22) map to a more distal portion of BTAX compared to BBUX. In addition, two pairs of LG2 markers exhibit inversions compared to BTAX (ILSTS017 and ATRX; XBM38 and PPEF1). Alternatively, when compared to the most recent bovine RH map (Bov-Gen 3000rads), BL1098 and BMS2227 from LG1 as well as PLS3 and BMS1820 from LG2 showed inverted positions on the BBUX map. These discrepancies in buffalo and cattle maps may reflect evolutionary divergence of the chromosomes or mapping errors in one of the two species. Although the set of mapped markers does not cover the entire X chromosome, this map is a starting point for the construction of a high-resolution map, which is necessary for characterization of small rearrangements that might have occurred between the Bubalus bubalis and Bos taurus X chromosomes.