Multi-breed genomic predictions and functional variants for fertility of tropical bulls.

Worldwide, most beef breeding herds are naturally mated. As such, the ability to identify and select fertile bulls is critically important for both productivity and genetic improvement. Here, we collected ten fertility-related phenotypes for 6,063 bulls from six tropically adapted breeds. Phenotypes were comprised of four bull conformation traits and six traits directly related to the quality of the bull's semen. We also generated high-density DNA genotypes for all the animals. In total, 680,758 single nucleotide polymorphism (SNP) genotypes were analyzed. The genomic correlation of the same trait observed in different breeds was positive for scrotal circumference and sheath score on most breed comparisons, but close to zero for the percentage of normal sperm, suggesting a divergent genetic background for this trait. We confirmed the importance of a breed being present in the reference population to the generation of accurate genomic estimated breeding values (GEBV) in an across-breed validation scenario. Average GEBV accuracies varied from 0.19 to 0.44 when the breed was not included in the reference population. The range improved to 0.28 to 0.59 when the breed was in the reference population. Variants associated with the gene HDAC4, six genes from the spermatogenesis-associated (SPATA) family of proteins, and 29 transcription factors were identified as candidate genes. Collectively these results enable very early in-life selection for bull fertility traits, supporting genetic improvement strategies currently taking place within tropical beef production systems. This study also improves our understanding of the molecular basis of male fertility in mammals.

Genomic correlation: harnessing the benefit of combining two unrelated populations for genomic selection.

BACKGROUND: The success of genomic selection in animal breeding hinges on the availability of a large reference population on which genomic-based predictions of additive genetic or breeding values are built. Here, we explore the benefit of combining two unrelated populations into a single reference population. METHODS: The datasets consisted of 1829 Brahman and 1973 Tropical Composite cattle with measurements on five phenotypes relevant to tropical adaptation and genotypes for 71,726 genome-wide single nucleotide polymorphisms (SNPs). The underlying genomic correlation for the same phenotype across the two breeds was explored on the basis of consistent linkage disequilibrium (LD) phase and marker effects in both breeds. RESULTS: The proportion of genetic variance explained by the entire set of SNPs ranged from 37.5 to 57.6 %. Estimated genomic correlations were drastically affected by the process used to select SNPs and went from near 0 to more than 0.80 for most traits when using the set of SNPs with significant effects and the same LD phase in the two breeds. We found that, by carefully selecting the subset of SNPs, the missing heritability can be largely recovered and accuracies in genomic predictions can be improved six-fold. However, the increases in accuracy might come at the expense of large biases. CONCLUSIONS: Our results offer hope for the effective implementation of genomic selection schemes in situations where the number of breeds is large, the sample size within any single breed is small and the breeding objective includes many phenotypes.

Non-synonymous mutations mapped to chromosome X associated with andrological and growth traits in beef cattle.

BACKGROUND: Previous genome-wide association analyses identified QTL regions in the X chromosome for percentage of normal sperm and scrotal circumference in Brahman and Tropical Composite cattle. These traits are important to be studied because they are indicators of male fertility and are correlated with female sexual precocity and reproductive longevity. The aim was to investigate candidate genes in these regions and to identify putative causative mutations that influence these traits. In addition, we tested the identified mutations for female fertility and growth traits. RESULTS: Using a combination of bioinformatics and molecular assay technology, twelve non-synonymous SNPs in eleven genes were genotyped in a cattle population. Three and nine SNPs explained more than 1% of the additive genetic variance for percentage of normal sperm and scrotal circumference, respectively. The SNPs that had a major influence in percentage of normal sperm were mapped to LOC100138021 and TAF7L genes; and in TEX11 and AR genes for scrotal circumference. One SNP in TEX11 was explained ~13% of the additive genetic variance for scrotal circumference at 12 months. The tested SNP were also associated with weight measurements, but not with female fertility traits. CONCLUSIONS: The strong association of SNPs located in X chromosome genes with male fertility traits validates the QTL. The implicated genes became good candidates to be used for genetic evaluation, without detrimentally influencing female fertility traits.

Transcriptome profiles of Penaeus (Marsupenaeus) japonicus animal and vegetal half-embryos: identification of sex determination, germ line, mesoderm, and other developmental genes.

There is virtually no knowledge of the molecular events controlling early embryogenesis in Penaeid shrimp. A combination of controlled spawning environment, shrimp embryo micro-dissection techniques, and next-generation sequencing was used to produce transcriptome EST datasets of Penaeus japonicus animal and vegetal half-embryos. Embryos were collected immediately after spawning, and then blastomeres were separated at the two-cell stage and allowed to develop to late gastrulation, then pooled for RNA isolation and cDNA synthesis. Ion Torrent sequencing of cDNA from approximately 500 pooled animal and vegetal half-embryos from multiple spawnings resulted in 560,516 and 493,703 reads, respectively. Reads from each library were assembled and Gene Ontogeny analysis produced 3479 annotated animal contigs and 4173 annotated vegetal contigs, with 159/139 hits for developmental processes in the animal/vegetal contigs, respectively. Contigs were subject to BLAST for selected developmental toolbox genes. Some of the genes found included the sex determination genes sex-lethal and transformer; the germ line genes argonaute 1, boule, germ cell-less, gustavus, maelstrom, mex-3, par-1, pumilio, SmB, staufen, and tudor; the mesoderm genes brachyury, mef2, snail, and twist; the axis determination/segmentation genes ß-catenin, deformed, distal-less, engrailed, giant, hairy, hunchback, kruppel, orthodenticle, patched, tailless, and wingless/wnt-8c; and a number of cell-cycle regulators. Animal and vegetal contigs were computationally subtracted from each other to produce sets unique to either half-embryo library. Genes expressed only in the animal half included bmp1, kruppel, maelstrom, and orthodenticle. Genes expressed only in the vegetal half included boule, brachyury, deformed, dorsal, engrailed, hunchback, spalt, twist, and wingless/wnt-8c.

Exploiting genomic data to identify proteins involved in abalone reproduction.

Aside from their critical role in reproduction, abalone gonads serve as an indicator of sexual maturity and energy balance, two key considerations for effective abalone culture. Temperate abalone farmers face issues with tank restocking with highly marketable abalone owing to inefficient spawning induction methods. The identification of key proteins in sexually mature abalone will serve as the foundation for a greater understanding of reproductive biology. Addressing this knowledge gap is the first step towards improving abalone aquaculture methods. Proteomic profiling of female and male gonads of greenlip abalone, Haliotis laevigata, was undertaken using liquid chromatography-mass spectrometry. Owing to the incomplete nature of abalone protein databases, in addition to searching against two publicly available databases, a custom database comprising genomic data was used. Overall, 162 and 110 proteins were identified in females and males respectively with 40 proteins common to both sexes. For proteins involved in sexual maturation, sperm and egg structure, motility, acrosomal reaction and fertilization, 23 were identified only in females, 18 only in males and 6 were common. Gene ontology analysis revealed clear differences between the female and male protein profiles reflecting a higher rate of protein synthesis in the ovary and higher metabolic activity in the testis. BIOLOGICAL SIGNIFICANCE: A comprehensive mass spectrometry-based analysis was performed to profile the abalone gonad proteome providing the foundation for future studies of reproduction in abalone. Key proteins involved in both reproduction and energy balance were identified. Genomic resources were utilised to build a database of molluscan proteins yielding >60% more protein identifications than in a standard workflow employing public protein databases.

Variation in genes involved in epigenetic processes offers insights into tropically adapted cattle diversity.

We evaluated the relevance of the BovineHD Illumina SNP chip with respect to genes involved in epigenetic processes. Genotypes for 729,068 SNP on two tropical cattle breeds of Australia were used: Brahman (n = 2112) and Tropical Composite (n = 2550). We used data mining approaches to compile a list of bovine protein-coding genes involved in epigenetic processes. These genes represent 9 functional categories that contain between one (histone demethylases) and 99 (chromatin remodeling factors) genes. A total of 3091 SNP mapped to positions within 3000 bp of the 193 coding regions of those genes, including 113 SNP in transcribed regions, 2738 in intronic regions and 240 in up- or down-stream regions. For all these SNP categories, we observed differences in the allelic frequencies between Brahman and Tropical Composite cattle. These differences were larger than those observed for the entire set of 729,068 SNP (P = 1.79 x 10(-5)). A multidimensional scaling analysis using only the 113 SNP in transcribed regions allowed for the separation of the two populations and this separation was comparable to the one obtained with a random set of 113 SNP (Principal Component 1 r (2) > 0.84). To further characterize the differences between the breeds we defined a gene-differentiation metric based on the average genotypic frequencies of SNP connected to each gene and compared both cattle populations. The 10% most differentiated genes were distributed across 10 chromosomes, with significant (P < 0.05) enrichment on BTA 3 and 10. The 10% most conserved genes were located in 12 chromosomes. We conclude that there is variation between cattle populations in genes connected to epigenetic processes, and this variation can be used to differentiate cattle breeds. More research is needed to fully characterize the use of these SNP and its potential as means to further our understanding of biological variation and epigenetic processes.

Using regulatory and epistatic networks to extend the findings of a genome scan: identifying the gene drivers of pigmentation in merino sheep.

Extending genome wide association analysis by the inclusion of gene expression data may assist in the dissection of complex traits. We examined piebald, a pigmentation phenotype in both human and Merino sheep, by analysing multiple data types using a systems approach. First, a case control analysis of 49,034 ovine SNP was performed which confirmed a multigenic basis for the condition. We combined these results with gene expression data from five tissue types analysed with a skin-specific microarray. Promoter sequence analysis of differentially expressed genes allowed us to reverse-engineer a regulatory network. Likewise, by testing two-loci models derived from all pair-wise comparisons across piebald-associated SNP, we generated an epistatic network. At the intersection of both networks, we identified thirteen genes with insulin-like growth factor binding protein 7 (IGFBP7), platelet-derived growth factor alpha (PDGFRA) and the tetraspanin platelet activator CD9 at the kernel of the intersection. Further, we report a number of differentially expressed genes in regions containing highly associated SNP including ATRN, DOCK7, FGFR1OP, GLI3, SILV and TBX15. The application of network theory facilitated co-analysis of genetic variation with gene expression, recapitulated aspects of the known molecular biology of skin pigmentation and provided insights into the transcription regulation and epistatic interactions involved in piebald Merino sheep.

A genomics-informed, SNP association study reveals FBLN1 and FABP4 as contributing to resistance to fleece rot in Australian Merino sheep.

BACKGROUND: Fleece rot (FR) and body-strike of Merino sheep by the sheep blowfly Lucilia cuprina are major problems for the Australian wool industry, causing significant losses as a result of increased management costs coupled with reduced wool productivity and quality. In addition to direct effects on fleece quality, fleece rot is a major predisposing factor to blowfly strike on the body of sheep. In order to investigate the genetic drivers of resistance to fleece rot, we constructed a combined ovine-bovine cDNA microarray of almost 12,000 probes including 6,125 skin expressed sequence tags and 5,760 anonymous clones obtained from skin subtracted libraries derived from fleece rot resistant and susceptible animals. This microarray platform was used to profile the gene expression changes between skin samples of six resistant and six susceptible animals taken immediately before, during and after FR induction. Mixed-model equations were employed to normalize the data and 155 genes were found to be differentially expressed (DE). Ten DE genes were selected for validation using real-time PCR on independent skin samples. The genomic regions of a further 5 DE genes were surveyed to identify single nucleotide polymorphisms (SNP) that were genotyped across three populations for their associations with fleece rot resistance. RESULTS: The majority of the DE genes originated from the fleece rot subtracted libraries and over-representing gene ontology terms included defense response to bacterium and epidermis development, indicating a role of these processes in modulating the sheep's response to fleece rot. We focused on genes that contribute to the physical barrier function of skin, including keratins, collagens, fibulin and lipid proteins, to identify SNPs that were associated to fleece rot scores. CONCLUSIONS: We identified FBLN1 (fibulin) and FABP4 (fatty acid binding protein 4) as key factors in sheep's resistance to fleece rot. Validation of these markers in other populations could lead to vital tests for marker assisted selection that will ultimately increase the natural fleece rot resistance of Merino sheep.

Analysis of the complement and molecular evolution of tRNA genes in cow.

BACKGROUND: Detailed information regarding the number and organization of transfer RNA (tRNA) genes at the genome level is becoming readily available with the increase of DNA sequencing of whole genomes. However the identification of functional tRNA genes is challenging for species that have large numbers of repetitive elements containing tRNA derived sequences, such as Bos taurus. Reliable identification and annotation of entire sets of tRNA genes allows the evolution of tRNA genes to be understood on a genomic scale. RESULTS: In this study, we explored the B. taurus genome using bioinformatics and comparative genomics approaches to catalogue and analyze cow tRNA genes. The initial analysis of the cow genome using tRNAscan-SE identified 31,868 putative tRNA genes and 189,183 pseudogenes, where 28,830 of the 31,868 predicted tRNA genes were classified as repetitive elements by the RepeatMasker program. We then used comparative genomics to further discriminate between functional tRNA genes and tRNA-derived sequences for the remaining set of 3,038 putative tRNA genes. For our analysis, we used the human, chimpanzee, mouse, rat, horse, dog, chicken and fugu genomes to predict that the number of active tRNA genes in cow lies in the vicinity of 439. Of this set, 150 tRNA genes were 100% identical in their sequences across all nine vertebrate genomes studied. Using clustering analyses, we identified a new tRNA-GlyCCC subfamily present in all analyzed mammalian genomes. We suggest that this subfamily originated from an ancestral tRNA-GlyGCC gene via a point mutation prior to the radiation of the mammalian lineages. Lastly, in a separate analysis we created phylogenetic profiles for each putative cow tRNA gene using a representative set of genomes to gain an overview of common evolutionary histories of tRNA genes. CONCLUSION: The use of a combination of bioinformatics and comparative genomics approaches has allowed the confident identification of a set of cow tRNA genes that will facilitate further studies in understanding the molecular evolution of cow tRNA genes.

A gene coexpression network for bovine skeletal muscle inferred from microarray data.

We present the application of large-scale multivariate mixed-model equations to the joint analysis of nine gene expression experiments in beef cattle muscle and fat tissues with a total of 147 hybridizations, and we explore 47 experimental conditions or treatments. Using a correlation-based method, we constructed a gene network for 822 genes. Modules of muscle structural proteins and enzymes, extracellular matrix, fat metabolism, and protein synthesis were clearly evident. Detailed analysis of the network identified groupings of proteins on the basis of physical association. For example, expression of three components of the z-disk, MYOZ1, TCAP, and PDLIM3, was significantly correlated. In contrast, expression of these z-disk proteins was not highly correlated with the expression of a cluster of thick (myosins) and thin (actin and tropomyosins) filament proteins or of titin, the third major filament system. However, expression of titin was itself not significantly correlated with the cluster of thick and thin filament proteins and enzymes. Correlation in expression of many fast-twitch muscle structural proteins and enzymes was observed, but slow-twitch-specific proteins were not correlated with the fast-twitch proteins or with each other. In addition, a number of significant associations between genes and transcription factors were also identified. Our results not only recapitulate the known biology of muscle but have also started to reveal some of the underlying associations between and within the structural components of skeletal muscle.

Gene expression profiling of bovine in vitro adipogenesis using a cDNA microarray.

The gene expression profile of bovine bone marrow stromal cells undergoing adipogenesis was established using a custom cDNA microarray. Cells that were treated with adipogenic stimulants and those that were not were collected at each of the six time points, and gene expression differences between the treated and untreated samples within each time point were compared using a microarray. Statistical analyses revealed that 158 genes showed a minimum fold change of 2 in at least one of the five post-differentiation time points. These genes are involved in various cellular pathways and functions, including lipogenesis, glycolysis, cytoskeleton remodelling, extracellular matrix, transcription as well as various signalling pathways such as insulin, calcium and wingless signalling. The experiment also identified 17 differentially expressed (DE) microarray elements with no assigned function. Quantitative real-time PCR was employed to validate eight DE genes, and the PCR data were found to reproduce the microarray data for these eight genes. Subsequent gene ontology annotation was able to provide a global overview of the molecular function of DE genes during adipogenesis. This analysis was able to indicate the importance of different gene categories at various stages of adipogenic conversion, thereby providing further insights into the molecular changes during bovine adipogenesis.

Transcriptional profiling of skeletal muscle tissue from two breeds of cattle.

We used a 9.6 K cattle muscle/fat cDNA microarray to study gene expression differences between the longuissimus dorsi (LD) muscle of Japanese Black (JB) and Holstein (HOL) cattle. JB cattle exhibit an unusual ability to accumulate intramuscular adipose tissue with fat melting points lower than that in other breeds. The LD biopsies from three JB (Tajima strain) and three HOL animals were used in this breed comparison. Seventeen genes were identified as preferentially expressed in LD samples from JB and seven genes were found to be expressed more highly in HOL. The expression of six selected differentially expressed genes was confirmed by quantitative real-time PCR. The genes more highly expressed in JB are associated with unsaturated fatty acid synthesis, fat deposition, and the thyroid hormone pathway. These results are consistent with the increased amounts and proportions of monounsaturated fatty acids observed in the muscle of JB animals. By discovering as yet uncharacterized genes that are differentially regulated in this comparison, the work may lead us to a better understanding of the regulatory pathways involved in the development of intramuscular adipose tissue.

An interactive bovine in silico SNP database (IBISS).

An interactive bovine in silico SNP (IBISS) database has been created through the clustering and aligning of bovine EST and mRNA sequences. Approximately 324,000 EST and mRNA sequences were clustered to produce 29,965 clusters (producing 48,679 consensus sequences) and 48,565 singletons. A SNP screening regime was placed on variations detected in the multiple sequence alignment files to determine which SNPs are more likely to be real rather than sequencing errors. A small subset of predicted SNPs was validated on a diverse set of bovine DNA samples using PCR amplification and sequencing. Fifty percent of the predicted SNPs in the "putative >1" category were polymorphic in the population sampled. The IBISS database represents more than just a SNP database; it is also a genomic database containing uniformly annotated predicted gene mRNA and protein sequences, gene structure, and genomic organization information.