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

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

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

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

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

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

Growth characteristics, reproductive performance, and evaluation of their associative relationships in Brangus cattle managed in a Chihuahuan Desert production system1.

Balancing growth and reproductive performance in beef cattle managed in desert environments is challenging. Our objectives were to 1) evaluate trends in growth and reproductive traits, and 2) assess associative relationships between growth characteristics and reproductive performance in a Brangus herd managed in a Chihuahuan Desert production system from 1972 to 2006. Data were from bull (n = 597) and heifer calves (n = 585; 1988 to 2006) and cows (n = 525; repeated records of cows, n = 2,611; 1972 to 2006). Variables describing the growth curve of each cow were estimated using a nonlinear logistic function (each cow needed 6 yr of data). Mixed-effect models and logistic regression were used to analyze trends across years in growth and reproductive traits (both continuous and categorical). For continuous traits of calves, a slight cubic response (P < 0.01) described the dynamics of birth weight, 205-d BW, and 365-d BW across years. For categorical traits of females, positive linear trends (P < 0.05) across years were observed in percent pregnant as yearlings, calved at 2 yr of age, and first-calf heifer rebreeding (slopes ranged from 0.007 to 0.014%/yr). Autumn cow BW increased gradually until 1997 (509 kg +/- 8.8) and then decreased gradually by 0.6 kg/yr, whereas pregnancy percentage decreased gradually until 1995 (78.4% +/- 1.0) and then increased slightly by 0.2%/yr. A quadratic effect best described the dynamics of these 2 variables across years (P < 0.01) as well as estimates describing the growth curve of each cow. Specifically, asymptotic BW and age increased (P < 0.05) from 1972 to 1983 and 1990, respectively. Asymptotic age then decreased by 27% from 1983 to 1996 (P < 0.05). The maturing rate index was negatively correlated with age at first calving and calving interval (r = -0.42 and -0.18, P < 0.01), which suggested that early-maturing cows had enhanced fertility in this environment and production system. In summary, minimal changes were observed in measures of growth in bulls and heifers in a Brangus herd managed in the Chihuahuan Desert. Opposing relationships were observed among measures of cow size and fertility; as growth curves shifted toward earlier maturity, measures of reproductive performance suggested that fertility improved.