Loci associated with conception rate in crossbred beef heifers.

The inability of beef cattle to maintain full term pregnancies has become an economic concern for the beef industry. Herd management and nutritional improvements have alleviated environmental impacts on embryonic and fetal loss, yet additional gains can be made through genomic selection. The objectives of this study were to identify loci and gene-sets in crossbred beef heifers associated with the number of services required to become pregnant (TBRD) and heifer conception rate at first service (HCR1). Heifers (n = 709) from a commercial beef operation underwent one round of artificial insemination, before exposure to bulls for natural service for 50 days. Pregnancy and time of conception was determined by ultrasound 35 days after the breeding season. Heifers were genotyped using the GeneSeek (Lincoln, NE) Bovine GGP50K BeadChip prior to genome-wide association analyses (GWAA) conducted using an EIGENSTRAT-like model to identify loci associated (P < 1 × 10-5) with TBRD and HCR1. One locus was associated (P = 8.97 × 10-6) with TBRD on BTA19 and included the positional candidate gene ASIC2, which is differentially expressed in the endometrium of fertility classified heifers, and the positional candidate gene, SPACA3. Gene-set enrichment analyses using SNP (GSEA-SNP) data, was performed and identified one gene-set, oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen as enriched (NES = 3.15) with TBRD and contained nine leading edge genes that contributed to the enrichment of the gene set. The enriched gene-set is involved in catalyzing oxidation-reduction reactions, which have been associated with oxidative stressors impacting pregnancy success. No loci were associated nor gene-sets enriched with HCR1. Identification of loci, positional candidate genes, gene-sets and leading edge genes enriched for fertility facilitate genomic selection that allows producers to select for reproductively superior cattle, reduce costs associated with infertility, and increase percent calf crop.

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.

Genome-wide association analysis and gene set enrichment analysis with SNP data identify genes associated with 305-day milk yield in Holstein dairy cows.

Milk production traits, such as 305-day milk yield (305MY), have been under direct selection to improve production in dairy cows. Over the past 50 years, the average milk yield has nearly doubled, and over 56% of the increase is attributable to genetic improvement. As such, additional improvements in milk yield are still possible as new loci are identified. The objectives of this study were to detect SNPs and gene sets associated with 305MY in order to identify new candidate genes contributing to variation in milk production. A population of 781 primiparous Holstein cows from six central Washington dairies with records of 305MY and energy corrected milk were used to perform a genome-wide association analysis (GWAA) using the Illumina BovineHD BeadChip (777 962 SNPs) to identify QTL associated with 305MY (P < 1.0 × 10-5 ). A gene set enrichment analysis with SNP data (GSEA-SNP) was performed to identify gene sets (normalized enrichment score > 3.0) and leading edge genes (LEGs) influencing 305MY. The GWAA identified three QTL comprising 34 SNPs and 30 positional candidate genes. In the GSEA-SNP, five gene sets with 58 unique and 24 shared LEGs contributed to 305MY. Identification of QTL and LEGs associated with 305MY can provide additional targets for genomic selection to continue to improve 305MY in dairy cattle.

Gene set enrichment analysis of SNP data in dairy and beef cattle with bovine respiratory disease.

Bovine respiratory disease (BRD) is a complex disease that is associated with infection by bacterial and viral pathogens when cattle fail to adequately respond to stress. The objective of this study was to use gene set enrichment analysis of SNP data (GSEA-SNP) and a network analysis (ingenuity pathway analysis) to identify gene sets, genes within gene sets (leading-edge genes) and upstream regulators associated with BRD in pre-weaned dairy calves and beef feedlot cattle. BRD cases and controls were diagnosed using the McGuirk health scoring system. Holstein calves were sampled from commercial calf-raising facilities in California (1003 cases and 1011 controls) and New Mexico (376 cases and 372 controls). Commercial feedlot cattle were sampled from Colorado (500 cases and 499 controls) and Washington (504 cases and 497 controls). There were 102 and 237 unique leading-edge genes identified in the dairy calf and beef cattle populations respectively. Six leading-edge genes (ADIPOQ, HTR2A, MIF, PDE6G, PRDX3 and SNCA) were associated with BRD in both dairy and beef cattle. Network analysis identified glucose as the most influential upstream regulator in dairy cattle, whereas in beef cattle, TNF was the most influential upstream regulator. The genes, gene sets and upstream regulators associated with BRD have common functions associated with immunity, inflammation and pulmonary disease and provide insights into the mechanisms that are critical to BRD susceptibility in cattle.

Identification of genes associated with susceptibility to Mycobacterium avium ssp. paratuberculosis (Map) tissue infection in Holstein cattle using gene set enrichment analysis-SNP.

Multiple genome-wide association analyses have investigated susceptibility to bovine paratuberculosis, but few loci have been identified across independent cattle populations. A SNP-based gene set enrichment analysis (GSEA-SNP) allows expanded identification of genes with moderate effects on a trait through the enrichment of gene sets instead of identifying only few loci with large effects. Therefore, the objective of this study was to identify genes that were moderately associated with Mycobacterium avium ssp. paratuberculosis (Map) tissue infection using GSEA-SNP in Holstein cattle from the Pacific Northwest (PNW; n = 205) and from the PNW and Northeast (PNW+NE; n = 245) which were previously genotyped with the Illumina BovineSNP50 BeadChip. The GSEA-SNP utilized 4389 gene sets from five databases. For each annotated gene in the UMD3.1 assembly (n = 19,723), the most significant SNP within each gene and its surrounding region (10 kb up- and downstream) was selected as a proxy for that gene. Any gene set with a normalized enrichment score > 2.5 was considered enriched. Thirteen gene sets (8 PNW GSEA-SNP; 5 PNW+NE) were enriched in these analyses and all have functions that relate to nuclear factor kappa beta. Nuclear factor kappa beta is critical to gut immune responses, implicated in host immune responses to other mycobacterial diseases, and has established roles in inflammation as well as cancer. Gene sets and genes moderately associated with Map infection could be used in genomic selection to allow producers to select for less susceptible cattle, lower the prevalence of the disease, and reduce economic losses.

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.

Identification of loci associated with susceptibility to subspecies () tissue infection in cattle.

Johne's disease is a contagious bacterial infection of cattle caused by ssp. (). A previous genome-wide association analysis (GWAA) in Holstein cattle identified QTL on BTA3 and BTA9 that were highly associated (P < 5 × 10) and on BTA1, BTA16, and BTA21 that were moderately associated (P < 5 × 10) with Map tissue infection. The objectives of this study were to validate previous GWAA results in Jersey cattle ( = 57), Holstein cattle from the Pacific Northwest (PNW, = 205) and a combined Holstein population from the PNW and the Northeast (PNW + NE, = 423), and also identify new loci associated with tissue infection. DNA was genotyped using the Illumina BovineSNP50 BeadChip, and the PNW + NE data was also imputed to whole genome sequence level using Run4 of the 1000 Bull Genomes project with Beagle v 4.1 and FImpute. Cases were ileocecal node positive and controls were negative for by quantitative PCR (qPCR). Individuals were removed for SNP call rate < 90%, and SNP were removed for genotype call rate < 90% or minor allele frequency < 1%. For the Jersey, PNW, and PNW + NE, GWAA were conducted using an allelic dosage model. For the PNW and the PNW + NE, an additional efficient mixed-model association eXpedited (EMMAX) analysis was performed using additive, dominance and recessive models. Seven QTL on BTA22 were identified in the Jersey population with the most significant ( = 4.45 × 10) located at 21.7 megabases (Mb). Six QTL were associated in the PNW and the PNW + NE analyses, including a QTL previously identified on BTA16 in the NE population. The most significant locus for the PNW was located on BTA21 at 61 Mb ( = 8.61 × 10) while the most significant locus for the PNW + NE was on BTA12 at 90 Mb ( = 2.33 × 10). No additional QTL were identified with the imputed GWAA. Putative positional candidate genes were identified within 50 kb 5' and 3' of each QTL. Two positional candidate genes were identified in Jersey cattle, 1 identified in the PNW and 8 in the PNW + NE populations. Many identified positional candidate genes are involved in signal transduction, have immunological functions, or have putative functional relevance in entry into host cells. This study supported 2 previously identified SNP within a QTL on BTA16 and identified 16 new QTL, including 2 found in the PNW and the PNW+NE, associated with tissue infection.

Estimation of inbreeding and effective population size of full-blood Wagyu cattle registered with the American Wagyu Cattle Association.

The objective of this research was to examine the population structure of full-blood (100%) Wagyu cattle registered in the United States with the American Wagyu Association, with the aim of estimating and comparing the levels of inbreeding from both pedigree and genotypic data. A total of 4132 full-blood Wagyu cattle pedigrees were assessed and used to compute the inbreeding coefficients (FIT and FST ) and the effective population size (Ne ) from pedigree data for the period 1994 to 2011. In addition to pedigree analysis, 47 full-blood Wagyu cattle representing eight prominent sire lines in the American Wagyu cattle population were genotyped using the Illumina BovineSNP50 BeadChip. Genotypic data were then used to estimate genomic inbreeding coefficients (FROH ) by calculating runs of homozygosity. The mean inbreeding coefficient based on the pedigree data was estimated at 4.80%. The effective population size averaged 17 between the years 1994 and 2011 with an increase of 42.9 in 2000 and a drop of 1.8 in 2011. Examination of the runs of homozygosity revealed that the 47 Wagyu cattle from the eight prominent sire lines had a mean genomic inbreeding coefficient (FROH ) estimated at 9.08% compared to a mean inbreeding coefficient based on pedigree data of 4.8%. These data suggest that the mean genotype inbreeding coefficient of full-blood Wagyu cattle exceeds the inbreeding coefficient identified by pedigree. Inbreeding has increased slowly at a rate of 0.03% per year over the past 17 years. Wagyu breeders should continue to utilize many sires from divergent lines and consider outcrossing to other breeds to enhance genetic diversity and minimize the adverse effects of inbreeding in Wagyu.

Invited review: Genomic analysis of data from physiological studies.

Physiology deals with the functions of living organisms and their systems, and its scientific endeavors can be viewed as having temporally occurred in 3 phases. The first phase of physiology studies focused on determining the functions of particular organs and tissues and their functional differences according to physiological status. The second phase of studies focused on characterizing differences in these functions according to the environment, or productivity. The third phase of studies focuses on determining the physiological causes of differences in productivity. Distinguishing cause from effect in physiological systems of inter-related processes is problematic, such that science has struggled to identify the root physiological mechanisms and their role in the network of genes leading to differences in productivity. Genomics is the study of the entire genome and provides powerful new tools that will accelerate third-phase discoveries of causal physiological processes. That research exploits information on DNA polymorphisms known as markers, complete DNA sequence, RNA sequence, and RNA expression in particular tissues at specific life stages. Physiologists can determine the genetic cause of mutant animals, identify genetic differences between cases and controls, and identify genes responsible for differences in performance between average and above-average animals. In some species, physiologists can leverage genomic data being used to predict genetic merit in elite seedstock populations, as a starting point to identify genes that will then motivate detailed physiological studies in the organs or tissues and stages of life in which those genes are expressed. Such work will increase our knowledge of biology and may lead to novel approaches to manipulate animal performance.

Alpharma Beef Cattle Nutrition Symposium: nutrition and the genome.

It has long been appreciated that animals fed the same diet may perform differently. This is due to the ability of nutrients to interact with and affect molecular pathways that result in differences in BW gain, production performance, or disease resistance. To understand these effects, studies are being undertaken to discover how the differential expression and function of genes occur with different diets. These studies are using new technologies, genomic resources, and analysis techniques that have recently become available for domestic animals. Nutrigenomics and nutrigenetics are new research approaches that strive to optimize health by looking beyond the diet to understand the effects of food at the genetic and epigenetic levels. Nutrigenomics is focused on the effects of diet on health through an understanding of how bioactive chemicals in foods and supplements alter gene expression or the structure of the genome of an animal. Nutrigenetics focuses on how the genetic composition (i.e., genetic variation) of an animal influences their response to a given diet. Results from these studies will aid in formulating nutritionally appropriate diets that may be optimized for animals based on their genomic underpinnings. Nutrigenomics and nutrigenetics unite many fields: nutrition, bioinformatics, molecular biology, genomics, functional genomics, epidemiology, and epigenomics. The use of multi-disciplinary tools promises new opportunities to investigate the complex interactions of the genome and the diet of an animal. Through these new approaches, the partnerships of the genome and nutrition will be revealed resulting in improved efficiency of diets, enhanced sustainability of animals as a protein source, and improved methods for preventing illnesses.

Short communication: refinement of genetic regions associated with Mycobacterium avium subspecies paratuberculosis tissue infection and tolerance to Johne's disease.

Johne's disease is a highly transmissible bacterial disease caused by Mycobacterium avium ssp. paratuberculosis (MAP). The objective of this study was to refine the locus associated with MAP tissue infection and the locus associated with tolerance to Johne's disease. Using a genome-wide association analysis, single nucleotide polymorphisms associated with MAP tissue infection and tolerance to Johne's disease on Bos taurus autosome (BTA)3 and BTA15, respectively, have previously been identified. A 235-kb region on BTA3 was evaluated with 42 single nucleotide polymorphisms, and a 193-kb region on BTA15 was evaluated with 54 single nucleotide polymorphisms in a group of 209 Holstein cows. Using a single marker association analysis and haplotype tests, we refined a region of 10.6 kb on BTA3 as being associated with MAP tissue infection and a region of 6.5 kb on BTA15 as being associated with tolerance to Johne's disease.

Identification of loci associated with tolerance to Johne's disease in Holstein cattle.

Johne's disease, caused by Mycobacterium avium subspecies paratuberculosis (Map), is a fatal disease in cattle. The objective of this study was to identify loci associated with tolerance in cows infected with Map. Tolerance was defined as a cow's fitness at a given level of Map infection intensity. Fitness was measured by Map faecal cultures, and Map infection intensity was measured by culturing four gut tissues. The quantitative phenotype of tolerance was defined by numerical indexes of cultures of peak (peak tolerance, PT) and average (average tolerance, AT) faecal and tissue Map from 245 Holstein cows. The categorical phenotype was defined as: ≥ 100 cfu Map tissue infection, and faecal shedding ≥ 75 cfu (intolerant) or <10 cfu (tolerant cows). In 94 cows, Map was identified in ≥ 1 tissue, including 44 cows with ≥ 100 Map tissue cfu and 36 with ≥ 1 faecal cfu. A genome-wide association analysis was performed after filtering, leaving genotypes for 45,789 SNPs in 90 animals for the quantitative phenotype and 16 cases and 25 controls for the categorical analysis of tolerance. rs41748405:A>C (BTA15) was associated with PT (P = 1.12 × 10(-7)) and AT (P = 2.17 × 10(-6)). Associations were identified with PT and adjacent SNPs ss61512613:A>G and ss61530518:A>G (BTA6) (P < 3.0 × 10(-5)), and with AT for ss61469568:A>G (BTA 2) (P = 3.3 × 10(-5)) and ss86284768:A>G (BTA1) (P = 3.31 × 10(-5)). For the categorical phenotype, an association was found with ss8632653:A>G (BTA6) (P < 5.0 × 10(-5)). This is the first study to identify loci associated with tolerance to Johne's disease.

Polymorphic markers for the arylsulfatase A gene reveal a greatly expanded meiotic map for the human 22q telomeric region.

Two microsatellite markers, D22S1743 and D22S1744, were developed for the arylsulfatase A (ARSA) region of chromosome 22q. Linkage analysis for 171 families, using nine reference markers covering all of 22q, placed these new markers 2.0 Kosambi cM distal to D22S526, making them more distal than any microsatellite markers currently on the Généthon or Marshfield linkage maps. Recombination between proximal markers D22S270/D22S683 and D22S446/D22S311 exhibited increased rates of female meiotic recombination compared to male recombination (P < 0.01). In contrast, the region encompassing sJCW16, D22S526, D22S1743, and D22S1744 exhibited relatively greater recombination in males (1.1 cM for females and 7.5 cM for males; chi(2); P < 0.005). These four distal markers lie in a region of hyperrecombination having a sex-averaged recombination ratio of between 8.3 (D22S1843/D22S1744) and 12 cM (sJCW16/D22S526) per megabase.

Differentiation of F18ab+ from F18ac+ Escherichia coli by single-strand conformational polymorphism analysis of the major fimbrial subunit gene (fedA).

Toxin-producing Escherichia coli expressing F18 fimbriae colonizes the small intestines of weaned pigs and causes diarrhea, edema disease, or both. The F18 family is composed of two antigenic variants, F18ab and F18ac. Because many strains do not express F18 fimbriae in vitro, identification and differentiation of these two variants are difficult. Single-strand conformational polymorphism (SSCP) analysis is a rapid method for identifying genetic mutations and polymorphisms. The F18 major fimbrial subunit genes (fedA) of 138 strains were amplified by PCR, and genetic differences were detected by SSCP analysis. The SSCP analysis of the fedA gene differentiated F18ab+ strains from F18ac+ strains. Most strains classified as F18ab+ by SSCP analysis contained Shiga toxin 2e and enterotoxin genes. Most strains classified as F18ac+ by SSCP analysis contained only enterotoxin genes. The SSCP analysis was a useful method for predicting the antigenicity of F18+ E. coli and could also be used for analysis of other virulence genes in E. coli and other pathogenic bacteria.

The natural resistance-associated macrophage protein gene is associated with Crohn's disease.

BACKGROUND: Identification of the genes causing inflammatory bowel disease (IBD) would enhance the understanding of and the treatment options for this disease. A hyperreactive immune response toward the intestinal flora has been implicated in the pathology of IBD. The natural resistance-associated macrophage protein (NRAMP) gene is believed to regulate macrophage function, especially the ability to fight intracellular pathogens. Genetic differences of NRAMP might, therefore, be associated with IBD. METHODS: Two DNA markers (D2S434 and D2S1323) near NRAMP were polymerase chain reaction (PCR) amplified and genotyped with DNA from 103 patients with Crohn's disease, 85 patients with ulcerative colitis, and 98 control subjects. Clinical data were obtained for all patients. Comparisons were made by chi-squared analysis. Disease association with significant haplotypes was expressed as odds ratio. RESULTS: Allele and genotype distributions were similar for both markers among all groups. Haplotype frequencies were different among Crohn's disease and control groups (p = 0.024). Two individual haplotypes of the patients with Crohn's disease were significant compared with control subjects: DA (p = 0.023; odds ratio, 0.5; 95% confidence interval, 0.3 to 0.9) and EA (p = 0.001; odds ratio, 3.5; 95% confidence interval, 1.6 to 3.2). The haplotype distribution was different within three age-of-onset groups of patients with Crohn's disease (p = 0.05). CONCLUSIONS: This study is the first to report an association between the NRAMP gene and Crohn's disease.

Crohn's disease and ulcerative colitis are associated with the DNA repair gene MLH1.

OBJECTIVE: The purpose of this study was to determine whether there is an association between Crohn's disease and ulcerative colitis with MLH1. SUMMARY BACKGROUND DATA: Identification of genes involved in the etiology of inflammatory bowel disease may lead to the development of markers that objectively can define disease and permit therapy. The treatment of Crohn's disease of the colon and ulcerative colitis also is complicated by difficulties in differentiating the two conditions. METHODS: The DNA and clinical data were obtained on 126 unrelated individuals (45 Crohn's disease, 36 ulcerative colitis, and 45 control subjects without intestinal disease). Polymerase chain reaction products were analyzed by single-strand conformation polymorphisms (MLH1 exons 9, 11, 14, 15, and 16) and polyacrylamide gel electrophoresis (markers D3S1611 and D3S1768). All comparisons were analyzed by chi square test. The association between single haplotypes and disease was expressed as relative odds. RESULTS: MLH1 exons 9, 11, 14, and 16 were monomorphic. Two, four, and six alleles were detected in MLH1 exon 15, D3S1611, and D3S1768, respectively. Significant associations were observed for MLH1 exon 15/D3S1611 haplotypes AB (OR = 5.5; p = 0.007) and BA (p = 0.002) with Crohn's disease and for haplotypes AB (OR = 4.0; p = 0.042), BA (p = 0.035), and BC (OR = 6.1; p = 0.016) with ulcerative colitis. Family history of inflammatory bowel disease was associated with D3S1768/D3S1611 (p = 0.05) and MLH1 exon 15/D3S1611 haplotypes (p = 0.03). D3S1611/D3S1768 haplotype CD (OR = 11.3; p = 0.03) was associated with disease, whereas MLH1 exon 15/D3S1611 haplotype AA (OR = 0.25; p = 0.02) was protective. Comparisons of MLH1 exon 15/D3S1611 haplotypes of Crohn's colitis and patients with ulcerative colitis were significant (p = 0.037). CONCLUSIONS: This study identifies a novel genetic and clinical association between MLH1 and inflammatory bowel disease.