71 CIRCULATING microRNAs AS POTENTIAL BIOMARKERS OF EARLY PREGNANCY IN HIGH-PRODUCING DAIRY COWS

Lactation induces changes in the metabolic status of postpartum dairy cows that negatively affects the likelihood of pregnancy establishment. At present, pregnancy diagnosis with confidence is only possible after the third week of gestation. Earlier diagnosis could facilitate earlier re-breeding, reduce calving intervals, and improve profits for the industry. MicroRNAs (miRNAs) released in body fluids have been identified as minimally invasive biomarkers of several diseases. In addition, distinct miRNA profiles have been directly related to specific stages of human pregnancy. The aim of this study was to profile circulating miRNAs in the blood of high-producing dairy cows in order to identify biomarkers of early pregnancy. In-calf primiparous Holstein-Friesian cows (n = 22) with a similar economic breeding index were used. At calving, cows were randomly assigned to one of two groups: (1) lactating (n = 11; milked twice per day) or (2) non-lactating (n = 11; dried off immediately). Around 65 to 75 days postpartum (dpp), oestrous cycles were synchronized and a single embryo from a superovulated Holstein-Friesian donor was transferred at Day 7 post-oestrus. Plasma samples were analysed at Day 13 (initiation of conceptus elongation) and at Day 19 (initiation of implantation). Pregnancy rate, established by the presence of conceptus at Day 19, was 5/11 (45%) for lactating and 8/11 (73%) for non-lactating cows, respectively. Circulating miRNA levels were profiled in 4 animals per group in non-lactating pregnant, and lactating pregnant and nonpregnant cows at selected timepoints using Illumina HiSEqn 2000 (Illumina Inc., San Diego, CA, USA) for smallRNA sequencing. Annotation and discovery of miRNAs were done using MirDeep2, and read counts were analysed using edgeR to identify differentially expressed miRNAs. Target genes analysis was run with miRWalk and pathways interactions were built using Cytoscape (P ≤ 0.05). Differentially abundant miRNAs between lactating and non-lactating cows were found at both time points (FDR ≤ 0.05). At Day 13, non-lactating cows had a distinct miRNA profile compared with lactating cows showing 8 differentially expressed miRNA (6 v. pregnant and 2 v. nonpregnant cows). At Day 19, no significant differences were found within pregnant cows, but 5 differentially expressed miRNAs were identified between pregnant and nonpregnant cows, regardless of metabolic status. Interestingly, one miRNA, bta-mir140, was up-regulated in non-lactating pregnant cows from Day 13 onwards compared with nonpregnant cows. Furthermore, the same miRNA was up-regulated in lactating pregnant v. nonpregnant at Day 19. Among bta-mir140 target genes, CD274, SLC44A4, CXCL12, and SIRPA were strictly associated with immune tolerance. In conclusion, the maternal plasma miRNome may represent an early indicator of pregnancy status. In particular, the up-regulation of bta-mir140 in pregnant cows suggests that this miRNA may be a good candidate as an early biomarker of fertility. Furthermore, the positive correlation between this miRNA and pathways involved in T-cell response may indicate a role of immune tolerance in preventing rejection at the initiation of implantation.

Expression of genes involved in the T cell signalling pathway in circulating immune cells of cattle 24 months following oral challenge with Bovine Amyloidotic Spongiform Encephalopathy (BASE)

Background

Bovine Amyloidotic Spongiform Encephalopathy (BASE) is a variant of classical BSE that affects cows and can be transmitted to primates and mice. BASE is biochemically different from BSE and shares some molecular and histo-pathological features with the MV2 sub-type of human sporadic Creutzfeld Jakob Disease (sCJD).

Results

The present work examined the effects of BASE on gene expression in circulating immune cells. Ontology analysis of genes differentially expressed between cattle orally challenged with brain homogenate from cattle following intracranial inoculation with BASE and control cattle identified three main pathways which were affected. Within the immune function pathway, the most affected genes were related to the T cell receptor-mediated T cell activation pathways. The differential expression of these genes in BASE challenged animals at 10,12 and 24 months following challenge, vs unchallenged controls, was investigated by real time PCR.

Conclusions

The results of this study show that the effects of prion diseases are not limited to the CNS, but involve the immune system and particularly T cell signalling during the early stage following challenge, before the appearance of clinical signs.

SNPchiMp v.3: integrating and standardizing single nucleotide polymorphism data for livestock species

BACKGROUND

In recent years, the use of genomic information in livestock species for genetic improvement, association studies and many other fields has become routine. In order to accommodate different market requirements in terms of genotyping cost, manufacturers of single nucleotide polymorphism (SNP) arrays, private companies and international consortia have developed a large number of arrays with different content and different SNP density. The number of currently available SNP arrays differs among species: ranging from one for goats to more than ten for cattle, and the number of arrays available is increasing rapidly. However, there is limited or no effort to standardize and integrate array- specific (e.g. SNP IDs, allele coding) and species-specific (i.e. past and current assemblies) SNP information.

RESULTS

Here we present SNPchiMp v.3, a solution to these issues for the six major livestock species (cow, pig, horse, sheep, goat and chicken). Original data was collected directly from SNP array producers and specific international genome consortia, and stored in a MySQL database. The database was then linked to an open-access web tool and to public databases. SNPchiMp v.3 ensures fast access to the database (retrieving within/across SNP array data) and the possibility of annotating SNP array data in a user-friendly fashion.

CONCLUSIONS

This platform allows easy integration and standardization, and it is aimed at both industry and research. It also enables users to easily link the information available from the array producer with data in public databases, without the need of additional bioinformatics tools or pipelines. In recognition of the open-access use of Ensembl resources, SNPchiMp v.3 was officially credited as an Ensembl E!mpowered tool. Availability at http://bioinformatics.tecnoparco.org/SNPchimp.

Community-driven development for computational biology at Sprints, Hackathons and Codefests

Background

Computational biology comprises a wide range of technologies and approaches. Multiple technologies can be combined to create more powerful workflows if the individuals contributing the data or providing tools for its interpretation can find mutual understanding and consensus. Much conversation and joint investigation are required in order to identify and implement the best approaches.

Traditionally, scientific conferences feature talks presenting novel technologies or insights, followed up by informal discussions during coffee breaks. In multi-institution collaborations, in order to reach agreement on implementation details or to transfer deeper insights in a technology and practical skills, a representative of one group typically visits the other. However, this does not scale well when the number of technologies or research groups is large.

Conferences have responded to this issue by introducing Birds-of-a-Feather (BoF) sessions, which offer an opportunity for individuals with common interests to intensify their interaction. However, parallel BoF sessions often make it hard for participants to join multiple BoFs and find common ground between the different technologies, and BoFs are generally too short to allow time for participants to program together.

Results

This report summarises our experience with computational biology Codefests, Hackathons and Sprints, which are interactive developer meetings. They are structured to reduce the limitations of traditional scientific meetings described above by strengthening the interaction among peers and letting the participants determine the schedule and topics. These meetings are commonly run as loosely scheduled "unconferences" (self-organized identification of participants and topics for meetings) over at least two days, with early introductory talks to welcome and organize contributors, followed by intensive collaborative coding sessions. We summarise some prominent achievements of those meetings and describe differences in how these are organised, how their audience is addressed, and their outreach to their respective communities.

Conclusions

Hackathons, Codefests and Sprints share a stimulating atmosphere that encourages participants to jointly brainstorm and tackle problems of shared interest in a self-driven proactive environment, as well as providing an opportunity for new participants to get involved in collaborative projects.

A unique mutation in a MYB gene cosegregates with the nectarine phenotype in peach

Nectarines play a key role in peach industry; the fuzzless skin has implications for consumer acceptance. The peach/nectarine (G/g) trait was described as monogenic and previously mapped on chromosome 5. Here, the position of the G locus was delimited within a 1.1 cM interval (635 kb) based on linkage analysis of an F₂ progeny from the cross ‘Contender’ (C, peach) x ‘Ambra’ (A, nectarine). Careful inspection of the genes annotated in the corresponding genomic sequence (Peach v1.0), coupled with variant discovery, led to the identification of MYB gene PpeMYB25 as a candidate for trichome formation on fruit skin. Analysis of genomic re-sequencing data from five peach/nectarine accessions pointed to the insertion of a LTR retroelement in exon 3 of the PpeMYB25 gene as the cause of the recessive glabrous phenotype. A functional marker (indelG) developed on the LTR insertion cosegregated with the trait in the CxA F₂ progeny and was validated on a broad panel of genotypes, including all known putative donors of the nectarine trait. This marker was shown to efficiently discriminate between peach and nectarine plants, indicating that a unique mutational event gave rise to the nectarine trait and providing a useful diagnostic tool for early seedling selection in peach breeding programs.

SNPchiMp: a database to disentangle the SNPchip jungle in bovine livestock

BACKGROUND

Currently, six commercial whole-genome SNP chips are available for cattle genotyping, produced by two different genotyping platforms. Technical issues need to be addressed to combine data that originates from the different platforms, or different versions of the same array generated by the manufacturer. For example: i) genome coordinates for SNPs may refer to different genome assemblies; ii) reference genome sequences are updated over time changing the positions, or even removing sequences which contain SNPs; iii) not all commercial SNP ID's are searchable within public databases; iv) SNPs can be coded using different formats and referencing different strands (e.g. A/B or A/C/T/G alleles, referencing forward/reverse, top/bottom or plus/minus strand); v) Due to new information being discovered, higher density chips do not necessarily include all the SNPs present in the lower density chips; and, vi) SNP IDs may not be consistent across chips and platforms. Most researchers and breed associations manage SNP data in real-time and thus require tools to standardise data in a user-friendly manner.

DESCRIPTION

Here we present SNPchiMp, a MySQL database linked to an open access web-based interface. Features of this interface include, but are not limited to, the following functions: 1) referencing the SNP mapping information to the latest genome assembly, 2) extraction of information contained in dbSNP for SNPs present in all commercially available bovine chips, and 3) identification of SNPs in common between two or more bovine chips (e.g. for SNP imputation from lower to higher density). In addition, SNPchiMp can retrieve this information on subsets of SNPs, accessing such data either via physical position on a supported assembly, or by a list of SNP IDs, rs or ss identifiers.

CONCLUSIONS

This tool combines many different sources of information, that otherwise are time consuming to obtain and difficult to integrate. The SNPchiMp not only provides the information in a user-friendly format, but also enables researchers to perform a large number of operations with a few clicks of the mouse. This significantly reduces the time needed to execute the large number of operations required to manage SNP data.

BioHackathon series in 2011 and 2012: penetration of ontology and linked data in life science domains

The application of semantic technologies to the integration of biological data and the interoperability of bioinformatics analysis and visualization tools has been the common theme of a series of annual BioHackathons hosted in Japan for the past five years. Here we provide a review of the activities and outcomes from the BioHackathons held in 2011 in Kyoto and 2012 in Toyama. In order to efficiently implement semantic technologies in the life sciences, participants formed various sub-groups and worked on the following topics: Resource Description Framework (RDF) models for specific domains, text mining of the literature, ontology development, essential metadata for biological databases, platforms to enable efficient Semantic Web technology development and interoperability, and the development of applications for Semantic Web data. In this review, we briefly introduce the themes covered by these sub-groups. The observations made, conclusions drawn, and software development projects that emerged from these activities are discussed.

Transcript profiling of common bean nodules subjected to oxidative stress

Several environmental stresses generate high amounts of reactive oxygen species (ROS) in plant cells, resulting in oxidative stress. Symbiotic nitrogen fixation (SNF) in the legume-rhizobia symbiosis is sensitive to damage from oxidative stress. Active nodules of the common bean (Phaseolus vulgaris) exposed to the herbicide paraquat (1,1'-dimethyl-4,4'-bipyridinium dichloride hydrate), which stimulates ROS accumulation, exhibited reduced nitrogenase activity and ureide content. We analyzed the global gene response of nodules subjected to oxidative stress using the Bean Custom Array 90K, which includes probes from 30,000 expressed sequence tags (ESTs). A total of 4280 ESTs were differentially expressed in stressed bean nodules; of these, 2218 were repressed. Based on Gene Ontology analysis, these genes were grouped into 42 different biological process categories. Analysis with the PathExpress bioinformatic tool, adapted for bean, identified five significantly repressed metabolic pathways related to carbon/nitrogen metabolism, which is crucial for nodule function. Quantitative reverse transcription (qRT)-PCR analysis of transcription factor (TF) gene expression showed that 67 TF genes were differentially expressed in nodules exposed to oxidative stress. Putative cis-elements recognized by highly responsive TF were detected in promoter regions of oxidative stress regulated genes. The expression of oxidative stress responsive genes and of genes important for SNF in bacteroids analyzed in stressed nodules revealed that these conditions elicited a transcriptional response.

De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis

De novo assembly of RNA-seq data enables researchers to study transcriptomes without the need for a genome sequence; this approach can be usefully applied, for instance, in research on 'non-model organisms' of ecological and evolutionary importance, cancer samples or the microbiome. In this protocol we describe the use of the Trinity platform for de novo transcriptome assembly from RNA-seq data in non-model organisms. We also present Trinity-supported companion utilities for downstream applications, including RSEM for transcript abundance estimation, R/Bioconductor packages for identifying differentially expressed transcripts across samples and approaches to identify protein-coding genes. In the procedure, we provide a workflow for genome-independent transcriptome analysis leveraging the Trinity platform. The software, documentation and demonstrations are freely available from http://trinityrnaseq.sourceforge.net. The run time of this protocol is highly dependent on the size and complexity of data to be analyzed. The example data set analyzed in the procedure detailed herein can be processed in less than 5 h.

De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis

De novo assembly of RNA-seq data enables researchers to study transcriptomes without the need for a genome sequence; this approach can be usefully applied, for instance, in research on 'non-model organisms' of ecological and evolutionary importance, cancer samples or the microbiome. In this protocol we describe the use of the Trinity platform for de novo transcriptome assembly from RNA-seq data in non-model organisms. We also present Trinity-supported companion utilities for downstream applications, including RSEM for transcript abundance estimation, R/Bioconductor packages for identifying differentially expressed transcripts across samples and approaches to identify protein-coding genes. In the procedure, we provide a workflow for genome-independent transcriptome analysis leveraging the Trinity platform. The software, documentation and demonstrations are freely available from http://trinityrnaseq.sourceforge.net. The run time of this protocol is highly dependent on the size and complexity of data to be analyzed. The example data set analyzed in the procedure detailed herein can be processed in less than 5 h.

De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis

De novo assembly of RNA-seq data enables researchers to study transcriptomes without the need for a genome sequence; this approach can be usefully applied, for instance, in research on 'non-model organisms' of ecological and evolutionary importance, cancer samples or the microbiome. In this protocol we describe the use of the Trinity platform for de novo transcriptome assembly from RNA-seq data in non-model organisms. We also present Trinity-supported companion utilities for downstream applications, including RSEM for transcript abundance estimation, R/Bioconductor packages for identifying differentially expressed transcripts across samples and approaches to identify protein-coding genes. In the procedure, we provide a workflow for genome-independent transcriptome analysis leveraging the Trinity platform. The software, documentation and demonstrations are freely available from http://trinityrnaseq.sourceforge.net. The run time of this protocol is highly dependent on the size and complexity of data to be analyzed. The example data set analyzed in the procedure detailed herein can be processed in less than 5 h.

Common bean (Phaseolus vulgaris L.) PvTIFY orchestrates global changes in transcript profile response to jasmonate and phosphorus deficiency

BACKGROUND

TIFY is a large plant-specific transcription factor gene family. A subgroup of TIFY genes named JAZ (Jasmonate-ZIM domain) has been identified as repressors of jasmonate (JA)-regulated transcription in Arabidopsis and other plants. JA signaling is involved in many aspects of plant growth/development and in defense responses to biotic and abiotic stresses. Here, we identified the TIFY genes (designated PvTIFY) from the legume common bean (Phaseolus vulgaris) and functionally characterized PvTIFY10C as a transcriptional regulator.

RESULTS

Nineteen genes from the PvTIFY gene family were identified through whole-genome sequence analysis. Most of these were induced upon methyl-JA elicitation. We selected PvTIFY10C as a representative JA-responsive PvTIFY gene for further functional analysis. Transcriptome analysis via microarray hybridization using the newly designed Bean Custom Array 90 K was performed on transgenic roots of composite plants with modulated (RNAi-silencing or over-expression) PvTIFY10C gene expression. Data were interpreted using Gene Ontology and MapMan adapted to common bean. Microarray differential gene expression data were validated by real-time qRT-PCR expression analysis. Comparative global gene expression analysis revealed opposite regulatory changes in processes such as RNA and protein regulation, stress responses and metabolism in PvTIFY10C silenced vs. over-expressing roots. These data point to transcript reprogramming (mainly repression) orchestrated by PvTIFY10C. In addition, we found that several PvTIFY genes, as well as genes from the JA biosynthetic pathway, responded to P-deficiency. Relevant P-responsive genes that participate in carbon metabolic pathways, cell wall synthesis, lipid metabolism, transport, DNA, RNA and protein regulation, and signaling were oppositely-regulated in control vs. PvTIFY10C-silenced roots of composite plants under P-stress. These data indicate that PvTIFY10C regulates, directly or indirectly, the expression of some P-responsive genes; this process could be mediated by JA-signaling.

CONCLUSION

Our work contributes to the functional characterization of PvTIFY transcriptional regulators in common bean, an agronomically important legume. Members from the large PvTIFY gene family are important global transcriptional regulators that could participate as repressors in the JA signaling pathway. In addition, we propose that the JA-signaling pathway involving PvTIFY genes might play a role in regulating the plant response/adaptation to P-starvation.

Response of the goat mammary gland to infection with Staphylococcus aureus revealed by gene expression profiling in milk somatic and white blood cells

Background

S. aureus is one of the main pathogens responsible for the intra-mammary infection in dairy ruminants. Although much work has been carried out to understand the complex physiological and cellular events that occur in the mammary gland in response to S. aureus, the protective mechanisms are still poorly understood. The objectives of the present study were to investigate gene expression during the early response of the goat mammary gland to an experimental challenge with S. aureus, in order to better understand the local and systemic response and to compare them in two divergent lines of goat selected for high and low milk somatic cell scores.

Results

No differences in gene expression were found between high and low SCS (Somatic Cells Score) selection lines. Analysing the two groups together, an expression of 300 genes were found to change from T0 before infection, and T4 at 24 hours and T5 at 30 hours following challenge. In blood derived white blood cells 8 genes showed increased expression between T0 and T5 and 1 gene has reduced expression. The genes showing the greatest increase in expression following challenge (5.65 to 3.16 fold change) play an important role in (i) immune and inflammatory response (NFKB1, TNFAIP6, BASP1, IRF1, PLEK, BATF3); (ii) the regulation of innate resistance to pathogens (PTX3); and (iii) the regulation of cell metabolism (CYTH4, SLC2A6, ARG2). The genes with reduced expression (−1.5 to −2.5 fold) included genes involved in (i) lipid metabolism (ABCG2, FASN), (ii) chemokine, cytokine and intracellular signalling (SPPI), and (iii) cell cytoskeleton and extracellular matrix (KRT19).

Conclusions

Analysis of genes with differential expression following infection showed an inverse relationship between immune response and lipid metabolism in the early response of the mammary gland to the S. aureus challenge. PTX3 showed a large change in expression in both milk and blood, and is therefore a candidate for further studies on immune response associated with mastitis.

Meta-analysis of two genome-wide association studies of bovine paratuberculosis

Background

Bovine paratuberculosis (ParaTB) also known as Johne's disease, is a contagious fatal disease resulting from infection by Mycobacterium avium subspecies paratuberculosis (MAP). Previous studies have identified loci associated with ParaTB using different measurements to define cases and controls. The objective of this study was to combine the data from two recent studies to identify genetic loci associated with MAP tissue infection and humoral immune response, defined by MAP ELISA-positive cattle, by comparing cases and control animals for one or both measures of infection.

Methodology/Principal Findings

The two populations used for the association analyses were a cohort of MAP tissue infected animals and control Holstein cows from the USA and the second cohort composed of ELISA-positive and ELISA-negative Holstein cows from Italy. Altogether 1190 cattle were genotyped with the Illumina BovineSNP50 BeadChip. SNP markers were removed if the minor allele frequency <0.01 or genotyping failure was >5%. Animals were removed with >5% genotyping failure. Whole genome association analyses were conducted with the GRAMMAR-CG method using two different definitions of control populations.

Conclusion/Significance

The analyses identified several loci (P<5 e-05) associated with ParaTB, defined by positive ELISA and presence of bacteria in tissue compared to ELISA and tissue negative animals, on chromosomes 1, 12 and 15 and one unassigned SNP. These results confirmed associations on chromosome 12 and the unassigned SNP with ParaTB which had been found in the Italian population alone. Furthermore, several additional genomic regions were found associated with ParaTB when ELISA and tissue positive animals were compared with tissue negative samples. These loci were on chromosomes 1, 6, 7, 13, 16, 21,23 and 25 (P<5 e-05). The results clearly indicate the importance of the phenotype definition when seeking to identify markers associated with different disease responses.

Biogem: an effective tool-based approach for scaling up open source software development in bioinformatics

SUMMARY

Biogem provides a software development environment for the Ruby programming language, which encourages community-based software development for bioinformatics while lowering the barrier to entry and encouraging best practices. Biogem, with its targeted modular and decentralized approach, software generator, tools and tight web integration, is an improved general model for scaling up collaborative open source software development in bioinformatics.

AVAILABILITY

Biogem and modules are free and are OSS. Biogem runs on all systems that support recent versions of Ruby, including Linux, Mac OS X and Windows. Further information at http://www.biogems.info. A tutorial is available at http://www.biogems.info/howto.html

CONTACT

bonnal@ingm.org.

Genome wide scan for somatic cell counts in holstein bulls

Background

Mastitis is the most costly disease for dairy production, and control of the disease is often difficult, due to its multi-factorial nature. Susceptibility to mastitis is under partial genetic control and the industry uses indirect selection for decreased concentrations of somatic cells in milk to reduce mastitis.

Methods

A genome-wide scan was performed to identify genomic regions associated with deregressed estimated breeding values (EBVs) for somatic cell counts (SCC) in Holstein bulls. In total 1183 proven bulls of the Italian of Holstein population, were genotyped with the BovineSNP50 BeadChip (Illumina, San Diego, CA) and a whole genome association analysis was performed using the R package GenABEL.

Results

Two chromosomal regions showed association with SCC, a region on chromosome 14 with high significance (P < 5x10^-6) and a region on chromosome 6 with moderate significance (P < 5x10^-5).

Conclusions

Two regions with effects on SCC have been identified with good statistical support. A further study of these candidate regions will be performed to verify the results and identify the causal mutations.

A Ruby API to query the Ensembl database for genomic features

The Ensembl database makes genomic features available via its Genome Browser. It is also possible to access the underlying data through a Perl API for advanced querying. We have developed a full-featured Ruby API to the Ensembl databases, providing the same functionality as the Perl interface with additional features. A single Ruby API is used to access different releases of the Ensembl databases and is also able to query multi-species databases.

AVAILABILITY AND IMPLEMENTATION

Most functionality of the API is provided using the ActiveRecord pattern. The library depends on introspection to make it release independent. The API is available through the Rubygem system and can be installed with the command gem install ruby-ensembl-api.

Analysis of gene expression in white blood cells of cattle orally challenged with bovine amyloidotic spongiform encephalopathy

Bovine amyloidotic spongiform encephalopathy (BASE) is one of the recently discovered atypical forms of BSE, which is transmissible to primates, and may be the bovine equivalent of sporadic Creutzfeldt-Jacob disease (CJD) in humans. Although it is transmissible, it is unknown whether BASE is acquired through infection or arises spontaneously. In the present study, the gene expression of white blood cells (WBCs) from 5 cattle at 1 yr after oral BASE challenge was compared with negative controls using a custom microarray containing 43,768 unique gene probes. In total, 56 genes were found to be differentially expressed between BASE and control animals with a log fold change of 2 or greater. Of these, 39 were upregulated in BASE animals, while 17 were downregulated. The majority of these genes are related to immune function. In particular, BASE animals appeared to have significantly modified expression of genes linked to T- and B-cell development and activation, and to inflammatory responses. The potential impacts of these gene expression changes are described.

Genetic loci involved in antibody response to Mycobacterium avium ssp. paratuberculosis in cattle

BACKGROUND

Mycobacterium avium subsp. paratuberculosis (MAP) causes chronic enteritis in a wide range of animal species. In cattle, MAP causes a chronic disease called Johne's disease, or paratuberculosis, that is not treatable and the efficacy of vaccine control is controversial. The clinical phase of the disease is characterised by diarrhoea, weight loss, drop in milk production and eventually death. Susceptibility to MAP infection is heritable with heritability estimates ranging from 0.06 to 0.10. There have been several studies over the last few years that have identified genetic loci putatively associated with MAP susceptibility, however, with the availability of genome-wide high density SNP maker panels it is now possible to carry out association studies that have higher precision.

METHODOLOGY/PRINCIPAL FINDINGS

The objective of the current study was to localize genes having an impact on Johne's disease susceptibility using the latest bovine genome information and a high density SNP panel (Illumina BovineSNP50 BeadChip) to perform a case/control, genome-wide association analysis. Samples from MAP case and negative controls were selected from field samples collected in 2007 and 2008 in the province of Lombardy, Italy. Cases were defined as animals serologically positive for MAP by ELISA. In total 966 samples were genotyped: 483 MAP ELISA positive and 483 ELISA negative. Samples were selected randomly among those collected from 119 farms which had at least one positive animal.

CONCLUSION/SIGNIFICANCE

THE ANALYSIS OF THE GENOTYPE DATA IDENTIFIED SEVERAL CHROMOSOMAL REGIONS ASSOCIATED WITH DISEASE STATUS: a region on chromosome 12 with high significance (P<5x10(-6)), while regions on chromosome 9, 11, and 12 had moderate significance (P<5x10(-5)). These results provide evidence for genetic loci involved in the humoral response to MAP. Knowledge of genetic variations related to susceptibility will facilitate the incorporation of this information into breeding programmes for the improvement of health status.

The other side of comparative genomics: genes with no orthologs between the cow and other mammalian species

Background

With the rapid growth in the availability of genome sequence data, the automated identification of orthologous genes between species (orthologs) is of fundamental importance to facilitate functional annotation and studies on comparative and evolutionary genomics. Genes with no apparent orthologs between the bovine and human genome may be responsible for major differences between the species, however, such genes are often neglected in functional genomics studies.

Results

A BLAST-based method was exploited to explore the current annotation and orthology predictions in Ensembl. Genes with no orthologs between the two genomes were classified into groups based on alignments, ontology, manual curation and publicly available information. Starting from a high quality and specific set of orthology predictions, as provided by Ensembl, hidden relationship between genes and genomes of different mammalian species were unveiled using a highly sensitive approach, based on sequence similarity and genomic comparison.

Conclusions

The analysis identified 3,801 bovine genes with no orthologs in human and 1010 human genes with no orthologs in cow, among which 411 and 43 genes, respectively, had no match at all in the other species. Most of the apparently non-orthologous genes may potentially have orthologs which were missed in the annotation process, despite having a high percentage of identity, because of differences in gene length and structure. The comparative analysis reported here identified gene variants, new genes and species-specific features and gave an overview of the other side of orthology which may help to improve the annotation of the bovine genome and the knowledge of structural differences between species.