Identification of novel alternative splicing associated with mastitis disease in Holstein dairy cows using large gap read mapping.

BACKGROUND: Mastitis is a very common disease in the dairy industry that producers encounter daily. Transcriptomics, using RNA-Sequencing (RNA - Seq) technology, can be used to study the functional aspect of mastitis resistance to identify animals that have a better immune response to mastitis. When the cow has mastitis, not only genes but also specific mRNA isoforms generated via alternative splicing (AS) could be differentially expressed (DE), leading to the phenotypic variation observed. Therefore, the objective of this study was to use large gap read mapping to identify mRNA isoforms DE between healthy and mastitic milk somatic cell samples (N = 12). These mRNA isoforms were then categorized based on being 1) annotated mRNA isoforms for gene name and length, 2) annotated mRNA isoforms with different transcript length and 3) novel mRNA isoforms of non - annotated genes. RESULTS: Analysis identified 333 DE transcripts (with at least 2 mRNA isoforms annotated, with at least one being DE) between healthy and mastitic samples corresponding to 303 unique genes. Of these 333 DE transcripts between healthy and mastitic samples, 68 mRNA isoforms are annotated in the bovine genome reference (ARS.UCD.1.2), 249 mRNA isoforms had novel transcript lengths of known genes and 16 were novel transcript lengths of non - annotated genes in the bovine genome reference (ARS.UCD.1.2). Functional analysis including gene ontology, gene network and metabolic pathway analysis was performed on the list of 288 annotated and unique DE mRNA isoforms. In total, 67 significant metabolic pathways were identified including positive regulation of cytokine secretion and immune response. Additionally, numerous DE novel mRNA isoforms showed potential involvement with the immune system or mastitis. Lastly, QTL annotation analysis was performed on coding regions of the DE mRNA isoforms, identifying overlapping QTLs associated with clinical mastitis and somatic cell score. CONCLUSION: This study identified novel mRNA isoforms generated via AS that could lead to differences in the immune response of Holstein dairy cows and be potentially implemented in future breeding programs.

Symposium review: Advances in sequencing technology herald a new frontier in cattle genomics and genome-enabled selection.

The cattle reference genome assembly has underpinned major innovations in beef and dairy genetics through genome-enabled selection, including removal of deleterious recessive variants and selection for favorable alleles affecting quantitative production traits. The initial reference assemblies, up to and including UMD3.1 and Btau4.1, were based on a combination of clone-by-clone sequencing of bacterial artificial chromosome clones generated from blood DNA of a Hereford bull and whole-genome shotgun sequencing of blood DNA from his inbred daughter/granddaughter named L1 Dominette 01449 (Dominette). The approach introduced assembly gaps, misassemblies, and errors, and it limited the ability to assemble regions that undergo rearrangement in blood cells, such as immune gene clusters. Nonetheless, the reference supported the creation of genotyping tools and provided a basis for many studies of gene expression. Recently, long-read sequencing technologies have emerged that facilitated a re-assembly of the reference genome, using lung tissue from Dominette to resolve many of the problems and providing a bridge to place historical studies in common context. The new reference, ARS-UCD1.2, successfully assembled germline immune gene clusters and improved overall continuity (i.e., reduction of gaps and inversions) by over 250-fold. This reference properly places nearly all of the legacy genetic markers used for over a decade in the industry. In this review, we discuss the improvements made to the cattle reference; remaining issues present in the assembly; tools developed to support genome-based studies in beef and dairy cattle; and the emergence of newer genome assembly methods that are producing even higher-quality assemblies for other breeds of cattle at a fraction of the cost. The new frontier for cattle genomics research will likely include a transition from the individual Hereford reference genome, to a "pan-genome" reference, representing all the DNA segments existing in commonly used cattle breeds, bringing the cattle reference into line with the current direction of human genome research.

Genetic mechanisms regulating the host response during mastitis.

Mastitis is a very costly and common disease in the dairy industry. The study of the transcriptome from healthy and mastitic milk somatic cell samples using RNA-Sequencing technology can provide measurements of transcript levels associated with the immune response to the infection. The objective of this study was to characterize the Holstein milk somatic cell transcriptome from 6 cows to determine host response to intramammary infections. RNA-Sequencing was performed on 2 samples from each cow from 2 separate quarters, one classified as healthy (n = 6) and one as mastitic (n = 6). In total, 449 genes were differentially expressed between the healthy and mastitic quarters (false discovery rate <0.05, fold change >±2). Among the differentially expressed genes, the most expressed genes based on reads per kilobase per million mapped reads (RPKM) in the healthy group were associated with milk components (CSN2 and CSN3), and in the mastitic group they were associated with immunity (B2M and CD74). In silico functional analysis was performed using the list of 449 differentially expressed genes, which identified 36 significantly enriched metabolic pathways (false discovery rate <0.01), some of which were associated with the immune system, such as cytokine-cytokine interaction and cell adhesion molecules. Seven functional candidate genes were selected, based on the criteria of being highly differentially expressed between healthy and mastitic groups and significantly enriched in metabolic pathways that are relevant to the inflammatory process (GLYCAM1, B2M, CD74, BoLA-DRA, FCER1G, SDS, and NFKBIA). Last, we identified the differentially expressed genes that are located in quantitative trait locus regions previously known to be associated with mastitis, specifically clinical mastitis, somatic cell count, and somatic cell score. It was concluded that multiple genes within quantitative trait locus regions could potentially affect host response to mastitis-causing agents, making some cows more susceptible to intramammary infections. The identification of potential candidate genes with functional, statistical, biological, and positional relevance associated with host defense to infection will contribute to a better understanding of the underlying genetic architecture associated with mastitis. This in turn will improve the sustainability of agricultural practices by facilitating the selection of cows with improved host defense leading to increased resistance to mastitis.

Implementation of Bayesian methods to identify SNP and haplotype regions with transmission ratio distortion across the whole genome: TRDscan v.1.0.

Realized deviations from the expected Mendelian inheritance of alleles from heterozygous parents have been previously reported in a broad range of organisms (i.e., transmission ratio distortion; TRD). Various biological mechanisms affecting gametes, embryos, fetuses, or even postnatal offspring can produce patterns of TRD. However, knowledge about its prevalence and potential causes in livestock species is still scarce. Specific Bayesian models have been recently developed for the analyses of TRD for biallelic loci, which accommodated a wide range of population structures, enabling TRD investigation in livestock populations. The parameterization of these models is flexible and allows the study of overall (parent-unspecific) TRD and sire- and dam-specific TRD. This research aimed at deriving Bayesian models for fitting TRD on the basis of haplotypes, testing the models for both haplotype- and SNP-based methods in simulated data and actual Holstein genotypes, and developing a specific software for TRD analyses. Results obtained on simulated data sets showed that the statistical power of the analysis increased with sample size of trios (n), proportion of heterozygous parents, and the magnitude of the TRD. On the other hand, the statistical power to detect TRD decreased with the number of alleles at each loci. Bayesian analyses showed a strong Pearson correlation coefficient (≥0.97) between simulated and estimated TRD that reached the significance level of Bayes factor ≥10 for both single-marker and haplotype analyses when n ≥ 25. Moreover, the accuracy in terms of the mean absolute error decreased with the increase of the sample size and increased with the number of alleles at each loci. Using real data (55,732 genotypes of Holstein trios), SNP- and haplotype-based distortions were detected with overall TRD, sire-TRD, or dam-TRD, showing different magnitudes of TRD and statistical relevance. Additionally, the haplotype-based method showed more ability to capture TRD compared with individual SNP. To discard possible random TRD in real data, an approximate empirical null distribution of TRD was developed. The program TRDscan v.1.0 was written in Fortran 2008 language and provides a powerful statistical tool to scan for TRD regions across the whole genome. This developed program is freely available at http://www.casellas.info/files/TRDscan.zip.

Pre- and post-puberty expression of genes and proteins in the uterus of Bos indicus heifers: the luteal phase effect post-puberty.

Progesterone signaling and uterine function are crucial in terms of pregnancy establishment. To investigate how the uterine tissue and its secretion changes in relation to puberty, we sampled tissue and uterine fluid from six pre- and six post-pubertal Brahman heifers. Post-pubertal heifers were sampled in the luteal phase. Gene expression of the uterine tissue was investigated with RNA-sequencing, whereas the uterine fluid was used for protein profiling with mass spectrometry. A total of 4034 genes were differentially expressed (DE) at a nominal P-value of 0.05, and 26 genes were significantly DE after Bonferroni correction (P < 3.1 × 10-6 ). We also identified 79 proteins (out of 230 proteins) that were DE (P < 1 × 10-5 ) in the uterine fluid. When we compared proteomics and transcriptome results, four DE proteins were identified as being encoded by DE genes: OVGP1, GRP, CAP1 and HBA. Except for CAP1, the other three had lower expression post-puberty. The function of these four genes hypothetically related to preparation of the uterus for a potential pregnancy is discussed in the context of puberty. All DE genes and proteins were also used in pathway and ontology enrichment analyses to investigate overall function. The DE genes were enriched for terms related to ribosomal activity. Transcription factors that were deemed key regulators of DE genes are also reported. Transcription factors ZNF567, ZNF775, RELA, PIAS2, LHX4, SOX2, MEF2C, ZNF354C, HMG20A, TCF7L2, ZNF420, HIC1, GTF3A and two novel genes had the highest regulatory impact factor scores. These data can help to understand how puberty influences uterine function.

Global differential gene expression in the pituitary gland and the ovaries of pre- and postpubertal Brahman heifers.

To understand genes, pathways, and networks related to puberty, we characterized the transcriptome of two tissues: the pituitary gland and ovaries. Samples were harvested from pre- and postpubertal Brahman heifers (same age group). Brahman heifers () are older at puberty compared with , a productivity issue. With RNA sequencing, we identified differentially expressed (DEx) genes and important transcription factors (TF) and predicted coexpression networks. The number of DEx genes detected in the pituitary gland was 284 ( < 0.05), and was the most DEx gene (fold change = 4.12, = 0.01). The gene promotes bone mineralization through transforming growth factor-ß (TGFß) signaling. Further studies of the link between bone mineralization and puberty could target . In ovaries, 3,871 genes were DEx ( < 0.05). Four highly DEx genes were noteworthy for their function: (a γ-aminobutyric acid [GABA] transporter), (), and () and its receptor . These genes had higher ovarian expression in postpubertal heifers. The GABA and its receptors and transporters were expressed in the ovaries of many mammals, suggesting a role for this pathway beyond the brain. The pathway has been known to influence the timing of puberty in rats, via modulation of GnRH. The effects of at the hypothalamus, pituitary gland, and ovaries have been documented. and its receptors are known factors in the release of GnRH, similar to and GABA, although their roles in ovarian tissue are less clear. Pathways previously related to puberty such as TGFß signaling ( = 6.71 × 10), Wnt signaling ( = 4.1 × 10), and peroxisome proliferator-activated receptor (PPAR) signaling ( = 4.84 × 10) were enriched in our data set. Seven genes were identified as key TF in both tissues: , , , , , , and a novel gene. An ovarian subnetwork created with TF and significant ovarian DEx genes revealed five zinc fingers as regulators: , , , , and . Recent work of hypothalamic gene expression also pointed to zinc fingers as TF for bovine puberty. Although some zinc fingers may be ubiquitously expressed, the identification of DEx genes in common across tissues points to key regulators of puberty. The hypothalamus and pituitary gland had eight DEx genes in common. The hypothalamus and ovaries had 89 DEx genes in common. The pituitary gland and ovaries had 48 DEx genes in common. Our study confirmed the complexity of puberty and suggested further investigation on genes that code zinc fingers.

SNP detection using RNA-sequences of candidate genes associated with puberty in cattle.

Fertility traits, such as heifer pregnancy, are economically important in cattle production systems, and are therefore, used in genetic selection programs. The aim of this study was to identify single nucleotide polymorphisms (SNPs) using RNA-sequencing (RNA-Seq) data from ovary, uterus, endometrium, pituitary gland, hypothalamus, liver, longissimus dorsi muscle, and adipose tissue in 62 candidate genes associated with heifer puberty in cattle. RNA-Seq reads were assembled to the bovine reference genome (UMD 3.1.1) and analyzed in five cattle breeds; Brangus, Brahman, Nellore, Angus, and Holstein. Two approaches used the Brangus data for SNP discovery 1) pooling all samples, and 2) within each individual sample. These approaches revealed 1157 SNPs. These were compared with those identified in the pooled samples of the other breeds. Overall, 172 SNPs within 13 genes (CPNE5, FAM19A4, FOXN4, KLF1, LOC777593, MGC157266, NEBL, NRXN3, PEPT-1, PPP3CA, SCG5, TSG101, and TSHR) were concordant in the five breeds. Using Ensembl's Variant Effector Predictor, we determined that 12% of SNPs were in exons (71% synonymous, 29% nonsynonymous), 1% were in untranslated regions (UTRs), 86% were in introns, and 1% were in intergenic regions. Since these SNPs were discovered in RNA, the variants were predicted to be within exons or UTRs. Overall, 160 novel transcripts in 42 candidate genes and five novel genes overlapping five candidate genes were observed. In conclusion, 1157 SNPs were identified in 62 candidate genes associated with puberty in Brangus cattle, of which, 172 were concordant in the five cattle breeds. Novel transcripts and genes were also identified.

Transcriptome analyses identify five transcription factors differentially expressed in the hypothalamus of post- versus prepubertal Brahman heifers.

Puberty onset is a developmental process influenced by genetic determinants, environment, and nutrition. Mutations and regulatory gene networks constitute the molecular basis for the genetic determinants of puberty onset. The emerging knowledge of these genetic determinants presents opportunities for innovation in the breeding of early pubertal cattle. This paper presents new data on hypothalamic gene expression related to puberty in (Brahman) in age- and weight-matched heifers. Six postpubertal heifers were compared with 6 prepubertal heifers using whole-genome RNA sequencing methodology for quantification of global gene expression in the hypothalamus. Five transcription factors (TF) with potential regulatory roles in the hypothalamus were identified in this experiment: , , , , and . These TF genes were significantly differentially expressed in the hypothalamus of postpubertal versus prepubertal heifers and were also identified as significant according to the applied regulatory impact factor metric ( < 0.05). Two of these 5 TF, and , were zinc fingers, belonging to a gene family previously reported to have a central regulatory role in mammalian puberty. The gene belongs to the family of homologues of Drosophila sine oculis () genes implicated in transcriptional regulation of gonadotrope gene expression. Tumor-related genes such as and are known to affect basic cellular processes that are relevant in both cancer and developmental processes. Mutations in were associated with puberty in humans. Mutations in these TF, together with other genetic determinants previously discovered, could be used in genomic selection to predict the genetic merit of cattle (i.e., the likelihood of the offspring presenting earlier than average puberty for Brahman). Knowledge of key mutations involved in genetic traits is an advantage for genomic prediction because it can increase its accuracy.

Bovine and murine tissue expression of insulin like growth factor-I.

The genomic architecture and expression of the Igf-1 gene are complex yielding multiple IGF-I transcript isoforms with putative functional contributions to growth and metabolism. Using RNA-seq on different tissues, physiological states, and species, the breadth of transcripts expressed was determined. Tissues from pre- and post-pubertal heifers and mature mice were collected and the transcript isoforms characterized. Three different IGF-I isoforms were detected in heifers with Class 1 transcripts most abundantly expressed. The pituitary reduced IGF-I expression post-pubertally whereas the uterus increased expression. Murine IGF-I transcript expression was more diverse utilizing multiple exons, start sites, and 3'UTRs. The RNA-seq methodology to characterize expression profiles permits assessment of the transcript isoforms yielding insight into functional roles of each transcript.

RNA sequencing to study gene expression and single nucleotide polymorphism variation associated with citrate content in cow milk.

The technological properties of milk have significant importance for the dairy industry. Citrate, a normal constituent of milk, forms one of the main buffer systems that regulate the equilibrium between Ca(2+) and H(+) ions. Higher-than-normal citrate content is associated with poor coagulation properties of milk. To identify the genes responsible for the variation of citrate content in milk in dairy cattle, the metabolic steps involved in citrate and fatty acid synthesis pathways in ruminant mammary tissue using RNA sequencing were studied. Genetic markers that could influence milk citrate content in Holstein cows were used in a marker-trait association study to establish the relationship between 74 single nucleotide polymorphisms (SNP) in 20 candidate genes and citrate content in 250 Holstein cows. This analysis revealed 6 SNP in key metabolic pathway genes [isocitrate dehydrogenase 1 (NADP+), soluble (IDH1); pyruvate dehydrogenase (lipoamide) ß (PDHB); pyruvate kinase (PKM2); and solute carrier family 25 (mitochondrial carrier; citrate transporter), member 1 (SLC25A1)] significantly associated with increased milk citrate content. The amount of the phenotypic variation explained by the 6 SNP ranged from 10.1 to 13.7%. Also, genotype-combination analysis revealed the highest phenotypic variation was explained combining IDH1_23211, PDHB_5562, and SLC25A1_4446 genotypes. This specific genotype combination explained 21.3% of the phenotypic variation. The largest citrate associated effect was in the 3' untranslated region of the SLC25A1 gene, which is responsible for the transport of citrate across the mitochondrial inner membrane. This study provides an approach using RNA sequencing, metabolic pathway analysis, and association studies to identify genetic variation in functional target genes determining complex trait phenotypes.

Genetic variants in a lipid regulatory pathway as potential tools for improving the nutritional quality of grass-fed beef.

The aim of this study was to evaluate the effect of genetic variants on candidate genes corresponding to the sterol recognition element-binding protein-1 (SREBP-1) signaling pathway and stearoyl-CoA desaturases (SCD1 and SCD5) on muscle fatty acid (FA) composition of Brangus steers fattened on grass. FA profiles were measured on Longissimus lumborum muscle samples using a gas chromatography-flame ionization detection technique. A total of 43 tag single-nucleotide polymorphisms on the SCD1, SCD5, SREBP-1, SCAP, INSIG1, INSIG2, MBTPS1, MBTPS2, and SRPR genes were genotyped on 246 steers to perform a marker-trait association study. To evaluate the influence of the Indicine breed in the composite breed, additional groups of 48 Angus, 18 Hereford, 75 Hereford x Angus, and 36 Limousin x Hereford-Angus steers were also genotyped. To perform the association analysis, FA data were grouped according to the number of carbon atoms and/or number of double bonds (i.e. SFA, MUFA, PUFA, etc.). In addition, different indexes that reflect the activity of FA desaturase and elongase enzymes were calculated. SCD1 markers significantly affected C14:1/(C14:0 + C14:1) and C18:1/(C18:0 + C18:1) indexes, whereas one SNP in SCD5 was correlated with the C16:1/(C16:0 + C16:1) index. Polymorphisms in the signal recognition particle receptor (SRPR) gene were associated with all the estimated desaturase indexes. Because the evaluated markers showed no effect on total lipid content of beef, this work supports the potential utilization of these markers for the improvement of grass-fed beef without undesirable side effects.

Gene network analyses of first service conception in Brangus heifers: use of genome and trait associations, hypothalamic-transcriptome information, and transcription factors.

Measures of heifer fertility are economically relevant traits for beef production systems and knowledge of candidate genes could be incorporated into future genomic selection strategies. Ten traits related to growth and fertility were measured in 890 Brangus heifers (3/8 Brahman × 5/8 Angus, from 67 sires). These traits were: BW and hip height adjusted to 205 and 365 d of age, postweaning ADG, yearling assessment of carcass traits (i.e., back fat thickness, intramuscular fat, and LM area), as well as heifer pregnancy and first service conception (FSC). These fertility traits were collected from controlled breeding seasons initiated with estrous synchronization and AI targeting heifers to calve by 24 mo of age. The BovineSNP50 BeadChip was used to ascertain 53,692 SNP genotypes for ∼802 heifers. Associations of genotypes and phenotypes were performed and SNP effects were estimated for each trait. Minimally associated SNP (P < 0.05) and their effects across the 10 traits formed the basis for an association weight matrix and its derived gene network related to FSC (57.3% success and heritability = 0.06 ± 0.05). These analyses yielded 1,555 important SNP, which inferred genes linked by 113,873 correlations within a network. Specifically, 1,386 SNP were nodes and the 5,132 strongest correlations (|r| ≥ 0.90) were edges. The network was filtered with genes queried from a transcriptome resource created from deep sequencing of RNA (i.e., RNA-Seq) from the hypothalamus of a prepubertal and a postpubertal Brangus heifer. The remaining hypothalamic-influenced network contained 978 genes connected by 2,560 edges or predicted gene interactions. This hypothalamic gene network was enriched with genes involved in axon guidance, which is a pathway known to influence pulsatile release of LHRH. There were 5 transcription factors with 21 or more connections: ZMAT3, STAT6, RFX4, PLAGL1, and NR6A1 for FSC. The SNP that identified these genes were intragenic and were on chromosomes 1, 5, 9, and 11. Chromosome 5 harbored both STAT6 and RFX4. The large number of interactions and genes observed with network analyses of multiple sources of genomic data (i.e., GWAS and RNA-Seq) support the concept of FSC being a polygenic trait.

Genetic markers of body composition and carcass quality in grazing Brangus steers.

The somatotropic axis is a major regulatory pathway of energy metabolism during postnatal growth in mammals. Genes involved in this pathway influence many economically important traits. The association of selected SNPs in these genes with carcass traits was examined in grazing Brangus steers. These traits included final live weight, ultrasound backfat thickness (UBFT), rib-eye area, kidney fat weight, hot carcass weight, and intramuscular fat percentage (%IMF). Genomic DNA (N = 246) was genotyped for a panel of 15 tag SNPs located in the growth hormone receptor (GHR), insulin-like growth factor I, insulin-like growth factor-binding protein 6, pro-melanin-concentrating hormone, suppressor of cytokine signaling 2, and signal transducer and activator of transcription 6 (STAT6) genes. Allelic and haplotype frequencies were compared with those of a sample of European breeds (N = 177 steers). Two tag SNPs in the GHR affected %IMF; one of them (ss86273136) was also strongly associated with UBFT (P < 0.003). The frequency of the most favorable GHR haplotype for %IMF was lower in Brangus steers. Moreover, the haplotype carrying two unfavorable alleles was present at a frequency of 31% in this group. Four tag SNPs on STAT6 had a significant effect on UBFT. One of these, SNP ss115492467, was also associated with %IMF. The STAT6 haplotype, including all the alleles favoring UBFT, was the most abundant variant (34%) in the European cattle, while it had a frequency of 14% in the Brangus steers. The four less favorable variants (absent in the European cattle) were found at a frequency of 38% in the Brangus steers. These results support the association of GHR and STAT6 SNP with carcass traits in composite breeds, such as Brangus, under grazing conditions.

Hot topic: performance of bovine high-density genotyping platforms in Holsteins and Jerseys.

Two high-density single nucleotide polymorphism (SNP) genotyping arrays have recently become available for bovine genomic analyses, the Illumina High-Density Bovine BeadChip Array (777,962 SNP) and the Affymetrix Axiom Genome-Wide BOS 1 Array (648,874 SNP). These products each have unique design and chemistry attributes, and the extent of marker overlap and their potential utility for quantitative trait loci fine mapping, detection of copy number variation, and multibreed genomic selection are of significant interest to the cattle community. This is the first study to compare the performance of these 2 arrays. Deoxyribonucleic acid samples from 16 dairy cattle (10 Holstein, 6 Jersey) were used for the comparison. An independent set of DNA samples taken from 46 Jersey cattle and 18 Holstein cattle were used to ascertain the amount of SNP variation accounted by the 16 experimental samples. Data were analyzed with SVS7 software (Golden Helix Inc., Bozeman, MT) to remove SNP having a call rate less than 90%, and linkage disequilibrium pruning was used to remove linked SNP (r² ≥ 0.9). Maximum, average, and median gaps were calculated for each analysis based on genomic position of SNP on the bovine UMD3.1 genome assembly. All samples were successfully genotyped (≥ 98% SNP genotyped) with both platforms. The average number of genotyped SNP in the Illumina platform was 775,681 and 637,249 for the Affymetrix platform. Based on genomic position, a total of 107,896 SNP were shared between the 2 platforms; however, based on genotype concordance, only 96,031 SNP had complete concordance at these loci. Both Affymetrix BOS 1 and Illumina BovineHD genotyping platforms are well designed and provide high-quality genotypes and similar coverage of informative SNP. Despite fewer total SNP on BOS 1, 19% more SNP remained after linkage disequilibrium pruning, resulting in a smaller gap size (5.2 vs. 6.9 kb) in Holstein and Jersey samples relative to BovineHD. However, only 224,115 Illumina and 241,038 Affymetrix SNP remained following removal of SNP with a minor allele frequency of zero in Holstein and Jersey samples, resulting in an average gap size of 11,887 bp and 11,018 bp, respectively. Combining the 354,348 informative (r² ≥ 0.9), polymorphic (minor allele frequency ≥ 0), unique SNP data from both platforms decreased the average gap size to 7,560 bp. Genome-wide copy number variant analyses were performed using intensity files from both platforms. The BovineHD platform provided an advantage to the copy number variant data compared with the BOS 1 because of the larger number of SNP, higher intensity signals, and lower background effects. The combined use of both platforms significantly improved coverage over either platform alone and decreased the gap size between SNP, providing a valuable tool for fine mapping quantitative trait loci and multibreed animal evaluation.

Genetic analyses involving microsatellite ETH10 genotypes on bovine chromosome 5 and performance trait measures in Angus- and Brahman-influenced cattle.

ETH10 is a dinucleotide microsatellite within the promoter of signal transducer and activator of transcription 6 (STAT6) gene on bovine chromosome 5. ETH10 is included in the panel of genetic markers used in parentage testing procedures of cattle breed associations. Allelic sizes of ETH10 PCR amplicons range from 199 to 225 bp. Objectives of this study were to use microsatellite data from beef cattle breed associations to investigate genetic distance and population stratification among Angus- and Brahman-influenced cattle and to use ETH10 genotypes and growth and ultrasound carcass data to investigate their statistical relationships. Three series of genotype to phenotype association analyses were conducted with 1) Angus data (n=5,094), 2) Brangus data (3/8 Brahman × 5/8 Angus; n=2,296), and 3) multibreed data (n=4,426) of Angus and Brangus cattle. Thirteen alleles and 38 genotypes were observed, but frequencies varied among breed groups. Tests of genetic identity and distance among 6 breed composition groups increasing in Brahman influence from 0 to 75% revealed that as Brahman-influence increased to ≥50%, genetic distance from Angus ranged from 18.3 to 43.5%. This was accomplished with 10 microsatellite loci. A mixed effects model involving genotype as a fixed effect and sire as a random source of variation suggested that Angus cattle with the 217/219 genotype tended to have 2.1% heavier (P=0.07) 205-d BW than other genotypes. In Brangus cattle, allele combinations were classified as small (≤215 bp) or large (≥217 bp). Brangus cattle with the small/large genotype had 2.0% heavier (P<0.05) birth weight, yet cattle with the large/large genotype had approximately 5.1% greater (P<0.05) percentage of fat within LM and more LM per BW than cattle with small/large or small/small genotypes. Genotype-to-phenotype relationships were not detected in multibreed analyses. The ETH10 locus appears to be associated with growth and carcass traits in Angus and Brangus cattle. Results from this study provide support for STAT6 as one of the candidate genes underlying cattle growth QTL on chromosome 5.

Variants in the pregnancy-associated plasma protein-A2 gene on Bos taurus autosome 16 are associated with daughter calving ease and productive life in Holstein cattle.

Reproductive disorders in dairy herds have a negative effect on farm profitability and sustainability of milk production. Given the substantial evidence of the role of the pregnancy-associated plasma protein (PAPP) gene family in the regulation of reproduction in humans and mice, its role in insulin-like growth factor metabolism, quantitative trait loci effects in the mouse, and location of a calving ease QTL on bovine chromosome 16, the PAPP-A2 gene was chosen as a candidate gene to perform an association study for reproductive health in cattle. Single nucleotide polymorphisms (SNP) were identified in coding and conserved noncoding regions of the PAPP-A2 gene in 3 dairy breeds. A total of 7 tag SNP were genotyped in 662 Holstein bulls (UCD-bulls) to perform marker trait association analysis. Three SNP (SNP 13, 15, and 16) were in strong linkage disequilibrium in Holsteins, showing significant positive associations with daughter calving ease, productive life, milk yield, and protein yield. These results were validated by genotyping SNP15 in a larger population of 992 bulls from the cooperative dairy DNA repository (CDDR-bulls). Our results demonstrate that the PAPP-A2 gene is associated with reproductive health in Holstein cattle and that the identified SNP can be used as genetic markers in dairy breeding due to their positive association with reproductive and productive traits. Functional studies need to be conducted to identify the mechanisms for the association of SNP with these traits.

Single nucleotide polymorphisms in the growth hormone-insulin-like growth factor axis in straightbred and crossbred Angus, Brahman, and Romosinuano heifers: population genetic analyses and association of genotypes with reproductive phenotypes.

The growth endocrine axis influences reproduction. The objectives of this study were to evaluate population genetic characteristics of SNP genotypes within genes of the GH-IGF axis in straightbred and crossbred Angus, Brahman, and Romosinuano heifers (n = 650) and to test the association of these genotypes with measures of reproduction. These objectives were achieved using 73 SNP within 7 genes on chromosome 5 and the pregnancy-associated plasma protein A2 (PAPP-A2) and GH-receptor genes, which map to chromosomes 16 and 20, respectively. The SNP were elucidated by resequencing conserved regions of each gene by using DNA from familial-unrelated cattle of a multibreed discovery population. A multiplex SNP assay yielded 59 biallelic SNP useful for evaluating genetic identity and distance. Specifically, the divergence of straightbred Brahman cattle was approximately 15.5% from 5 Bos taurus-influenced breed groups. In the straightbred groups used as a validation population, only 3 SNP had minor allele frequencies >10%. These SNP were in the genes PAPP-A2 (ss115492449-A/C and ss115492450-G/T within intron 10) and signal transducers and activators of transcription 2 (STAT2; ss252841035-A/G within the 5' untranslated region), and they met the conditions of Hardy-Weinberg equilibrium (P > 0.31). The other 56 SNP were useful for assigning each animal into ancestral clusters (n = 3 proportions) to account for population stratification in genotype to phenotype association analyses. The 2 SNP in the PAPP-A2 gene influenced (P < 0.05) traits indicative of first-calf heifer rebreeding (i.e., calving interval, days to calving, and pregnancy rate). A STAT2 SNP genotype (i.e., GG) × primary ancestral cluster interaction (P < 0.05) suggested heifers primarily of B. taurus ancestry had a reduction of approximately 16.4 ± 0.1% in calving interval and days to calving relative to heifers clustering primarily as Bos indicus ancestry. Even though additional research is needed to delineate the allelic variation attributed to genes of the GH-IGF axis, results of this study provide support for STAT2 and PAPP-A2 as potential candidate genes associated with first-calf heifer rebreeding traits.

Dairy Cattle Breeding Simulation Program: a simulation program to teach animal breeding principles and practices.

A dairy cattle breeding simulation program (DCBSP v.4.9) has been developed to teach undergraduate and graduate students animal breeding principles associated with selection for multiple traits in dairy cattle. The current version of the program was written in FORTRAN 90, and a web-based interface was developed for the students to interact with the program in the teaching environment. This software simulates a population of dairy cattle herds and artificial insemination bulls through several generations by integrating students' decisions about mating, culling, and selection of new heifers and bulls based on a multivariate animal mixed model evaluation and marker-assisted selection. All simulation parameters (e.g., number of herds and cows per herd, variance components, effect of genetic markers) can be defined by the administrator of the program in relation to the animal breeding course. During each running period, the program simulates the composition of each herd during a virtual year, generating new calves and new productive records and performing a genetic evaluation for all productive traits. A herd-specific productive summary of all demographic, productive, and genetic data is provided to the students at the end of each simulation period. After several running periods, the genetic trend can be evaluated, providing a realistic experience for the development of animal breeding skills that will be relevant to students with a basic knowledge of animal breeding. Earlier versions of this program have been used at several universities where it has proven to be a very useful teaching tool to illustrate the theoretical basis of animal breeding in livestock.

Polymorphisms in the STAT6 gene and their association with carcass traits in feedlot cattle.

Identification of the genes and polymorphisms underlying quantitative traits, and understanding how these genes and polymorphisms affect economic traits, are important for successful marker-assisted selection and more efficient management strategies in commercial cattle populations. Signal transducer and activator of transcription 6 (STAT6) gene is tightly connected to IL-4 and IL-13 signalling and plays a key role in T(H)2 polarization of the immune system. In addition, STAT6 acts as a mediator of leptin signalling and has been associated with body weight regulation. The objective of this study was to determine if SNPs within the bovine STAT6 gene are associated with economically important traits in feedlot cattle. The approach consisted of resequencing STAT6 using a panel of DNA from unrelated animals of different beef breeds. Specifically, 16 kb of STAT6 was resequenced in 47 animals and the process revealed 39 SNPs. From the 39 SNPs, a panel of 15 tag SNPs was genotyped in 1500 beef cattle samples with phenotypes to perform a marker-trait association analysis. Among the 15 tag SNPs, five and six were polymorphic in Bos taurus and Bos indicus respectively. An association analysis was performed between the 15 tag SNPs and 14 performance and production traits. SNP ss115492459:C > A, ss115492461:A > G and ss115492458:G > C were significantly associated with back fat, calculated yield grade, cutability, hot carcass weight, dry matter intake, days on feed, back fat rate and average daily gain. These three SNPs were present in all Bos taurus beef breeds examined. Our results provide evidence that polymorphisms in STAT6 are associated with carcass and growth efficiency traits, and may be used for marker-assisted selection and management in feedlot cattle.

Fine mapping and association analysis of a quantitative trait locus for milk production traits on Bos taurus autosome 4.

To fine map a quantitative trait locus (QTL) affecting milk production traits previously associated with microsatellite RM188, we implemented an interval mapping analysis by using microsatellite markers in a large Israeli Holstein half-sib sire family, and linkage disequilibrium (LD) mapping in a large set of US Holstein bulls. Interval mapping located the target QTL to the near vicinity of RM188. For the LD mapping, we identified 42 single nucleotide polymorphisms (SNP) in 15 genes in a 12-Mb region on bovine chromosome 4. A total of 24 tag SNP were genotyped in 882 bulls belonging to the University of California Davis archival collection of Holstein bull DNA samples with predicted transmitted ability phenotypes. Marker-to-marker LD analysis revealed 2 LD blocks, with intrablock r(2) values of 0.10 and 0.46, respectively; outside the blocks, r(2) values ranged from 0.002 to 0.23. A standard additive/dominance model using the generalized linear model procedure of SAS and the regression module of HelixTree software were used to test marker-trait associations. Single nucleotide polymorphism 9 on ARL4A, SNP10 on XR_027435.1, SNP12 on ETV1, SNP21 on SNX13, and SNP24 were significantly associated with milk production traits. We propose the interval encompassing ARL4A and SNX13 genes as a candidate region in bovine chromosome 4 for a concordant QTL related to milk protein traits in dairy cattle. Functional studies are needed to confirm this result.