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.

A single nucleotide polymorphism-derived regulatory gene network underlying puberty in 2 tropical breeds of beef cattle.

Harsh tropical environments impose serious challenges on poorly adapted species. In beef cattle, tropical adaptation in the form of temperature and disease resistance, coupled with acclimatization to seasonal and limited forage, comes at a cost to production efficiency. Prominent among these costs is delayed onset of puberty, a challenging phenotype to manipulate through traditional breeding mechanisms. Recently, system biology approaches, including gene networks, have been applied to the genetic dissection of complex phenotypes. We aimed at developing and studying gene networks underlying cattle puberty. Our starting material comprises the association results of ~50,000 SNP on 22 traits, including age at puberty, and 2 cattle breed populations: Brahman (n = 843) and Tropical Composite (n = 866). We defined age at puberty as the age at first corpus luteum (AGECL). By capturing the genes harboring mutations minimally associated (P < 0.05) to AGECL or to a set of traits related with AGECL, we derived a gene network for each breed separately and a third network for the combined data set. At the intersection of the 3 networks, we identified candidate genes and pathways that were common to both breeds. Resulting from these analyses, we identified an enrichment of genes involved in axon guidance, cell adhesion, ErbB signaling, and glutamate activity, pathways that are known to affect pulsatile release of GnRH, which is necessary for the onset of puberty. Furthermore, we employed network connectivity and centrality parameters along with a regulatory impact factor metric to identify the key transcription factors (TF) responsible for the molecular regulation of puberty. As a novel finding, we report 5 TF (HIVEP3, TOX, EYA1, NCOA2, and ZFHX4) located in the network intersecting both breeds and interacting with other TF, forming a regulatory network that harmonizes with the recent literature of puberty. Finally, we support our network predictions with evidence derived from gene expression in hypothalamic tissue of adult cows.

Duplication and divergence: the evolution of nematode globins.

In common with many other groups, nematodes express globins with unknown functions. Nematode globin-like genes can be divided into class 1 globins, similar to vertebrate myoglobins, and a wide range of additional classes. Here we show that class 1 nematode globins possess a huge amount of diversity in gene sequence and structure. There is evidence for multiple events of gene duplication, intron insertion and loss between species, and for allelic variation effecting both synonymous and non-synonymous sites within species. We have also examined gene expression patterns in class I globins from a variety of species. The results show variation in the degree of gene expression, but the tissue specificity and temporal specificity of expression may be more conserved in the phylum. Because the structure-function relationships for the binding and transport of oxygen by globins are well understood, the consequences of genetic variation causing amino acid changes are explored. The gene family shows great promise for discovering unique insights into both structure-function relationships of globins and their physiologial roles.

Predicting genetic merit for mastitis and fertility in dairy cattle using genome wide selection and high density SNP screens.

Two novel methods for genome wide selection (GWS) were examined for predicting the genetic merit of animals using SNP information alone. A panel of 1,546 dairy bulls with reliable EBVs was genotyped for 15,380 SNPs that spanned the whole bovine genome. Two complexity reduction methods were used, partial least squares (PLS) and regression using a genetic algorithm (GAR), to find optimal solutions of EBVs against SNP information. Extensive internal cross-validation was used tofind the best predictive models followed by external validation (without direct use of the pedigree or SNP location). Both PLS and GAR provided both accurate fit to the training data set for somatic cell count (SCC) (max r = 0.83) and fertility (max r = 0.88) and showed an accuracy of prediction of r = 0.47 for SCC, and r = 0.72 for fertility. This is the first empirical demonstration that genome wide selection can account for a very high proportion of additive genetic variation in fitness traits whilst exploiting only a small percentage of available SNP information, without use of pedigree or QTL mapping. PLS was computationally more efficient than GAR.

A validated whole-genome association study of efficient food conversion in cattle.

The genetic factors that contribute to efficient food conversion are largely unknown. Several physiological systems are likely to be important, including basal metabolic rate, the generation of ATP, the regulation of growth and development, and the homeostatic control of body mass. Using whole-genome association, we found that DNA variants in or near proteins contributing to the background use of energy of the cell were 10 times as common as those affecting appetite and body-mass homeostasis. In addition, there was a genic contribution from the extracellular matrix and tissue structure, suggesting a trade-off between efficiency and tissue construction. Nevertheless, the largest group consisted of those involved in gene regulation or control of the phenotype. We found that the distribution of micro-RNA motifs was significantly different for the genetic variants associated with residual feed intake than for the genetic variants in total, although the distribution of promoter sequence motifs was not different. This suggests that certain subsets of micro-RNA are more important for the regulation of this trait. Successful validation depended on the sign of the allelic association in different populations rather than on the strength of the initial association or its size of effect.

Skeletal muscle extracellular matrix remodelling after aestivation in the green striped burrowing frog, Cyclorana alboguttata.

Connective tissue has recently been found to play a role in mediating mammalian skeletal muscle atrophy. We investigated connective tissue remodelling in the skeletal muscle of a species of the Australian burrowing frog, Cyclorana alboguttata. Despite being inactive whilst aestivating, the frog shows an inhibition of muscle atrophy. Connective tissue size and distribution was measured in histological sections of the cruralis muscle of control and aestivating C. alboguttata. Using a custom written software application we could detect no significant difference in any connective tissue morphological parameter between the two treatment groups. Biochemical measurements of gelatinase activity showed 2-fold higher activity in aestivating gastrocnemius muscle than in controls (p<0.001). We measured the messenger RNA transcript levels for C. alboguttata metalloproteinase 2 (MMP2) and tissue inhibitor of metalloproteinase 2 (TIMP2) in cruralis skeletal muscle using quantitative real-time PCR. The trend of reduced expression of the two genes in the aestivators did not meet statistical significance. This work indicates that aestivation in C. alboguttata leads to subtle and specific changes in some extracellular matrix remodelling factors. Their main impact is to maintain proportional representation of extracellular matrix components of skeletal muscle and therefore preserve the active frog phenotype.

A rapid method for computationally inferring transcriptome coverage and microarray sensitivity.

MOTIVATION: There are many different gene expression technologies, including cDNA and oligo-based microarrays, SAGE and MPSS. For each organism of interest, coverage of the transcriptome and the genome will be different. We address the question of what level of coverage is required to exploit the sensitivity of the different technologies, and what is the sensitivity of the different approaches in the experimental study. RESULTS: We estimate the transcriptome coverage by randomly sampling transcripts from a pre-defined tag-to-gene mapping function. For a given microarray experiment, we locate the thresholds in intensities that define the distribution of transcript abundance. These values are compared against the distribution obtained by applying the same thresholds to the intensities from differentially expressed genes. The ratio of these two distributions meets at the equilibrium defining sensitivity. We conclude that a collection of approximately 340,000 sequences is adequate for microarrays, but not large enough for maximum utilization of tag-based technologies. In the absence of large-scale sequencing, the majority of the tags detected by the latter approaches will remain unidentified until the genome sequence is available.