NormalizeMets: assessing,selecting and implementing statistical methods for normalizing metabolomics data

Introduction

In metabolomics studies, unwanted variation inevitably arises from various sources. Normalization, that is the removal of unwanted variation, is an essential step in the statistical analysis of metabolomics data. However, metabolomics normalization is often considered an imprecise science due to the diverse sources of variation and the availability of a number of alternative strategies that may be implemented.

Objectives

We highlight the need for comparative evaluation of different normalization methods and present software strategies to help ease this task for both data-oriented and biological researchers.

Methods

We present NormalizeMets—a joint graphical user interface within the familiar Microsoft Excel and freely-available R software for comparative evaluation of different normalization methods. The NormalizeMets R package along with the vignette describing the workflow can be downloaded from https://cran.r-project.org/web/packages/NormalizeMets/. The Excel Interface and the Excel user guide are available on https://metabolomicstats.github.io/ExNormalizeMets.

Results

NormalizeMets allows for comparative evaluation of normalization methods using criteria that depend on the given dataset and the ultimate research question. Hence it guides researchers to assess, select and implement a suitable normalization method using either the familiar Microsoft Excel and/or freely-available R software. In addition, the package can be used for visualisation of metabolomics data using interactive graphical displays and to obtain end statistical results for clustering, classification, biomarker identification adjusting for confounding variables, and correlation analysis.

Conclusion

NormalizeMets is designed for comparative evaluation of normalization methods, and can also be used to obtain end statistical results. The use of freely-available R software offers an attractive proposition for programming-oriented researchers, and the Excel interface offers a familiar alternative to most biological researchers. The package handles the data locally in the user’s own computer allowing for reproducible code to be stored locally.

     

MM2S: personalized diagnosis of medulloblastoma patients and model systems

Background

Medulloblastoma (MB) is a highly malignant and heterogeneous brain tumour that is the most common cause of cancer-related deaths in children. Increasing availability of genomic data over the last decade had resulted in improvement of human subtype classification methods, and the parallel development of MB mouse models towards identification of subtype-specific disease origins and signaling pathways. Despite these advances, MB classification schemes remained inadequate for personalized prediction of MB subtypes for individual patient samples and across model systems. To address this issue, we developed the Medullo-Model to Subtypes ( MM2S ) classifier, a new method enabling classification of individual gene expression profiles from MB samples (patient samples, mouse models, and cell lines) against well-established molecular subtypes [Genomics 106:96-106, 2015]. We demonstrated the accuracy and flexibility of MM2S in the largest meta-analysis of human patients and mouse models to date. Here, we present a new functional package that provides an easy-to-use and fully documented implementation of the MM2S method, with additional functionalities that allow users to obtain graphical and tabular summaries of MB subtype predictions for single samples and across sample replicates. The flexibility of the MM2S package promotes incorporation of MB predictions into large Medulloblastoma-driven analysis pipelines, making this tool suitable for use by researchers.

Results

The MM2S package is applied in two case studies involving human primary patient samples, as well as sample replicates of the GTML mouse model. We highlight functions that are of use for species-specific MB classification, across individual samples and sample replicates. We emphasize on the range of functions that can be used to derive both singular and meta-centric views of MB predictions, across samples and across MB subtypes.

Conclusions

Our MM2S package can be used to generate predictions without having to rely on an external web server or additional sources. Our open-source package facilitates and extends the MM2S algorithm in diverse computational and bioinformatics contexts. The package is available on CRAN, at the following URL: https://cran.r-project.org/web/packages/MM2S/, as well as on Github at the following URLs: https://github.com/DGendoo and https://github.com/bhklab.

     

UBCG: Up-to-date bacterial core gene set and pipeline for phylogenomic tree reconstruction

Genome-based phylogeny plays a central role in the future taxonomy and phylogenetics of Bacteria and Archaea by replacing 16S rRNA gene phylogeny. The concatenated core gene alignments are frequently used for such a purpose. The bacterial core genes are defined as single-copy, homologous genes that are present in most of the known bacterial species. There have been several studies describing such a gene set, but the number of species considered was rather small. Here we present the up-to-date bacterial core gene set, named UBCG, and software suites to accommodate necessary steps to generate and evaluate phylogenetic trees. The method was successfully used to infer phylogenomic relationship of Escherichia and related taxa and can be used for the set of genomes at any taxonomic ranks of Bacteria. The UBCG pipeline and file viewer are freely available at https://www.ezbiocloud.net/tools/ubcg and https://www.ezbiocloud.net/tools/ubcg_viewer, respectively.     

popRange: a highly flexible spatially and temporally explicit Wright-Fisher simulator

Background

Sequencing and genotyping technology advancements have led to massive, growing repositories of spatially explicit genetic data and increasing quantities of temporal data (i.e., ancient DNA). These data will allow more complex and fine-scale inferences about population history than ever before; however, new methods are needed to test complex hypotheses.

Results

This article presents popRange, a forward genetic simulator, which incorporates large-scale genetic data with stochastic spatially and temporally explicit demographic and selective models. Features such as spatially and temporally variable selection coefficients and demography are incorporated in a highly flexible manner. popRange is implemented as an R package and presented with an example simulation exploring a selected allele’s trajectory in multiple subpopulations.

Conclusions

popRange allows researchers to evaluate and test complex scenarios by simulating large-scale data with complicated demographic and selective features. popRange is available for download at http://cran.r-project.org/web/packages/popRange/index.html.

     

Interspecific responses of wild African carnivores to odour of 3-mercapto-3-methylbutanol,a component of wildcat and leopard urine

The scent of 3-mercapto-3-methylbutanol (3-M-3-MB), a volatile component of leopard (Panthera pardus) and domestic cat (Felis silvestris catus) urine, released at about 10 ng/s from slow-release dispensers, elicited scent-marking from African civet (Civettictis civetta), small-spotted genet (Genetta genetta) and slender mongoose (Galerella sanguinea), as well as African wildcat (F. s. cafra). A female leopard was apparently repelled by the scent. The scent-marking and scent-rubbing by species other than African wildcats and leopards were unexpected and have important implications for the design of studies to investigate chemical communication between wild mammals and the use of camera traps to estimate animal numbers. Videos showing the behaviours referred to in this article are available at; http://www.momo-p.com/showdetail-e.php?movieid=momo161223fs01a; http://www.momo-p.com/showdetail-e.php?movieid=momo161223gs01a; http://www.momo-p.com/showdetail-e.php?movieid=momo161223gg01a.     

The scent of your enemy is my friend? The acquisition of large carnivore scent by a smaller carnivore

Scent marking is critical to intraspecific communication in many mammal species, but little is known regarding its role in communication among different species. We used 4 years of motion-triggered video to document the use of scent marking areas—termed “community scrapes”—by pumas (Puma concolor) (http://www.momo-p.com/showdetail-e.php?movieid=momo160812pc01a) and other carnivore species. We found that gray foxes (Urocyon cinereoargenteus) routinely rubbed their cheeks on puma scrapes (http://www.momo-p.com/showdetail-e.php?movieid=momo160812uc01a), and tested a series of hypotheses to determine its function. We found that gray foxes selected puma scrapes over other objects, and cheek rubbing by foxes was also correlated with how recently a puma had visited the scrape, suggesting that foxes were intent upon accumulating fresh puma scent. Cheek rubbing by foxes was not correlated with their breeding season or with how recently another fox had visited the site. Finally we found a cascading pattern in the occurrence of pumas, coyotes (Canis latrans) and gray foxes at community scrapes, suggesting that gray foxes may use puma scent to deter predation. This is the first published study to find evidence of a subordinate species using the scent of a dominant species to communicate with heterospecifics. The behavioral cascade we found in scent marking patterns also suggests that scent marking could be a mechanism that impacts the distribution and abundance of species. Additional videos pertaining to this article include http://www.momo-p.com/showdetail-e.php?movieid=momo160812uc02a, and http://www.momo-p.com/showdetail-e.php?movieid=momo160812uc03a.     

Mating and aggressive behaviour of <Emphasis Type="Italic">Brachyrhaphis olomina</Emphasis> (Cyprinodontiformes: Poeciliidae)

Despite the increasing interest in the use of intromittent male genitalia and coercive mating behaviour in poeciliids, detailed studies of the mating behaviour of most species in this family are lacking. We describe here the mating and aggressive behaviours of Brachyrhaphis olomina, and correlate them with the condition of the female’s ovum and embryos (immature, mature and pregnant). B. olomina performed a wide range of aggressive (sidle spread, tail beating, coordinate) and mating behaviours (approximation, touch, lateral display, touch-lateral display). Some behaviours (e.g. tail beating) are shared with other poeciliids, but two sexual behaviours (touch and lateral display) and one aggressive (coordinate) behaviour may be unique to B. olomina and were not reported in a previous study. Differences in male behaviour when paired with a female with mature ovum (more mating displays, no agonistic movements) suggest that males detect the female’s reproductive condition from some distance. The distinctive nature of mating behaviour in B. olomina highlights the importance of studying different species to have a better understanding of the evolution of mating and aggressive behaviours in poeciliids. Digital video images related to the article are available at http://www.momo-p.com/showdetail-e.php?movieid=momo170720bo01a, http://www.momo-p.com/showdetail-e.php?movieid=momo170720bo02a and http://www.momo-p.com/showdetail-e.php?movieid=momo170720bo03a.     

Ethological description of a fixed action pattern in a kissing bug (Triatominae): vision,gustation, proboscis extension and drinking of water and guava

Triatomines (Heteroptera, Reduviidae) are vectors of Trypanosoma cruzi, the etiological agent of Chagas disease in America. These true bugs have traditionally been considered to be blood suckers, although some species have been catalogued as being entomophagous. By using their highly specialized mouthparts, these insects have evolved a stereotyped habit which includes lifting up the proboscis, piercing and sucking, when the occasion arises. Most triatomines bite their sleeping and unaware vertebrate or invertebrate hosts, but they can also search for other targets, guided, in part, by visual and chemical stimuli. In this study, we observed that triatomines apparently visually identify a drop of water in the distance, then taste it with their legs, upon which proboscis extension and sucking ensues. This invariant behavior or fixed action pattern, observed in several triatomine species (Rhodnius prolixus, Triatoma infestans and Panstrongylus geniculatus), was also elicited by a dummy drop of water and guava fruit. We discuss evolutionary and ethological aspects of this innate behavior. Digital video images related to this article are available at http://www.momo-p.com/showdetail-e.php?movieid=momo180314rp01a and http://www.momo-p.com/showdetail-e.php?movieid=momo180314rp03a.     

Larval host records of butterflies in Japan

Using Japanese literature, we created a consolidated list of host records of butterflies in Japan. The list used the host records described in eight major illustrated reference books, two checklists, and 14 other pieces of literature. The presence of larvae on plants, the observation of larvae eating plants or insects in the field were considered as host records. We collected all species recorded in Japan. Scientific, family, and Japanese names of butterflies were consolidated using the BINRAN database (http://binran.lepimages.jp/). Scientific and Japanese names of host plants were based on the YList database (http://ylist.info/). If scientific names of host plants were not found in YList, we used scientific names based on The Plant List (http://www.theplantlist.org/). Family names of host plants were based on the Catalogue of Life database (http://www.catalogueoflife.org/). Scientific, family, and Japanese names of host insects were based on the MOKUROKU database (http://konchudb.agr.agr.kyushu-u.ac.jp/mokuroku/) for Hymenoptera and the catalogue of the Paraneoptera of Japan published by the Entomological Society of Japan for Hemiptera. We also provided the references of each host record and the original names described in the referred literature. Two datasets, HostDB and ReferenceDB, were created to include 3600 records of butterfly larval hosts in Japan, along with scientific and Japanese names of each species and a literature list. These datasets will be useful for basic and applied biological studies of butterflies. Data files are stored in the Ecological Research Data Archives (http://db.cger.nies.go.jp/JaLTER/ER_DataPapers/) and available from http://hostbj.lepumus.net/. These datasets are published under the Creative Commons License Attribution-ShareAlike 4.0 (CC BY-SA, https://creativecommons.org/licenses/by-sa/4.0/).     

Analyses of density-dependent effects are needed to understand how and when <Emphasis Type="Italic">Wolbachia</Emphasis> can control dengue vectors

Releases of Wolbachia-infected mosquitoes have been shown to be an effective method of controlling Aedes aegypti, the main vector of dengue fever, in Australia. A study in BMC Biology from Penelope Hancock and others shows that incorporation of density-dependent effects into population models can provide major improvements in understanding how and when the infected populations can become established.See research article: https://bmcbiol.biomedcentral.com/articles/10.1186/s12915-016-0319-5.     

Structural influence of gene networks on their inference: analysis of C3NET

Background

The availability of large-scale high-throughput data possesses considerable challenges toward their functional analysis. For this reason gene network inference methods gained considerable interest. However, our current knowledge, especially about the influence of the structure of a gene network on its inference, is limited.

Results

In this paper we present a comprehensive investigation of the structural influence of gene networks on the inferential characteristics of C3NET - a recently introduced gene network inference algorithm. We employ local as well as global performance metrics in combination with an ensemble approach. The results from our numerical study for various biological and synthetic network structures and simulation conditions, also comparing C3NET with other inference algorithms, lead a multitude of theoretical and practical insights into the working behavior of C3NET. In addition, in order to facilitate the practical usage of C3NET we provide an user-friendly R package, called c3net, and describe its functionality. It is available from https://r-forge.r-project.org/projects/c3net and from the CRAN package repository.

Conclusions

The availability of gene network inference algorithms with known inferential properties opens a new era of large-scale screening experiments that could be equally beneficial for basic biological and biomedical research with auspicious prospects. The availability of our easy to use software package c3net may contribute to the popularization of such methods.

Reviewers

This article was reviewed by Lev Klebanov, Joel Bader and Yuriy Gusev.

     

Paternal care behaviors of Japanese giant salamander <Emphasis Type="Italic">Andrias japonicus</Emphasis> in natural populations

Parental care among salamanders is typically provided by females. A rare case of parental care by male salamanders appears to occur in Cryptobranchidae. Yet, paternal behaviors have rarely been reported from natural populations of any Cryptobranchid salamanders, and their adaptive significance is poorly understood. The present study aimed to examine paternal care behaviors in a fully aquatic Japanese giant salamander (Andrias japonicus) in situ. At the beginning of the summer breeding season, large males, called den-masters, occupy burrows along stream banks for breeding and nesting. We videotaped post-breeding behaviors of two den-masters that stayed with the eggs, one in a natural and the other in an artificial nest in natural streams. We identified three behaviors, tail fanning, agitating and egg eating, to be parental care. Tail fanning provides oxygenated water for the eggs. We found that the den-master in the artificial nest, where dissolved oxygen level was lower, displayed tail fanning more frequently. Agitating the eggs with its head and body likely prevents yolk adhesions. The den-masters selectively ate whiter eggs that appeared to be dead or infected with water mold. This behavior, which we termed hygienic filial cannibalism, likely prevents water mold from spreading over healthy eggs. Digital video images relating to this article are available at http://www.momo-p.com/showdetail-e.php?movieid=momo140906aj01a, http://www.momo-p.com/showdetail-e.php?movieid=momo140906aj02a, http://www.momo-p.com/showdetail-e.php?movieid=momo140906aj03a and http://www.momo-p.com/showdetail-e.php?movieid=momo140906aj04a.     

Frugal fat or munificent muscle: genomic imprinting and metabolism

Variation in body composition is a popular obsession. The culturally ‘ideal’ body type is light on fat and heavy on muscle but the human population is collectively laying on fat. A new study finds antagonistic effects of two imprinted genes, Grb10 and Dlk1, on body composition in mice. These findings pose the question whether there is an evolutionary conflict between genes of maternal and paternal origin over the optimal proportions of body fat and lean muscle mass.See research article: http://www.biomedcentral.com/1741-7007/12/771     

Hsp90 depletion goes wild

Hsp90 reveals phenotypic variation in the laboratory, but is Hsp90 depletion important in the wild? Recent work from Chen and Wagner in BMC Evolutionary Biology has discovered a naturally occurring Drosophila allele that downregulates Hsp90, creating sensitivity to cryptic genetic variation. Laboratory studies suggest that the exact magnitude of Hsp90 downregulation is important. Extreme Hsp90 depletion might reactivate transposable elements and/or induce aneuploidy, in addition to revealing cryptic genetic variation.See research article http://wwww.biomedcentral.com/1471-2148/12/25     

The Year of the Wisent

Delving into European prehistory, two recent studies analyze ancient DNA from bison species depicted by our ancestors on the walls of their caves. The DNA tells a story of migrations driven by climate change but leaves some mystery clouding the genetic descent and climate preference of the still-extant wisent, otherwise known as the European bison.See research articles: https://bmcbiol.biomedcentral.com/articles/10.1186/s12915-016-0317-7 http://www.nature.com/articles/ncomms13158     

Lateral gene transfer in eukaryotes: tip of the iceberg or of the ice cube?

Lateral gene transfer (LGT) is the transmission of genes, sometimes across species barriers, outwith the classic vertical inheritance from parent to offspring. LGT is recognized as an important phenomenon that has shaped the genomes and biology of prokaryotes. Whether LGT in eukaryotes is important and widespread remains controversial. A study in BMC Biology concludes that LGT in eukaryotes is neither continuous nor prevalent and suggests a rule of thumb for judging when apparent LGT may reflect contamination.See research article: http://bmcbiol.biomedcentral.com/articles/10.1186/s12915-016-0315-9.     

<Emphasis Type="Italic">snpReady</Emphasis>: a tool to assist breeders in genomic analysis

The snpReady R package is a new instrument developed to help breeders in genomic projects such as genomic prediction and association studies. This package offers three different methods to build the genomic relationship matrix, a new imputation method for missing markers based on Wright’s theory, and a population genetic overview. Therefore, we implemented three functions (raw.data, G.matrix, and popgen). Hence, this tool allows the raw data to be transformed from different genotyping platforms to numeric matrices and performs quality control (missing data and allele frequency). Moreover, the package generates and exports four different relationship matrices (proposed by Yang et al. (N 569:565–569, 2010), VanRaden (JDS 91:4414–23, 2008), and the Gaussian kernel) depending on the purpose and software to be used in further analysis. Finally, based on the genotypic matrix, the package estimates the genetic variability, effective population size, and endogamy, among other population genetic parameters. Empirical comparisons between the method of imputation proposed and other well-known approaches have shown a lower accuracy of imputation, however, with no significant impact on the genome prediction accuracies when a lower amount of missing data is allowed. The functions and arguments were designed to carry out the preparation of genomic datasets in a straightforward, fast, and more computationally efficient way.The package and its details are available at CRAN or http://www.github.com/italo-granato/snpReady.