PEtab—Interoperable specification of parameter estimation problems in systems biology

Reproducibility and reusability of the results of data-based modeling studies are essential. Yet, there has been—so far—no broadly supported format for the specification of parameter estimation problems in systems biology. Here, we introduce PEtab, a format which facilitates the specification of parameter estimation problems using Systems Biology Markup Language (SBML) models and a set of tab-separated value files describing the observation model and experimental data as well as parameters to be estimated. We already implemented PEtab support into eight well-established model simulation and parameter estimation toolboxes with hundreds of users in total. We provide a Python library for validation and modification of a PEtab problem and currently 20 example parameter estimation problems based on recent studies.     

Unconstrained gene flow between populations of a widespread epiphytic lichen Usnea subfloridana (Parmeliaceae,Ascomycota) in Estonia

Few studies have investigated the genetic diversity of populations of common and widespread lichenized fungi using microsatellite markers, especially the relationships between different measures of genetic diversity and environmental heterogeneity. The main aim of our study was to investigate the population genetics of a widespread and mainly clonally reproducing Usnea subfloridana at the landscape scale, focusing on the comparison of lichen populations within hemiboreal forest stands. Particular attention has been paid to the genetic differentiation of lichen populations in two geographically distinct regions in Estonia and the relationships between forest characteristics and measures of genetic diversity. We genotyped 578 Usnea thalli from eleven lichen populations using seven specific fungal microsatellite markers. Measures of genetic diversity (allelic richness, Shannon's information index, Nei's unbiased genetic diversity, clonal diversity, the number of multilocus genotypes, the number of private alleles, and the minimum number of colonization events) were calculated and compared between Usnea populations. Shared haplotypes, gene flow and AMOVA analyses suggest that unconstrained gene flow and exchange of multilocus genotypes exist between the two geographically remote regions in Estonia. Stand age, mean circumference of the host tree, size of forest site and tree species composition did not show any significant influence on allelic richness, Shannon's information index, Nei's unbiased genetic diversity, clonal diversity, the number of private alleles, and the minimum number of colonization events of U. subfloridana populations. Therefore it was concluded that other factors of habitat heterogeneity could probably have a more significant effect on population genetics of U. subfloridana populations.     

Endoplasmic reticulum stress in wake-active neurons progresses with aging

Fragmentation of wakefulness and sleep are expected outcomes of advanced aging. We hypothesize that wake neurons develop endoplasmic reticulum dyshomeostasis with aging, in parallel with impaired wakefulness. In this series of experiments, we sought to more fully characterize age-related changes in wakefulness and then, in relevant wake neuronal populations, explore functionality and endoplasmic reticulum homeostasis. We report that old mice show greater sleep/wake transitions in the active period with markedly shortened wake periods, shortened latencies to sleep, and less wake time in the subjective day in response to a novel social encounter. Consistent with sleep/wake instability and reduced social encounter wakefulness, orexinergic and noradrenergic wake neurons in aged mice show reduced c-fos response to wakefulness and endoplasmic reticulum dyshomeostasis with increased nuclear translocation of CHOP and GADD34. We have identified an age-related unfolded protein response injury to and dysfunction of wake neurons. It is anticipated that these changes contribute to sleep/wake fragmentation and cognitive impairment in aging.     

Volatile secondary metabolites of Micromeria dalmatica Benth. (Lamiaceae): biosynthetical and chemotaxonomical aspects

Analysis by GC and GC/MS of the essential oil obtained from above-ground parts of Micromeria dalmatica Benth. allowed the identification of 116 components, comprising 93.6% of the total oil composition. The major compounds are 3-oxygenated p-menthane monoterpenes and were identified as pulegone (29.6%), menthone (11.7%), and piperitenone (10.8%). The chemical composition of this and additional 30 oils obtained from selected Micromeria Benth. taxa were compared by using multivariate statistical analysis (agglomerative hierarchical cluster analysis and principal component analysis (PCA)). The results of statistical analyses, as well as the domination of different concurrent p-menthane-skeleton-type monoterpene biosynthetical sub-branches in the compared M. dalmatica samples, implied the occurrence of at least two different chemotypes of the mentioned species.     

Translation termination factor eRF1 of the ciliate <Emphasis Type="Italic">Blepharisma japonicum</Emphasis> recognizes all three stop codons

In species with variant genetic codes, one or two stop codons encode amino acid residues and are not recognized by the intrinsic class I translation termination factor (eRF1). Ciliata include a large number of species with variant genetic codes. The stop codon specificity of the Blepharisma japonicum translation termination factor eRF1 was determined in an in vitro eukaryotic translation system and in an in vivo assay (a dual reporter system). It was shown that eRF1 of B. japonicum retained specificity to all three stop codons, although the efficiency of peptydyl-tRNA hydrolysis in the presence of UGA was reduced in the in vitro assay. Since Heterotrichea (including B. japonicum) are the earliest diverged lineage in the phylogenetic tree of ciliates, B. japonicum probably possesses a universal genetic code similar to the putative ciliate ancestor group.     

Spatial-Orientation Priming Impedes Rather than Facilitates the Spontaneous Control of Hand-Retraction Speeds in Patients with Parkinson’s Disease

Background

Often in Parkinson’s disease (PD) motor-related problems overshadow latent non-motor deficits as it is difficult to dissociate one from the other with commonly used observational inventories. Here we ask if the variability patterns of hand speed and acceleration would be revealing of deficits in spatial-orientation related decisions as patients performed a familiar reach-to-grasp task. To this end we use spatial-orientation priming which normally facilitates motor-program selection and asked whether in PD spatial-orientation priming helps or hinders performance.

Methods

To dissociate spatial-orientation- and motor-related deficits participants performed two versions of the task. The biomechanical version (DEFAULT) required the same postural- and hand-paths as the orientation-priming version (primed-UP). Any differences in the patients here could not be due to motor issues as the tasks were biomechanically identical. The other priming version (primed-DOWN) however required additional spatial and postural processing. We assessed in all three cases both the forward segment deliberately aimed towards the spatial-target and the retracting segment, spontaneously bringing the hand to rest without an instructed goal.

Results and Conclusions

We found that forward and retracting segments belonged in two different statistical classes according to the fluctuations of speed and acceleration maxima. Further inspection revealed conservation of the forward (voluntary) control of speed but in PD a discontinuity of this control emerged during the uninstructed retractions which was absent in NC. Two PD groups self-emerged: one group in which priming always affected the retractions and the other in which only the more challenging primed-DOWN condition was affected. These PD-groups self-formed according to the speed variability patterns, which systematically changed along a gradient that depended on the priming, thus dissociating motor from spatial-orientation issues. Priming did not facilitate the motor task in PD but it did reveal a breakdown in the spatial-orientation decision that was independent of the motor-postural path.     

Structural analyses of the ankyrin repeat domain of TRPV6 and related TRPV ion channels

Transient receptor potential (TRP) proteins are cation channels composed of a transmembrane domain flanked by large N- and C-terminal cytoplasmic domains. All members of the vanilloid family of TRP channels (TRPV) possess an N-terminal ankyrin repeat domain (ARD). The ARD of mammalian TRPV6, an important regulator of calcium uptake and homeostasis, is essential for channel assembly and regulation. The 1.7 A crystal structure of the TRPV6-ARD reveals conserved structural elements unique to the ARDs of TRPV proteins. First, a large twist between the fourth and fifth repeats is induced by residues conserved in all TRPV ARDs. Second, the third finger loop is the most variable region in sequence, length and conformation. In TRPV6, a number of putative regulatory phosphorylation sites map to the base of this third finger. Size exclusion chromatography and crystal packing indicate that the TRPV6-ARD does not assemble as a tetramer and is monomeric in solution. Adenosine triphosphate-agarose and calmodulin-agarose pull-down assays show that the TRPV6-ARD does not interact with either ligand, indicating a different functional role for the TRPV6-ARD than in the paralogous thermosensitive TRPV1 channel. Similar biochemical findings are also presented for the highly homologous mammalian TRPV5-ARD. The implications of the structural and biochemical data on the role of the ankyrin repeats in different TRPV channels are discussed.