Using provenance to manage knowledge of in silico experiments

This article offers a briefing in one of the knowledge management issues of in silico experimentation in bioinformatics. Recording of the provenance of an experiment-what was done; where, how and why, etc. is an important aspect of scientific best practice that should be extended to in silico experimentation. We will do this in the context of eScience which has been part of the move of bioinformatics towards an industrial setting. Despite the computational nature of bioinformatics, these analyses are scientific and thus necessitate their own versions of typical scientific rigour. Just as recording who, what, why, when, where and how of an experiment is central to the scientific process in laboratory science, so it should be in silico science. The generation and recording of these aspects, or provenance, of an experiment are necessary knowledge management goals if we are to introduce scientific rigour into routine bioinformatics. In Silico experimental protocols should themselves be a form of managing the knowledge of how to perform bioinformatics analyses. Several systems now exist that offer support for the generation and collection of provenance information about how a particular in silico experiment was run, what results were generated, how they were generated, etc. In reviewing provenance support, we will review one of the important knowledge management issues in bioinformatics.     

Quantifying uncertainty in the phylogenetics of Australian numeral systems

Researchers have long been interested in the evolution of culture and the ways in which change in cultural systems can be reconstructed and tracked. Within the realm of language, these questions are increasingly investigated with Bayesian phylogenetic methods. However, such work in cultural phylogenetics could be improved by more explicit quantification of reconstruction and transition probabilities. We apply such methods to numerals in the languages of Australia. As a large phylogeny with almost universal ‘low-limit'' systems, Australian languages are ideal for investigating numeral change over time. We reconstruct the most likely extent of the system at the root and use that information to explore the ways numerals evolve. We show that these systems do not increment serially, but most commonly vary their upper limits between 3 and 5. While there is evidence for rapid system elaboration beyond the lower limits, languages lose numerals as well as gain them. We investigate the ways larger numerals build on smaller bases, and show that there is a general tendency to both gain and replace 4 by combining 2 + 2 (rather than inventing a new unanalysable word ‘four''). We develop a series of methods for quantifying and visualizing the results.     

Locally Epistatic Genomic Relationship Matrices for Genomic Association and Prediction

In plant and animal breeding studies a distinction is made between the genetic value (additive plus epistatic genetic effects) and the breeding value (additive genetic effects) of an individual since it is expected that some of the epistatic genetic effects will be lost due to recombination. In this article, we argue that the breeder can take advantage of the epistatic marker effects in regions of low recombination. The models introduced here aim to estimate local epistatic line heritability by using genetic map information and combining local additive and epistatic effects. To this end, we have used semiparametric mixed models with multiple local genomic relationship matrices with hierarchical designs. Elastic-net postprocessing was used to introduce sparsity. Our models produce good predictive performance along with useful explanatory information.     

Please don't misuse the museum: ‘declines’ may be statistical

Detecting declines in populations at broad spatial scales takes enormous effort, and long‐term data are often more sparse than is desired for estimating trends, identifying drivers for population changes, framing conservation decisions, or taking management actions. Museum records and historic data can be available at large scales across multiple decades, and are therefore an attractive source of information on the comparative status of populations. However, changes in populations may be real (e.g. in response to environmental covariates) or resulting from variation in our ability to observe the true population response (also possibly related to environmental covariates). This is a (statistical) nuisance in understanding the true status of a population. Evaluating statistical hypotheses alongside more interesting ecological ones is important in the appropriate use of museum data. Two statistical considerations are generally applicable to use of museum records: first without initial random sampling, comparison with contemporary results cannot provide inference to the entire range of a species, and second the availability of only some individuals in a population may respond to environmental changes. Changes in the availability of individuals may reduce the proportion of the population that is present and able to be counted on a given survey event, resulting in an apparent decline even when population size is stable.     

Arctic warming: nonlinear impacts of sea‐ice and glacier melt on seabird foraging

Arctic climate change has profound impacts on the cryosphere, notably via shrinking sea‐ice cover and retreating glaciers, and it is essential to evaluate and forecast the ecological consequences of such changes. We studied zooplankton‐feeding little auks (Alle alle), a key sentinel species of the Arctic, at their northernmost breeding site in Franz‐Josef Land (80°N), Russian Arctic. We tested the hypothesis that little auks still benefit from pristine arctic environmental conditions in this remote area. To this end, we analysed remote sensing data on sea‐ice and coastal glacier dynamics collected in our study area across 1979–2013. Further, we recorded little auk foraging behaviour using miniature electronic tags attached to the birds in the summer of 2013, and compared it with similar data collected at three localities across the Atlantic Arctic. We also compared current and historical data on Franz‐Josef Land little auk diet, morphometrics and chick growth curves. Our analyses reveal that summer sea‐ice retreated markedly during the last decade, leaving the Franz‐Josef Land archipelago virtually sea‐ice free each summer since 2005. This had a profound impact on little auk foraging, which lost their sea‐ice‐associated prey. Concomitantly, large coastal glaciers retreated rapidly, releasing large volumes of melt water. Zooplankton is stunned by cold and osmotic shock at the boundary between glacier melt and coastal waters, creating new foraging hotspots for little auks. Birds therefore switched from foraging at distant ice‐edge localities, to highly profitable feeding at glacier melt‐water fronts within <5 km of their breeding site. Through this behavioural plasticity, little auks maintained their chick growth rates, but showed a 4% decrease in adult body mass. Our study demonstrates that arctic cryosphere changes may have antagonistic ecological consequences on coastal trophic flow. Such nonlinear responses complicate modelling exercises of current and future polar ecosystem dynamics.     

The power of statistical tests using field trial count data of nontarget organisms in environmental risk assessment of genetically modified plants

相似文献   

Linked genetic variants on chromosome 10 control ear morphology and body mass among dog breeds

Background

The domestic dog is a rich resource for mapping the genetic components of phenotypic variation due to its unique population history involving strong artificial selection. Genome-wide association studies have revealed a number of chromosomal regions where genetic variation associates with morphological characters that typify dog breeds. A region on chromosome 10 is among those with the highest levels of genetic differentiation between dog breeds and is associated with body mass and ear morphology, a common motif of animal domestication. We characterised variation in this region to uncover haplotype structure and identify candidate functional variants.

Results

We first identified SNPs that strongly associate with body mass and ear type by comparing sequence variation in a 3 Mb region between 19 breeds with a variety of phenotypes. We next genotyped a subset of 123 candidate SNPs in 288 samples from 46 breeds to identify the variants most highly associated with phenotype and infer haplotype structure. A cluster of SNPs that associate strongly with the drop ear phenotype is located within a narrow interval downstream of the gene MSRB3, which is involved in human hearing. These SNPs are in strong genetic linkage with another set of variants that correlate with body mass within the gene HMGA2, which affects human height. In addition we find evidence that this region has been under selection during dog domestication, and identify a cluster of SNPs within MSRB3 that are highly differentiated between dogs and wolves.

Conclusions

We characterise genetically linked variants that potentially influence ear type and body mass in dog breeds, both key traits that have been modified by selective breeding that may also be important for domestication. The finding that variants on long haplotypes have effects on more than one trait suggests that genetic linkage can be an important determinant of the phenotypic response to selection in domestic animals.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1702-2) contains supplementary material, which is available to authorized users.     

ATGL and HSL are not coordinately regulated in response to fuel partitioning in fasted rats

Prolonged fasting is characterized by lipid mobilization (Phase 2), followed by protein breakdown (Phase 3). Knowing that body lipids are not exhausted in Phase 3, we investigated whether changes in the metabolic status of prolonged fasted rats are associated with differences in the expression of epididymal adipose tissue proteins involved in lipid mobilization. The final body mass, body lipid content, locomotor activity and metabolite and hormone plasma levels differed between groups. Compared with fed rats, adiposity and epididymal fat mass decreased in Phase 2 (approximately two- to threefold) and Phase 3 (∼4.5-14-fold). Plasma nonesterified fatty acids (NEFA) concentrations were increased in Phase 2 (approximately twofold) and decreased in Phase 3 (approximately twofold). Daily locomotor activity was markedly increased in Phase 3 (∼11-fold). Compared with the fed state, expressions of adipose triglyceride lipase (ATGL; mRNA and protein), hormone-sensitive lipase (HSL; mRNA) and phosphorylated HSL at residue Ser660 (HSL Ser⁶⁶⁰) were increased during Phase 2 (∼1.5-2-fold). HSL (mRNA and protein) and HSL Ser⁶⁶⁰ levels were lowered during Phase 3 (∼3-12-fold). Unlike HSL and HSL Ser⁶⁶⁰, ATGL expression did not correlate with circulating NEFA, mostly due to data from animals in Phase 3. At this stage, ATGL could play an essential role for maintaining a low mobilization rate of NEFA, possibly to sustain muscle performance and hence increased locomotor activity. We conclude that ATGL and HSL are not coordinately regulated in response to changes in fuel partitioning during prolonged food deprivation, ATGL appearing as the major lipase in late fasting.     

The value of DNA sequence data for studying landscape genetics

In a recent Opinion article in Molecular Ecology, Wang (2010) emphasizes the fact that current patterns of genetic differentiation among populations reflect processes that have acted over temporal scales ranging from contemporary to ancient. He draws a sharp distinction between the fields of phylogeography (as the study of historical processes) and landscape genetics (which he restricts to very recent processes). Wang characterizes DNA sequence data as being inappropriate for the study of contemporary population processes and further states that studies which only include mitochondrial DNA or chloroplast DNA data cannot be considered part of landscape genetics. In this response, we clarify the generally accepted view that DNA sequence data can be analysed with methods that separate contemporary and historical processes. To illustrate this point, we summarize the study of Vandergast et al. (2007), which Wang mischaracterizes as being confused in terms of temporal scale. Although additional focus should be placed on the important issue of correct data interpretation, we disagree strongly with the implication that contemporary and historic processes cannot be separated in the analyses of DNA sequence data.