Analysis of nonsynonymous SNPs in candidate genes that influence bovine temperament and evaluation of their effect in Brahman cattle

BACKGROUND

Temperament is an important production trait in cattle and multiple strategies had been developed to generate molecular markers to assist animal selection. As nonsynonymous single nucleotide polymorphisms are markers with the potential to affect gene functions, they could be useful to predict phenotypic effects. Genetic selection of less stress-responsive, temperamental animals is desirable from an economic and welfare point of view.

METHODS AND RESULTS

Two nonsynonymous single nucleotide polymorphisms identified in HTR1B and SLC18A2 candidate genes for temperament were analyzed in silico to determine their effects on protein structure. Those nsSNPs allowing changes in proteins were selected for a temperament association analysis in a Brahman population. Transversion effects on protein structure were evaluated in silico for each amino acid change model, revealing structural changes in the proteins of the HTR1B and SLC18A2 genes. The selected nsSNPs were genotyped in a Brahman population (n = 138), and their genotypic effects on three temperament traits were analyzed: exit velocity, pen score, and temperament score. Only the SNP rs209984404-HTR1B (C/A) showed a significant association (P = 0.0144) with pen score. The heterozygous genotype showed a pen score value 1.17 points lower than that of the homozygous CC genotype.

CONCLUSION

The results showed that in silico analysis could direct the selection of nsSNPs with the potential to change the protein. Non-synonymous single nucleotide polymorphisms causing structural changes and reduced protein stability were identified. Only rs209984404-HTR1B shows that the allele affecting protein stability was associated with the genotype linked to docility in cattle.

Genetic evaluation of crossbred Bos indicus cow temperament at parturition

Cow temperament at parturition may be mostly a measure of aggressiveness. The heritability of cow temperament at parturition in Bos taurus cows has been reported to be low. The objectives of this study were to estimate the heritability of cow temperament at parturition, conduct a genome-wide association analysis of cow temperament at the time of parturition, and estimate the correspondence of cow temperament at the time of parturition with cow productive performance and early-life temperament traits in Bos indicus crossbreds. Cow temperament was assessed from 1 to 5 indicating increasing levels of aggressiveness of cows (937 cows and 4,337 parturitions) from 2005 to 2022. Estimates of heritability and repeatability were 0.12 ± 0.024 and 0.24 ± 0.018. The estimates of proportion of phenotypic variance were 0.13 ± 0.019 and 0.02 ± 0.011 for permanent and maternal permanent environmental components, respectively. Estimates of heritability for maximum lifetime temperament score and proportions of temperament scores >1 were 0.18 ± 0.07 and 0.13 ± 0.072. Within cycles (generations), 2-yr-old cows had lower temperament score means than cows in most other age categories. There were low to moderate positive estimates of unadjusted correlation coefficients (r = 0.22 to 0.29; P < 0.05) of unadjusted temperament score with temperament measured on the same females when they were 8 mo old. There were low to moderate positive estimates of correlation coefficients (r = 0.09 to 0.37; P < 0.05) of unadjusted temperament score with calving rate, weaning rate, weaning weight per cow exposed, and weaning weight per 454 kg cow weight at weaning. Cows with the lowest temperament score had lower (P < 0.05) calving and weaning rate than cows in other temperament categories. Within 3 of 5 cycles, cows with the lowest temperament score (totally docile) had lower (P < 0.05) weaning weight per cow exposed than cows in other temperament categories. There were 2 SNP on BTA 4 associated with maximum lifetime temperament score (FDR < 0.05). The non-genetic influence of a cow's mother was documented in her own temperament measured at the time of calving; this may be a consequence of learned behavior. Less aggressiveness displayed by cows at the time of calving may be accompanied by lower reproductive and maternal performance.

Differential Expression of Circadian Clock Genes in the Bovine Neuroendocrine Adrenal System

Knowledge of circadian rhythm clock gene expression outside the suprachiasmatic nucleus is increasing. The purpose of this study was to determine whether expression of circadian clock genes differed within or among the bovine stress axis tissues (e.g., amygdala, hypothalamus, pituitary, adrenal cortex, and adrenal medulla). Tissues were obtained at an abattoir from eight mature nonpregnant Brahman cows that had been maintained in the same pasture and nutritional conditions. Sample tissues were stored in RNase-free sterile cryovials at -80 °C until the total RNA was extracted, quantified, assessed, and sequenced (NovaSeq 6000 system; paired-end 150 bp cycles). The trimmed reads were then mapped to a Bos taurus (B. taurus) reference genome (Umd3.1). Further analysis used the edgeR package. Raw gene count tables were read into RStudio, and low-expression genes were filtered out using the criteria of three minimum reads per gene in at least five samples. Normalization factors were then calculated using the trimmed mean of M values method to produce normalized gene counts within each sample tissue. The normalized gene counts important for a circadian rhythm were analyzed within and between each tissue of the stress axis using the GLM and CORR procedures of the Statistical Analysis System (SAS). The relative expression profiles of circadian clock genes differed (p < 0.01) within each tissue, with neuronal PAS domain protein 2 (NPAS2) having greater expression in the amygdala (p < 0.01) and period circadian regulator (PER1) having greater expression in all other tissues (p < 0.01). The expression among tissues also differed (p < 0.01) for individual circadian clock genes, with circadian locomotor output cycles protein kaput (CLOCK) expression being greater within the adrenal tissues and nuclear receptor subfamily 1 group D member 1 (NR1D1) expression being greater within the other tissues (p < 0.01). Overall, the results indicate that within each tissue, the various circadian clock genes were differentially expressed, in addition to being differentially expressed among the stress tissues of mature Brahman cows. Future use of these findings may assist in improving livestock husbandry and welfare by understanding interactions of the environment, stress responsiveness, and peripheral circadian rhythms.

Telomere Dynamics in Livestock

Telomeres are repeated sequences of nucleotides at the end of chromosomes. They deteriorate across mitotic divisions of a cell. In Homo sapiens this process of lifetime reduction has been shown to correspond with aspects of organismal aging and exposure to stress or other insults. The early impetus to characterize telomere dynamics in livestock related to the concern that aged donor DNA would result in earlier cell senescence and overall aging in cloned animals. Telomere length investigations in dairy cows included breed effects, estimates of additive genetic control (heritability 0.12 to 0.46), and effects of external stressors on telomere degradation across animal life. Evaluation of telomeres with respect to aging has also been conducted in pigs and horses, and there are fewer reports of telomere biology in beef cattle, sheep, and goats. There were minimal associations of telomere length with animal productivity measures. Most, but not all, work in livestock has documented an inverse relationship between peripheral blood cell telomere length and age; that is, a longer telomere length was associated with younger age. Because livestock longevity affects productivity and profitability, the role of tissue-specific telomere attrition in aging may present alternative improvement strategies for genetic improvement while also providing translational biomedical knowledge.

Effect of ewe birth litter size and estimation of genetic parameters on ewe reproductive life traits

Ewe lifetime productivity has economic implications for producers because shorter lifetime productivity results in less profit. Productive years of ewes from extensive, range-based systems of the United States West are generally less than ewes from more temperate regions of the United States. Accordingly, ewes from range-based systems, especially those employing shed-lambing strategies, have been selected for increased litter size to offset decreased lifetime productivity. However, the relationship of the ewe's birth litter size (ELSB) has not been considered a potential contributor to lifetime productivity. Longevity (number of productive years, n = 1 per population) and stayability (probability to survive to the next age; ages 2-7 years, n = 6 per population) were investigated to understand ELSB effects on productive life. Columbia, Polypay, Rambouillet, and Targhee breeds were used in this study. Across-breed (n = 11 550) and within-breed (Columbia, n = 4 398; Polypay, n = 4 534; Rambouillet, n = 5 922; Targhee, n = 6 482) analyses were used. Depending on the population, records spanned from 1950 to 2008, where ewe's birth year was included as a fixed effect in the animal model using restricted maximum likelihood estimation procedures. Fixed effects investigated included ELSB (single-, twin-, or triplet-born) and ewe breed (across-breed analyses only). Regardless of trait or population used, heritability ranged from 0.06 ± 0.02 to 0.34 ± 0.03, where stayability at younger ages had the highest estimates. The breed effect was significant in all across-breed analyses (0.0001 ≤ P ≤ 0.038; n = 7), where Polypay, a breed selected for accelerated lambing and increased fertility, averaged shorter productive life or had a lower probability of survival to the next age compared with other breeds (longevity: 0.009 ≤ P ≤ 0.223; stayability: 0.000 ≤ P ≤ 0.842). The ELSB was significant in 60% (n = 5) and 37% (n = 30) of longevity and stayability analyses, respectively. Except for Targhee, all analyses showed ewes born in smaller litter sizes were associated with longer productive lives or higher probability of surviving to the next age, particularly in across-breed analyses (e.g., longevity: single- vs twin-born ewes, P = 0.004; vs triplet-born ewes, P = 0.003). This study provides evidence that increasing prolificacy in ewes from extensive, range-based production systems may impact productive life. Due to the low heritability of these traits, additional investigation into modeling these traits with dominance effects and litter size needs to be conducted.

Precision and Accuracy of Field Versus Laboratory Bioassay Insecticide Efficacy for the Control of Immature Bemisia tabaci

Ecotoxicological studies often result in reports on the limitation and sometime failures of biological assay data to predict field response to similar treatments. Nevertheless, it is widely accepted that controlled bioassays can better quantify the specific mortality response of a target pest species to a specific toxin. To quantify the relationship between whitefly bioassay and field response data, we evaluated a controlled laboratory bioassay and a concurrent cucurbit field trial method to assess insecticide efficacy for controlling the sweetpotato whitefly, Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae). This was based on oviposition and nymphal development. We specifically tested the assumptions that a maximum dose bioassay would more precisely measure insecticide efficacy as compared with a comparable field spray test evaluation, and the response would be equal between the bioassay and the field as a measure of control accuracy for both adult oviposition and development of nymphal stages. To make a direct comparison, we tested the same whitefly population subsamples from 352 plots in eight cucurbit field experiments in Georgia, USA, in 2021 and 2022. The bioassays provide significantly precision for estimating proportional whitefly response. As expected, treatment-specific nonequivalence in immature whitefly counts between the bioassay and field, i.e., a lack of accuracy, only occurred with insecticides that were not highly toxic to all growth stages of whiteflies.

Comparison of Telomere Length in Age-Matched Primiparous and Multiparous Brahman Cows

Physiological and psychological stressors have been associated with the attrition of telomeres, which are the protective caps of chromosomes. This study compares the telomere length (TL) in 4-year-old Brahman cows grouped by the first parity (n = 8) and the second parity (n = 11). The cows were bled via jugular venipuncture, weighed, and had their body condition scores recorded at Day -28 prior to calving and at Day + 7 and Day + 28 post-calving. The duration of labor (Dlabor) and parturition ease were recorded. The peripheral leukocytes were isolated, the leukocyte blood count with differential was recorded, and the genomic DNA was extracted. The relative quantity of telomere products, which is proportional to the average TL, was determined via multiplex quantitative PCR using the ratio (T/S ratio) of bovine telomere and β-globulin DNA. Standards of the bovine telomere (10¹²-10⁷ dilution series) and β-globulin (10⁹-10⁴ dilution series) genes were utilized to produce relative copy numbers. The samples were assayed in triplicate and were included if the triplicate C_q difference was less than 0.25 cycles. The parity was the fixed effect, and the random effects included the sire and day repeated with the cow as the subject. Statistical significance was not observed in the leukocyte number or type (p > 0.1). A reduction in the TL of approximately 9225 telomeric copies was found between Parity 1 and Parity 2 (p = 0.02). A trend was found between the TL and Dlabor (p = 0.06). The stress of parturition and raising the first calf of a cow's life may be responsible for TL attenuation. Parity may be considered a stressor of cow longevity.

Determining Field Insecticide Efficacy on Whiteflies with Maximum Dose Bioassays

We conducted a rapid bioassay method to assess insecticide efficacy for controlling adult sweetpotato whitefly Bemisia tabaci in squash and cucumber crops before insecticide applications. The study aimed to evaluate the accuracy of a 24-hour laboratory bioassay in determining maximum dose insecticide efficacy in the field. Ten insecticides were evaluated using leaf-dip bioassays, and their effectiveness was tested across eight cucurbit field experiments in Georgia, USA, during the 2021 and 2022 field seasons. The maximum dose, defined as the highest labeled rate of an insecticide diluted in the equivalent of 935 L ha^-1 of water, was used for all bioassays. Adult survival observed in the bioassay was compared to adult field count-based survival 24 h after treatment. A low concentration (1/10 rate) was used for imidacloprid, flupyradifurone, pyriproxyfen, and cyantraniliprole to assess insecticide tolerance in the whitefly population. Overall, significant positive correlation between laboratory bioassay and field efficacy was reported, explaining 50-91% of the observed variation. The addition of the low dosage was helpful, indicating that no rate response was consistent with susceptibility to the tested insecticide, while a rate response was associated with a loss of susceptibility between 2021 and 2022.

Genetic Variants and Their Putative Effects on microRNA-Seed Sites: Characterization of the 3' Untranslated Region of Genes Associated with Temperament

The 3' untranslated region has an important role in gene regulation through microRNAs, and it has been estimated that microRNAs regulate up to 50% of coding genes in mammals. With the aim of allelic variant identification of 3' untranslated region microRNA seed sites, the 3' untranslated region was searched for seed sites of four temperament-associated genes (CACNG4, EXOC4, NRXN3, and SLC9A4). The microRNA seed sites were predicted in the four genes, and the CACNG4 gene had the greatest number with 12 predictions. To search for variants affecting the predicted microRNA seed sites, the four 3' untranslated regions were re-sequenced in a Brahman cattle population. Eleven single nucleotide polymorphisms were identified in the CACNG4, and eleven in the SLC9A4. Rs522648682:T>G of the CACNG4 gene was located at the predicted seed site for bta-miR-191. Rs522648682:T>G evidenced an association with both exit velocity (p = 0.0054) and temperament score (p = 0.0097). The genotype TT had a lower mean exit velocity (2.93 ± 0.4 m/s) compared with the TG and GG genotypes (3.91 ± 0.46 m/s and 3.67 ± 0.46 m/s, respectively). The allele associated with the temperamental phenotype antagonizes the seed site, disrupting the bta-miR-191 recognition. The G allele of CACNG4-rs522648682 has the potential to influence bovine temperament through a mechanism associated with unspecific recognition of bta-miR-191.

A Review of Successes and Impeding Challenges of IoT-Based Insect Pest Detection Systems for Estimating Agroecosystem Health and Productivity of Cotton

Using artificial intelligence (AI) and the IoT (Internet of Things) is a primary focus of applied engineering research to improve agricultural efficiency. This review paper summarizes the engagement of artificial intelligence models and IoT techniques in detecting, classifying, and counting cotton insect pests and corresponding beneficial insects. The effectiveness and limitations of AI and IoT techniques in various cotton agricultural settings were comprehensively reviewed. This review indicates that insects can be detected with an accuracy of between 70 and 98% using camera/microphone sensors and enhanced deep learning algorithms. However, despite the numerous pests and beneficial insects, only a few species were targeted for detection and classification by AI and IoT systems. Not surprisingly, due to the challenges of identifying immature and predatory insects, few studies have designed systems to detect and characterize them. The location of the insects, sufficient data size, concentrated insects on the image, and similarity in species appearance are major obstacles when implementing AI. Similarly, IoT is constrained by a lack of effective field distance between sensors when targeting insects according to their estimated population size. Based on this study, the number of pest species monitored by AI and IoT technologies should be increased while improving the system's detection accuracy.

Inter-Individual Variation in DNA Methylation Patterns across Two Tissues and Leukocytes in Mature Brahman Cattle

Quantifying the natural inter-individual variation in DNA methylation patterns is important for identifying its contribution to phenotypic variation, but also for understanding how the environment affects variability, and for incorporation into statistical analyses. The inter-individual variation in DNA methylation patterns in female cattle and the effect that a prenatal stressor has on such variability have yet to be quantified. Thus, the objective of this study was to utilize methylation data from mature Brahman females to quantify the inter-individual variation in DNA methylation. Pregnant Brahman cows were transported for 2 h durations at days 60 ± 5; 80 ± 5; 100 ± 5; 120 ± 5; and 140 ± 5 of gestation. A non-transport group was maintained as a control. Leukocytes, amygdala, and anterior pituitary glands were harvested from eight cows born from the non-transport group (Control) and six from the transport group (PNS) at 5 years of age. The DNA harvested from the anterior pituitary contained the greatest variability in DNA methylation of cytosine-phosphate-guanine (mCpG) sites from both the PNS and Control groups, and the amygdala had the least. Numerous variable mCpG sites were associated with retrotransposable elements and highly repetitive regions of the genome. Some of the genomic features that had high variation in DNA methylation are involved in immune responses, signaling, responses to stimuli, and metabolic processes. The small overlap of highly variable CpG sites and features between tissues and leukocytes supports the role of variable DNA methylation in regulating tissue-specific gene expression. Many of the CpG sites that exhibited high variability in DNA methylation were common between the PNS and Control groups within a tissue, but there was little overlap in genomic features with high variability. The interaction between the prenatal environment and the genome could be responsible for the differences in location of the variable DNA methylation.

DNA methylation patterns and gene expression from amygdala tissue of mature Brahman cows exposed to prenatal stress

Prenatal stress can alter postnatal performance and temperament of cattle. These phenotypic effects may result from changes in gene expression caused by stress-induced epigenetic alterations. Specifically, shifts in gene expression caused by DNA methylation within the brain’s amygdala can result in altered behavior because it regulates fear, stress response and aggression in mammals Thus, the objective of this experiment was to identify DNA methylation and gene expression differences in the amygdala tissue of 5-year-old prenatally stressed (PNS) Brahman cows compared to control cows. Pregnant Brahman cows (n = 48) were transported for 2-h periods at 60 ± 5, 80 ± 5, 100 ± 5, 120 ± 5, and 140 ± 5 days of gestation. A non-transported group (n = 48) were controls (Control). Amygdala tissue was harvested from 6 PNS and 8 Control cows at 5 years of age. Overall methylation of gene body regions, promoter regions, and cytosine-phosphate-guanine (CpG) islands were compared between the two groups. In total, 202 genes, 134 promoter regions, and 133 CpG islands exhibited differential methylation (FDR ≤ 0.15). Following comparison of gene expression in the amygdala between the PNS and Control cows, 2 differentially expressed genes were identified (FDR ≤ 0.15). The minimal differences observed could be the result of natural changes of DNA methylation and gene expression as an animal ages, or because this degree of transportation stress was not severe enough to cause lasting effects on the offspring. A younger age may be a more appropriate time to assess methylation and gene expression differences produced by prenatal stress.

9 Evaluating the Utility of Individual Cattle Brush Use Behavior as a Novel Behavioral Phenotype Regarding Productivity, Temperament, and Feeding Behavior of Bos Indicus Steers and Heifers Housed in Dry Lots

Little is understood about the relationships among environmental enrichment, temperament, productivity, and feeding behavior in cattle. The objective of this study was to examine the effect of divergent behavioral phenotypes for brush usage on temperament, productivity, and feeding behavior patterns in beef cattle. Four pens of Brahman cattle (n = 2 pens with 15 heifers/pen, and n = 2 pens with 19 steers/pen) were housed in dry lots each equipped with four electronic feedbunks (GrowSafe Systems). Animals were video recorded throughout an 85d commercial gain test. Video recordings were decoded on d 1, 2, 4, 8, 16, 32, and 64 using continuous observations for 15 minutes every 30 minutes from 08:00 to 20:00 to determine the total daily brush usage for focal animals that were selected based on d -19 exit velocities. Exit velocity and body weight were measured on d -19, 0, 1, 25, 50, and 70 and feeding behavior was monitored daily using the GrowSafe 4000E system. Animals were classified within gender into one of three phenotypes (high, medium, and low brush usage) based on ±0.5SD from the brush usage mean. A General Linear Mixed Model (PROC GLIMMIX) with brush use phenotype, sex and the interaction of these as fixed effects and pen as random effect was used for data analysis. For high and medium brush use phenotypes, steers had greater final BW than heifers (P = 0.0005), while heifer and steer final body weight did not differ when brush usage was low. High brush use steers had greater ADG (P = 0.0002) than high and medium brush use heifers, and low brush use steers and heifers. This same pattern was observed for DMI (P = 0.0003). Heifers were more temperamental than steers (P = 0.0008), and high and medium brush use heifers had faster exit velocities than low brush use steers. These results suggest that monitoring animal’s pleasurable behavior (brush use), rather than fear response alone, may be a useful behavioral indicator of cattle productivity.

A Target Site Mutation Associated With Diamide Insecticide Resistance in the Diamondback Moth Plutella xylostella (Lepidoptera: Plutellidae) is Widespread in South Georgia and Florida Populations

Laboratory colonies of diamondback moth (DBM) larvae were established from larvae collected from four sites in Georgia and Florida where diamide, specifically chlorantraniliprole, insecticide resistance was recently documented. Based on dose-response experiments, these colonies exhibited 109- to 4,298-fold resistance to chlorantraniliprole, compared to a commercially available susceptible control colony. Colonies exhibited 50- to 107-fold resistance to another diamide, cyantraniliprole, based on similar dose-response experiments. All colonies were screened for the presence of four known mutations in the ryanodine receptor (RyR), the target of diamide insecticides, previously associated with resistance in Asian DBM populations. One mutation, G4946E, was identified in colonies from all four field sites, but not the susceptible control colony. Three additional RyR target site mutations, E1338D, Q4594L, and I4790M, were not identified in any of the screened samples. The estimated allele frequency of the G4946E mutation in these colonies ranged from 32 to 90%. These data are consistent with recently reported chlorantraniliprole control failures in Georgia and Florida. It is likely that the G4946E mutation is currently an important contributing factor to chlorantraniliprole resistance in Georgia and Florida DBM populations.

Skeletal Muscle Expression of Actinin-3 (ACTN3) in Relation to Feed Efficiency Phenotype of F2Bos indicus - Bos taurus Steers

In this study, actinin-3 (ACTN3) gene expression was investigated in relation to the feed efficiency phenotype in Bos indicus - Bos taurus crossbred steers. A measure of relative feed efficiency based on residual feed intake relative to predictions from the NRC beef cattle model was analyzed by the use of a mixed linear model that included sire and family nested within sire as fixed effects and age, animal type, sex, condition, and breed as random effects for 173 F₂ Nellore-Angus steers. Based on these residual intake observations, individuals were ranked from most efficient to least efficient. Skeletal muscle samples were analyzed from 54 steers in three groups of 18 (high efficiency, low efficiency, and a statistically average group). ACTN3, which encodes a muscle-specific structural protein, was previously identified as a candidate gene from a microarray analysis of RNA extracted from muscle samples obtained from a subset of steers from each of these three efficiency groups. The expression of ACTN3 was evaluated by quantitative reverse transcriptase PCR analysis. The expression of ACTN3 in skeletal muscle was 1.6-fold greater in the inefficient steer group than in the efficient group (p = 0.007). In addition to expression measurements, blocks of SNP haplotypes were assessed for breed or parent of origin effects. A maternal effect was observed for ACTN3 inheritance, indicating that a maternal B. indicus block conferred improved residual feed efficiency relative to the B. taurus copy (p = 0.03). A SNP haplotype analysis was also conducted for m-calpain (CAPN2) and fibronectin 1 (FN1), and a significant breed effect was observed for both genes, with B. indicus and B. taurus alleles each conferring favorable efficiency when inherited maternally (p = 0.03 and p = 0.04). Because the ACTN3 structural protein is specific to fast-twitch (type II) muscle fibers and not present in slow-twitch muscle fibers (type I), muscle samples used for expression analysis were also assayed for fiber type ratio (type II/type I). Inefficient animals had a fast fiber type ratio 1.8-fold greater than the efficient animals (p = 0.027). Because these fiber-types exhibit different metabolic profiles, we hypothesize that animals with a greater proportion of fast-twitch muscle fibers are also less feed efficient.

First report of turnip yellows virus infecting cabbage (Brassica oleracea var. capitata) in the USA

During the spring of 2021, cabbage (Brassica oleracea var. capitata) planted in the research farm at the University of Georgia, Tifton, exhibited leaf distortion, yellow and purple discoloration at the leaf margin of older leaves, and severe stunting. Symptoms were present on nearly 30% of the plants in the field. To identify the potential agents associated, leaf tissues from two symptomatic plants were sent for high throughput sequencing (HTS) of small RNA (sRNA; DNB sequencing, SE read 1x75bp) to Beijing Genomics Institute, China. From each sample, ~ 18 million raw reads were generated. The reads with poor quality and adapter sequences were removed using CLC Genomics Workbench 21.2 (Qiagen, Germantown, MD). Of the total reads, 2,093 and 3,889 reads aligned to the genome of turnip yellows virus (TuYV) in samples one and two, respectively. Reads of turnip mosaic virus (TuMV) were also detected (data not shown). Partial sequences of TuYV assembled from samples one and two showed 89.5% and 89.9% match and 86% and 93% coverage, respectively, with the genome of the type isolate of TuYV (NC_003743) from the United Kingdom. To confirm the presence of TuYV in the samples collected from the same location, specific primers were designed targeting the P0 region (FP- 5'ACAAAAGAAACCAG- GAGGGAATCC3'; RP-5'GCCTTTTCATACAAACATTTCGGTG3') and coat protein (CP) region (FP-5'GTTAATGAATACGGTCGTGGGTAG3'; RP-5'ATTCTGAAAGAACCAGCT- ATCGATG3') of the virus. Eight of 20 (40%) symptomatic samples were determined to be infected with TuYV based on the amplification of expected size products of the P0 (786 nt) and the CP gene (581 nt) in reverse transcription-PCR (RT-PCR). All samples were also tested for the presence of TuMV by RT-PCR as in Sanchez et al. (2003), but none tested positive despite being identified in HTS. Symptoms on samples from which eithervirus could not be detected indicates the involvement of other factors and would require further studies. The partial P0 and CP gene amplicons of TuYV from two samples each were Sanger sequenced bi-directionally at Genewiz (South Plainfield, NJ) and confirmed as TuYV using BLASTn. The partial CP gene sequences from two samples shared 98.7% nucleotide sequence identity with each other and 88.0% (OK349421) and 87.1% (OK349422) identity with the type isolate. The partial P0 gene sequences (OK349423 and OK349424) shared 99.6% nucleotide sequence identity with each other and 92.2% identity with the type isolate. TuYV, formerly known as beet western yellows virus (BWYV) (Mayo, 2002), genus Palerovirus, family Solemoviridae (Walker et al., 2021), is transmitted persistently by aphids (Stevens et al., 2008), and is distributed throughout temperate regions of the world (Kawakubo et al., 2021). TuYV has a wide host range, including brassica, vegetables and weeds (Stevens et al., 2008). However, losses have been reported primarily on canola (B. napus) in Australia (Jones, 2007) and Europe (Stevens et al., 2008). On cabbage, TuYV infections have been reported from China (Zhang et al., 2016), Serbia (Milošević et al., 2020) and the Philippines (Buxton-Kirk et al, 2020). TuYV (BWYV) has been found infecting shepherd's purse (Capsella bursa-pastoris) in California (Falk and Duffus, 1984), but there are no reports of the virus from any cultivated crops in the USA. To our knowledge, this is the first report of TuYV in cabbage in the USA. More studies are needed to understand its occurrence and impact on cabbage crops in Georgia as well as other regions in the USA.

30 Calf Performance and Sex Ratios in Second Generation Reciprocal Nellore-angus Crosses

The objectives of this study were to analyze calf weight and sex distribution differences in reciprocal F2 Nellore-Angus calves (n = 559) born 2010–2015 and calves out of these F2 females through 4 yr age (n = 469) born 2012–2018. The F2 animals were produced by 4 F1 parent type combinations from Nellore-sired (NA) vs. Angus-sired (AN) sires and dams. Weight traits were analyzed using mixed models, and sex distribution was analyzed through PROC GLIMMIX in SAS. Sire type and dam type affected birth weight (P &lt; 0.05) but not weaning weight of the F2 calves. However, the interaction between sire type and dam type accounted for weaning weight differences (P = 0.038). Birth weight averaged 3.5 kg heavier from NA vs. AN sires, and 2.4 kg heavier from NA vs. AN dams. Weaning weight displayed a different pattern than observed in birth weight with the heaviest calves being NAAN (224.5 ± 3.7 kg) and lightest being ANAN (215.9 ± 3.9 kg). Sire type accounted for important variation (P = 0.036) in the calf sex distribution at birth among F2 calves. Percentage males at birth exceeded the expected ratio for NA-sired calves (66.6%, n = 168) but not for AN-sired calves (49.9%, n = 105). The sire and dam type main effects in the F2 females did not influence weight traits or sex distribution in their calves (from predominantly Angus sires). However, the interaction of sire type x dam type x sex showed a trend in birth weight (P = 0.080) and importance for weaning weight (P = 0.010) where only from ANNA and NANA dams were males heavier than females. These results suggest that some of the calf weight reciprocal differences typically observed among F1 Bos indicus-Bos taurus crosses can persist in later generations.

91 Assessing Variability of Herd Sire Economic Value for Commercial Operations

Number of calves (n = 596) and weaning calf value (n = 574) from natural service, multi-sire breeding groups were evaluated in a crossbred research herd. Sires and dams that were F1 crosses of Angus (A) and Nellore (N) were mated annually during 60 to 90-day breeding seasons to produce spring-born calves in 2009–2015. Numbers of sires used annually varied from 5 to 9; numbers of females exposed per bull ranged from 14 to 22 across years. Bulls were pastured together throughout the year as well as breeding seasons. Sires were identified based on calf DNA genotyping. Calf number, birth date, birth weight, weaning weight, and economic value were determined per bull annually. Calf value was based on weaning weight and reported prices for respective weaning dates. Calf prices were calculated separately for steers and heifers relative to 22.7 kg (50 lb) increments for corresponding weight class and year. Mixed model analyses were conducted that included fixed categorical effects of calf birth year, type of F1 sire (A-sired vs. N-sired), sire nested within type, calf sex, and the interaction of F1 sire type with calf sex. Covariates of Julian birth date and calf weaning age were included for birth weight and weaning weight, respectively. Large differences in calf numbers and performance (P &lt; 0.05) were observed. Weaned calves produced annually per sire ranged from 0 to 48. Average annual calf performance per sire ranged from 28.4 to 50.8 kg for birth weight and 146.5 to 249.0 kg for weaning weight. Annual economic value per sire ranged from $0 to $30,870 when considering half of each calf’s value as attributed to the sire. Commercial producers should consider potential sire variability for calf numbers, birth date distribution, and ratio of female-to-male calves in combination with calf performance for economic assessments.

Analysis of Polycerate Mutants Reveals the Evolutionary Co-option of HOXD1 for Horn Patterning in Bovidae

In the course of evolution, pecorans (i.e., higher ruminants) developed a remarkable diversity of osseous cranial appendages, collectively referred to as “headgear,” which likely share the same origin and genetic basis. However, the nature and function of the genetic determinants underlying their number and position remain elusive. Jacob and other rare populations of sheep and goats are characterized by polyceraty, the presence of more than two horns. Here, we characterize distinct POLYCERATE alleles in each species, both associated with defective HOXD1 function. We show that haploinsufficiency at this locus results in the splitting of horn bud primordia, likely following the abnormal extension of an initial morphogenetic field. These results highlight the key role played by this gene in headgear patterning and illustrate the evolutionary co-option of a gene involved in the early development of bilateria to properly fix the position and number of these distinctive organs of Bovidae.

Improving Whitefly Management

Whiteflies (Hemiptera: Aleyrodidae), especially the sweetpotato or cotton whitefly, Bemisia tabaci (Gennadius), are among the most destructive and difficult to manage polyphagous insect pests around the globe [...].

Evaluation of Prenatal Transportation Stress on the Number of Pituitary Corticotrophs in Mature Brahman Cows

Prenatal transportation stress (PNS) results in calves that are more temperamental and have greater circulating concentrations of cortisol compared to control calves. The objective of this experiment was to evaluate whether PNS alters the number of pituitary corticotrophs in mature Brahman cows. We hypothesized that the increased circulating cortisol concentrations previously characterized in this bovine model is associated with developmental changes in the anterior pituitary leading to an increased number of corticotrophs. Pregnant Brahman cows (n = 48) were transported in trailers for 2-hour periods at 60±5, 80±5, 100±5, 120±5, and 140±5 days of gestation. Non-transported pregnant cows (n = 48) were designated as the Control group. Control and PNS offspring heifers were managed together under the same environmental conditions. At approximately 5 yr of age, randomly selected non-pregnant cows (Control, n = 8; PNS, n = 6) were humanely harvested and the whole pituitaries were collected. Pituitaries were weighed, fixed in paraformaldehyde, serially dehydrated with graded ethanol, embedded in paraffin blocks, and cut into 5-μm sections. Immunohistochemistry was performed to detect cells expressing adrenocorticotropic hormone (ACTH) as a marker for corticotrophs. Three comparable sections from the midsagittal plane from each animal were processed using an ovine ACTH-specific antibody (Dr. A.F. Parlow, NIDDK). Five fields of view were analyzed per section (15 fields per animal). Anterior pituitary gland weight did not differ (P &gt; 0.10) between groups (Control = 2.11 ± 0.12 g; PNS = 2.10 ± 0.15 g). The mean number of ACTH-positive cells between control (531±56 cells/section) and PNS cows (477±49 cells/section) also did not differ (P &gt; 0.10). In conclusion, the number of pituitary corticotrophs in mature Brahman cows was not affected by prenatal transportation stress and is likely not a mechanism mediating the increased circulating cortisol concentrations seen in this bovine model of fetal programming.

Evaluation of Peripheral Blood Leukocytes as a Surrogate for RNA-Seq Detection of Gene Expression in Stress Axis Tissues of Mature Brahman Cows

Longitudinal studies of gene expression where sampling certain organ systems is precluded may become more feasible if a more accessible, representative biological sample could be identified. The objective of this study was to determine whether peripheral blood leukocytes (PBL) could be an effective surrogate for assessing gene expression differences within neuroendocrine tissues that regulate stress responses. Specifically, leukocyte samples and stress axis tissues, including the paraventricular region of the hypothalamus, anterior pituitary gland, adrenal cortex, and adrenal medulla were isolated from humanely slaughtered 5-yr-old Brahman cows (n = 8). RNA was extracted from each sample and gene expressions were calculated using bioinformatic approaches of RNA-Seq analysis. Differentially expressed genes (DEG) between adrenal axis tissue RNA and PBL RNA were identified using the edgeR Package from Bioconductor, R. More than 3,000 DEGs were identified in each tissue relative to the PBLs. Overall, higher gene expression levels were observed in the hypothalamus vs PBL, while the remaining tissues had reduced expression compared to PBL (Table 1). Pearson correlation values were calculated from the mean normalized gene count in the PBL and tissues for genes (n = 17) important for stress response and metabolism to assess PBL as a possible surrogate for gene expression analyses (Table 2). Only two genes, oxytocin and 11β-hydroxylase, were positively correlated (P &lt; 0.05) between PBL and specific tissue gene expression. As typically observed among different cell types, gene expression differences exist between the PBL and various tissues of the HPA axis. Therefore, it remains to be determined if gene expression in the surrogate PBL will shift in a direction similar to that of specific adrenal axis tissue of mature Brahman cows exposed to stressors during key developmental periods.

The Potential Use of DNA Methylation Patterns from Peripheral Blood Leukocytes as a Surrogate for Stress Axis Tissues in Mature Brahman Cows

DNA methylation (DNAm) patterns are tissue specific and aid in tissue specific gene expression changes. The use of DNAm patterns from peripheral blood leukocytes (PBL) as a surrogate for patterns in other tissues is common, especially in longitudinal studies when sampling of tissues is not plausible. Thus, the objective of this study was to investigate the suitability of using DNAm patterns of PBL as a surrogate for the DNAm patterns in neuroendocrine tissues responsible for stress responses and energy metabolism. Samples from the paraventricular region of the hypothalamus, anterior pituitary gland, adrenal cortex, and the adrenal medulla were harvested from 5-yr-old Brahman cows (n = 8) and DNA was extracted from each sample. Methylation was assessed using reduced representation sodium bisulfite sequencing and differentially methylated regions (DMR) between the PBL DNA and tissue DNA were identified using EdgeR from Bioconductor, R. Analysis revealed over 15,000 DMRs located within promoter regions of genes in each tissue, with the majority of the sites having increased methylation in the PBL (Table 1). To further evaluate the use of PBL DNA as a surrogate, Pearson correlation values were calculated for genes (n = 20) pertinent to each respective tissue using the mean methylation of the specific gene in the PBL and in the tissue (Table 2). Three correlations were significant (P ≤ 0.05), two of which were negative. The sizable differences indicate that DNA methylation patterns from PBL do not compare well to patterns from hypothalamic, pituitary, adrenal cortex, and adrenal medulla tissues from 5-yr-old Brahman cows. This is especially the case for the majority of the specific genes examined in this study. Whether DNAm in the surrogate PBL will shift in a direction similar to that of specific tissues of Brahman cows exposed to stressful stimuli during developmental periods remains to be determined.

Differential Expression of Circadian Clock Genes in the Bovine Neuroendocrine Adrenal System

The primary objective of this investigation was to determine whether circadian clock genes were differentially expressed within or among bovine hypothalamic paraventricular nucleus (PVN), anterior pituitary gland (AP), adrenocortical (AC) and adrenomedullary (AM) tissues. The PVN, AP, AC, and AM were isolated from 5-yr-old Brahman cows (n = 8) harvested humanely at an abattoir between 0800-1100 h. Expression of target genes in each sample was evaluated via RNA-sequencing analyses. Gene counts were normalized using the trimmed mean of M values (TMM) method in the edgeR Package from Bioconductor, R. The normalized gene counts of genes important for circadian rhythm were statistically analyzed using the GLM Procedure of SAS. The genes analyzed were circadian locomotor output cycles protein kaput (CLOCK), cryptochrome circadian regulator 1 and 2 (CRY1 and CRY2), aryl hydrocarbon receptor nuclear translocator like (ARNTL), period circadian regulator 1 and 2 (PER1 and PER2), neuronal PAS domain protein 2 (NPAS2), and nuclear receptor subfamily 1 group D member 1 (NR1D1). Overall, relative expression profiles of clock genes differed (P &lt; 0.01) within each tissue with PER1 having greater expression in all tissues (P &lt; 0.01). Within the PVN expression of CLOCK, CRY1, ARNTL, and PER2 was less than that of CRY2, NPAS2, and NR1D1 (P &lt; 0.01). In the AP, with the exception of PER1, no other clock gene differed in degree of expression. In the AC, expression of CLOCK and NPAS2 was greater than CRY1, ARNTL, PER2, and NR1D1 (P &lt; 0.05), whereas CRY2 expression exceeded only CRY1 (P &lt; 0.05). Within the AM, CLOCK and CRY2 expression was greater than CRY1 and ARNTL (P &lt; 0.05). Overall, clock gene expression among tissues differed (P &lt; 0.01) for each individual clock gene. The AC and AM had similar clock gene expression, except expression of CRY2 and PER2 was greater in AM (P &lt; 0.05). The AC and AM had greater expression of CLOCK than the PVN and AP (P &lt; 0.01), with PVN having greater expression than AP (P &lt; 0.01). The AP had greater expression of NPAS2, followed by PVN, with the least expression in the AC and AM (P &lt; 0.01). Both PVN and AP had greater CRY1 and NR1D1 expression than AC or AM (P &lt; 0.01). The AP had greater PER1 expression than PVN, AC, and AM (P &lt; 0.01), whereas PVN, AC, and AM had greater ARNTL expression than AP (P &lt; 0.05). Both AP and AM had greater expression of PER2 than PVN or AC (P &lt; 0.01). The PVN had greater expression of CRY2 than the AP, AC, and AM (P &lt; 0.01). These results indicated that within each tissue the various clock genes were expressed in different quantities. Also, the clock genes were expressed differentially among the tissues of the bovine neuroendocrine adrenal system. Temporal relationships of these genes with the primary endocrine products of these tissues should be investigated to define the roles of peripheral clock genes in regulation of metabolism and health.

Distributional characterizations and testing for differences of relatedness and inbreeding of a subpopulation of American Hereford bulls

Beta distributions are characterized by two determining parameters and a parameter space from 0 to 1, and may be useful for examining population genetic parameters such as the relationship or inbreeding coefficients. Often subpopulations exist within breeds that are congregated around particular lineages of cattle or ancestors that breeders value. These subpopulations are more related to each other than to the majority of other animals; they may have higher inbreeding as well. Value may be added to these subpopulations because of their relatedness with important or renowned ancestors. The objectives of this work were to compare the relatedness and inbreeding of a group of 26 modern bulls from a subpopulation of the American Hereford breed relative to 1) 30 males with the most descendants present in the pedigree, 2) 15 renowned American Hereford bulls considered important individuals in the breed's history, and 3) 19 prominent subpopulation male ancestors. Conformance of the mean relationship coefficients of the bulls with the three groups and the mean inbreeding coefficient with all pedigree animals to beta distributions was assessed by 1) visually determining the parameters of the beta distributions based on the entire pedigree, 2) testing the mean relationship coefficient or inbreeding coefficient of the group of subpopulation bulls for its positional inclusion in those distributions, and 3) bootstrap sampling methodology. The mean relationship coefficients of the 26 Trask bulls with the 30 bulls with the most descendants, the 15 renowned ancestors, and the 19 Trask male ancestors were 0.15, 0.132, and 0.208, respectively. Testing of these means in beta distributions indicated that the group of 26 Trask bulls were no more related to the three groups of bulls than all of the animals in the pedigree (0.06 < P < 0.25). Bootstrap sampling indicated that the 26 bulls were more related to the three groups of male ancestors than the remainder of the animals in the pedigree (P < 0.0001). The mean inbreeding coefficient of the 26 bulls (0.13) did not differ from the overall inbreeding coefficient (0.056) when tested using a beta distribution; however, bootstrap sampling indicated otherwise (P < 0.0001). Results may indicate the inadequacy of visually parameterizing a beta distribution. Quantification of pedigree relatedness of a group of animals to key ancestors, especially with no DNA available, may add value to that group and individuals.

145 President Oral Presentation Pick: Prenatal stress increases skeletal muscle mitochondrial volume density and function in yearling Brahman calves

We previously determined that prenatal stress (PNS) differentially affected methylation of DNA from leukocytes of 28-d-old calves. Specifically, COX14 (cytochrome c oxidase (COX) assembly factor) and CKMT1B (mitochondrial creatine kinase U-type) were hypomethylated and COA5 (COX assembly factor 5), COX5A (COX subunit 5A), NRF1 (nuclear respiratory factor 1), and GSST1 (glutathione S-transferase theta-1) were hypermethylated in PNS compared to non-PNS calves (P ≤ 0.05). Our current objective was to test the hypothesis that PNS exhibit impaired mitochondrial function and greater oxidative stress than non-PNS calves. Blood and longissimus dorsi muscle samples were collected from yearling Brahman calves whose mothers were stressed by 2 h transportation at 60, 80, 100, 120, and 140 days of gestation (PNS; 8 bulls, 6 heifers) and non-PNS calves (4 bulls, 6 heifers). Serum was evaluated for the stress hormone, cortisol, and muscle damage marker, creatine kinase; muscle was analyzed for mitochondrial volume density and function by citrate synthase (CS) and COX activities, respectively, concentration of malondialdehyde, a lipid peroxidation marker, and activity of the antioxidant, superoxide dismutase (SOD). Data were analyzed using mixed linear models with treatment and sex as fixed effects. Serum cortisol was numerically higher in PNS than non-PNS calves but was not statistically different. Muscle CS and COX activities relative to protein were greater in PNS than non-PNS calves (P ≤ 0.03), but COX relative to CS activity was similar between groups. Activity of COX was greater in bulls than heifers (P = 0.03), but no other measure was affected by sex. All other measures were unaffected by PNS. Prenatal stress did not affect markers of muscle damage and oxidative stress in yearling Brahman calves at rest but mitochondrial volume density and function were greater in PNS calves. Acute stressors induce oxidative stress, so implications of differences in mitochondria in PNS calves following a stressor should be investigated.

PSV-12 Evaluator effect on the prediction of genetic merit using subjective measures of beef cattle temperament

The objective of this study was to determine evaluator effect on genetic predictions of docility score (DS), temperament score (TS), and qualitative behavior attributes (QBA; 6 positive and 6 negative QBA). Weaning age calves (n = 1,542) were evaluated using DS, TS, and QBA scores over 4-year period. Fixed effects included evaluator (n = 11 total), primary breed (n = 2), sex (n = 2), and random effect of calf using repeated measure design (across-evaluator model only). Variance components, heritability, and breeding values (EBV) were estimated using pedigree in ASReml 4.2. Evaluator effect on EBV was based on 1) Spearman rank correlation coefficients (rs) and 2) 3-quartiles change in rankings of calves among evaluators per trait. Results showed evaluator scoring was different in DS, TS, and 12 QBA (P &lt; 0.002). Correlations (rs) on EBV across and within evaluators ranged from 0.52 to 0.84 (DS) and 0.53 to 0.87 (TS) (P &lt; 0.0001). For QBA, all rs on EBV were significant (P &lt; 0.05) and ranged from -0.26 to 0.90 (positive) and 0.11 to 0.82 (negative), except on positively occupied and distressed QBAs (P = 0.46 to 0.65). When comparing evaluator specific EBV to across evaluator model, DS and TS had 1.10% and 1.56% 3-quartiles change, whereas comparison of evaluator to evaluator EBV showed 4.54% and 3.76%. For positive and negative QBA attributes, 3-quartiles change ranged from 0% to 35.13% and 0.52% to 3.76% (evaluator to across-evaluator model); and 1.30% to 21.09% and 3.31% to 13.81% (within evaluators). This illustrates that DS and TS had less variability in EBV prediction compared to QBA. Within QBA, negative QBA had less variability compared to positive QBA. In conclusion, evaluators scored differently using these methods. However, in predicting genetic merit, evaluator has negligible effect for scoring systems already implemented by breed associations.

84 Evaluation of hair coat shedding ability as an adaptive trait in Angus cattle in the southern U.S

The objective of this study was to evaluate winter hair coat shedding ability and the association to cow production performance. This work was a collaborative effort from the S-1064 multi-state research project. Data were collected on spring and fall-calving purebred and commercial Angus females with calves (n = 1192) from March until July from 2008 to 2018 in Mississippi, Arkansas and Texas. Dams were observed once monthly by two trained technicians for winter hair coat shedding and give a visual hair shedding score of 1 to 5 with 1 indicating 100% shed, 2 = 75% shed, 3 = 50% shed, 4 = 25% shed, and 5 indicating 0% shedding of the winter hair coat. Month of first shedding (MFS) was determined once a female reached an average hair shedding score of ≤ 3.5 for any given month. Calf birth and weaning weights were considered as performance of the dam. For cow performance, the response variables included AdjBW and d205wt of the calf. Data were analyzed using the MIXED procedure of SAS with significance declared at P ≤ 0.05. The model included MFS as a fixed effect with calf sex as a covariate and calf sire as a random effect. Females that obtained a shedding score of ≤ 3.5 by May had heavier calves at birth (P ≤ 0.01) when compared with those that shed in June and July, with July weights being the lightest. Females that scored ≤ 3.5 by May has significantly heavier d205wt (P ≤ 0.01) when compared to those that shed by June and July. When comparing locations, females in Arkansas shed later (P ≤ 0.01) when compared with both Mississippi and Texas, which were similar. Shedding ability could be an important economic trait to be considered in cow-calf operations.

86 Evaluation of the effect of prenatal transportation stress on endocrine and immune tissues of neonatal Brahman calves

The objective of this experiment was to evaluate whether prenatal transportation stress (PNS) affects the weight of endocrine and immune tissues of calves. Mature Brahman cows inseminated to a single Brahman sire in 2018 were assigned to be either Control (n = 35; not transported) or PNS (n = 37; 2 h of transportation at 60, 80, 100, 120, and 140 ± 5 d of gestation). Of the calves born in 2019, 16 Control (8 bulls and 8 heifers) and 16 PNS (8 bulls and 8 heifers) calves were studied. Pen score, body weight, and blood samples were obtained from calves at 25 ± 2 d of age. At that time, calves were euthanized by barbiturate overdose in order to collect tissues (brain, pituitary and adrenal glands, spleen and thymus). Tissues were trimmed and weighed. Serum cortisol was determined by RIA. Data were analyzed using ANOVA, GLM, and CORR procedures of SAS with body weight at sample collection as a covariate. Whole pituitary weight tended (P = 0.08) to be greater in females compared to males, whereas anterior pituitary weight tended (P = 0.06) to be greater in PNS. Whole brain weight was greater (P &lt; 0.01) in males. The interaction of treatment and sex on whole brain weight of PNS tended (P = 0.09) to fall between that of the Control females and Control males. Whole pituitary weight was positively correlated with total adrenal weight (r = 0.32; P = 0.08). Total thymus weight positively correlated with whole brain and total adrenal weight (r = 0.43 and 0.41, respectively; P &lt; 0.05). Serum cortisol was negatively correlated with spleen weight and total immune tissue weight (r = -0.37 and -0.38, respectively; P &lt; 0.05). Results suggest next steps should include a closer look at function of immune tissues by studying PNS effect on thymus response to immunization.

Novel genes involved in the genetic architecture of temperament in Brahman cattle

Cattle temperament is a complex and economically relevant trait. The objective of this study was to identify genomic regions and genes associated with cattle temperament. From a Brahman cattle population of 1,370 animals evaluated for temperament traits (Exit velocity-EV, Pen Score-PS, Temperament Score-TS), two groups of temperament-contrasting animals were identified based on their EV-average values ±1/2 standard deviation (SD). To be considered in the calm group, the EV of females ranged between 0.16–1.82 m/s (n = 50) and the EV of males ranged between 0.4–1.56 m/s (n = 48). Females were classified as temperamental if their EV ranged between 3.13–7.66 m/s (n = 46) and males were classified as temperamental if their EV ranged between 3.05–10.83 m/s (n = 45). Selected animals were genotyped using a total of 139,376 SNPs (GGP-HD-150K), evaluated for their association with EV. The Genome-Wide Association analysis (GWAS) identified fourteen SNPs: rs135340276, rs134895560, rs110190635, rs42949831, rs135982573, rs109393235, rs109531929, rs135087545, rs41839733, rs42486577, rs136661522, rs110882543, rs110864071, rs109722627, (P<8.1E-05), nine of them were located on intergenic regions, harboring seventeen genes, of which only ACER3, VRK2, FANCL and SLCO3A1 were considered candidate associated with bovine temperament due to their reported biological functions. Five SNPs were located at introns of the NRXN3, EXOC4, CACNG4 and SLC9A4 genes. The indicated candidate genes are implicated in a wide range of behavioural phenotypes and complex cognitive functions. The association of the fourteen SNPs on bovine temperament traits (EV, PS and TS) was evaluated; all these SNPs were significant for EV; only some were associated with PS and TS. Fourteen SNPs were associated with EV which allowed the identification of twenty-one candidate genes for Brahman temperament. From a functional point of view, the five intronic SNPs identified in this study, are candidates to address control of bovine temperament, further investigation will probe their role in expression of this trait.

Genetic characteristics of Korean Jeju Black cattle with high density single nucleotide polymorphisms

Objective

Conservation and genetic improvement of cattle breeds require information about genetic diversity and population structure of the cattle. In this study, we investigated the genetic diversity and population structure of the three cattle breeds in the Korean peninsula.

Methods

Jeju Black, Hanwoo, Holstein cattle in Korea, together with six foreign breeds were examined. Genetic diversity within the cattle breeds was analyzed with minor allele frequency (MAF), observed and expected heterozygosity (H_O and H_E), inbreeding coefficient (F_IS) and past effective population size. Molecular variance and population structure between the nine breeds were analyzed using a model-based clustering method. Genetic distances between breeds were evaluated with Nei’s genetic distance and Weir and Cockerham’s F_ST.

Results

Our results revealed that Jeju Black cattle had lowest level of heterozygosity (H_E = 0.21) among the studied taurine breeds, and an average MAF of 0.16. The level of inbreeding was −0.076 for Jeju Black, while −0.018 to −0.118 for the other breeds. Principle component analysis and neighbor-joining tree showed a clear separation of Jeju Black cattle from other local (Hanwoo and Japanese cattle) and taurine/indicine cattle breeds in evolutionary process, and a distinct pattern of admixture of Jeju Black cattle having no clustering with other studied populations. The F_ST value between Jeju Black cattle and Hanwoo was 0.106, which was lowest across the pair of breeds ranging from 0.161 to 0.274, indicating some degree of genetic closeness of Jeju Black cattle with Hanwoo. The past effective population size of Jeju Black cattle was very small, i.e. 38 in 13 generation ago, whereas 209 for Hanwoo.

Conclusion

This study indicates genetic uniqueness of Jeju Black cattle. However, a small effective population size of Jeju Black cattle indicates the requirement for an implementation of a sustainable breeding policy to increase the population for genetic improvement and future conservation.

Deciphering Cattle Temperament Measures Derived From a Four-Platform Standing Scale Using Genetic Factor Analytic Modeling

The animal's reaction to human handling (i.e., temperament) is critical for work safety, productivity, and welfare. Subjective phenotyping methods have been traditionally used in beef cattle production. Even so, subjective scales rely on the evaluator's knowledge and interpretation of temperament, which may require substantial experience. Selection based on such subjective scores may not precisely change temperament preferences in cattle. The objectives of this study were to investigate the underlying genetic interrelationships among temperament measurements using genetic factor analytic modeling and validate a movement-based objective method (four-platform standing scale, FPSS) as a measure of temperament. Relationships among subjective methods of docility score (DS), temperament score (TS), 12 qualitative behavior assessment (QBA) attributes and objective FPSS including the standard deviation of total weight on FPSS over time (SSD) and coefficient of variation of SSD (CVSSD) were investigated using 1,528 calves at weaning age. An exploratory factor analysis (EFA) identified two latent variables account for TS and 12 QBA attributes, termed difficult and easy from their characteristics. Inclusion of DS in EFA was not a good fit because it was evaluated under restraint and other measures were not. A Bayesian confirmatory factor analysis inferred the difficult and easy scores discovered in EFA. This was followed by fitting a pedigree-based Bayesian multi-trait model to characterize the genetic interrelationships among difficult, easy, DS, SSD, and CVSSD. Estimates of heritability ranged from 0.18 to 0.4 with the posterior standard deviation averaging 0.06. The factors of difficult and easy exhibited a large negative genetic correlation of -0.92. Moderate genetic correlation was found between DS and difficult (0.36), easy (-0.31), SSD (0.42), and CVSSD (0.34) as well as FPSS with difficult (CVSSD: 0.35; SSD: 0.42) and easy (CVSSD: -0.35; SSD: -0.4). Correlation coefficients indicate selection could be performed with either and have similar outcomes. We contend that genetic factor analytic modeling provided a new approach to unravel the complexity of animal behaviors and FPSS-like measures could increase the efficiency of genetic selection by providing automatic, objective, and consistent phenotyping measures that could be an alternative of DS, which has been widely used in beef production.

Cattle adapted to tropical and subtropical environments: genetic and reproductive considerations

Efforts to understand biological functions and develop management schemes specific to Bos indicus-influenced cattle raised in tropical and subtropical environments are critical to meet the increasing global demand for protein. In the United States, B. indicus breeds are mostly used to generate B. indicus × B. taurus crosses with increased thermal and parasite tolerance, while retaining some productive characteristics of B. taurus cattle. Although crossbreeding represents a proven strategy to improve cattle adaptation almost immediately, research has also attempted to identify B. taurus genetics that can withstand subtropical and tropical climates. Reduced milk production and delayed reproductive maturation appear to be related with tropical adaptation of B. taurus breeds, as a means to conserve energy under stressful conditions and limited nutrition. Moreover, longevity may be the ultimate adaptation response to unfavorable environments, and retention of bulls and heifers from proven cows is the recommended strategy to improve longevity in B. indicus-influenced herds. Besides selection for longevity, other aspects should be considered when planning reproductive management in tropical and subtropical regions. Bos indicus and B. taurus breeds have multiple differences pertaining to reproductive function, including age at puberty, ovarian dynamics, and pregnancy development. Nutritional strategies such as the stair-step regimen, and use of exogenous progesterone (P4) inserts are options to hasten puberty attainment of late-maturing B. indicus-influenced heifers. Yet, limited pharmacological alternatives are available for reproductive management of B. indicus-influenced females in the United States, which rely on GnRH-based protocols not specifically designed to the reproductive function of B. indicus breeds. In contrast, hormonal protocols based on exogenous P4, estradiol esters, and equine chorionic gonadotropin are available for use in B. indicus females in South America. These include protocols tailored to prepubertal heifers, anestrous cows, and cycling nulliparous or parous females, which often yield pregnancy rates of 50% to fixed-time artificial insemination. The global dairy industry also faces similar challenges in increasing demand and production as the beef industry. Selection of cows capable of sustaining optimal milk yield, reproductive success, and health status in hot and humid conditions is essential for optimal dairy production in subtropical and tropical regions.

Reproductive Biology and Evidence of Diapause in the Cowpea Curculio (Coleoptera: Curculionidae)

Chalcodermus aeneus Boheman (Coleoptera: Curculionidae) has been the most destructive insect pest of black-eyed peas or cowpeas, Vigna unguiculata L. (Fabales: Fabaceae), over the last century in the southeastern United States. The historical distribution of this semitropical pest suggests the likelihood that diapause plays a key role in the overwintering success in parts of the United States. However, this report is the first to document biological evidence for diapause in C. aeneus. Our study assessed larval emergence from cowpea pods in the summer to fall growing seasons, egg development in female adults over the first (summer) and second (fall) generations, and adult emergence from infested soil after the first and second generations. There was a clear reduction in larval emergence from summer to fall. Egg and follicle development in female C. aeneus dropped off dramatically by September of each year. There was an extended emergence pattern of weevil adults from the soil in the fall as compared to the summer generation. Any future regional management of cowpea curculio will have to take into account the ability of this insect to diapause, thereby increasing its capacity to overwinter in regions where the cowpea crop, a warm-season, semitropical plant, is terminated with winter freezing temperatures.

Shedding and characterization of gastrointestinal nematodes of growing beef heifers in Central Texas

•

Bos indicus crossbred females had higher fecal eggs as calves than as yearlings.

•

More Bos indicus background had higher eggs per gram of feces through 1 year of age.

•

Cooperia and Haemonchus species predominated through 1 year of heifer age.

•

Ostertagia and Trichostrongylus species were detected mostly after 1 year of age.

•

Eggs per gram of feces correlated with few production traits of heifers.

Gastrointestinal nematodes (GIN) do not detrimentally affect cattle to the extent of small ruminants. However, they are developing resistance to drugs used to treat them. Genetic strategies to control the nematodes and/or their detrimental effects could be a sustainable alternative to treatment with drugs. An essential first step in development of such a strategy is characterization of nematode populations in commonly used breed types of cattle in local conditions. Fecal egg counts (FEC) were obtained every two months on a cohort of 53 crossbred Nellore-Angus heifers grazing Central Texas pastures from an average heifer age of 3 months to approximately 2 years of age. For 10 of those 12 sets of samples, coprocultures were set up to characterize gastrointestinal nematode species present. Heifers were ½ Nellore ½ Angus (n = 18) or ¾ Angus ¼ Nellore (n = 35). They were born in the spring of 2014 to cows that were from 3–5 years old. They were maintained as a group throughout weaning, postweaning, exposure to bulls as yearlings, and as pregnant cows through the birth of their first calves. An interaction of breed group with sampling time (P < 0.0001) highlighted favorable FEC of ¼ Nellore heifers as compared to ½ Nellore in all but two sampling times. Fecal egg count means were, in general, higher for heifers of both groups in sampling times up to one year of age. Season effects on FEC may be important, but the effect of age may have obscured their detection. There were few significant correlation coefficients for FEC traits with a variety of production traits of these females. Average FEC residuals were positively correlated (r = 0.28 and 0.41; P < 0.05) with winter coat shedding score evaluated at approximately 17 and 24 months of age. Residual correlations of average FEC with calf weaning weight and incidence of shedding with calf age at weaning (r> 0.3) may be indicative of the increased susceptibility of females that lactate heavily or longer to internal parasite infection. Proportions of GIN genera by sampling day differed from χ2 expectation (P < 0.0001). Cooperia and Haemonchus species were detected in large proportions in sampling dates that corresponded to heifers less than one year of age. Ostertagia and Trichostrongylus species predominated in sampling dates after heifers reached one year of age.

Random regression of Hereford percentage intramuscular fat on geographical coordinates

Accounting for genotype-environment interactions may improve genetic prediction and parameter estimation. The objective was to use random regression analyses to estimate variances and thereby heritability for intramuscular fat (IMF) across longitude and latitude coordinates within the continental United States. Records from the American Hereford Association (n = 169,440) were used. Analyses were first conducted using the continental United States in its entirety, and then as subdivided into two or four regions. Data were analyzed with an animal model, and linear and quadratic random regressions of additive genetic merit on longitude or latitude as covariate (separately). Subdivided data were analyzed with linear random regressions unique to regions. Regions were North and South separated at 40°N latitude, or West and East separated at 99°W longitude using longitude or latitude as covariate, respectively. Further subdivision to four regions included additional boundaries of 44.46° and 36.46°N latitude and 104.55° and 92.22°W longitude. The estimated heritability of IMF from the traditional model was 0.19 ± 0.004. Without regional subdivision of data, quadratic random regression had the best fit for the data based on likelihood ratio tests using longitude or latitude as covariate (P < 0.01). Estimates of heritability from quadratic random regression on latitude ranged from 0.12 in the South to a high of 0.27 at the extreme Northern latitude. Estimates of heritability from quadratic random regression on longitude ranged from 0.17 in the middle of the parameter space (corresponding to the central United States) to 0.37; higher estimates were noted at the extremes, that is, the far West and East longitudes. Random regression analyses of data divided into regions were conducted with a linear coefficient, as increasing to a quadratic polynomial was never accomplished. Results from random regression on latitude in the East region were similar to results from analyses without regions (h2 ranged from 0.09 to 0.32); however, estimates of heritability in the West region had a lower range from South to North (0.14 to 0.27). Estimates of heritability from random regression on longitude with data divided into two regions were similar to those from analyses that did not include region. Estimates in the South region were somewhat lower and had a lower range (0.15 to 0.31) than those from the North region (0.19 to 0.47). When data were further subdivided, estimation of only a subset of covariances among random regression coefficients was possible, that is, within-region covariances of intercept and linear terms (latitude); those and covariances between all linear random regression coefficients were estimated when longitude was the covariate. Results from random regression analyses of data with four regions modeled produced very high estimates of heritability in low latitudes in the furthest West and high latitudes in the furthest East region, with approximate difference of 0.3 and 0.2 between estimates in the two West regions and the two East regions, respectively. Results from random regression on longitude indicated higher estimates of heritability in North region, especially at the furthest East longitudes of the most Northern region. There appeared to be substantial additive genetic variance differences, as well as estimates of heritability, that correspond to different geographical environments as modeled by random regressions on within-region latitude or longitude coordinates.

Influence of prenatal transportation stress-induced differential DNA methylation on the physiological control of behavior and stress response in suckling Brahman bull calves

The objective of this experiment was to examine potential differential methylation of DNA as a mechanism for altered behavioral and stress responses in prenatally stressed (PNS) compared with nonprenatally stressed (Control) young bull calves. Mature Brahman cows (n = 48) were transported for 2-h periods at 60 ± 5, 80 ± 5, 100 ± 5, 120 ± 5, and 140 ± 5 d of gestation (Transported group) or maintained as nontransported Controls (n = 48). From the offspring born to Transported and Control cows, a subset of 28-d-old intact bulls (n = 7 PNS; n = 7 Control) were evaluated for methylation of DNA of behavior and stress response-associated genes. Methylation of DNA from white blood cells was assessed via reduced representation bisulfite sequencing methods. Because increased methylation of DNA within gene promoter regions has been associated with decreased transcriptional activity of the corresponding gene, differentially methylated (P ≤ 0.05) CG sites (cytosine followed by a guanine nucleotide) located within promoter regions (n = 1,205) were used to predict (using Ingenuity Pathway Analysis software) alterations to canonical pathways in PNS compared with Control bull calves. Among differentially methylated genes (P ≤ 0.05) related to behavior and the stress response were OPRK1, OPRM1, PENK, POMC, NR3C2, TH, DRD1, DRD5, COMT, HTR6, HTR5A, GABRA4, GABRQ, and GAD2. Among altered (P < 0.05) signaling pathways related to behavior and the stress response were Opioid Signaling, Corticotropin-Releasing Hormone Signaling, Dopamine Receptor Signaling, Dopamine-DARPP32 Feedback in cAMP Signaling, Serotonin Receptor Signaling, and GABA Receptor Signaling. Alterations to behavior and stress response-related genes and canonical pathways supported previously observed elevations in temperament score and serum cortisol through weaning in the larger population of PNS calves from which bulls in this study were derived. Differential methylation of DNA and predicted alterations to behavior and stress response-related pathways in PNS compared with Control bull calves suggest epigenetic programming of behavior and the stress response in utero.

303 Identification of candidate genes related to temperament in Brahman cattle

The objective of this study was to identify genomic regions and genes associated with beef cattle temperament. Temperament, measured as exit velocity (EV; m/s), was recorded in 1,370 Brahman cattle from Texas A&M AgriLife Research at Overton, TX. We identified two groups of temperament-contrasting animals. Cows were calm if their EV of 0.16–3.41 m/s and bulls if their EV was 0.4–3.12 m/s (n-119). Cows were temperamental if their EV was 3.55–7.66 m/s and bulls if their EV was 3.13–10.83 m/s (n = 79). The 198 animals were genotyped using the GGP-HD-150K chip. 139,376 SNPs were evaluated for association with temperament. 13 SNP′s were associated with EV (P &lt; 4.0E-05). The SNPs GABRG2-26484, NRXN3-26436 and TBX20-191081 are located in introns of the GABRG2, NRXN3 and TBX20 genes, respectively. The GABRG2 gene encodes a GABA receptor, the major inhibitory neurotransmitter in the mammalian brain. The NRXN3 gene encodes receptor proteins related to chemical transmission at synapses. TBX20 is a member of the T-box transcription factor family expressed in the developing stages of heart, limbs, eye and ventral neural tube. To test the effect of these 3 SNP′s on EV, Pen-Score and Temperament-Score, a general linear model was fitted including the fixed effects of sex of calf and year of birth, and the individual effect of the 3 SNPs. The marker TBX20-191081 was associated with the three traits evaluated (P &lt; 0.01), where the GG genotype was associated with the calmest temperament. The GG genotype had a significant effect on EV (P &lt; 0.0001) that was 1.35 and 1.95 m/s slower than AG and AA, respectively. For TS, the GG genotype had a TS that was 1.41 and 1.24 DS less than those of the AA and GA genotypes. Our study indicates that genetic control of cattle temperament has a wide network of genes with divergent functions and genetic background specificity.

304 Transgenerational influence of prenatal stress on temperament, birth weight, and weaning weight of Brahman calves

The objective was to evaluate the transgenerational influence of prenatal transportation stress (PNS) on temperament, birth weight, and weaning weight in Brahman calves. Dams were derived from multiparous cows that were transported for 2 h at 60, 80, 100, 120, and 140 d of gestation (n = 48) or maintained as Controls (n = 48). Heifers born to Transported or Control cows were followed through maturity and their calves (n = 79; 35 calves born to Control dams and 44 to PNS dams) studied across multiple years. Calf birth and weaning weights were recorded. Calf temperament at weaning was evaluated utilizing pen score (PS; 1=calm and 5=excitable), exit velocity (EV; m/sec), and temperament score (TS; (PS+EV)/2). Data were analyzed using mixed linear models. Fixed effects included dam treatment, calf sex, and birth season nested within birth year (season nested within year was identified as an efficient parameterization). Calf sire was a random effect. Birth weight was not affected by dam treatment, calf sex, or season nested within year (P &gt; 0.1). Weaning weight was greater in males than females (P &lt; 0.01). Weaning weight was affected by season nested within year (P &lt; 0.01). DOA at weaning was a covariate for weaning weight (0.74±0.09 kg; P &lt; 0.01), but not PS, EV, or TS (P &lt; 0.1). Weaning weight was not affected by dam treatment (P &gt; 0.1). TS was affected by season nested within year (P = 0.04) but not dam treatment or calf sex (P &gt; 0.1). EV was not affected by dam treatment, calf sex, or season nested within year (P &gt; 0.1). PS was affected by season nested within year (P &lt; 0.01) but not calf sex (P &gt; 0.1). Calves born to PNS dams had greater (P = 0.02) PS than calves born to Controls (3.17±0.20 and 2.58±0.22, respectively). There was a transgenerational influence of PNS on temperament but not birth or weaning weight.

37 Genetic aspects of livestock adaptation

Livestock adaptation to less than favorable ambient has a genetic basis. Estimates of additive genetic variance and narrow sense heritability for % intramuscular fat in Hereford varied across United States geography quantifications evaluated in random regression analyses. Shedding and regrowth of winter coats of Angus cows in subtropical areas may indicate differential adaptation in cattle not adapted to the subtropics. Acquired adaptation to local conditions (multiple generations across more than 50 years) may was evidenced by performance of Florida Angus relative to prominent U.S. Angus in subtropical Florida; later sexual maturation may be a prominent component of such adaptation, and increasing milk production may be antagonistic to adaptation in natural conditions. Cattle temperament may be indicative of adaptation and is highly heritable; however, results from random regression analyses suggest that the additive genetic component appears to decrease in importance and the permanent environmental component of phenotype appears to become more important as calves age. Crossbreeding represents a proven strategy to improve adaptation almost immediately. Heterosis influences cattle body temperature maintenance, reproduction, survival, and, to a lesser extent, temperament in subtropical or other stressful environmental conditions (for example, in toxic fescue). Prenatal stress alters patterns of methylation (and likely other epigenetic mechanisms) and thereby encourages or inhibits gene expression to promote postnatal fitness. Brahman exposed to prenatal stress exhibited substantially different patterns of methylation across the genome in lymphocytes in both male and female calves; those patterns differed by sex. Female longevity may be the ultimate adaptation trait, as annual compliance to reproductive standards may be an appropriate assessment of a combination of attributes that represent adaptation. Longevity has documented heterotic influence; the additive genetic component is less well characterized but real. A simple, effective way to improve longevity may be to select bulls from aged, proven cows.

PSVIII-40 Effect of prenatal transportation stress on DNA methylation in Brahman heifers

Prenatal stress in beef cattle has the potential to result in unexpected, long term phenotypic changes to an animal’s performance. One possible explanation for this phenomenon is stress induced epigenetic changes to the genome (e.g., methylation or histone modifications). The objective of this experiment was to identify differential DNA methylation patterns in prenatally stressed (PNS) Brahman females. Multiparous Brahman cows (n = 48) were transported for a 2hr periods on 60 ± 5, 80 ± 5, 100 ± 5, 120 ± 5, and 140 ± 5 d of gestation. A non-transported group (n = 48) was maintained as a control. DNA was extracted from the white blood cells of their 28-d old heifer calves (Control n = 8, PNS n = 6) and methylation was assessed using reduced representation sodium bisulfite sequencing. Mean methylation ratios at locations across the genome in PNS and Control females were compared using a student’s t-test. Analysis revealed 16,377 CpG sites, 311 CHG, and 612 CHH (C= cytosine; G= guanine; H = adenine, cytosine or thymine) sites that were differentially methylated in PNS females compared to the Control females (P ≤ 0.05). The top 2,000 differentially methylated regions (DMRs) were classified by functional location in introns, exons, CpG islands, promoters, and 5’UTR regions. Methylation differences in these regions are known to cause gene expression differences. The identified regions were analyzed using Ingenuity Pathway Analysis (IPA) in order to identify involvement in canonical pathways. Certain regions identified from the hypermethylated DMRs showed the strongest overlap (P = 1.93E-6) with the complex pathway of axonal guidance signaling (important for neural development). These results suggest that prenatal transportation stress does alter the epigenome of Brahman heifers and can influence biological pathways, potentially impacting animal performance.

PSVIII-39 Genome-wide DNA methylation alteration in prenatally stressed Brahman heifer calves with the advancement of age

This study investigated whether DNA methylation patterns changed over the first five yr of life within prenatally stressed (PNS) heifer calves compared to change within a Control group. Prenatal stress was induced by the transportation of pregnant Brahman cows for 2-hr periods at 60±5, 80± 5, 100±5, 120±5, and140±5d of gestation. White blood cells were sampled from the same 6 PNS heifer calves and 8 Control heifer calves at 28 d and 5 yr of age. The DNA methylation data were generated through Reduced Representation Bisulfite Sequencing. Based on results of mapping and bioinformatics analyses, 73,758 hypermethylated and 73,367 hypomethylated CpG sites, 375 hypermethylated and 377 hypomethylated CHG sites, 735 hypermethylated and 842 hypomethylated CHH (C = cytosine; G = guanine; H = either adenine, thymine, or cytosine) sites were obtained from 28-d-old PNS calves compared to when they had matured into 5-yr-old PNS cows (P ≤ 0.05). The 28-d-old Control heifer calves contained 53,005 hypermethylated and 57,103 hypomethylated CpG sites, 200 hypermethylated and 202 hypomethylated CHG sites, 439 hypermethylated and 535 hypomethylated CHH sites compared to when they matured into 5-yr-old Control cows (P ≤ 0.05). As DNA methylation of gene promoter regions is associated with reduced transcription activity, strongly hypermethylated and hypomethylated CpG sites located in promoter regions underwent Ingenuity Pathway Analysis. The top canonical pathways altered by strongly hypermethylated and hypomethylated CpG sites between 28-d-old and 5-yr-old PNS cows were 4-1BB Signaling in T Lymphocytes (P = 0.00169) and Transcriptional Regulatory Network in Embryonic Stem Cells (P = 0.000744). Mineralocorticoid Biosynthesis (P = 0.00901) and Transcriptional Regulatory Network in Embryonic Stem Cells (P = 0.000804) were the other top canonical pathways altered between 28-d-old and 5-yr-old Control cows. PNS calves appeared to develop an altered epigenome compared to Control group calves during the first five yr from birth.

376 Awardee Talk - Prenatal transportation stress alters physiology of suckling Brahman calves as mediated by changes in DNA methylation

Mechanisms by which prenatal programming alters the developmental trajectory of biological systems in utero is of relevance to animal agriculture. We hypothesized that prenatally stressed (PNS) calves would exhibit altered postnatal phenotype and differential methylation of DNA. Mature Brahman cows (n = 48) were transported for 2-h periods at 60, 80, 100, 120, and 140 (± 5) d of gestation or maintained as non-transported Controls (n = 48). PNS calves exhibited more excitable temperaments and increased circulating concentrations of cortisol. In response to a glucose tolerance test administered to a subset of PNS and Control bulls, PNS calves returned to basal concentrations of insulin more rapidly and exhibited a lower total area under the insulin response curve. In response to a lipopolysaccharide challenge administered to a subset of PNS and Control bulls, PNS calves exhibited greater circulating concentrations of cortisol and a greater change from baseline for IFN-γ and monocytes. Reduced representation bisulfite sequencing was used to assess methylation of DNA from white blood cells in a subset of 28-d-old bull calves. Because increased methylation of DNA in a gene promoter region typically results in decreased transcriptional activity, differentially methylated CG sites located within promoter regions (n = 1,205) were used to predict alterations to canonical pathways using Ingenuity Pathway Analysis software. There were 113 pathways altered (P ≤ 0.05) in PNS calves. Among these were pathways related to behavior, stress response, metabolism, and immune function. In support of our hypotheses, differential methylation in gene regions related to temperament, stress response, metabolism, and immune response corresponded to phenotypic differences observed in PNS calves, suggesting DNA methylation as a mechanistic basis for prenatal programming. The role of the prenatal epigenome in shaping postnatal outcomes provides novel opportunities to improve animal health and production traits.

PSXIII-13 Comparison of telomere length in age-matched primiparous and multiparous Brahman cows

Telomeres are comprised of G-rich nucleotide sequences (5’-TTAGGG-3’) at the chromosome termini that are responsible for protecting chromosomes; however, attrition of these sequences has been observed in conditions of physiological and psychological stress. The purpose of this study was to compare telomere length (TL) in 4-yr-old Brahman cows grouped by first parity (n = 8) and second parity (n = 11). Cows were bled by jugular venipuncture, weighed, and had BCS recorded d+28 prior to calving and d-7 and d-28 post calving. Cows were observed for duration of labor (Tlabor) and calving ease (CE) at the time of parturition. Calf birth weight (CBW) and gender (CG) were recorded. Peripheral leukocytes were isolated, complete blood counts (CBC) were recorded, and genomic DNA was extracted utilizing silicone membrane spin column kits. The relative quantity of telomere products, which is proportional to the average TL, was determined by multiplex quantitative PCR analysis using the ratio of bovine telomere and β-globulin DNA. An absolute standard of bovine telomere (1012–107 dilution series) and β-globulin (109-104 dilution series) genes was utilized to produce relative copy number. All samples were run in triplicate and samples were included if triplicate Cq difference was less than 0.25 cycles. Parity was the fixed effect of interest and random effects included sire and day repeated with cow as the subject. No differences in CBC were seen. Tlabor, CE, CG, and CBW did not impact TL (P &gt; 0.1). A trend was observed for day-parity interaction (P = 0.0918). TL between parity differed most on d-28 (P = 0.1046; parity one 127292 ± 6483; parity two 111045 ± 5376). The stress of parturition and raising the first calf of a cow’s life may be responsible for slight attenuation in TL.

Integrated Management of Tomato Yellow Leaf Curl Virus and its Whitefly Vector in Tomato

Whitefly-transmitted Tomato yellow leaf curl virus (Family Geminiviridae; Genus Begomovirus) severely restricts tomato production in the Southeastern United States. Whitefly and tomato yellow leaf curl virus management studies typically investigate control tactics individually, but successful management of this pest complex more often relies on a combination of tactics. This study examined the individual and combined effects of tomato yellow leaf curl virus-resistant cultivars, insecticides, and metallic reflective mulch on whiteflies, tomato yellow leaf curl virus disease incidence, and marketable tomato yields using split-split plot trials over 3 yr. Reflective mulch significantly reduced whitefly adults and nymphs and tomato yellow leaf curl virus symptom severity in all 3 yr of the study. Reflective mulch treatments also provided greater marketable tomato yield in 2 out of 3 yr. Imidacloprid and cyantraniliprole treatments reduced whitefly adults and nymphs' establishment and marginally increased yields, but there was no significant insecticide effect on tomato yellow leaf curl virus incidence/symptom severity compared with the non-treated check. Virus-resistant tomato cultivars did not influence whitefly populations, but provided consistent reduction in virus disease incidence. Interactions between host plant resistance and insecticide treatments ranged from strongly additive in the standard white plastic mulch treatment to only marginally additive in the reflective mulch treatments in terms of enhancing tomato yields. tomato yellow leaf curl virus-resistant tomato cultivars and reflective mulch provided the bulk of the protection against tomato yellow leaf curl virus disease incidence. However, it was the combination of all the best tactics (reflective mulch, cyantraniliprole, 'Security' hyb. tomato yellow leaf curl virus-resistant) that provided the maximum increase in marketable tomato yield (2.8-fold) over the least effective combination (white mulch, no whitefly insecticide, 'FL47' hyb. tomato yellow leaf curl virus-susceptible).

114 The relationship of cell mediated immune response of Brahman cows on calving interval

The relationships between cell-mediated immune response (CMIR) of cows and calving interval and calf growth were evaluated. Calving interval was evaluated for each of the initial 7 life-time opportunities that females had to produce a calf. The number of days from when the cow had a calf until she calved again determined the cow’s calving interval for each year. The Spring breeding seasons were divided into 45-d artificial insemination and 45-d natural service periods. The Fall breeding season was natural service each year. Multiparous Brahman cows (n = 435) had records for breeding and calving dates and calf weaning weights. The CMIR (hypersensitivity response to Candida albicans measured using tail-fold thickness) was determined once for each female in the Fall and early Winter 2015. Cows were classified into response groups based on mean and SD of CMIR. High responders were those with CMIR ½ SD ≥ the mean (≥ 2.8 mm), Intermediate cows were within ½ SD of the mean (2.7-1.8 mm), and cows ½ SD ≤ the mean were Low (≤ 1.7 mm). Data were analyzed using repeated measures, mixed, and Chi Square procedures (SAS 9.4). Variables included calving interval, cow age, CMIR response group, pregnancy status, calf count, calf sex, cow sire, calf sire, and calf adjusted 180-d weaning weight. A greater proportion of High CMIR cows made it to their 6th calving interval. Chi Square analysis indicated that High CMIR cows were more likely to remain in the herd longer and produce a calf each year. Cows with High CMIR had a greater stayability (P = 0.1). Selection for High CMIR can result in selection of cows which are more likely to remain in the herd long enough to be profitable as they have a greater stayability and shorter calving intervals (about 25 days) than their Low CMIR herdmates.

154 Evaluation of repeated exposure to lipopolysaccharide on vaginal temperature and complete blood count in two-year-old Brahman heifers

The effects of escalating lipopolysaccharide (LPS) doses on vaginal temperature (VT) and complete blood count (CBC) were evaluated to develop a model of low-level inflammatory response in cattle. Brahman heifers (2-yr-old; 326 kg BW) paired by birthdate, sire, and weaning temperament score were fitted with self-contained, indwelling vaginal temperature probes 7 d before initiation of treatment (d 0) and randomly assigned to Control (n = 6) or LPS (n = 6) treatment. Heifers were maintained as a group with free choice access to Coastal bermudagrass hay and fed 3.6 kg of a 3:1 corn:corn gluten grain mix per head per day. Heifers were weighed, BCS recorded, and blood samples were collected weekly by jugular venapuncture before delivery of treatments. Treated heifers received LPS (from Escherichia coli) via subcutaneous injection in the neck (d 0: 0.25 μg / kg BW; d7: 0.50 μg /kg BW; d 14: 0.75 μg /kg BW) and Control heifers received saline. Sickness behavior score (SBS) was monitored at 30-min intervals for 12 hr after LPS. Treatment was the fixed effect of interest. Random effects included sire and day repeated with heifer as the subject. No dose-treatment interactions were detected. The SBS did not change after subcutaneous LPS or saline (P > 0.1). Neutrophil and lymphocyte numbers 7 d after each injection were not affected by LPS (P > 0.1). A febrile response began within 1 hr and persisted up to 14 hr after LPS injections. The febrile effect of the 0.25 μg LPS dose was less pronounced (P < 0.05) than that of the 0.50 μg and 0.75 μg LPS doses which were similar to one another. Lipopolysaccharide treatment increased maximum VT (P < 0.0001), hastened time to attain maximum VT (P = 0.0429), and increased the change in VT (P < 0.0001). Weekly subcutaneous injection of LPS facilitates the study of the effect of sub-clinical illness in beef cattle.