StORF-Reporter: finding genes between genes.

Large regions of prokaryotic genomes are currently without any annotation, in part due to well-established limitations of annotation tools. For example, it is routine for genes using alternative start codons to be misreported or completely omitted. Therefore, we present StORF-Reporter, a tool that takes an annotated genome and returns regions that may contain missing CDS genes from unannotated regions. StORF-Reporter consists of two parts. The first begins with the extraction of unannotated regions from an annotated genome. Next, Stop-ORFs (StORFs) are identified in these unannotated regions. StORFs are open reading frames that are delimited by stop codons and thus can capture those genes most often missing in genome annotations. We show this methodology recovers genes missing from canonical genome annotations. We inspect the results of the genomes of model organisms, the pangenome of Escherichia coli, and a set of 5109 prokaryotic genomes of 247 genera from the Ensembl Bacteria database. StORF-Reporter extended the core, soft-core and accessory gene collections, identified novel gene families and extended families into additional genera. The high levels of sequence conservation observed between genera suggest that many of these StORFs are likely to be functional genes that should now be considered for inclusion in canonical annotations.

Improving Orthologous Signal and Model Fit in Datasets Addressing the Root of the Animal Phylogeny.

There is conflicting evidence as to whether Porifera (sponges) or Ctenophora (comb jellies) comprise the root of the animal phylogeny. Support for either a Porifera-sister or Ctenophore-sister tree has been extensively examined in the context of model selection, taxon sampling, and outgroup selection. The influence of dataset construction is comparatively understudied. We re-examine five animal phylogeny datasets that have supported either root hypothesis using an approach designed to enrich orthologous signal in phylogenomic datasets. We find that many component orthogroups in animal datasets fail to recover major lineages as monophyletic with the exception of Ctenophora, regardless of the supported root. Enriching these datasets to retain orthogroups recovering ≥3 major lineages reduces dataset size by up to 50% while retaining underlying phylogenetic information and taxon sampling. Site-heterogeneous phylogenomic analysis of these enriched datasets recovers both Porifera-sister and Ctenophora-sister positions, even with additional constraints on outgroup sampling. Two datasets which previously supported Ctenophora-sister support Porifera-sister upon enrichment. All enriched datasets display improved model fitness under posterior predictive analysis. While not conclusively rooting animals at either Porifera or Ctenophora, we do see an increase in signal for Porifera-sister and a decrease in signal for Ctenophore-sister when data are filtered for orthologous signal. Our results indicate that dataset size and construction as well as model fit influence animal root inference.

No one tool to rule them all: prokaryotic gene prediction tool annotations are highly dependent on the organism of study.

MOTIVATION: The biases in CoDing Sequence (CDS) prediction tools, which have been based on historic genomic annotations from model organisms, impact our understanding of novel genomes and metagenomes. This hinders the discovery of new genomic information as it results in predictions being biased towards existing knowledge. To date, users have lacked a systematic and replicable approach to identify the strengths and weaknesses of any CDS prediction tool and allow them to choose the right tool for their analysis. RESULTS: We present an evaluation framework (ORForise) based on a comprehensive set of 12 primary and 60 secondary metrics that facilitate the assessment of the performance of CDS prediction tools. This makes it possible to identify which performs better for specific use-cases. We use this to assess 15 ab initio- and model-based tools representing those most widely used (historically and currently) to generate the knowledge in genomic databases. We find that the performance of any tool is dependent on the genome being analysed, and no individual tool ranked as the most accurate across all genomes or metrics analysed. Even the top-ranked tools produced conflicting gene collections, which could not be resolved by aggregation. The ORForise evaluation framework provides users with a replicable, data-led approach to make informed tool choices for novel genome annotations and for refining historical annotations. AVAILABILITY AND IMPLEMENTATION: Code and datasets for reproduction and customisation are available at https://github.com/NickJD/ORForise. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

On the complexity of haplotyping a microbial community.

MOTIVATION: Population-level genetic variation enables competitiveness and niche specialization in microbial communities. Despite the difficulty in culturing many microbes from an environment, we can still study these communities by isolating and sequencing DNA directly from an environment (metagenomics). Recovering the genomic sequences of all isoforms of a given gene across all organisms in a metagenomic sample would aid evolutionary and ecological insights into microbial ecosystems with potential benefits for medicine and biotechnology. A significant obstacle to this goal arises from the lack of a computationally tractable solution that can recover these sequences from sequenced read fragments. This poses a problem analogous to reconstructing the two sequences that make up the genome of a diploid organism (i.e. haplotypes) but for an unknown number of individuals and haplotypes. RESULTS: The problem of single individual haplotyping was first formalized by Lancia et al. in 2001. Now, nearly two decades later, we discuss the complexity of 'haplotyping' metagenomic samples, with a new formalization of Lancia et al.'s data structure that allows us to effectively extend the single individual haplotype problem to microbial communities. This work describes and formalizes the problem of recovering genes (and other genomic subsequences) from all individuals within a complex community sample, which we term the metagenomic individual haplotyping problem. We also provide software implementations for a pairwise single nucleotide variant (SNV) co-occurrence matrix and greedy graph traversal algorithm. AVAILABILITY AND IMPLEMENTATION: Our reference implementation of the described pairwise SNV matrix (Hansel) and greedy haplotype path traversal algorithm (Gretel) is open source, MIT licensed and freely available online at github.com/samstudio8/hansel and github.com/samstudio8/gretel, respectively.

The rumen eukaryotome is a source of novel antimicrobial peptides with therapeutic potential.

BACKGROUND: The rise of microbial antibiotic resistance is a leading threat to the health of the human population. As such, finding new approaches to tackle these microbes, including development of novel antibiotics is vital. RESULTS: In this study, we mined a rumen eukaryotic metatranscriptomic library for novel Antimicrobial peptides (AMPs) using computational approaches and thereafter characterised the therapeutic potential of the AMPs. We identified a total of 208 potentially novel AMPs from the ruminal eukaryotome, and characterised one of those, namely Lubelisin. Lubelisin (GIVAWFWRLAR) is an α-helical peptide, 11 amino acid long with theoretical molecular weight of 1373.76 D. In the presence of Lubelisin, strains of methicillin-resistant Staphylococcus aureus (MRSA) USA300 and EMRSA-15 were killed within 30 min of exposure with ≥103 and 104 CFU/mL reduction in viable cells respectively. Cytotoxicity of Lubelisin against both human and sheep erythrocytes was low resulting in a therapeutic index of 0.43. Membrane permeabilisation assays using propidium iodide alongside transmission electron microscopy revealed that cytoplasmic membrane damage may contribute to the antimicrobial activities of Lubelisin. CONCLUSIONS: We demonstrate that the rumen eukaryotome is a viable source for the discovery of antimicrobial molecules for the treatment of bacterial infections and further development of these may provide part of the potential solution to the ongoing problem of antimicrobial resistance. The role of these AMPs in the ecological warfare within the rumen is also currently unknown.

Insights into the skin of caecilian amphibians from gene expression profiles.

BACKGROUND: Gene expression profiles can provide insights into the molecular machinery behind tissue functions and, in turn, can further our understanding of environmental responses, and developmental and evolutionary processes. During vertebrate evolution, the skin has played a crucial role, displaying a wide diversity of essential functions. To unravel the molecular basis of skin specialisations and adaptations, we compared gene expression in the skin with eight other tissues in a phylogenetically and ecologically diverse species sample of one of the most neglected vertebrate groups, the caecilian amphibians (order Gymnophiona). RESULTS: The skin of the five studied caecilian species showed a distinct gene expression profile reflecting its developmental origin and showing similarities to other epithelial tissues. We identified 59 sequences with conserved enhanced expression in the skin that might be associated with caecilian dermal specialisations. Some of the up-regulated genes shared expression patterns with human skin and potentially are involved in skin functions across vertebrates. Variation trends in gene expression were detected between mid and posterior body skin suggesting different functions between body regions. Several candidate biologically active peptides were also annotated. CONCLUSIONS: Our study provides the first atlas of differentially expressed sequences in caecilian tissues and a baseline to explore the molecular basis of the skin functions in caecilian amphibians, and more broadly in vertebrates.

Inadvertent Paralog Inclusion Drives Artifactual Topologies and Timetree Estimates in Phylogenomics.

Increasingly, large phylogenomic data sets include transcriptomic data from nonmodel organisms. This not only has allowed controversial and unexplored evolutionary relationships in the tree of life to be addressed but also increases the risk of inadvertent inclusion of paralogs in the analysis. Although this may be expected to result in decreased phylogenetic support, it is not clear if it could also drive highly supported artifactual relationships. Many groups, including the hyperdiverse Lissamphibia, are especially susceptible to these issues due to ancient gene duplication events and small numbers of sequenced genomes and because transcriptomes are increasingly applied to resolve historically conflicting taxonomic hypotheses. We tested the potential impact of paralog inclusion on the topologies and timetree estimates of the Lissamphibia using published and de novo sequencing data including 18 amphibian species, from which 2,656 single-copy gene families were identified. A novel paralog filtering approach resulted in four differently curated data sets, which were used for phylogenetic reconstructions using Bayesian inference, maximum likelihood, and quartet-based supertrees. We found that paralogs drive strongly supported conflicting hypotheses within the Lissamphibia (Batrachia and Procera) and older divergence time estimates even within groups where no variation in topology was observed. All investigated methods, except Bayesian inference with the CAT-GTR model, were found to be sensitive to paralogs, but with filtering convergence to the same answer (Batrachia) was observed. This is the first large-scale study to address the impact of orthology selection using transcriptomic data and emphasizes the importance of quality over quantity particularly for understanding relationships of poorly sampled taxa.

What lies beneath? Molecular evolution during the radiation of caecilian amphibians.

BACKGROUND: Evolution leaves an imprint in species through genetic change. At the molecular level, evolutionary changes can be explored by studying ratios of nucleotide substitutions. The interplay among molecular evolution, derived phenotypes, and ecological ranges can provide insights into adaptive radiations. Caecilians (order Gymnophiona), probably the least known of the major lineages of vertebrates, are limbless tropical amphibians, with adults of most species burrowing in soils (fossoriality). This enigmatic order of amphibians are very distinct phenotypically from other extant amphibians and likely from the ancestor of Lissamphibia, but little to nothing is known about the molecular changes underpinning their radiation. We hypothesised that colonization of various depths of tropical soils and of freshwater habitats presented new ecological opportunities to caecilians. RESULTS: A total of 8540 candidate groups of orthologous genes from transcriptomic data of five species of caecilian amphibians and the genome of the frog Xenopus tropicalis were analysed in order to investigate the genetic machinery behind caecilian diversification. We found a total of 168 protein-coding genes with signatures of positive selection at different evolutionary times during the radiation of caecilians. The majority of these genes were related to functional elements of the cell membrane and extracellular matrix with expression in several different tissues. The first colonization of the tropical soils was connected to the largest number of protein-coding genes under positive selection in our analysis. From the results of our study, we highlighted molecular changes in genes involved in perception, reduction-oxidation processes, and aging that likely were involved in the adaptation to different soil strata. CONCLUSIONS: The genes inferred to have been under positive selection provide valuable insights into caecilian evolution, potentially underpin adaptations of caecilians to their extreme environments, and contribute to a better understanding of fossorial adaptations and molecular evolution in vertebrates.

Improvement of Feed Efficiency in Pigs through Microbial Modulation via Fecal Microbiota Transplantation in Sows and Dietary Supplementation of Inulin in Offspring.

As previous studies have demonstrated a link between the porcine intestinal microbiome and feed efficiency (FE), microbiota manipulation may offer a means of improving FE in pigs. A fecal microbiota transplantation procedure (FMTp), using fecal extracts from highly feed-efficient pigs, was performed in pregnant sows (n = 11), with a control group (n = 11) receiving no FMTp. At weaning, offspring were allocated, within sow treatment, to (i) control (n = 67; no dietary supplement) or (ii) inulin (n = 65; 6-week dietary inulin supplementation) treatments. The sow FMTp, alone or in combination with inulin supplementation in offspring, reduced offspring body weight by 8.1 to 10.6 kg at ∼140 days of age, but there was no effect on feed intake. It resulted in better FE, greater bacterial diversity, and higher relative abundances of potentially beneficial bacterial taxa (Fibrobacter and Prevotella) in offspring. Due to the FMTp and/or inulin supplementation, relative abundances of potential pathogens (Chlamydia and Treponema) in the ileum and cecal concentrations of butyric acid were significantly lower. The maternal FMTp led to a greater number of jejunal goblet cells in offspring. Inulin supplementation alone did not affect growth or FE but upregulated duodenal genes linked to glucose and volatile fatty acid homeostasis and increased the mean platelet volume but reduced ileal propionic acid concentrations, granulocyte counts, and serum urea concentrations. Overall, the FMTp in pregnant sows, with or without dietary inulin supplementation in offspring, beneficially modulated offspring intestinal microbiota (albeit mostly low-relative-abundance taxa) and associated physiological parameters. Although FE was improved, the detrimental effect on growth limits the application of this FMTp-inulin strategy in commercial pig production.IMPORTANCE As previous research suggests a link between microbiota and FE, modulation of the intestinal microbiome may be effective in improving FE in pigs. The FMTp in gestating sows, alone or in combination with postweaning dietary inulin supplementation in offspring, achieved improvements in FE and resulted in a higher relative abundance of intestinal bacteria associated with fiber degradation and a lower relative abundance of potential pathogens. However, there was a detrimental effect on growth, although this may not be wholly attributable to microbiota transplantation, as antibiotic and other interventions were also part of the FMT regimen. Therefore, further work with additional control groups is needed to disentangle the effects of each component of the FMTp in order to develop a regimen with practical applications in pig production. Additional research based on findings from this study may also identify specific dietary supplements for the promotion/maintenance of the microbiota transferred via the maternal FMTp, thereby optimizing pig growth and FE.

Spherical: an iterative workflow for assembling metagenomic datasets.

BACKGROUND: The consensus emerging from the study of microbiomes is that they are far more complex than previously thought, requiring better assemblies and increasingly deeper sequencing. However, current metagenomic assembly techniques regularly fail to incorporate all, or even the majority in some cases, of the sequence information generated for many microbiomes, negating this effort. This can especially bias the information gathered and the perceived importance of the minor taxa in a microbiome. RESULTS: We propose a simple but effective approach, implemented in Python, to address this problem. Based on an iterative methodology, our workflow (called Spherical) carries out successive rounds of assemblies with the sequencing reads not yet utilised. This approach also allows the user to reduce the resources required for very large datasets, by assembling random subsets of the whole in a "divide and conquer" manner. CONCLUSIONS: We demonstrate the accuracy of Spherical using simulated data based on completely sequenced genomes and the effectiveness of the workflow at retrieving lost information for taxa in three published metagenomics studies of varying sizes. Our results show that Spherical increased the amount of reads utilized in the assembly by up to 109% compared to the base assembly. The additional contigs assembled by the Spherical workflow resulted in a significant (P < 0.05) changes in the predicted taxonomic profile of all datasets analysed. Spherical is implemented in Python 2.7 and freely available for use under the MIT license. Source code and documentation is hosted publically at: https://github.com/thh32/Spherical .

Effect of short term diet restriction on gene expression in the bovine hypothalamus using next generation RNA sequencing technology.

BACKGROUND: Negative energy balance (NEB) is an imbalance between energy intake and energy requirements for lactation and body maintenance affecting high-yielding dairy cows and is of considerable economic importance due to its negative impact on fertility and health in dairy herds. It is anticipated that the cow hypothalamus experiences extensive biochemical changes during the early post partum period in an effort to re-establish metabolic homeostasis. However, there is variation in the tolerance to NEB between individual cows. In order to understand the genomic regulation of ovulation in hypothalamic tissue during NEB, mRNA transcriptional patterns between tolerant and sensitive animals were examined. A short term dietary restriction heifer model was developed which induced abrupt onset of anoestrus in some animals (Restricted Anovulatory; RA) while others maintained oestrous cyclicity (Restricted Ovulatory; RO). A third control group (C) received a higher level of normal feeding. RESULTS: A total of 15,295 genes were expressed in hypothalamic tissue. Between RA and C groups 137 genes were differentially expressed, whereas between RO and C, 32 genes were differentially expressed. Differentially expressed genes were involved in the immune response and cellular motility in RA and RO groups, respectively, compared to C group. The largest difference between groups was observed in the comparison between RA and RO heifers, with 1094 genes shown to be significantly differentially expressed (SDE). Pathway analysis showed that these SDE genes were associated with 6 canonical pathways (P < 0.01), of which neuroactive ligand-receptor interaction was the most significant. Within the comparisons the main over-represented pathway functions were immune response including neuroprotection (CXCL10, Q1KLR3, IFIH1, IL1 and IL8; RA v C and RA v RO); energy homeostasis (AgRP and NPY; RA v RO); cell motility (CADH1, DSP and TSP4; RO v C) and prevention of GnRH release (NTSR1 IL1α, IL1ß, NPY and PACA; RA v RO). CONCLUSIONS: This information will assist in understanding the genomic factors regulating the influence of diet restriction on fertility and may assist in optimising nutritional and management systems for the improvement in reproductive performance.

Differential gene expression in the endometrium reveals cytoskeletal and immunological genes in lactating dairy cows genetically divergent for fertility traits.

Profitable milk production in dairy cows requires good reproductive performance. Calving interval is a trait used to measure reproductive efficiency. Herein we used a novel lactating Holstein cow model of fertility that displayed genetic and phenotypic divergence in calving interval, a trait used to define reproductive performance using a national breeding index in Ireland. Cows had similar genetic merit for milk production traits, but either very good genetic merit for fertility (Fert+; n=7) or very poor genetic merit for fertility (Fert-; n=6). We tested the hypothesis that Fert+ cows would have a corresponding detectable difference in endometrial gene expression compared with the Fert- cows. To do this, we sequenced the transcriptome of endometrial biopsies collected on Day 7 of the oestrous cycle (non-pregnant). This is an important stage for uterine remodelling and initiation of histotroph secretion. Significant differential expression (false discovery rate-adjusted P<0.1) of 403 genes between Fert+ and Fert- cows was found. A novel network-based functional analysis highlighted 123 genes from three physiologically relevant networks of the endometrium: (1) actin and cytoskeletal components; (2) immune function; and (3) ion transportation. In particular, our results indicate an overall downregulation of inflammation-related genes and an upregulation of multiple ion transporters and gated-voltage channels and cytoskeletal genes in Fert+ cows. These three topics, which are discussed in terms of the uterus and in the context of fertility, provide molecular evidence for an association between gene expression in the uterine environment and genetic merit for fertility in dairy cows.

Transcriptomics of liver and muscle in Holstein cows genetically divergent for fertility highlight differences in nutrient partitioning and inflammation processes.

BACKGROUND: The transition between pregnancy and lactation is a major physiological change for dairy cows. Complex systemic and local processes involving regulation of energy balance, galactopoiesis, utilisation of body reserves, insulin resistance, resumption of oestrous cyclicity and involution of the uterus can affect animal productivity and hence farm profitability. Here we used an established Holstein dairy cow model of fertility that displayed genetic and phenotypic divergence in calving interval. Cows had similar genetic merit for milk production traits, but either very good genetic merit for fertility traits ('Fert+'; n = 8) or very poor genetic merit for fertility traits ('Fert-'; n = 8). We used RNA sequencing to investigate gene expression profiles in both liver and muscle tissue biopsies at three distinct time-points: late pregnancy, early lactation and mid lactation (-18, 1 and 147 days relative to parturition, respectively). RESULTS: We found 807 and 815 unique genes to be differentially expressed in at least one time-point in liver and muscle respectively, of which 79 % and 83 % were only found in a single time-point; 40 and 41 genes were found differentially expressed at every time-point indicating possible systemic or chronic dysregulation. Functional annotation of all differentially expressed genes highlighted two physiological processes that were impacted at every time-point in the study, These were immune and inflammation, and metabolic, lipid and carbohydrate-binding. CONCLUSION: These pathways have previously been identified by other researchers. We show that several specific genes which are differentially regulated, including IGF-1, might impact dairy fertility. We postulate that an increased burden of reactive oxidation species, coupled with a chronic inflammatory state, might reduce dairy cow fertility in our model.

An ultra-high density genetic linkage map of perennial ryegrass (Lolium perenne) using genotyping by sequencing (GBS) based on a reference shotgun genome assembly.

BACKGROUND AND AIMS: High density genetic linkage maps that are extensively anchored to assembled genome sequences of the organism in question are extremely useful in gene discovery. To facilitate this process in perennial ryegrass (Lolium perenne L.), a high density single nucleotide polymorphism (SNP)- and presence/absence variant (PAV)-based genetic linkage map has been developed in an F2 mapping population that has been used as a reference population in numerous studies. To provide a reference sequence to which to align genotyping by sequencing (GBS) reads, a shotgun assembly of one of the grandparents of the population, a tenth-generation inbred line, was created using Illumina-based sequencing. METHODS: The assembly was based on paired-end Illumina reads, scaffolded by mate pair and long jumping distance reads in the range of 3-40 kb, with >200-fold initial genome coverage. A total of 169 individuals from an F2 mapping population were used to construct PstI-based GBS libraries tagged with unique 4-9 nucleotide barcodes, resulting in 284 million reads, with approx. 1·6 million reads per individual. A bioinformatics pipeline was employed to identify both SNPs and PAVs. A core genetic map was generated using high confidence SNPs, to which lower confidence SNPs and PAVs were subsequently fitted in a straightforward binning approach. KEY RESULTS: The assembly comprises 424 750 scaffolds, covering 1·11 Gbp of the 2·5 Gbp perennial ryegrass genome, with a scaffold N50 of 25 212 bp and a contig N50 of 3790 bp. It is available for download, and access to a genome browser has been provided. Comparison of the assembly with available transcript and gene model data sets for perennial ryegrass indicates that approx. 570 Mbp of the gene-rich portion of the genome has been captured. An ultra-high density genetic linkage map with 3092 SNPs and 7260 PAVs was developed, anchoring just over 200 Mb of the reference assembly. CONCLUSIONS: The combined genetic map and assembly, combined with another recently released genome assembly, represent a significant resource for the perennial ryegrass genetics community.

L.U.St: a tool for approximated maximum likelihood supertree reconstruction.

BACKGROUND: Supertrees combine disparate, partially overlapping trees to generate a synthesis that provides a high level perspective that cannot be attained from the inspection of individual phylogenies. Supertrees can be seen as meta-analytical tools that can be used to make inferences based on results of previous scientific studies. Their meta-analytical application has increased in popularity since it was realised that the power of statistical tests for the study of evolutionary trends critically depends on the use of taxon-dense phylogenies. Further to that, supertrees have found applications in phylogenomics where they are used to combine gene trees and recover species phylogenies based on genome-scale data sets. RESULTS: Here, we present the L.U.St package, a python tool for approximate maximum likelihood supertree inference and illustrate its application using a genomic data set for the placental mammals. L.U.St allows the calculation of the approximate likelihood of a supertree, given a set of input trees, performs heuristic searches to look for the supertree of highest likelihood, and performs statistical tests of two or more supertrees. To this end, L.U.St implements a winning sites test allowing ranking of a collection of a-priori selected hypotheses, given as a collection of input supertree topologies. It also outputs a file of input-tree-wise likelihood scores that can be used as input to CONSEL for calculation of standard tests of two trees (e.g. Kishino-Hasegawa, Shimidoara-Hasegawa and Approximately Unbiased tests). CONCLUSION: This is the first fully parametric implementation of a supertree method, it has clearly understood properties, and provides several advantages over currently available supertree approaches. It is easy to implement and works on any platform that has python installed. AVAILABILITY: bitBucket page - https://afro-juju@bitbucket.org/afro-juju/l.u.st.git. CONTACT: Davide.Pisani@bristol.ac.uk.

Whole genome association study identifies regions of the bovine genome and biological pathways involved in carcass trait performance in Holstein-Friesian cattle.

BACKGROUND: Four traits related to carcass performance have been identified as economically important in beef production: carcass weight, carcass fat, carcass conformation of progeny and cull cow carcass weight. Although Holstein-Friesian cattle are primarily utilized for milk production, they are also an important source of meat for beef production and export. Because of this, there is great interest in understanding the underlying genomic structure influencing these traits. Several genome-wide association studies have identified regions of the bovine genome associated with growth or carcass traits, however, little is known about the mechanisms or underlying biological pathways involved. This study aims to detect regions of the bovine genome associated with carcass performance traits (employing a panel of 54,001 SNPs) using measures of genetic merit (as predicted transmitting abilities) for 5,705 Irish Holstein-Friesian animals. Candidate genes and biological pathways were then identified for each trait under investigation. RESULTS: Following adjustment for false discovery (q-value < 0.05), 479 quantitative trait loci (QTL) were associated with at least one of the four carcass traits using a single SNP regression approach. Using a Bayesian approach, 46 QTL were associated (posterior probability > 0.5) with at least one of the four traits. In total, 557 unique bovine genes, which mapped to 426 human orthologs, were within 500kbs of QTL found associated with a trait using the Bayesian approach. Using this information, 24 significantly over-represented pathways were identified across all traits. The most significantly over-represented biological pathway was the peroxisome proliferator-activated receptor (PPAR) signaling pathway. CONCLUSIONS: A large number of genomic regions putatively associated with bovine carcass traits were detected using two different statistical approaches. Notably, several significant associations were detected in close proximity to genes with a known role in animal growth such as glucagon and leptin. Several biological pathways, including PPAR signaling, were shown to be involved in various aspects of bovine carcass performance. These core genes and biological processes may form the foundation for further investigation to identify causative mutations involved in each trait. Results reported here support previous findings suggesting conservation of key biological processes involved in growth and metabolism.

Genome sequence of Ensifer adhaerens OV14 provides insights into its ability as a novel vector for the genetic transformation of plant genomes.

BACKGROUND: Recently it has been shown that Ensifer adhaerens can be used as a plant transformation technology, transferring genes into several plant genomes when equipped with a Ti plasmid. For this study, we have sequenced the genome of Ensifer adhaerens OV14 (OV14) and compared it with those of Agrobacterium tumefaciens C58 (C58) and Sinorhizobium meliloti 1021 (1021); the latter of which has also demonstrated a capacity to genetically transform crop genomes, albeit at significantly reduced frequencies. RESULTS: The 7.7 Mb OV14 genome comprises two chromosomes and two plasmids. All protein coding regions in the OV14 genome were functionally grouped based on an eggNOG database. No genes homologous to the A. tumefaciens Ti plasmid vir genes appeared to be present in the OV14 genome. Unexpectedly, OV14 and 1021 were found to possess homologs to chromosomal based genes cited as essential to A. tumefaciens T-DNA transfer. Of significance, genes that are non-essential but exert a positive influence on virulence and the ability to genetically transform host genomes were identified in OV14 but were absent from the 1021 genome. CONCLUSIONS: This study reveals the presence of homologs to chromosomally based Agrobacterium genes that support T-DNA transfer within the genome of OV14 and other alphaproteobacteria. The sequencing and analysis of the OV14 genome increases our understanding of T-DNA transfer by non-Agrobacterium species and creates a platform for the continued improvement of Ensifer-mediated transformation (EMT).

Metasecretome-selective phage display approach for mining the functional potential of a rumen microbial community.

BACKGROUND: In silico, secretome proteins can be predicted from completely sequenced genomes using various available algorithms that identify membrane-targeting sequences. For metasecretome (collection of surface, secreted and transmembrane proteins from environmental microbial communities) this approach is impractical, considering that the metasecretome open reading frames (ORFs) comprise only 10% to 30% of total metagenome, and are poorly represented in the dataset due to overall low coverage of metagenomic gene pool, even in large-scale projects. RESULTS: By combining secretome-selective phage display and next-generation sequencing, we focused the sequence analysis of complex rumen microbial community on the metasecretome component of the metagenome. This approach achieved high enrichment (29 fold) of secreted fibrolytic enzymes from the plant-adherent microbial community of the bovine rumen. In particular, we identified hundreds of heretofore rare modules belonging to cellulosomes, cell-surface complexes specialised for recognition and degradation of the plant fibre. CONCLUSIONS: As a method, metasecretome phage display combined with next-generation sequencing has a power to sample the diversity of low-abundance surface and secreted proteins that would otherwise require exceptionally large metagenomic sequencing projects. As a resource, metasecretome display library backed by the dataset obtained by next-generation sequencing is ready for i) affinity selection by standard phage display methodology and ii) easy purification of displayed proteins as part of the virion for individual functional analysis.

Mitochondrial data are not suitable for resolving placental mammal phylogeny.

Mitochondrial data have traditionally been used in reconstructing a variety of species phylogenies. The low rates of recombination and thorough characterization of mitochondrial data across vertebrate species make it a particularly attractive phylogenetic marker. The relatively low number of fully sequenced mammal genomes and the lack of extensive sampling within Superorders have posed a serious problem for reaching agreement on the placement mammal species. The use of mitochondrial data sequences from large numbers of mammals could serve to circumvent the taxon-sampling deficit. Here we assess the suitability of mitochondrial data as a phylogenetic marker in mammal phylogenetics. MtDNA datasets of mammal origin have been filtered as follows: (i) we have sampled sparsely across the phylogenetic tree, (ii) we have constrained our sampling to genes with high taxon coverage, (iii) we have categorised rates across sites in a phylogeny independent manner and have removed fast evolving sites, and (iv), we have sampled from very shallow divergence times to reduce phylogenetic conflict. However, topologies obtained using these filters are not consistent with previous studies and are discordant across different genes. Individual mitochondrial genes, and indeed all mitochondrial genes analysed as a supermatrix, resulted in poor resolution of the species phylogeny. Overall, our study highlights the limitations of mitochondrial data, not only for resolving deep divergences and but also for shallow divergences in the mammal phylogeny.

Rumen methanogenic genotypes differ in abundance according to host residual feed intake phenotype and diet type.

Methane is an undesirable end product of rumen fermentative activity because of associated environmental impacts and reduced host feed efficiency. Our study characterized the rumen microbial methanogenic community in beef cattle divergently selected for phenotypic residual feed intake (RFI) while offered a high-forage (HF) diet followed by a low-forage (LF) diet. Rumen fluid was collected from 14 high-RFI (HRFI) and 14 low-RFI (LRFI) animals at the end of both dietary periods. 16S rRNA gene clone libraries were used, and methanogen-specific tag-encoded pyrosequencing was carried out on the samples. We found that Methanobrevibacter spp. are the dominant methanogens in the rumen, with Methanobrevibacter smithii being the most abundant species. Differences in the abundance of Methanobrevibacter smithii and Methanosphaera stadtmanae genotypes were detected in the rumen of animals offered the LF compared to the HF diet while the abundance of Methanobrevibacter smithii genotypes was different between HRFI and LRFI animals irrespective of diet. Our results demonstrate that while a core group of methanogen operational taxonomic units (OTUs) exist across diet and phenotype, significant differences were observed in the distribution of genotypes within those OTUs. These changes in genotype abundance may contribute to the observed differences in methane emissions between efficient and inefficient animals.