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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08430-x.

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.

3 Embryo knockout efficiency improved when targeting ovine suppressor of cytokine signalling 2 with 2 small guide RNA

Suppressor of cytokine signalling (SOCS2) is a key regulator of growth. Growth hormone (GH) and insulin growth factor-1 induce SOCS2 expression, which acts as a suppressor of GH and insulin growth factor-1 signalling by preventing the activation of the signal transducer and activator of transcription proteins. In mice, SOCS2 knockout results in a gigantic phenotype with elevated GH levels, longer bones, enlargement of internal organs, and 30 to 40% heavier body weight in comparison with wild-type animals. High growth was also noted in naturally occurring SOCS2 mutant pigs. The aim of this project was to produce high-growth sheep by introducing sense mutations to SOCS2 by direct injection of CRISPR/Cas9 in ovine embryos. In this study, we compared the editing efficiency of 2 small guide RNA (sgRNA), which were designed to target the first exon of SOCS2. Ovine ovaries were collected from a local slaughterhouse and cumulus-oocyte complexes were aspirated from 2- to 6-mm antral follicles, matured invitro for 22h, and then denuded from the surrounding cumulus cells. Approximately 6 pL of CRISPR/Cas9 solution (80ng µL−1 of sgRNA, 80ng µL−1 of Cas9 mRNA, and 160ng µL−1 of Cas9 protein in low-TE buffer) was microinjected into the cytoplasm of metaphase II oocytes using a laser-assisted microinjection technique (Bogliotti et al. 2016 J. Vis. Exp. 116, e54465). Groups of 50 oocytes were injected with the CRISPR/Cas9 mix containing sgRNA1 (group G1), sgRNA2 (group G2), sgRNA1 and sgRNA2 (group G1/2), or low-TE buffer alone (control). After microinjection, IVF was carried out using fresh semen, and oocytes were co-incubated with the sperm for 16h. Embryos were cultured in groups of 25 in 50-µL drops of KSOM (Evolve, Zenith Technologies, Cork, Ireland) with 4mg mL−1 of BSA under oil at 38.5°C, 5% CO2, 5% O2, and 90% N2. Blastocysts were collected on Day 7 post-fertilization to assess mutation efficiency. The CRISPR/Cas9 target sequence was amplified using nested PCR from single-blastocyst whole lysates and analysed by Sanger sequencing. Mutations in each embryo were identified by visual inspection of the sequence traces aided by trace deconvolution software (SnapGene, GSL Biotech, Chicago, IL, USA). The experiment was replicated 3 times. Blastocyst rates were 22.6, 21.3, 21.3, and 29.3% for G1, G2, G1/2, and control, respectively. Targeting the exon with 2 sgRNA (G1/2) resulted in the highest mutation rate (87.3%), with an average indel size of 55.6bp, leading to 100% nonfunctional SOCS2, and with 87.6% of the sequenced embryos presenting bi-allelic mutations. Using a single sgRNA caused mutations in 47.3 and 38.6% of embryos for G1 and G2, respectively, with averages indel size of 8.2bp in G1 and 13.8bp in G2. In G1, 88.6% of the mutations were biallelic and 90% would result in nonfunctional SOCS2. In G2, 50% of the mutations were biallelic and 40% translated to nonfunctional SOCS2. Based on these results, the percentage of SOCS2 knockout embryos was 37.7, 7.7, and 76.4% in G1, G2, and G1/2, respectively. In conclusion, the knockout efficiency improved after targeting ovine SOCS2 with 2 sgRNA.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Promoter region of the bovine growth hormone receptor gene: single nucleotide polymorphism discovery in cattle and association with performance in Brangus bulls

Expression of the GH receptor (GHR) gene and its binding with GH is essential for growth and fat metabolism. A GT microsatellite exists in the promoter of bovine GHR segregating short (11 bp) and long (16 to 20 bp) allele sequences. To detect SNP and complete an association study of genotype to phenotype, we resequenced a 1,195-bp fragment of DNA including the GT microsatellite and exon 1A. Resequencing was completed in 48 familialy unrelated Holstein, Jersey, Brown Swiss, Simmental, Angus, Brahman, and Brangus cattle. Nine SNP were identified. Phylogeny analyses revealed minor distance (i.e., <5%) in DNA sequence among the 5 Bos taurus breeds; however, sequence from Brahman cattle averaged 27.4 +/- 0.07% divergence from the Bos taurus breeds, whereas divergence of Brangus was intermediate. An association study of genotype to phenotype was completed with data from growing Brangus bulls (n = 553 from 96 sires) and data from 4 of the SNP flanking the GT microsatellite. These SNP were found to be in Hardy-Weinberg equilibrium and in phase based on linkage disequilibrium analyses (r(2) = 0.84 and D'= 0.92). An A/G tag SNP was identified (ss86273136) and was located in exon 1A, which began 88 bp downstream from the GT microsatellite. Minor allele frequency of the tag SNP was greater than 10%, and Mendelian segregation was verified in 3 generation pedigrees. The A allele was derived from Brahman, and the G allele was derived from Angus. This tag SNP genotype was a significant effect in analyses of rib fat data collected with ultrasound when bulls were ~365 d of age. Specifically, bulls of the GG genotype had 6.1% more (P = 0.0204) rib fat than bulls of the AA and AG genotypes, respectively. Tag SNP (ss86273136), located in the promoter of GHR, appears to be associated with a measure of corporal fat in Bos taurus x Bos indicus composite cattle.

Birth of large calves that developed from in vitro-derived bovine embryos

High birth weights were observed in calves that developed from bovine embryos produced by in vitro maturation (IVM) and in vitro fertilization (IVF) procedures. After IVM and IVF, embryos were either co-cultured in vitro with oviductal epithelial cells or transferred into the sheep oviduct for development to the blastocyst stage. Blastocysts were transferred to the reproductive tracts of recipient heifers and cows for development to term. Birth weights and gestation periods were compared between calves that developed from in vitro-derived embryos and calves born after artificial insemination (AI) of cows in the herd from which recipient females were selected. Gestation periods were not different among the groups (P > 0.05), but calves that developed from IVM/IVF-derived embryos co-cultured in vitro were larger at birth than calves born from IVM/IVF-derived embryos that developed into blastocysts in the sheep oviduct and calves born from AI (P < 0.001). Dystocia and calf mortality were associated with large calf size at birth. These data were collected from an experiment designed for other purposes, and confounding variables and small sample size could have influenced the observed differences in birth weights. Nevertheless, the extreme birth weights of some calves suggest that abnormal prenatal growth occurs in some IVM/IVF-derived bovine embryos and that conditions for co-culture to the blastocyst stage may exacerbate the problem.

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.

Identification of sex chromosome molecular markers using RAPDs and fluorescent in situ hybridization in rainbow trout

The goal of this work is to identify molecular markers associated with the sex chromosomes in rainbow trout to study the mode of sex determination mechanisms in this species. Using the RAPD assay and bulked segregant analysis, two markers were identified that generated polymorphic bands amplifying preferentially in males of the Mount Lassen and Scottish strains of rainbow trout. Chromosomal localization using fluorescent in situ hybridization of a 900 bp probe developed from one of these markers revealed a brightly defined signal on a chromosome that could morphologically be classified as the Y chromosome.

Hot Topic: Endogenous Production of n-3 and n-6 Fatty Acids in Mammalian Cells

Polyunsaturated fatty acids (PUFA) are important components of mammalian diets, and the beneficial effects of n-3 PUFA on human development and cardiovascular health have been well documented. Caenorhabditis elegans is one of the few animals known to be able to produce linoleic (LA, 18:2n-6) and alpha-linolenic (ALA, 18:3n-3) essential fatty acids. These essential PUFA are generated by the action of desaturases that successively direct the conversion of monounsaturated fatty acids (MUFA) to PUFA. The cDNA coding sequences of the C. elegans Delta(12) and n-3 fatty acid desaturases were each placed under the control of separate constitutive eukaryotic promoters and simultaneously introduced into HC11 mouse mammary epithelial cells by adenoviral transduction. Phospholipids from transduced cells showed a significant decrease in the ratios of both MUFA:PUFA and n-6:n-3 fatty acids relative to control cultures. The fatty acid profile of transduced cellular phospholipids revealed significant decreases in MUFA and arachidonic acid (20:4n-6), and increases in LA, ALA, and eicosapentaenoic acid (20:5n-3). The fatty acid composition of triacylglycerols derived from transduced cells was similarly, but less dramatically, affected. These results demonstrate the functionality of C. elegans fatty acid desaturase enzymes in mammalian cells. Expression of these desaturases in livestock might act to counterbalance the saturating effect that rumen microbial biohydrogenation has on the fatty acid profile of ruminant products, and allow for the development of novel, land-based dietary sources of n-3 PUFA.

Identification of putative homology between horse microsatellite flanking sequences and cross-species ESTs, mRNAs and genomic sequences

In this study the flanking sequences of 1534 horse microsatellites were used in a BLAST search to identify putative human-horse homologies. BLAST searches revealed 129 flanking sequences with significant blastn matches [alignment scores (S) > or = 60 and sum probability values (E) < or = 3.0E-6], also, 25 of these produced significant blastx matches. To provide a reference point in the human genome the flanking sequences with matches were subjected to a BLAT search of the University of California Santa Cruz (UCSC) human genome assembly (July 2003 freeze). Eighty-three of the flanking sequences showed high similarity to sequence of known or putative human genes and the remaining 46 demonstrated high similarity to human intragenic regions. Interestingly, 87 of the microsatellites showed conservation of the tandem repeat in addition to flanking regions. Overall, 41 of the microsatellites had been mapped in the horse and of these 37 localized to the expected syntenic location. The other four did not and represent new putative regions of human-horse synteny. The results of this study contribute 79 new putative human-horse homologies, increasing the density of markers on the human-horse comparative map.

Putative in silico mapping of DNA sequences to livestock genome maps using SSLP flanking sequences

In this study, an in silico approach was developed to identify homologies existing between livestock microsatellite flanking sequences and GenBank nucleotide sequences. Initially, 1955 bovine, 1570 porcine and 1121 chicken microsatellites were downloaded and the flanking sequences were compared with the nr and dbEST databases of GenBank. A total of 74 bovine, 44 porcine and 37 chicken microsatellite flanking sequences passed our criteria and had at least one significant match to human genomic sequence, genes/expressed sequence tags (ESTs) or both. GenBank annotation and BLAT searches of the UCSC human genome assembly revealed that 38 bovine, 13 porcine and 17 chicken microsatellite flanking sequences were highly similar to known human genes. Map locations were available for 67 bovine, 44 porcine and 21 chicken microsatellite flanking sequences, providing useful links in the comparative maps of humans and livestock. In support of our approach, 112 alignments with both microsatellite and match mapping information were located in the expected chromosomal regions based on previously reported syntenic relationships. The development of this in silico mapping approach has significantly increased the number of genes and EST sequences anchored to the bovine, porcine and chicken genome maps and the number of links in various human-livestock comparative maps.

Characterization of sex chromosomes in rainbow trout and coho salmon using fluorescence in situ hybridization (FISH)

With the aim of characterizing the sex chromosomes of rainbow trout (Oncorhynchus mykiss) and to identify the sex chromosomes of coho salmon (O. kisutch), we used molecular markers OmyP9, 5S rDNA, and a growth hormone gene fragment (GH2), as FISH probes. Metaphase chromosomes were obtained from lymphocyte cultures from farm specimens of rainbow trout and coho salmon. Rainbow trout sex marker OmyP9 hybridizes on the sex chromosomes of rainbow trout, while in coho salmon, fluorescent signals were localized in the medial region of the long arm of one subtelocentric chromosome pair. This hybridization pattern together with the hybridization of a GH2 intron probe on a chromosome pair having the same morphology, suggests that a subtelocentric pair could be the sex chromosomes in this species. We confirm that in rainbow trout, one of the two loci for 5S rDNA genes is on the X chromosome. In males of this species that lack a heteromorphic sex pair (XX males), the 5S rDNA probe hybridized to both subtelocentrics. This finding is discussed in relation to the hypothesis of intraspecific polymorphism of sex chromosomes in rainbow trout.

Cloning and localization of the bovine and ovine lysophosphatidic acid acyltransferase (LPAAT) genes that codes for an enzyme involved in triglyceride biosynthesis

Lysophosphatidic acid acyltransferase (LPAAT) catalyzes the addition of fatty acyl moieties to the sn-2 position of the glycerol backbone of lysophosphatidic acid in triglyceride biosynthesis. In this study, we have cloned, sequenced, and characterized the bovine and ovine LPAAT cDNA. Both encode proteins of 287 amino acids with molecular masses of 32 and 31.9 kDa, respectively, differing only by a single amino acid residue. The bovine and ovine LPAAT are predicted to be transmembrane enzymes localized to the endoplasmic reticulum. We also characterized the sequence and genomic organization of the bovine LPAAT gene. The gene consists of seven exons and six introns, spanning a 7.5-kb distance. With the use of a whole genome radiation hybrid panel, we localized the bovine LPAAT to the central region of chromosome 23.

Inbred lines of mice derived from long-term growth selected lines: unique resources for mapping growth genes

Lines of mice selected for many generations for high or low growth in several laboratories around the world have been collected, and from these, inbred lines are being developed by recurrent full-sib mating in Edinburgh. There are seven high selected lines and four low lines (each low line is from the same base population as one of the high lines), and the histories of each are summarized. Mean body weight of males at 70 days of age in the Edinburgh laboratory in the heaviest inbred line (77 g) is 4.8-fold higher than in the lightest line (16 g), and 1.9-fold higher than in the least extreme high line (41 g). Litter size, food intake, and fat content also differ substantially. These inbred extreme selected lines are a uniquely valuable resource for QTL or gene mapping, candidate gene identification, and elucidation of epistatic effects.

Development and characterization of DNA sequence OmyP9 associated with the sex chromosomes in rainbow trout

This work describes the construction and characterization of a sequence characterized amplified DNA region (SCAR DNA), designated OmyP9, that was derived from a RAPD marker associated with the sex chromosomes in rainbow trout. A RsaI restriction fragment length polymorphism in OmyP9 identifies variants A, B and C. We found six OmyP9 variant phenotypes - A, B, C, AB, BC and ABC, in 186 individuals of seven different rainbow trout strains. The patterns of inheritance of OmyP9 in 139 fingerlings from 10 crosses of three strains of rainbow trout were studied. The males had a greater representation of the A variant (93.3%) suggesting an association with the Y chromosome. All male fingerlings analysed inherited the A variant from their male parents. These results support the hypothesis that OmyP9 is located on the sex chromosomes of rainbow trout, and that for the males studied the A variant is located on the Y chromosome in a region close to sex determinants and/or in a sector where the genetic recombination between X and Y is restricted. The present evidence also supports our previous hypothesis that OmyP9 is organized as a tandem repeated sequence in the sex chromosomes of rainbow trout. We feel that the OmyP9 RsaI marker can be used for sex identification in crosses where it is possible to determine the phenotype of the parents.

Quantitative trait loci affecting growth in high growth (hg) mice

A genome-wide scan was performed in order to identify Quantitative Trait Loci (QTL) associated with growth in a population segregating high growth (hg), a partially recessive mutation that enhances growth rate and body size in the mouse. A sample of 262 hg/hg mice was selected from a C57BL/6J-hg/hg x CAST/EiJ F2 cross and typed with 79 SSLP markers distributed across the genome. Eight significant loci were identified through interval mapping. Loci on Chromosomes (Chrs) 2 and 8 affected the growth rate of F2 mice. Loci on Chr 2 and 11 affected growth rate and carcass lean mass (protein and ash). A locus on Chr 9 modified femur length and another one in Chr 17 affected both carcass lean mass and femur length, but none of these had significant effects on growth rate. Loci on Chrs 5 and 9 modified carcass fat content. Additive effects were positive for C57BL/6J alleles, except for the two loci affecting carcass fatness. Typing of selected markers in 274 +/+ F2 mice revealed significant interactions between hg and other growth QTL, which were detected as changes in gene action (additive or dominant) and in allele substitution effects. Knowledge about interactions between loci, especially when major genes are involved, will help in the identification of positional candidate genes and in the understanding of the complex genetic regulation of growth rate and body size in mammals.

Lack of Socs2 expression causes the high-growth phenotype in mice

Characterizing causal molecular defects in mouse models of overgrowth or dwarfism helps to identify the key genes and pathways that regulate the growth process. We report here the molecular basis for high growth (hg), a spontaneous mutation that causes a 30-50% increase in postnatal growth. We conclude that hg is an allele of the suppressor of cytokine signaling 2 (Socs2), a member of a family of regulators of cytokine signal transduction. We demonstrate mapping of Socs2 to the hg region, lack of Socs2 mRNA expression, a disruption of the Socs2 locus in high-growth (HG) mice, and a similarity of phenotypes of HG mice and Socs2(-/-) mice generated by gene targeting. Characteristics of the HG phenotype suggest that Socs2 deficiency affects growth prenatally and postnatally most likely through deregulating the growth hormone (GH)/insulin-like growth factor I (IGF1). These results demonstrate a critical role for Socs2 in controlling growth.