Using Bioconductor Package BiGGR for Metabolic Flux Estimation Based on Gene Expression Changes in Brain

Predicting the distribution of metabolic fluxes in biochemical networks is of major interest in systems biology. Several databases provide metabolic reconstructions for different organisms. Software to analyze flux distributions exists, among others for the proprietary MATLAB environment. Given the large user community for the R computing environment, a simple implementation of flux analysis in R appears desirable and will facilitate easy interaction with computational tools to handle gene expression data. We extended the R software package BiGGR, an implementation of metabolic flux analysis in R. BiGGR makes use of public metabolic reconstruction databases, and contains the BiGG database and the reconstruction of human metabolism Recon2 as Systems Biology Markup Language (SBML) objects. Models can be assembled by querying the databases for pathways, genes or reactions of interest. Fluxes can then be estimated by maximization or minimization of an objective function using linear inverse modeling algorithms. Furthermore, BiGGR provides functionality to quantify the uncertainty in flux estimates by sampling the constrained multidimensional flux space. As a result, ensembles of possible flux configurations are constructed that agree with measured data within precision limits. BiGGR also features automatic visualization of selected parts of metabolic networks using hypergraphs, with hyperedge widths proportional to estimated flux values. BiGGR supports import and export of models encoded in SBML and is therefore interoperable with different modeling and analysis tools. As an application example, we calculated the flux distribution in healthy human brain using a model of central carbon metabolism. We introduce a new algorithm termed Least-squares with equalities and inequalities Flux Balance Analysis (Lsei-FBA) to predict flux changes from gene expression changes, for instance during disease. Our estimates of brain metabolic flux pattern with Lsei-FBA for Alzheimer’s disease agree with independent measurements of cerebral metabolism in patients. This second version of BiGGR is available from Bioconductor.     

Fast and Efficient XML Data Access for Next-Generation Mass Spectrometry

Motivation

In mass spectrometry-based proteomics, XML formats such as mzML and mzXML provide an open and standardized way to store and exchange the raw data (spectra and chromatograms) of mass spectrometric experiments. These file formats are being used by a multitude of open-source and cross-platform tools which allow the proteomics community to access algorithms in a vendor-independent fashion and perform transparent and reproducible data analysis. Recent improvements in mass spectrometry instrumentation have increased the data size produced in a single LC-MS/MS measurement and put substantial strain on open-source tools, particularly those that are not equipped to deal with XML data files that reach dozens of gigabytes in size.

Results

Here we present a fast and versatile parsing library for mass spectrometric XML formats available in C++ and Python, based on the mature OpenMS software framework. Our library implements an API for obtaining spectra and chromatograms under memory constraints using random access or sequential access functions, allowing users to process datasets that are much larger than system memory. For fast access to the raw data structures, small XML files can also be completely loaded into memory. In addition, we have improved the parsing speed of the core mzML module by over 4-fold (compared to OpenMS 1.11), making our library suitable for a wide variety of algorithms that need fast access to dozens of gigabytes of raw mass spectrometric data.

Availability

Our C++ and Python implementations are available for the Linux, Mac, and Windows operating systems. All proposed modifications to the OpenMS code have been merged into the OpenMS mainline codebase and are available to the community at https://github.com/OpenMS/OpenMS.     

Early Life Experiences and Exercise Associate with Canine Anxieties

Personality and anxiety disorders across species are affected by genetic and environmental factors. Shyness-boldness personality continuum exists across species, including the domestic dog, with a large within- and across-breed variation. Domestic dogs are also diagnosed for several anxiety-related behavioral conditions, such as generalized anxiety disorders, phobias, and separation anxiety. Genetic and environmental factors contributing to personality and anxiety are largely unknown. We collected questionnaire data from a Finnish family dog population (N = 3264) in order to study the associating environmental factors for canine fearfulness, noise sensitivity, and separation anxiety. Early life experiences and exercise were found to associate with anxiety prevalence. We found that fearful dogs had less socialization experiences (p = 0.002) and lower quality of maternal care (p < 0.0001) during puppyhood. Surprisingly, the largest environmental factor associating with noise sensitivity (p < 0.0001) and separation anxiety (p = 0.007) was the amount of daily exercise; dogs with noise sensitivity and separation anxiety had less daily exercise. Our findings suggest that dogs share many of the same environmental factors that contribute to anxiety in other species as well, such as humans and rodents. Our study highlights the importance of early life experiences, especially the quality of maternal care and daily exercise for the welfare and management of the dogs, and reveals important confounding factors to be considered in the genetic characterization of canine anxiety.     

Functional Significance of Labellum Pattern Variation in a Sexually Deceptive Orchid (Ophrys heldreichii): Evidence of Individual Signature Learning Effects

Mimicking female insects to attract male pollinators is an important strategy in sexually deceptive orchids of the genus Ophrys, and some species possess flowers with conspicuous labellum patterns. The function of the variation of the patterns remains unresolved, with suggestions that these enhance pollinator communication. We investigated the possible function of the labellum pattern in Ophrys heldreichii, an orchid species in which the conspicuous and complex labellum pattern contrasts with a dark background. The orchid is pollinated exclusively by males of the solitary bee, Eucera berlandi. Comparisons of labellum patterns revealed that patterns within inflorescences are more similar than those of other conspecific plants. Field observations showed that the males approach at a great speed and directly land on flowers, but after an unsuccessful copulation attempt, bees hover close and visually scan the labellum pattern for up to a minute. Learning experiments conducted with honeybees as an accessible model of bee vision demonstrated that labellum patterns of different plants can be reliably learnt; in contrast, patterns of flowers from the same inflorescence could not be discriminated. These results support the hypothesis that variable labellum patterns in O. heldreichii are involved in flower-pollinator communication which would likely help these plants to avoid geitonogamy.     

The CATP-8/P5A-type ATPase functions in multiple pathways during neuronal patterning

The assembly of neuronal circuits involves the migrations of neurons from their place of birth to their final location in the nervous system, as well as the coordinated growth and patterning of axons and dendrites. In screens for genes required for patterning of the nervous system, we identified the catp-8/P5A-ATPase as an important regulator of neural patterning. P5A-ATPases are part of the P-type ATPases, a family of proteins known to serve a conserved function as transporters of ions, lipids and polyamines in unicellular eukaryotes, plants, and humans. While the function of many P-type ATPases is relatively well understood, the function of P5A-ATPases in metazoans remained elusive. We show here, that the Caenorhabditis elegans ortholog catp-8/P5A-ATPase is required for defined aspects of nervous system development. Specifically, the catp-8/P5A-ATPase serves functions in shaping the elaborately sculpted dendritic trees of somatosensory PVD neurons. Moreover, catp-8/P5A-ATPase is required for axonal guidance and repulsion at the midline, as well as embryonic and postembryonic neuronal migrations. Interestingly, not all axons at the midline require catp-8/P5A-ATPase, although the axons run in the same fascicles and navigate the same space. Similarly, not all neuronal migrations require catp-8/P5A-ATPase. A CATP-8/P5A-ATPase reporter is localized to the ER in most, if not all, tissues and catp-8/P5A-ATPase can function both cell-autonomously and non-autonomously to regulate neuronal development. Genetic analyses establish that catp-8/P5A-ATPase can function in multiple pathways, including the Menorin pathway, previously shown to control dendritic patterning in PVD, and Wnt signaling, which functions to control neuronal migrations. Lastly, we show that catp-8/P5A-ATPase is required for localizing select transmembrane proteins necessary for dendrite morphogenesis. Collectively, our studies suggest that catp-8/P5A-ATPase serves diverse, yet specific, roles in different genetic pathways and may be involved in the regulation or localization of transmembrane and secreted proteins to specific subcellular compartments.     

The Arabidopsis thaliana cysteine-rich receptor-like kinases CRK6 and CRK7 protect against apoplastic oxidative stress

Receptor-like kinases are important regulators of many different processes in plants. Despite their large number only a few have been functionally characterized. One of the largest subgroups of receptor-like kinases in Arabidopsis is the cysteine-rich receptor like kinases (CRKs). High sequence similarity among the CRKs has been suggested as major cause for functional redundancy. The genomic localization of CRK genes in back-to-back repeats has made their characterization through mutant analysis unpractical. Expression profiling has linked the CRKs with reactive oxygen species, important signaling molecules in plants. Here we have investigated the role of two CRKs, CRK6 and CRK7, and analyzed their role in extracellular ROS signaling. CRK6 and CRK7 are active protein kinases with differential preference for divalent cations. Our results suggest that CRK7 is involved in mediating the responses to extracellular but not chloroplastic ROS production.     

Establishing genomic tools and resources for Guizotia abyssinica (L.f.) Cass.-the development of a library of expressed sequence tags, microsatellite loci, and the sequencing of its chloroplast genome

We present an EST library, chloroplast genome sequence, and nuclear microsatellite markers that were developed for the semi-domesticated oilseed crop noug (Guizotia abyssinica) from Ethiopia. The EST library consists of 25 711 Sanger reads, assembled into 17 538 contigs and singletons, of which 4781 were functionally annotated using the Arabidopsis Information Resource (TAIR). The age distribution of duplicated genes in the EST library shows evidence of two paleopolyploidizations—a pattern that noug shares with several other species in the Heliantheae tribe (Compositae family). From the EST library, we selected 43 microsatellites and then designed and tested primers for their amplification. The number of microsatellite alleles varied between 2 and 10 (average 4.67), and the average observed and expected heterozygosities were 0.49 and 0.54, respectively. The chloroplast genome was sequenced de novo using Illumina’s sequencing technology and completed with traditional Sanger sequencing. No large re-arrangements were found between the noug and sunflower chloroplast genomes, but 1.4% of sites have indels and 1.8% show sequence divergence between the two species. We identified 34 tRNAs, 4 rRNA sequences, and 80 coding sequences, including one region (trnH-psbA) with 15% sequence divergence between noug and sunflower that may be particularly useful for phylogeographic studies in noug and its wild relatives.     

APC/C-Cdh1: From cell cycle to cellular differentiation and genomic integrity

Anaphase-promoting complex/cyclosome (APC/C) is a multifunctional ubiquitin-protein ligase that targets various substrates for proteolysis inside and outside of the cell cycle. The activation of APC/C is dependent on two WD-40 domain proteins, Cdc20 and Cdh1. While APC/Cdc20 principally regulates mitotic progression, APC/Cdh1 shows a broad spectrum of substrates in and beyond cell cycle. In the past several years, numerous biochemical and mouse genetic studies have greatly attracted our attention to the emerging role of APC/Cdh1 in genomic integrity, cellular differentiation and human diseases. This review will aim to summarize the recently expanded understanding of APC/Cdh1 in regulating biological function and how its dysfunction may lead to diseases.Key words: APC/C, Cdh1, proteolysis, genomic integrity, signal transduction, differentiation, tumorigenesis