Feed restriction, but not l-carnitine infusion, alters the liver transcriptome by inhibiting sterol synthesis and mitochondrial oxidative phosphorylation and increasing gluconeogenesis in mid-lactation dairy cows.

Abomasal carnitine infusion during acute feed restriction increases hepatic fatty acid oxidation and decreases liver lipid in dairy cows. Eight mid-lactation Holstein cows were used in a replicated 4×4 Latin square design with 14-d periods. A 2×2 factorial arrangement was used to determine the effects of water infusion+ad libitum dry matter intake (DMI), water infusion+restricted DMI (50% of previous 5-d average), l-carnitine infusion (20 g/d)+ad libitum DMI, or l-carnitine infusion+restricted DMI. Liver RNA from 7 healthy cows was used for transcriptome profiling using a bovine microarray. An ANOVA with a false discovery rate was used to identify treatment and interaction effects. A substantial transcriptome change was observed only with DMI restriction, resulting in 312 (155 downregulated, 157 upregulated) differentially expressed genes. Quantitative PCR was performed to verify microarray data and measure expression of additional genes not present on the microarray. The quantitative PCR data confirmed the effect of feed restriction but not of l-carnitine treatment. Feed restriction increased expression of GPX3 and of genes associated with gluconeogenesis (PC, PDK4), inflammation (SAA3), and signaling (ADIPOR2). In contrast, feed restriction downregulated BBOX, a key for l-carnitine biosynthesis, and the transcription factor HNF4A. The bioinformatics functional analysis of genes affected by DMI restriction uncovered biosynthesis of cholesterol and energy generation by mitochondrial respiration as the most relevant and inhibited functions. The data also indicated an increase of flux toward gluconeogenesis. We interpreted those results as a likely response of the liver to spare energy and provide glucose for the lactating mammary gland during feed deprivation.

Level of nutrient intake affects mammary gland gene expression profiles in preweaned Holstein heifers.

Bovine mammary parenchyma (PAR) and fat pad (MFP) development are responsive to preweaning level of nutrient intake. We studied transcriptome alterations in PAR and MFP from Holstein heifer calves (n=6/treatment) fed different nutrient intakes from birth to ca. 65 d age. Conventional nutrient intake received 441 g of dry matter (DM)/d of a control milk replacer (MR) [CON; 20% crude protein (CP), 20% fat, DM basis]. Calves in the accelerated nutrition groups received 951 g/d of high-protein/low-fat MR (HPLF; 28% CP, 20% fat, DM basis), 951 g/d of high-protein/high-fat MR (HPHF; 28% CP, 28% fat, DM basis), or 1,431 g/d of HPHF (HPHF+) MR. Out of 13,000 genes evaluated, over 1,500 differentially expressed genes (DEG) were affected (false discovery rate <0.10) by level of nutrient intake in PAR or MFP. Feeding HPLF versus CON resulted in the most dramatic changes in gene expression, with 278 and 588 DEG having ≥1.5-fold change in PAR and MFP. In PAR, the most-altered molecular functions were associated with metabolism of the cell (molecular transport and lipid metabolism) with most of the genes downregulated in HPLF versus CON. In MFP, DEG also were primarily associated with metabolism but changes also occurred in genes linked to cell morphology, cell-to-cell signaling, and immune response. Compared with CON, feeding HPHF or HPHF+ did not result in substantial additional effects on DEG beyond those observed with HPLF. The pentose phosphate, mitochondrial dysfunction, and ubiquinone biosynthesis pathways were among the most enriched due to HPLF versus CON in PAR and were inhibited, whereas glycosphingolipid biosynthesis, arachidonic acid metabolism, and eicosanoid synthesis pathways were among the most enriched due to HPLF versus CON in MFP and were inhibited. These responses suggest that, in PAR, doubling nutrient intake from standard feeding rates inhibited energy metabolism and activity of oxidative pathways that partly serve to protect cells against oxidative stress. The MFP in those heifers appeared to decrease production of lipid-derived metabolites that may play roles in signaling pathways within the adipocyte. Overall, results indicated that prepubertal/preweaned mammary transcriptome is responsive to long-term enhanced nutrient supply to achieve greater growth rates before weaning. The biological significance of these results to future milk production remains to be elucidated.

Comparative analysis of the early transcriptome of Brucella abortus--infected monocyte-derived macrophages from cattle naturally resistant or susceptible to brucellosis.

Brucellosis is a worldwide zoonotic infectious disease that has a significant economic impact on animal production and human public health. We characterized the gene expression profile of B. abortus-infected monocyte-derived macrophages (MDMs) from naïve cattle naturally resistant (R) or susceptible (S) to brucellosis using a cDNA microarray technology. Our data indicate that (1) B. abortus induced a slightly increased genome activation in R MDMs and a down-regulated transcriptome in S MDMs, during the onset of infection, (2) R MDMs had the ability to mount a type 1 immune response against B. abortus infection which was impaired in S cells, and (3) the host cell activity was not altered after 12 h post-B. abortus infection in R MDMs while the cell cycle was largely arrested in infected S MDMs at 12 h p.i. These results contribute to an improved understanding of how host responses may be manipulated to prevent infection by brucellae.

Discovery, validation and characterization of 1039 cattle single nucleotide polymorphisms.

We identified approximately 13 000 putative single nucleotide polymorphisms (SNPs) by comparison of repeat-masked BAC-end sequences from the cattle RPCI-42 BAC library with whole-genome shotgun contigs of cattle genome assembly Btau 1.0. Genotyping of a subset of these SNPs was performed on a panel containing 186 DNA samples from 18 cattle breeds including 43 trios. Of 1039 SNPs confirmed as polymorphic in the panel, 998 had minor allele frequency > or =0.25 among unrelated individuals of at least one breed. When Btau 4.0 became available, 974 of these validated SNPs were assigned in silico to known cattle chromosomes, while 41 SNPs were mapped to unassigned sequence scaffolds, yielding one SNP every approximately 3 Mbp on average. Twenty-four SNPs identified in Btau 1.0 were not mapped to Btau 4.0. Of the 1015 SNPs mapped to Btau 4.0, 959 SNPs had nucleotide bases identical in Btau 4.0 and Btau 1.0 contigs, whereas 56 bases were changed, resulting in the loss of the in silico SNP in Btau 4.0. Because these 1039 SNPs were all directly confirmed by genotyping on the multi-breed panel, it is likely that the original polymorphisms were correctly identified. The 1039 validated SNPs identified in this study represent a new and useful resource for genome-wide association studies and applications in animal breeding.

Biological interpretations of transcriptomic profiles in mammalian oocytes and embryos.

The characterization of gene-expression profiles in oocytes and embryos is critical to understand the influence of genetic and environmental factors on preimplantation and fetal development. Numerous gene-expression microarray studies using different platforms and species are offering insights into the biological processes extensively represented among the genes exhibiting differential expression. Major advances on understanding the direct relationship between gene expression and developmental competence are being reported. Integration of information across studies using meta-analysis techniques can increase the precision and accuracy to identify expression profiles associated with embryo development. Gene network and pathway analyses are offering insights into gene interactions and expression profiles of embryos. All these advances are cementing the way toward a comparative and systems approach to understanding the complex processes underlying vertebrate development.

Comparative mapping of cattle chromosome 19: cytogenetic localization of 19 BAC clones.

Here we present the results of fluorescent in situ hybridization (FISH) mapping of a set of cattle BAC clones preselected for assignment on cattle chromosome 19 (BTA19). The BAC clones were anchored to human chromosome 17 (HSA17) sequences by BLASTn similarity search of cattle BAC-ends against the human genome sequence (NCBI build 33). Five blocks of homologous synteny were defined in the comparative map of BTA19 and HSA17 built with FISH data and the human genome coordinates. The positions for four evolutionary breakpoints in the bovine and human chromosomes were identified. Comparison of the FISH comparative map with previously published comparative RH, physical, and cytogenetic maps of BTA19 did not reveal major conflicts and allowed for the extension of the boundaries of homology between BTA19 and HSA17. Comparative analysis of HSA17, BTA19, and mouse chromosome 11 (MMU11) demonstrates that most likely mice retain the ancestral organization of the synteny group, and both cattle and human chromosomes underwent several major internal rearrangements after the divergence of Primates, Rodentia, and Cetartiodactyla.

Detection of quantitative trait loci influencing conformation traits and calving ease in Holstein-Friesian cattle.

An extension of our previous genome scan of a North American Holstein-Friesian population was conducted to identify quantitative trait loci (QTL) affecting conformation traits. Resource families consisted of 1404 sons of 10 elite sires. Genome coverage was estimated to be 2713.5 cM (90%) for 406 markers using a granddaughter design. Regression interval mapping was used to detect QTL affecting 22 conformation traits, including body, udder, feet and legs, and dairy conformation as well as calving ease. Analysis of the families jointly identified 41 chromosome-wise significant QTL influencing conformation traits and 3 significant QTL influencing calving ease on 20 chromosomes. The false discovery rate method was used to account for multiple testing and 3/4 of the suggestive and 5/6 of significant QTL should be real effects. Fourteen of the 44 QTL were significant at the genome-wise level. Comparison of these results with other published reports identifies common QTL affecting conformation traits. Regions on 10 chromosomes appear to affect multiple traits, including conformation, milk production, and somatic cell score, within these particular US Holstein families. Additional work is needed to determine the precise locations of the QTL and select positional candidate genes influencing these traits.

Detection of quantitative trait loci affecting milk production, health, and reproductive traits in Holstein cattle.

We report putative quantitative trait loci affecting female fertility and milk production traits using the merged data from two research groups that conducted independent genome scans in Dairy Bull DNA Repository grandsire families to identify quantitative trait loci (QTL) affecting economically important traits. Six families used by both groups had been genotyped for 367 microsatellite markers covering 2713.5 cM of the cattle genome (90%), with an average spacing of 7.4 cM. Phenotypic traits included PTA for pregnancy rate and daughter deviations for milk, protein and fat yields, protein and fat percentages, somatic cell score, and productive life. Analysis of the merged dataset identified putative quantitative trait loci that were not detected in the separate studies, and the pregnancy rate PTA estimates that recently became available allowed detection of pregnancy rate QTL for the first time. Sixty-one putative significant marker effects were identified within families, and 13 were identified across families. Highly significant effects were found on chromosome 3 affecting fat percentage and protein yield, on chromosome 6 affecting protein and fat percentages, on chromosome 14 affecting fat percentage, on chromosome 18 affecting pregnancy rate, and on chromosome 20 affecting protein percentage. Within-family analysis detected putative QTL associated with pregnancy rate on six chromosomes, with the effect on chromosome 18 being the most significant statistically. These findings may help identify the most useful markers available for QTL detection and, eventually, for marker-assisted selection for improvement of these economically important traits.

Consensus and comprehensive linkage maps of bovine chromosome 24.

This study describes development of a consensus genetic linkage map of bovine chromosome 24 (BTA24). Eight participating laboratories contributed data for 58 unique markers including a total of 25 409 meioses. Eighteen markers, which were typed in more than one reference population, were used as potential anchors to generate a consensus framework map. The framework map contained 16 loci ordered with odds greater than 1000:1 and spanned 79.3 cM. Remaining markers were included in a comprehensive map relative to these anchors. The resulting BTA24 comprehensive map was 98.3 cM in length. Average marker intervals were 6.1 and 2.5 cM for framework and comprehensive maps, respectively. Marker order was generally consistent with previously reported BTA24 linkage maps. Only one discrepancy was found when comparing the comprehensive map with the published USDA-MARC linkage map. Integration of genetic information from different maps provides a high-resolution BTA24 linkage map.

Interval and composite interval mapping of somatic cell score, yield, and components of milk in dairy cattle.

Single-marker, interval-mapping (IM) and composite interval mapping (CIM) were used to detect quantitative trait loci (QTL) controlling milk, fat and protein yields, and somatic cell score (SCS). A granddaughter design was used to combine molecular genetic information with predicted transmitting abilities (PTA) and estimated daughter yield deviations (DYD) from eight Dairy Bull DNA Repository Holstein families. Models that included and excluded weights accounting for the uncertainty of the response variable were evaluated in each trait, family and phenotype (DYD and PTA) combination. The genotypic information consisted of 174 microsatellite markers along 29 Bos taurus autosomes. The average number of informative markers per autosome was three and the number of informative sons per family and marker varied between 21 and 173. Within-family results from the least squares single-marker analyses were used in expectation-maximization likelihood IM and CIM implemented with QTL Cartographer. Different CIM model specifications, offering complementary control on the background QTL outside the interval under study, were evaluated. Permutation techniques were used to calculate the genome-wide threshold test statistic values based on 1,000 samples. Results from the DYD and PTA analyses were highly consistent across traits and families. The minor differences in the estimates from the models that accounted for or ignored the uncertainty of the DYD (variance) and PTA (inverse of reliability) may be associated to the elevated and consistent precision of the DYD and PTA among sons. The CIM model best supported by the data had 10 markers controlling for background effects. On autosome (BTA) three, a QTL at 32 cM influencing protein yield was located in family five and a QTL at 74 cM for fat yield was located in family eight. Two map positions associated with SCS were detected on BTA 21, one at 33 cM in family one and the other at 84 cM in family three. A QTL for protein yield was detected between 26 and 36 cM on BTA six, family six, and a QTL for milk yield was detected at 116 cM on BTA seven in family three. The IM and CIM approaches detected a QTL at 3 cM on BTA 14 influencing fat yield in family four. Two map positions on BTA 29 were associated with significant variation of milk (0 cM) and fat yield (14 cM) in family seven. These results suggest the presence of one QTL with pleiotropic effects on multiple traits or multiple QTL within the marker interval. Findings from this study could be used in subsequent fine-mapping work and applied to marker-assisted selection of dairy production and health traits.

Detection of quantitative trait loci influencing dairy traits using a model for longitudinal data.

A longitudinal-linkage analysis approach was developed and applied to an outbred population. Nonlinear mixed-effects models were used to describe the lactation patterns and were extended to include marker information following single-marker and interval mapping models. Quantitative trait loci (QTL) affecting the shape and scale of lactation curves for production and health traits in dairy cattle were mapped in three U.S. Holstein families (Dairy Bull DNA Repository families one, four, and five) using the granddaughter design. Information on 81 informative markers on six Bos taurus autosomes (BTA) was combined with milk yield, fat, and protein percentage and somatic cell score (SCS) test-day records. Six percent of the single-marker tests surpassed the experiment-wise significance threshold. Marker BL41 on BTA3 was associated with decrease in milk yield during mid-lactation in family one. The scale and shape of the protein percentage lactation curve in family four varied with BMC4203 (BTA6) allele that the son received from the grandsire. Some map locations were associated with variation in the lactation pattern of multiple traits. In family four, the marker HUJI177 (BTA3) was associated with changes in the milk yield and protein percentage curves suggesting a QTL with pleiotropic effects or multiple QTL in the region. The interval mapping model uncovered a QTL on BTA7 associated with variation in milk-yield pattern in family four and a QTL on BTA21 affecting SCS in family five. The developed approach can be extended to random regressions, covariance functions, spline, gametic and variance component models. The results from the longitudinal-QTL approach will help to understand the genetic factors acting at different stages of lactation and will assist in positional candidate gene research. Identified positions can be incorporated into marker-assisted selection decisions to alter the persistency and peak production or the fluctuation of SCS during a lactation.

Homogeneity of recombination rate within a conserved region on BTA3 that contains QTL.

The D3S20-D3S34-D3S3 region on BTA3 contains quantitative trait loci (QTL) controlling milk production traits. This region also displays extensive conservation of synteny among several species including cattle, humans, mice and sheep. In this study, we evaluated the adjacent intervals D3S20-D3S34 and D3S34-D3S3 for differences in recombination rate (theta) among bulls in order to assess the suitability of population-based estimates of theta for marker assisted selection and to explore the relationship between variation in theta and chromosome breakpoints associated with mammalian evolution. Using sperm typing, thetaD3S20-D3S34 and theta D3S34-D3S3 were estimated for six triply heterozygous bulls. Recombination frequency ranged from 6.2 to 12.5% and from 9.7 to 19.2% for the D3S20-D3S34 and D3S34-D3S3 intervals, respectively. However, significant variation in theta was not detected between bulls for either interval (D3S20-D3S34 chi(2)5 d.f.=2.59, P < 0.90; D3S34-D3S3 chi(2)5 d.f.=3.72, P < 0.75). The observed differences in theta were most readily attributed to differences in allele-specific amplification efficiencies among bulls. Our results suggest that the positions of QTL in this region can be reliably determined from population data and therefore accurate marker-assisted selection can be performed for desirable alleles without concern for variation in theta. Furthermore, when considered with results of earlier studies, these findings support a correlation between the existence of evolutionary breakpoints or chromosome rearrangements and variation in theta.

CD4(+) T lymphocytes from calves immunized with Anaplasma marginale major surface protein 1 (MSP1), a heteromeric complex of MSP1a and MSP1b, preferentially recognize the MSP1a carboxyl terminus that is conserved among strains.

Native major surface protein 1 (MSP1) of the ehrlichial pathogen Anaplasma marginale induces protective immunity in calves challenged with homologous and heterologous strains. MSP1 is a heteromeric complex of a single MSP1a protein covalently associated with MSP1b polypeptides, of which at least two (designated MSP1F1 and MSP1F3) in the Florida strain are expressed. Immunization with recombinant MSP1a and MSP1b alone or in combination fails to provide protection. The protective immunity in calves immunized with native MSP1 is associated with the development of opsonizing and neutralizing antibodies, but CD4(+) T-lymphocyte responses have not been evaluated. CD4(+) T lymphocytes participate in protective immunity to ehrlichial pathogens through production of gamma interferon (IFN-gamma), which promotes switching to high-affinity immunoglobulin G (IgG) and activation of phagocytic cells to produce nitric oxide. Thus, an effective vaccine for A. marginale and related organisms should contain both T- and B-lymphocyte epitopes that induce a strong memory response that can be recalled upon challenge with homologous and heterologous strains. This study was designed to determine the relative contributions of MSP1a and MSP1b proteins, which contain both variant and conserved amino acid sequences, in stimulating memory CD4(+) T-lymphocyte responses in calves immunized with native MSP1. Peripheral blood mononuclear cells and CD4(+) T-cell lines from MSP1-immunized calves proliferated vigorously in response to the immunizing strain (Florida) and heterologous strains of A. marginale. The conserved MSP1-specific response was preferentially directed to the carboxyl-terminal region of MSP1a, which stimulated high levels of IFN-gamma production by CD4(+) T cells. In contrast, there was either weak or no recognition of MSP1b proteins. Paradoxically, all calves developed high titers of IgG antibodies to both MSP1a and MSP1b polypeptides. These findings suggest that in calves immunized with MSP1 heteromeric complex, MSP1a-specific T lymphocytes may provide help to MSP1b-specific B lymphocytes. The data provide a basis for determining whether selected MSP1a CD4(+) T-lymphocyte epitopes and selected MSP1a and MSP1b B-lymphocyte epitopes presented on the same molecule can stimulate a protective immune response.

Development and mapping of microsatellites from a microdissected BTA 11-specific DNA library.

A chromosome-specific library was developed for Bos taurus autosome 11 by chromosome microdissection and microcloning using a bovine primary fibroblast culture, obtained from a t(X;23) heifer, that spontaneously developed a translocation chromosome involving bovine chromosome 11. The library was screened using (AC)12 oligos, positive clones selected, sequenced and primers developed to generate bovine chromosome 11-specific microsatellite markers. This study suggests that chromosome-specific libraries have great potential for development of microsatellite markers for the construction of marker-saturated linkage maps for each chromosome.

Comprehensive linkage map of bovine chromosome 11.

The results of genotypic data contributed to the International Society of Animal Genetics (ISAG) Bovine Chromosome 11 (BTA11) Workshop are presented. Six laboratories contributed a total of 26 199 informative meioses from 80 loci. Thirty-six loci were typed by at least two independent laboratories and were used to construct a consensus linkage map of the chromosome. The remaining loci were subsequently incorporated into a comprehensive map. The sex-averaged consensus map covered 128.9 cM. The female consensus map was 101.2 cM, while the male consensus map was 129.8 cM. The comprehensive sex-averaged map was 134.2 cM and the average genetic distance between loci was 1.72 cM.

Highly conserved regions of the immunodominant major surface protein 2 of the genogroup II ehrlichial pathogen Anaplasma marginale are rich in naturally derived CD4+ T lymphocyte epitopes that elicit strong recall responses.

Genogroup II ehrlichia, including the agent of human granulocytic ehrlichiosis, Ehrlichia phagocytophila, and the bovine pathogen Anaplasma marginale, express a markedly immunodominant outer membrane protein designated major surface protein 2 (MSP2). MSP2 is encoded by a multigene family, resulting in the expression of variant B cell epitopes. MSP2 variants are sequentially expressed in the repeated cycles of rickettsemia that characterize persistent A. marginale infection and control of each rickettsemic cycle is associated with development of a variant-specific IgG response. Importantly, these persistent rickettsemic cycles are controlled at levels 100-1000 times lower than those responsible for clinical disease during acute infection. Control of rickettsemia during persistence could result from an anamnestic Th lymphocyte response to conserved regions of MSP2 that enhances the primary Ab response against newly emergent variants. Comparison of MSP2 variants reveals conserved N and C termini flanking the central, surface-exposed hypervariable region that represents the variant B lymphocyte epitopes. We demonstrate MSP2-specific CD4(+) T lymphocyte recognition of epitopes common to several strains of A. marginale and the related pathogen A. ovis. Furthermore, T lymphocyte lines from three individuals identified six to nine overlapping peptides representing a minimum of four to seven dominant or subdominant epitopes in these conserved N and C termini. Immunodominant peptides induced high levels of IFN-gamma, a cytokine associated with protection against ehrlichia and needed for rapid generation of variant-specific IgG2. The presented data support the potential importance of a strong Th lymphocyte response to invariant MSP2 epitopes in controlling rickettsemia during persistent infection to subclinical levels.