Pre- and post-puberty expression of genes and proteins in the uterus of Bos indicus heifers: the luteal phase effect post-puberty.

Progesterone signaling and uterine function are crucial in terms of pregnancy establishment. To investigate how the uterine tissue and its secretion changes in relation to puberty, we sampled tissue and uterine fluid from six pre- and six post-pubertal Brahman heifers. Post-pubertal heifers were sampled in the luteal phase. Gene expression of the uterine tissue was investigated with RNA-sequencing, whereas the uterine fluid was used for protein profiling with mass spectrometry. A total of 4034 genes were differentially expressed (DE) at a nominal P-value of 0.05, and 26 genes were significantly DE after Bonferroni correction (P < 3.1 × 10-6 ). We also identified 79 proteins (out of 230 proteins) that were DE (P < 1 × 10-5 ) in the uterine fluid. When we compared proteomics and transcriptome results, four DE proteins were identified as being encoded by DE genes: OVGP1, GRP, CAP1 and HBA. Except for CAP1, the other three had lower expression post-puberty. The function of these four genes hypothetically related to preparation of the uterus for a potential pregnancy is discussed in the context of puberty. All DE genes and proteins were also used in pathway and ontology enrichment analyses to investigate overall function. The DE genes were enriched for terms related to ribosomal activity. Transcription factors that were deemed key regulators of DE genes are also reported. Transcription factors ZNF567, ZNF775, RELA, PIAS2, LHX4, SOX2, MEF2C, ZNF354C, HMG20A, TCF7L2, ZNF420, HIC1, GTF3A and two novel genes had the highest regulatory impact factor scores. These data can help to understand how puberty influences uterine function.

Prospecting major genes in dairy buffaloes.

BACKGROUND: Asian buffaloes (Bubalus bubalis) have an important socio-economic role. The majority of the population is situated in developing countries. Due to the scarce resources in these countries, very few species-specific biotechnology tools exist and a lot of cattle-derived technologies are applied to buffaloes. However, the application of cattle genomic tools to buffaloes is not straightforward and, as results suggested, despite genome sequences similarity the genetic polymorphisms are different. RESULTS: The first SNP chip genotyping platform designed specifically for buffaloes has recently become available. Herein, a genome-wide association study (GWAS) and gene network analysis carried out in buffaloes is presented. Target phenotypes were six milk production and four reproductive traits. GWAS identified SNP with significant associations and suggested candidate genes that were specific to each trait and also genes with pleiotropic effect, associated to multiple traits. CONCLUSIONS: Network predictions of interactions between these candidate genes may guide further molecular analyses in search of disruptive mutations, help select genes for functional experiments and evidence metabolism differences in comparison to cattle. The cattle SNP chip does not offer an optimal coverage of buffalo genome, thereafter the development of new buffalo-specific genetic technologies is warranted. An annotated reference genome would greatly facilitate genetic research, with potential impact to buffalo-based dairy production.

Evidence for pleiotropism and recent selection in the PLAG1 region in Australian Beef cattle.

A putative functional mutation (rs109231213) near PLAG1 (BTA14) associated with stature was studied in beef cattle. Data from 8199 Bos taurus, Bos indicus and Tropical Composite cattle were used to test the associations between rs109231213 and various phenotypes. Further, 23 496 SNPs located on BTA14 were tested for association with these phenotypes, both independently and fitted together with rs109231213. The C allele of rs109231213 significantly increased hip height, weight, net food intake, age at puberty in males and females and decreased IGF-I concentration in blood and fat depth. When rs109231213 was fitted as a fixed effect in the model, there was an overall reduction in associations between other SNPs and these traits but some SNPs remained associated (P < 10(-4) ). Frequency of the mutant C allele of rs109231213 differed among B. indicus (0.52), B. taurus (0.96) and Tropical Composite (0.68). Most chromosomes carrying the C allele had the same surrounding 10 SNP haplotype, probably because the C allele was introgressed into Brahman from B. taurus cattle. A region of reduced heterozygosity surrounds the C allele; this is small in B. taurus but 20 Mb long in Brahmans, indicating recent and strong selection for the mutant allele. Thus, the C allele appears to mark a mutation that has been selected almost to fixation in the B. taurus breeds studied here and introduced into Brahman cattle during grading up and selected to a frequency of 0.52 despite its negative effects on fertility.

Net protein contribution and enteric methane production of pasture and grain-finished beef cattle supply chains.

Ruminant red meat production systems around the world often include a grain feeding phase. The role of red meat in the food system is therefore often discussed in terms of the food vs feed debate, as well as invoking the comparatively poor feed conversion efficiency of ruminants and climate impacts from enteric methane. The concept of net protein contribution (NPC) incorporates the quality attributes of protein produced by livestock systems into estimates of the efficiency of production systems. We applied the NPC method to two Australian beef supply chains, i) Grass-fed and ii) Grain-finished beef, using an established model of ruminant grazing systems (GrassGro®) and these are reflective of beef production systems in other countries. The beef supply chains evaluated did not compete with humans for protein. The Grain-finished beef supply chain, while positively contributing to human protein requirements (NPC value 1.96), had markedly lower NPC values than the Grass-fed system (NPC value 1 597). However, Grass-fed beef production systems have a higher methane intensity than the Grain-finished supply chain. The two examples of pasture-based beef production systems examined provide a positive net protein contribution to human food supply, even with extended periods of finishing on grain-based diets. This is achieved by ruminant grazing on pastures converting low-quality forage into high value human edible protein. The efficiency of protein production varies according to the system design, and other considerations such as land use and enteric methane production are elements that need consideration in the overall assessment of the production footprint.

Evaluation of non-additive genetic variation in feed-related traits of broiler chickens.

Genome-wide association mapping and genomic predictions of phenotype of individuals in livestock are predominately based on the detection and estimation of additive genetic effects. Non-additive genetic effects are largely ignored. Studies in animals, plants, and humans to assess the impact of non-additive genetic effects in genetic analyses have led to differing conclusions. In this paper, we examined the consequences of including non-additive genetic effects in genome-wide association mapping and genomic prediction of total genetic values in a commercial population of 5,658 broiler chickens genotyped for 45,176 single nucleotide polymorphism (SNP) markers. We employed mixed-model equations and restricted maximum likelihood to analyze 7 feed related traits (TRT1 - TRT7). Dominance variance accounted for a significant proportion of the total genetic variance in all 7 traits, ranging from 29.5% for TRT1 to 58.4% for TRT7. Using a 5-fold cross-validation schema, we found that in spite of the large dominance component, including the estimated dominance effects in the prediction of total genetic values did not improve the accuracy of the predictions for any of the phenotypes. We offer some possible explanations for this counter-intuitive result including the possible confounding of dominance deviations with common environmental effects such as hatch, different directional effects of SNP additive and dominance variations, and the gene-gene interactions' failure to contribute to the level of variance.

Genomic inbreeding depression for climatic adaptation of tropical beef cattle.

Inbreeding has the potential to negatively impact animal performance. Strategies to monitor and mitigate inbreeding depression require that it can be accurately estimated. Here, we used genomewide SNP data to explore 3 alternative measures of genomic inbreeding: the diagonal elements of the genomic relationship matrix (FGRM), the proportion of homozygous SNP (FHOM), and the proportion of the genome covered by runs of homozygosity (FROH). We used 2,111 Brahman (BR) and 2,550 Tropical Composite (TC) cattle with phenotypes recorded for 10 traits of relevance to tropical adaptation. We further explored 3 marker densities ranging from a high-density chip (729,068 SNP), a medium-density chip (71,726 SNP) specifically designed for cattle, and a low-density chip (18,860 SNP) associated with the measures of inbreeding. Measures of FGRM were highly correlated across the 3 SNP densities and negatively correlated with FHOM and FROH in the BR population. In both populations, there was a strong positive correlation for each measure of inbreeding across the 3 SNP panels. We found significant ( < 0.01) inbreeding depression for various traits, particularly when using the highest-density SNP chip in the BR population, where inbreeding was negatively associated with coat color and coat type such that inbred animals presented shorter, slicker, and lighter coats. Based on FGRM using the medium-density chip, we found that a 1% increase in inbreeding in the BR and TC populations was associated with a decrease of 0.514 and 0.579 kg BW, respectively, in yearlings. In the TC population, a 1% increase in FHOM was associated with a decrease in BCS of -0.636% ( < 0.001). The low-density chip, comprising SNP associated with inbreeding, captured genes, and regions with pleiotropic effects ( < 0.001). However, it did not improve our ability to identify inbreeding depression, relative to the use of higher-density panels. We conclude that where heterogeneous populations are present, such as in tropical environments where composite animals abound, measures of inbreeding that do not depend on allele frequencies, such as FHOM and FROH, are preferable for estimating genomic inbreeding. Finally, the sustainable intensification of livestock systems in tropical regions will rely on genetic safeguards to ensure that productivity is improved while also adapting animals to cope with climate change. The results of this study are a step toward achieving that goal.

Global differential gene expression in the pituitary gland and the ovaries of pre- and postpubertal Brahman heifers.

To understand genes, pathways, and networks related to puberty, we characterized the transcriptome of two tissues: the pituitary gland and ovaries. Samples were harvested from pre- and postpubertal Brahman heifers (same age group). Brahman heifers () are older at puberty compared with , a productivity issue. With RNA sequencing, we identified differentially expressed (DEx) genes and important transcription factors (TF) and predicted coexpression networks. The number of DEx genes detected in the pituitary gland was 284 ( < 0.05), and was the most DEx gene (fold change = 4.12, = 0.01). The gene promotes bone mineralization through transforming growth factor-ß (TGFß) signaling. Further studies of the link between bone mineralization and puberty could target . In ovaries, 3,871 genes were DEx ( < 0.05). Four highly DEx genes were noteworthy for their function: (a γ-aminobutyric acid [GABA] transporter), (), and () and its receptor . These genes had higher ovarian expression in postpubertal heifers. The GABA and its receptors and transporters were expressed in the ovaries of many mammals, suggesting a role for this pathway beyond the brain. The pathway has been known to influence the timing of puberty in rats, via modulation of GnRH. The effects of at the hypothalamus, pituitary gland, and ovaries have been documented. and its receptors are known factors in the release of GnRH, similar to and GABA, although their roles in ovarian tissue are less clear. Pathways previously related to puberty such as TGFß signaling ( = 6.71 × 10), Wnt signaling ( = 4.1 × 10), and peroxisome proliferator-activated receptor (PPAR) signaling ( = 4.84 × 10) were enriched in our data set. Seven genes were identified as key TF in both tissues: , , , , , , and a novel gene. An ovarian subnetwork created with TF and significant ovarian DEx genes revealed five zinc fingers as regulators: , , , , and . Recent work of hypothalamic gene expression also pointed to zinc fingers as TF for bovine puberty. Although some zinc fingers may be ubiquitously expressed, the identification of DEx genes in common across tissues points to key regulators of puberty. The hypothalamus and pituitary gland had eight DEx genes in common. The hypothalamus and ovaries had 89 DEx genes in common. The pituitary gland and ovaries had 48 DEx genes in common. Our study confirmed the complexity of puberty and suggested further investigation on genes that code zinc fingers.

Genomic analyses of tropical beef cattle fertility based on genotyping pools of Brahman cows with unknown pedigree.

We introduce an innovative approach to lowering the overall cost of obtaining genomic EBV (GEBV) and encourage their use in commercial extensive herds of Brahman beef cattle. In our approach, the DNA genotyping of cow herds from 2 independent properties was performed using a high-density bovine SNP chip on DNA from pooled blood samples, grouped according to the result of a pregnancy test following their first and second joining opportunities. For the DNA pooling strategy, 15 to 28 blood samples from the same phenotype and contemporary group were allocated to pools. Across the 2 properties, a total of 183 pools were created representing 4,164 cows. In addition, blood samples from 309 bulls from the same properties were also taken. After genotyping and quality control, 74,584 remaining SNP were used for analyses. Pools and individual DNA samples were related by means of a "hybrid" genomic relationship matrix. The pooled genotyping analysis of 2 large and independent commercial populations of tropical beef cattle was able to recover significant and plausible associations between SNP and pregnancy test outcome. We discuss 24 SNP with significant association ( < 1.0 × 10) and mapped within 40 kb of an annotated gene. We have established a method to estimate the GEBV in young herd bulls for a trait that is currently unable to be predicted at all. In summary, our novel approach allowed us to conduct genomic analyses of fertility in 2 large commercial Brahman herds managed under extensive pastoral conditions.

Transcriptome analyses identify five transcription factors differentially expressed in the hypothalamus of post- versus prepubertal Brahman heifers.

Puberty onset is a developmental process influenced by genetic determinants, environment, and nutrition. Mutations and regulatory gene networks constitute the molecular basis for the genetic determinants of puberty onset. The emerging knowledge of these genetic determinants presents opportunities for innovation in the breeding of early pubertal cattle. This paper presents new data on hypothalamic gene expression related to puberty in (Brahman) in age- and weight-matched heifers. Six postpubertal heifers were compared with 6 prepubertal heifers using whole-genome RNA sequencing methodology for quantification of global gene expression in the hypothalamus. Five transcription factors (TF) with potential regulatory roles in the hypothalamus were identified in this experiment: , , , , and . These TF genes were significantly differentially expressed in the hypothalamus of postpubertal versus prepubertal heifers and were also identified as significant according to the applied regulatory impact factor metric ( < 0.05). Two of these 5 TF, and , were zinc fingers, belonging to a gene family previously reported to have a central regulatory role in mammalian puberty. The gene belongs to the family of homologues of Drosophila sine oculis () genes implicated in transcriptional regulation of gonadotrope gene expression. Tumor-related genes such as and are known to affect basic cellular processes that are relevant in both cancer and developmental processes. Mutations in were associated with puberty in humans. Mutations in these TF, together with other genetic determinants previously discovered, could be used in genomic selection to predict the genetic merit of cattle (i.e., the likelihood of the offspring presenting earlier than average puberty for Brahman). Knowledge of key mutations involved in genetic traits is an advantage for genomic prediction because it can increase its accuracy.

Genome-wide association for the outcome of fixed-time artificial insemination of Brahman heifers in northern Australia.

Fixed-time AI (FTAI) is a powerful tool for genetic improvement of extensively managed beef cattle. A genomewide association study (GWAS) was conducted to investigate genes and genetic markers associated with the outcome (pregnant or not pregnant) of FTAI in 614 commercial Brahman heifers genotyped for 18,895 SNP and imputed to 51,588 SNP. The likelihood of Brahman heifers becoming pregnant after hormonal treatment to synchronize ovulation followed by FTAI was influenced by the content of their genomes, as determined by a principal component analysis. The principal component analysis involved comparisons between the studied heifers and populations of known and ancestry. The heritability of FTAI outcome was = 0.18, which is higher than for most other reproductive outcome traits. The number of SNP associated with FTAI outcome was 101 ( < 0.001, false discovery rate = 0.53). Compared with all SNP tested, associated SNP had a tendency for highly divergent allelic frequencies between and . Associated SNP were located in nearly all chromosomes, a result that shows a complex genetic architecture that is typical of highly complex traits with low heritability. Considering this and previous GWAS that examined Brahman heifer puberty and postpartum anestrus interval, 3 genomic regions emerge as important for overall Brahman heifer fertility, which mapped to chromosomes 1, 7, and 9. Further analyses, including improved genome annotation, are required to elucidate the link between these regions and heifer fertility. Additional studies are needed to confirm SNP and gene associations reported herein and further elucidate the genetics of FTAI outcome. Future GWAS should target other Braham populations and additional cattle breeds with FTAI records, including breeds with higher ancestry.

An engineered PGK promoter and lac operator-repressor system for the regulation of gene expression in mammalian cells.

Previous reports have demonstrated that the Escherichia coli lac repressor can operate effectively in mammalian cells to repress expression of genes driven by modified viral or metallothionein (MT) promoters. We have developed a more general expression system using the promoter from the PGK1 gene (encoding murine 3-phosphoglycerate kinase) which is widely expressed in almost all cell types, including early embryonic and ES (embryonic stem) cells. Firstly, we engineered the lac repressor to include a nuclear localisation signal and placed it under control of the PGK1 promoter. Efficient nuclear localisation of the repressor was demonstrated by mobility-shift assays and immunofluorescence detection. For the target vectors, we modified the wild-type (wt) PGK1 promoter to include lac operator (lacO) sites for binding of the lac repressor and compared a number of different lacO positions and arrangements based on proximity to the native start points for transcription (tsp) and translation. In the absence of repressor, we observed reduced expression of the neo reporter gene for some placements of the lacO, but wt expression for placements near the tsp. When both target and repressor were present in the cells, we observed that the expression of neo could be strongly suppressed and reversibly regulated by induction with IPTG. In particular, for a promoter which contained two spaced lacO replacing native sequence around the major tsp, we observed 90-95% repression by the lac repressor for the neo reporter gene and up to 98% repression for the cat reporter gene. Efficient derepression by IPTG was observed in both cases.(ABSTRACT TRUNCATED AT 250 WORDS)

Isolation of a cDNA encoding a putative cellulase in the red claw crayfish Cherax quadricarinatus.

Amino acid sequences of cellulases have been determined in insects, nematodes, plants, slime moulds and bacteria but not in crustaceans. However, cellulase activity has been demonstrated in the hepatopancreas of the red claw crayfish, Cherax quadricarinatus. In order to obtain information on the nature of this cellulase, a C. quadricarinatus hepatopancreas cDNA library was screened with a PCR product generated using degenerate oligonucleotide primers derived from conserved regions of known cellulases. Two identical 1.56kb cDNAs with sequence similarities to known cellulases, particularly the termite endoglucanases, were identified and sequenced. The clones contain the complete cDNA open reading frame for an endo-1, 4-beta-glucanase of 469 amino acids termed Cherax quadricarinatus endoglucanase (CqEG). The endogenous origin of the gene was confirmed by PCR amplification and sequencing of a 1012bp PCR product from genomic DNA. This fragment contains four exon sequences identical to the cDNA and is interrupted by three introns of 371, 102, 194bp respectively, with one intron exhibiting typical eukaryotic splice sites. The isolation of an endo-1,4-beta-glucanase encoding cDNA from the crayfish C. quadricarinatus provides the first endogenous cellulase sequence in a crustacean species.

A mixed-model approach for the analysis of cDNA microarray gene expression data from extreme-performing pigs after infection with Actinobacillus pleuropneumoniae.

We proposed a novel statistical approach for the analysis of cDNA experiments based on mixed-model methodology combined with mixtures of distributions. Our objective was to detect genes that may be involved in conferring heritable differences in susceptibility to common infections in intensive pig production. We employed a microarray expression profiling strategy and a mixed-model approach to the analysis of the expression data. A cDNA microarray of pig with 6,420 probes from immune tissues and cells was used to compare gene expression in peripheral blood leukocytes of two pigs showing extreme performance in their response to infection with Actinobacillus pleuropneumoniae. Principal components analyses were used to identify the two most extreme-performing pigs after infection (i.e., pigs whose measured responses to infection fell at the extremes). Blood samples and expression profiles from 0 to 24 h after infection were compared using a bivariate, mixed-model approach, in which the effect gene x immunological status interaction was treated as a random effect. Bayesian model-based clustering via mixtures of normal distributions of the resulting BLUP of the random interaction was approached and resulted in a list of 307 differentially expressed genes, of which 179 were down-regulated in the susceptible pig. The majority of the differentially expressed genes were derived from a cDNA library of leukocytes of A. pleuropneumoniae-challenged pigs that were subtracted against leukocytes before the challenge. These results provide evidence that the proposed statistical approach was useful in enhancing the knowledge of the mechanisms involved in the genetics of the immune response.

A mixture model-based cluster analysis of DNA microarray gene expression data on Brahman and Brahman composite steers fed high-, medium-, and low-quality diets.

The objective of this study is to explore aspects of the statistical analysis of gene expression response at the muscle tissue level to varying levels of energy and protein in the diet. Eleven Brahman and Brahman composite steers (weighing 302 +/- 9.8 kg, on average) were allocated randomly into high- (HIGH), medium- (MED), and low- (LOW) quality forage diets for 27 d. After this period, a biopsy of the longissimus dorsi muscle was taken from each animal and total RNA was extracted to generate the labeled target for microarray experimentation. These targets were hybridized to a complementary DNA (cDNA) microarray of 9,274 probes from cattle muscle and subcutaneous fat cDNA libraries. After edits, 151,904 expression intensity levels of 4,747 genes were analyzed. Emphasis was given to the choice of power transformation of the intensity channel readings and to the consistency of readings within each diet quality group. The statistical approach to isolate differentially expressed genes was based on model-based clustering via a mixture of normal distributions estimated through maximal likelihood. The base-2 logarithm was found to be the optimal power transformation to normalize gene intensity levels. A two-sample t-statistic was defined as a measure of possible differential expression. For each of the three diet contrasts, HIGH vs. LOW, HIGH vs. MED, and MED vs. LOW, three clusters were found, two of which contained more than 94% genes with almost no altered gene expression levels, whereas the third cluster contained the remaining genes with a differential expression. Results from the HIGH vs. LOW contrast identified 27 genes with a greater than 95% posterior probability of belonging to the cluster of differentially expressed genes.

Numerical analysis of intensity signals resulting from genotyping pooled DNA samples in beef cattle and broiler chicken.

Pooled genomic DNA has been proposed as a cost-effective approach in genomewide association studies (GWAS). However, algorithms for genotype calling of biallelic SNP are not adequate with pooled DNA samples because they assume the presence of 2 fluorescent signals, 1 for each allele, and operate under the expectation that at most 2 copies of the variant allele can be found for any given SNP and DNA sample. We adapt analytical methodology from 2-channel gene expression microarray technology to SNP genotyping of pooled DNA samples. Using 5 datasets from beef cattle and broiler chicken of varying degrees of complexity in terms of design and phenotype, continuous and dichotomous, we show that both differential hybridization (M = green minus red intensity signal) and abundance (A = average of red and green intensities) provide useful information in the prediction of SNP allele frequencies. This is predominantly true when making inference about extreme SNP that are either nearly fixed or highly polymorphic. We propose the use of model-based clustering via mixtures of bivariate normal distributions as an optimal framework to capture the relationship between hybridization intensity and allele frequency from pooled DNA samples. The range of M and A values observed here are in agreement with those reported within the context of gene expression microarray and also with those from SNP array data within the context of analytical methodology for the identification of copy number variants. In particular, we confirm that highly polymorphic SNP yield a strong signal from both channels (red and green) while lowly or nonpolymorphic SNP yield a strong signal from 1 channel only. We further confirm that when the SNP allele frequencies are known, either because the individuals in the pools or from a closely related population are themselves genotyped, a multiple regression model with linear and quadratic components can be developed with high prediction accuracy. We conclude that when these approaches are applied to the estimation of allele frequencies, the resulting estimates allow for the development of cost-effective and reliable GWAS.

A marker-derived gene network reveals the regulatory role of PPARGC1A, HNF4G, and FOXP3 in intramuscular fat deposition of beef cattle.

High intramuscular fat (IMF) awards price premiums to beef producers and is associated with meat quality and flavor. Studying gene interactions and pathways that affect IMF might unveil causative physiological mechanisms and inform genomic selection, leading to increased accuracy of predictions of breeding value. To study gene interactions and pathways, a gene network was derived from genetic markers associated with direct measures of IMF, other fat phenotypes, feedlot performance, and a number of meat quality traits relating to body conformation, development, and metabolism that might be plausibly expected to interact with IMF biology. Marker associations were inferred from genomewide association studies (GWAS) based on high density genotypes and 29 traits measured on 10,181 beef cattle animals from 3 breed types. For the network inference, SNP pairs were assessed according to the strength of the correlation between their additive association effects across the 29 traits. The co-association inferred network was formed by 2,434 genes connected by 28,283 edges. Topological network parameters suggested a highly cohesive network, in which the genes are strongly functionally interconnected. Pathway and network analyses pointed towards a trio of transcription factors (TF) as key regulators of carcass IMF: PPARGC1A, HNF4G, and FOXP3. Importantly, none of these genes would have been deemed as significantly associated with IMF from the GWAS. Instead, a total of 313 network genes show significant co-association with the 3 TF. These genes belong to a wide variety of biological functions, canonical pathways, and genetic networks linked to IMF-related phenotypes. In summary, our GWAS and network predictions are supported by the current literature and suggest a cooperative role for the 3 TF and other interacting genes including CAPN6, STC2, MAP2K4, EYA1, COPS5, XKR4, NR2E1, TOX, ATF1, ASPH, TGS1, and TTPA as modulators of carcass and meat quality traits in beef cattle.

Genome-wide association study for inhibin, luteinizing hormone, insulin-like growth factor 1, testicular size and semen traits in bovine species.

The fertility of young bulls impacts on reproduction rates, farm profit and the rate of genetic progress in beef herds. Cattle researchers and industry therefore routinely collect data on the reproductive performance of bulls. Genome-wide association studies were carried out to identify genomic regions and genes associated with reproductive traits measured during the pubertal development of Tropical Composite bulls, from 4 to 24 months of age. Data from 1 085 bulls were collected for seven traits: blood hormone levels of inhibin at 4 months (IN), luteinizing hormone following a gonadotropin releasing hormone challenge at 4 months (LH), insulin-like growth factor 1 at 6 months (IGF1), scrotal circumference at 12 months (SC), sperm motility at 18 months (MOT), percentage of normal spermatozoa at 24 months (PNS) and age at a scrotal circumference of 26 cm (AGE26, or pubertal age). Data from 729 068 single-nucleotide polymorphisms were used in the association analysis. Significant polymorphism associations were discovered for IN, IGF1, SC, AGE26 and PNS. Based on these associations, INHBE, INHBC and HELB are proposed as candidate genes for IN regulation. Polymorphisms associated with IGF1 mapped to the PLAG1 gene region, validating a reported quantitative trait locus on chromosome 14 for IGF1. The X chromosome contained most of the significant associations found for SC, AGE26 and PNS. These findings will contribute to the identification of diagnostic genetic markers and informed genomic selection strategies to assist breeding of cattle with improved fertility. Furthermore, this work provides evidence contributing to gene function annotation in the context of male fertility.

Transplanted germ cells persist long-term in irradiated ram testes.

Testicular germ cell transplantation provides a tool to study transgenesis, spermatogenesis and to increase production efficiency in livestock industries. Isolated testicular germ cells can be transplanted into testes of livestock breeds to generate sperm of donor origin. In sheep, methods have been developed previously to isolate cell populations from ram testes and transplant these into irradiated testes of recipient rams. This has resulted in rams producing sperm derived from the donor cells and a number of the recipient animals have produced donor-derived offspring from the introduced spermatogonial cells. Microsatellite genotyping data presented here demonstrates that these rams continue to produce sperm of donor origin for at least 5 years post-transplantation. This research provides new evidence of the stability of transplanted germ cells in a commercially important species, and with further refinements to cell isolation, transplantation and recipient preparation, this technology should find use in breeding systems to increase livestock production efficiency.

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.

Genome-wide association studies of female reproduction in tropically adapted beef cattle.

The genetics of reproduction is poorly understood because the heritabilities of traits currently recorded are low. To elucidate the genetics underlying reproduction in beef cattle, we performed a genome-wide association study using the bovine SNP50 chip in 2 tropically adapted beef cattle breeds, Brahman and Tropical Composite. Here we present the results for 3 female reproduction traits: 1) age at puberty, defined as age in days at first observed corpus luteum (CL) after frequent ovarian ultrasound scans (AGECL); 2) the postpartum anestrous interval, measured as the number of days from calving to first ovulation postpartum (first rebreeding interval, PPAI); and 3) the occurrence of the first postpartum ovulation before weaning in the first rebreeding period (PW), defined from PPAI. In addition, correlated traits such as BW, height, serum IGF1 concentration, condition score, and fatness were also examined. In the Brahman and Tropical Composite cattle, 169 [false positive rate (FPR) = 0.262] and 84 (FPR = 0.581) SNP, respectively, were significant (P < 0.001) for AGECL. In Brahman, 41% of these significant markers mapped to a single chromosomal region on BTA14. In Tropical Composites, 16% of these significant markers were located on BTA5. For PPAI, 66 (FPR = 0.67) and 113 (FPR = 0.432) SNP were significant (P < 0.001) in Brahman and Tropical Composite, respectively, whereas for PW, 68 (FPR = 0.64) and 113 (FPR = 0.432) SNP were significant (P < 0.01). In Tropical Composites, the largest concentration of PPAI markers were located on BTA5 [19% (PPAI) and 23% (PW)], and BTA16 [17% (PPAI) and 18% (PW)]. In Brahman cattle, the largest concentration of markers for postpartum anestrus was located on BTA3 (14% for PPAI and PW) and BTA14 (17% PPAI). Very few of the significant markers for female reproduction traits for the Brahman and Tropical Composite breeds were located in the same chromosomal regions. However, fatness and BW traits as well as serum IGF1 concentration were found to be associated with similar genome regions within and between breeds. Clusters of SNP associated with multiple traits were located on BTA14 in Brahman and BTA5 in Tropical Composites.