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.

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.

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.

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.

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.

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

Extensive long-range and nonsyntenic linkage disequilibrium in livestock populations: deconstruction of a conundrum.

Great interest was aroused by reports, based on microsatellite markers, of high levels of statistically significant long-range and nonsyntenic linkage disequilibrium (LD) in livestock. Simulation studies showed that this could result from population family structure. In contrast, recent SNP-based studies of livestock populations report much lower levels of LD. In this study we show, on the basis of microsatellite data from four cattle populations, that high levels of long-range LD are indeed obtained when using the multi-allelic D' measure of LD. Long-range and nonsyntenic LD are exceedingly low, however, when evaluated by the standardized chi-square measure of LD, which stands in relation to the predictive ability of LD. Furthermore, specially constructed study populations provided no evidence for appreciable LD resulting from family structure at the grandparent level. We propose that the high statistical significance and family structure effects observed in the earlier studies are due to the use of large sample sizes, which accord high statistical significance to even slight deviations from asymptotic expectations under the null hypothesis. Nevertheless, even after taking sample size into account, our results indicate that microsatellites testify to the presence of usable LD at considerably wider separation distances than SNPs, suggesting that use of SNP haplotypes may considerably increase the usefulness of a given fixed SNP array.

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.

Valuation of milk composition and genotype in cheddar cheese production using an optimization model of cheese and whey production.

A mass balance optimization model was developed to determine the value of the kappa-casein genotype and milk composition in Cheddar cheese and whey production. Inputs were milk, nonfat dry milk, cream, condensed skim milk, and starter and salt. The products produced were Cheddar cheese, fat-reduced whey, cream, whey cream, casein fines, demineralized whey, 34% dried whey protein, 80% dried whey protein, lactose powder, and cow feed. The costs and prices used were based on market data from March 2004 and affected the results. Inputs were separated into components consisting of whey protein, ash, casein, fat, water, and lactose and were then distributed to products through specific constraints and retention equations. A unique 2-step optimization procedure was developed to ensure that the final composition of fat-reduced whey was correct. The model was evaluated for milk compositions ranging from 1.62 to 3.59% casein, 0.41 to 1.14% whey protein, 1.89 to 5.97% fat, and 4.06 to 5.64% lactose. The kappa casein genotype was represented by different retentions of milk components in Cheddar cheese and ranged from 0.715 to 0.7411 kg of casein in cheese/kg of casein in milk and from 0.7795 to 0.9210 kg of fat in cheese/kg of fat in milk. Milk composition had a greater effect on Cheddar cheese production and profit than did genotype. Cheese production was significantly different and ranged from 9,846 kg with a high-casein milk composition to 6,834 kg with a high-fat milk composition per 100,000 kg of milk. Profit (per 100,000 kg of milk) was significantly different, ranging from $70,586 for a high-fat milk composition to $16,490 for a low-fat milk composition. However, cheese production was not significantly different, and profit was significant only for the lowest profit ($40,602) with the kappa-casein genotype. Results from this model analysis showed that the optimization model is useful for determining costs and prices for cheese plant inputs and products, and that it can be used to evaluate the economic value of milk components to optimize cheese plant profits.

Characterization of variation in the canine suppressor of cytokine signaling-2 (SOCS2) gene.

Suppressor of cytokine signaling 2 (SOCS2) is a negative regulator of growth hormone signaling. The deletion of SOCS2 in mice results in a 30-50% increase in post-natal growth. In an effort to identify polymorphisms in the SOCS2 gene that may be associated with body size in dogs, we characterized the canine SOCS2 gene and analyzed its genetic diversity among small and large dog breeds. The study was carried out on a total of 520 dogs from 66 different breeds. Dogs were classified as large or small based on height and weight as determined by their respective American Kennel Club breed standards. The SH2 and SOCS domains of the canine SOCS2 gene were sequenced in 32 dogs from different breeds. Only one non-synonymous sequence variant (DQ415457:g.326G>T) was detected which corresponds to an amino acid change (Asp127Tyr). All samples were genotyped by PCR/RFLP and the allele frequencies were determined for each dog breed. The T allele was distributed primarily among European large dog breeds with a gene frequency ranging from 0.72 to 0.04. The nature of the nucleotide change and the effect on the protein together with the finding of a QTL related to body size in the same CFA15 region by other researchers suggest canine SOCS2 as a potential candidate gene for body size in dogs. Future studies will be needed to clarify the role of the 326G>T polymorphism and its interaction with genes like growth hormone and insulin-like growth factor 1.

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.

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.

Growth Rate, Body Composition, Cellular Growth and Enzyme Activities in Lines of TRIBOLIUM CASTANEUM Selected for 21-Day Pupa Weight.

Growth rate, body composition, cell number, cell size, and the activity of four dehydrogenase enzymes were studied from 10 to 25 days of age in one control (1C) and three lines (3, 9, 10) of Tribolium castaneum that had been subjected to long term selection for large 21-day pupae weight.-Selected lines were two- to three-fold larger in size than the control line throughout development. No major differences in percent of protein were detected among the lines but at any particular age, the selected lines were found to have a higher fat content than the control line. The differences in fat content were closely correlated with development such that all the lines reached very similar levels of percent of fat just prior to pupation. Water content showed an inverse relation to percent of fat.-Selection was observed to have caused major changes in the cellular response to growth. The selected lines had an average of from 17% to 48% larger cells (measured as protein/DNA) and were found to have from 37 to 62% more cells (measured as total DNA) than the control line at all ages from 10 to 19 days of age. In addition, the selected lines had a higher RNA content at all ages studied and higher RNA:DNA ratios at the young ages. In contrast the enzyme activities of ICDH and LDH were 60% lower. The results are interpreted as indicating that a more efficient metabolic machinery had evolved in the rapidly growing selected lines.