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.

Rapid Communication: Subclinical bovine respiratory disease - loci and pathogens associated with lung lesions in feedlot cattle.

Bovine respiratory disease (BRD) is an economically important disease of feedlot cattle that is caused by viral and bacterial pathogen members of the BRD complex. Many cases of subclinical BRD go untreated and are not detected until slaughter, when lung lesions are identified. The objectives of this study were to identify which BRD pathogens were associated with the presence of lung lesions at harvest and to identify genomic loci that were associated with susceptibility to lung lesions as defined by consolidation of the lung and/or the presence of fibrin tissue. Steers from a Colorado feedlot ( = 920) were tested for the presence of viral and bacterial pathogens using deep pharyngeal and mid-nasal swabs collected on entry into the study. Pathogen profiles were compared between cattle with or without lung consolidation (LC), fibrin tissue in the lung (FT), a combination of LC and FT in the same lung (lung lesions [LL]), and hyperinflated lungs (HIF) at harvest. Genotyping was conducted using the Illumina BovineHD BeadChip. Genomewide association analyses (GWAA) were conducted using EMMAX (efficient mixed-model association eXpedited), and pseudoheritabilities were estimated. The pathogen profile comparisons revealed that LC ( = 0.01, odds ratio [OR] = 3.37) and LL cattle ( = 0.04, OR = 4.58) were more likely to be infected with bovine herpes virus-1 and that HIF cattle were more likely to be infected with spp. ( = 0.04, OR = 4.33). Pseudoheritability estimates were 0.25 for LC, 0.00 for FT, 0.28 for LL, and 0.13 for HIF. Because pseudoheritability for FT was estimated to be 0, GWAA results for FT were not reported. There were 4 QTL that were moderately associated ( < 1 × 10) with only LC, 2 that were associated with only LL, and 1 that was associated with LC and LL. Loci associated with HIF included 12 that were moderately associated and 3 that were strongly associated (uncorrected P < 5 × 10-7). A 24-kb region surrounding significant lead SNP was investigated to identify positional candidate genes. Many positional candidate genes underlying or flanking the detected QTL have been associated with signal transduction, cell adhesion, or gap junctions, which have functional relevance to the maintenance of lung health. The identification of pathogens and QTL associated with the presence of lung abnormalities in cattle exhibiting subclinical BRD allows the identification of loci that may not be detected through manifestation of clinical disease alone.

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.

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.

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.

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.

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.

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 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.

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.

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.

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.

Polymorphism identification, RH mapping and association of placental lactogen gene with milk production traits of dairy cows.

Bovine placental lactogen (bPL) is structurally related to prolactin (PRL) and growth hormone (GH). In synergism with steroid and thyroid hormones, bPL is crucial in stimulating the development of the mammary gland, mammary cell differentiation and function. To further explore whether bPL gene is associated with milk production traits, we herein analyzed single-nucleotide polymorphisms (SNPs) within eight regions of bPL gene, which are potentially associated with five milk production traits on 1028 Chinese Holstein cows. Among these, two SNPs, NT7409(T-C) and Nt11246(G-A), were identified. The former is within exon 2; it induces an alteration of amino acid from Val to Ala. The later is within exon 4. It is a synonymous mutation. We found that there were significant associations between NT7409(T-C) and milk and protein yield. Cows of the AA genotype yielded less milk (P = 0.001) and less protein (P = 0.003) than those of genotypes AB and BB. However, on the NT11246(G-A) locus, no significant association was observed in the five milk production traits studied. In addition, bPL has been localized near markers RM185 and CC549051 with a distance of 23.2 cR on BTA 23. It is at the same position as the region including quantitative trait loci (QTLs) affecting milk and protein yields by previous linkage analysis. In summary, our findings demonstrated that the SNP within exon 2 of bPL (NT7409(T-C)) is associated with two milk production traits, and this provided further evidence that bPL could be a major gene-controlling milk production trait in Holstein dairy cattle.

Mapping MHC genes in river buffalo.

The major histocompatibility complex (MHC) contains a set of genes necessary for antigen presentation in the immune system. This gene dense and polymorphic region of the mammalian genome is of considerable interest due to the role of MHC genes in immune function and animal health. Previous cytogenetic studies have indicated that the MHC in river buffalo resides on the short arm of chromosome 2 (BBU2). A 5000-rad radiation hybrid mapping panel was recently generated to enable construction of a whole genome map of river buffalo. To this end, the aims of this project were to elucidate the general organization of the MHC on BBU2, and to compare gene order within this region to the MHC in cattle. PCR primers were selected from the bovine gene map and used with the BBURH5000 panel to map a set of ten MHC class II genes in river buffalo. Analysis indicates that these genes fall into two linkage groups, consistent with organization of the MHC in cattle. This comparison of buffalo and bovine MHC gene order provides the first insight into the organization of the MHC on river buffalo chromosome 2.

Mapping the BDA20 gene to BTAX.

The lipocalin family is a large group of proteins that exhibits great structural and functional variation both within and among species, including a significant number of animal-derived aeroallergens, such as the bovine BDA20 (major cow dander allergen). This protein is classified as an occupational allergen causing asthma and other work-related allergic disorders among dairy farmers. Using a somatic cell panel the BDA20 gene was assigned to the bovine X chromosome (BTAX) with a significant concordant value of 97% to the previously mapped reference marker MAF45. A radiation hybrid (RH) mapping approach confirmed the assignment of BDA20 to BTAX. Two-point LOD scores showed that BDA20 is linked to XBM451 with a LOD score of 22.1 for a theta value of 0.03.

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).

A high-resolution comparative radiation hybrid map of ovine chromosomal regions that are homologous to human chromosome 6 (HSA6).

In this study, we constructed high-resolution radiation hybrid (RH) and comparative maps of ovine chromosomes or chromosomal segments that are homologous to human chromosome 6 (HSA6). A total of 251 markers were successfully genotyped across the recently developed USUoRH5000 whole-genome panel; 208 of these markers were assigned to five RH linkage groups distributed on three ovine chromosomes (OAR8, 9 and 20). The RH maps have good correspondence with previous chromosome painting data, although a small centromeric region on OAR9 that is homologous to HSA6 had not been previously detected using human chromosome paints on ovine chromosomal spreads. High percentages of the ovine markers were identified as orthologues in the bovine (86.3%), dog (85.8%), horse (69.3%) and human (88.7%) genomes. These maps contribute to investigations in mammalian chromosome evolution and the search for economic trait loci in sheep.

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.

A radiation hybrid map of river buffalo (Bubalus bubalis) chromosome 1 (BBU1).

The largest chromosome in the river buffalo karyotype, BBU1, is a submetacentric chromosome with reported homology between BBU1q and bovine chromosome 1 and between BBU1p and BTA27. We present the first radiation hybrid map of this chromosome containing 69 cattle derived markers including 48 coding genes, 17 microsatellites and four ESTs distributed in two linkage groups spanning a total length of 1330.1 cR(5000). The RH map was constructed based on analysis of a recently developed river buffalo-hamster whole genome radiation hybrid (BBURH(5000)) panel. The retention frequency of individual markers across the panel ranged from 17.8 to 52.2%. With few exceptions, the order of markers within linkage groups is identical to the order established for corresponding cattle RH maps. The BBU1 map provides a starting point for comparison of gene order rearrangements between river buffalo chromosome 1 and its bovine homologs.