Bovine respiratory disease (BRD) is a multifactorial disease complex and the leading infectious disease in post-weaned beef cattle. Clinical manifestations of BRD are recognized in beef calves within a high-risk setting, commonly associated with weaning, shipping, and novel feeding and housing environments. However, the understanding of complex host immune interactions and genomic mechanisms involved in BRD susceptibility remain elusive. Utilizing high-throughput RNA-sequencing, we contrasted the at-arrival blood transcriptomes of 6 beef cattle that ultimately developed BRD against 5 beef cattle that remained healthy within the same herd, differentiating BRD diagnosis from production metadata and treatment records. We identified 135 differentially expressed genes (DEGs) using the differential gene expression tools edgeR and DESeq2. Thirty-six of the DEGs shared between these two analysis platforms were prioritized for investigation of their relevance to infectious disease resistance using WebGestalt, STRING, and Reactome. Biological processes related to inflammatory response, immunological defense, lipoxin metabolism, and macrophage function were identified. Production of specialized pro-resolvin mediators (SPMs) and endogenous metabolism of angiotensinogen were increased in animals that resisted BRD. Protein-protein interaction modeling of gene products with significantly higher expression in cattle that naturally acquire BRD identified molecular processes involving microbial killing. Accordingly, identification of DEGs in whole blood at arrival revealed a clear distinction between calves that went on to develop BRD and those that resisted BRD. These results provide novel insight into host immune factors that are present at the time of arrival that confer protection from BRD.

BackgroundMicroRNAs are small non-coding RNAs that regulate gene expression and thus play important roles in mammalian development. However, the comprehensive lists of microRNAs, as well as, molecular mechanisms by which microRNAs regulate gene expression during gamete and embryo development are poorly defined. The objectives of this study were to determine microRNAs in bull sperm and predict their functions.MethodsTo accomplish our objectives we isolated miRNAs from sperm of high and low fertility bulls, conducted microRNA microarray experiments and validated expression of a panel of microRNAs using real time RT-PCR. Bioinformatic approaches were carried out to identify regulated targets.ResultsWe demonstrated that an abundance of microRNAs were present in bovine spermatozoa, however, only seven were differentially expressed; hsa-aga-3155, -8197, -6727, -11796, -14189, -6125, -13659. The abundance of miRNAs in the spermatozoa and the differential expression in sperm from high vs. low fertility bulls suggests that the miRNAs possibly play important functions in the regulating mechanisms of bovine spermatozoa.ConclusionIdentification of specific microRNAs expressed in spermatozoa of bulls with different fertility phenotypes will help better understand mammalian gametogenesis and early development.

BackgroundAflatoxin is a potent carcinogen that can contaminate grain infected with the fungus Aspergillus flavus. However, resistance to aflatoxin accumulation in maize is a complex trait with low heritability. Here, two complementary analyses were performed to better understand the mechanisms involved. The first coupled results of a genome-wide association study (GWAS) that accounted for linkage disequilibrium among single nucleotide polymorphisms (SNPs) with gene-set enrichment for a pathway-based approach. The rationale was that the cumulative effects of genes in a pathway would give insight into genetic differences that distinguish resistant from susceptible lines of maize. The second involved finding non-pathway genes close to the most significant SNP-trait associations with the greatest effect on reducing aflatoxin in multiple environments. Unlike conventional GWAS, the latter analysis emphasized multiple aspects of SNP-trait associations rather than just significance and was performed because of the high genotype x environment variability exhibited by this trait.ResultsThe most significant metabolic pathway identified was jasmonic acid (JA) biosynthesis. Specifically, there was at least one allelic variant for each step in the JA biosynthesis pathway that conferred an incremental decrease to the level of aflatoxin observed among the inbred lines in the GWAS panel. Several non-pathway genes were also consistently associated with lowered aflatoxin levels. Those with predicted functions related to defense were: leucine-rich repeat protein kinase, expansin B3, reversion-to-ethylene sensitivity1, adaptor protein complex2, and a multidrug and toxic compound extrusion protein.ConclusionsOur genetic analysis provided strong evidence for several genes that were associated with aflatoxin resistance. Inbred lines that exhibited lower levels of aflatoxin accumulation tended to share similar haplotypes for genes specifically in the pathway of JA biosynthesis, along with several non-pathway genes with putative defense-related functions. Knowledge gained from these two complementary analyses has improved our understanding of population differences in aflatoxin resistance.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-1874-9) contains supplementary material, which is available to authorized users.

High throughput biological data need to be processed, analyzed, and interpreted to address problems in life sciences. Bioinformatics, computational biology, and systems biology deal with biological problems using computational methods. Clustering is one of the methods used to gain insight into biological processes, particularly at the genomics level. Clearly, clustering can be used in many areas of biological data analysis. However, this paper presents a review of the current clustering algorithms designed especially for analyzing gene expression data. It is also intended to introduce one of the main problems in bioinformatics - clustering gene expression data - to the operations research community.

The liver is the primary site for the metabolism of nutrients, drugs, and chemical agents. Although metabolic pathways are complex and tightly regulated, genetic variation among individuals, reflected in variations in gene expression levels, introduces complexity into research on liver disease. This study dissected genetic networks that control liver gene expression through the combination of large-scale quantitative mRNA expression analysis with genetic mapping in a reference population of BXD recombinant inbred mouse strains for which extensive singlenucleotide polymorphism, haplotype, and phenotypic data are publicly available. We profiled gene expression in livers of naive mice of both sexes from C57BL/6J, DBA/2J, B6D2F1, and 37 BXD strains using Agilent oligonucleotide microarrays. These data were used to map quantitative trait loci (QTLs) responsible for variations in the expression of about 19,000 transcripts. We identified polymorphic local and distant QTLs, including several loci that control the expression of large numbers of genes in liver, by comparing the physical transcript position with the location of the controlling QTL. Conclusion: The data are available through a public web-based resource (www.genenetwork.org) that allows custom data mining, identification of coregulated transcripts and correlated phenotypes, cross-tissue, and cross-species comparisons, as well as testing of a broad array of hypotheses. (HEPATOLOGY 2007;46:548-557.) T he maturation of gene expression technology has opened the door to the exploration of the genetics of gene expression. 1-3 Microarrays allow for the concurrent measurement of thousands of transcripts, with the resultant genomic data being increasingly used to improve the biological interpretation of data from mechanistic research. Phenotypic anchoring of observed phenotypes to gene expression changes has proven useful in uncovering the molecular mechanisms that lead to liver injury. 4,5 Such experiments connect the variation in the transcript expression to phenotypes. However, they do not lead to detailed gene expression networks in which the expression of 1 gene is found to control the expression of another.Recombinant inbred (RI) mice are created by the crossing of 2 parental strains followed by sib-mating for over 20 generations. 6 Strains created in this way have the advantage of being homozygous at almost every location along the genome. Each representative of an RI line will have limited phenotypic variation within that line, but the variation between lines is usually vast. RI panels are widely used to determine genotype-phenotype associations with quantitative trait locus (QTL) mapping techniques. 7 The relationships between the phenotypes and genotypes are calculated with a likelihood ratio statistic (LRS), which is a measure of the probability that a given genetic marker explains the variation in the phenotype. When mRNA levels are used as the phenotype, regions of the genome with a high LRS are likely to contain genes that control the expression of t...