Purification and characterization of cytochrome <Emphasis Type="Italic">c</Emphasis><Subscript><Emphasis Type="Italic">6</Emphasis></Subscript> from <Emphasis Type="Italic">Acaryochloris marina</Emphasis>

Cytochrome c ₆ , (cyt c ₆) a soluble monoheme electron transport protein, was isolated and characterized from the chlorophyll d-containing cyanobacterium Acaryochoris marina, the type strain MBIC11017. The protein was purified using ammonium sulfate precipitation, ion exchange and gel filtration column chromatography, and fast performance liquid chromatography. Its molecular mass and pI have been determined to be 8.87 kDa and less than 4.2, respectively, by mass spectrometry and isoelectrofocusing (IEF). The protein has an alpha helical structure as indicated by CD (circular dichroism) spectroscopy and a reduction midpoint potential (E _m) of +327 mV versus the normal hydrogen electrode (NHE) as determined by redox potentiometry. Its potential role in electron transfer processes is discussed.     

Natively oxidized amino acid residues in the spinach cytochrome <Emphasis Type="Italic">b</Emphasis><Subscript><Emphasis Type="Italic">6</Emphasis></Subscript><Emphasis Type="Italic">f</Emphasis> complex

The cytochrome b ₆ f complex of oxygenic photosynthesis produces substantial levels of reactive oxygen species (ROS). It has been observed that the ROS production rate by b ₆ f is 10–20 fold higher than that observed for the analogous respiratory cytochrome bc₁ complex. The types of ROS produced (O₂^{??, 1}O₂, and, possibly, H₂O₂) and the site(s) of ROS production within the b ₆ f complex have been the subject of some debate. Proposed sources of ROS have included the heme b _p , PQ _p ^?? (possible sources for O₂^??), the Rieske iron–sulfur cluster (possible source of O₂^?? and/or ¹O₂), Chl a (possible source of ¹O₂), and heme c _n (possible source of O₂^?? and/or H₂O₂). Our working hypothesis is that amino acid residues proximal to the ROS production sites will be more susceptible to oxidative modification than distant residues. In the current study, we have identified natively oxidized amino acid residues in the subunits of the spinach cytochrome b ₆ f complex. The oxidized residues were identified by tandem mass spectrometry using the MassMatrix Program. Our results indicate that numerous residues, principally localized near p-side cofactors and Chl a, were oxidatively modified. We hypothesize that these sites are sources for ROS generation in the spinach cytochrome b ₆ f complex.     

<Emphasis Type="Italic">Galleria</Emphasis><Emphasis Type="Italic">mellonella</Emphasis> are Resistant to <Emphasis Type="Italic">Pneumocystis</Emphasis><Emphasis Type="Italic">murina</Emphasis> Infection

Studying Pneumocystis has proven to be a challenge from the perspective of propagating a significant amount of the pathogen in a facile manner. The study of several fungal pathogens has been aided by the use of invertebrate model hosts. Our efforts to infect the invertebrate larvae Galleria mellonella with Pneumocystis proved futile since P. murina neither caused disease nor was able to proliferate within G. mellonella. It did, however, show that the pathogen could be rapidly cleared from the host.     

<Emphasis Type="Italic">Ty</Emphasis>1<Emphasis Type="Italic">/copia</Emphasis>- and <Emphasis Type="Italic">Ty</Emphasis>3<Emphasis Type="Italic">/gypsy</Emphasis>-like DNA sequences in <Emphasis Type="Italic">Helianthus </Emphasis>species

Two repeated DNA sequences isolated from a partial genomic DNA library of Helianthus annuus, p HaS13 and p HaS211, were shown to represent portions of the int gene of a Ty3 /gypsy retroelement and of the RNase-Hgene of a Ty1 /copia retroelement, respectively. Southern blotting patterns obtained by hybridizing the two probes to BglII- or DraI-digested genomic DNA from different Helianthus species showed p HaS13 and p HaS211 were parts of dispersed repeats at least 8 and 7 kb in length, respectively, that were conserved in all species studied. Comparable hybridization patterns were obtained in all species with p HaS13. By contrast, the patterns obtained by hybridizing p HaS211 clearly differentiated annual species from perennials. The frequencies of p HaS13- and p HaS211-related sequences in different species were 4.3x10(4)-1.3x10(5) copies and 9.9x10(2)-8.1x10(3) copies per picogram of DNA, respectively. The frequency of p HaS13-related sequences varied widely within annual species, while no significant difference was observed among perennial species. Conversely, the frequency variation of p HaS211-related sequences was as large within annual species as within perennials. Sequences of both families were found to be dispersed along the length of all chromosomes in all species studied. However, Ty3 /gypsy-like sequences were localized preferentially at the centromeric regions, whereas Ty1/ copia-like sequences were less represented or absent around the centromeres and plentiful at the chromosome ends. These findings suggest that the two sequence families played a role in Helianthusgenome evolution and species divergence, evolved independently in the same genomic backgrounds and in annual or perennial species, and acquired different possible functions in the host genomes.     

Over-expressing a UDP-glucosyltransferase gene (<Emphasis Type="Italic">Ta-UGT</Emphasis><Subscript><Emphasis Type="Italic">3</Emphasis></Subscript>) enhances <Emphasis Type="Italic">Fusarium</Emphasis> Head Blight resistance of wheat

Wheat Fusarium Head Blight (FHB), mainly caused by Fusarium graminearum (F.g), is a destructive fungal disease worldwide. FHB can not only cause considerable reduction in yield, but more seriously, can contaminate grain by trichothecene toxins released by the fungus. Here, we report new insights into the function and underlying mechanisms of a UDP-glycosyltransferase gene, Ta-UGT ₃ , that is involved in FHB resistance in wheat. In our previous study, Ta-UGT ₃ was found to enhance host tolerance against deoxynivalenol (DON) in Arabidopsis. In this study, four transgenic lines over-expressing Ta-UGT ₃ in a FHB highly susceptible wheat variety, Alondra’s, were obtained and characterized. 3 years of assays using single floret inoculation with F.g indicated that all four transgenic lines exhibited significantly enhanced type II resistance to FHB and less DON accumulation in the grains compared to the untransformed control. Histological observation using GFP labelled F.g was in agreement with the above test results since over-expression of Ta-UGT ₃ dramatically inhibited expansion of F.g. To explore the putative mechanism of resistance mediated by Ta-UGT ₃ , microarray analysis, qRT-PCR and hormone measurements were performed. Microarray analysis showed that DON up-regulated genes, such as TaNPR1, in the susceptible control, and down-regulated genes in F.g inoculated transgenic lines, while qRT-PCR showed that some defence related genes were up-regulated in F.g inoculated transgenic lines. Ta-UGT ₃ over-expression also changed the contents of the endogenous hormones SA and JA in the spikes. These data suggest that Ta-UGT ₃ positively regulates the defence responses to F.g, perhaps by regulating defence-related and DON-induced downstream genes.     

Double deletion of <Emphasis Type="Italic">dtsR1</Emphasis> and <Emphasis Type="Italic">pyc</Emphasis> induce efficient <Emphasis Type="SmallCaps">l</Emphasis>-glutamate overproduction in <Emphasis Type="Italic">Corynebacterium glutamicum</Emphasis>

Corynebacterium glutamicum strains are used for the fermentative production of l-glutamate. Five C. glutamicum deletion mutants were isolated by two rounds of selection for homologous recombination and identified by Southern blot analysis. The growth, glucose consumption and glutamate production of the mutants were analyzed and compared with the wild-type ATCC 13032 strain. Double disruption of dtsR1 (encoding a subunit of acetyl-CoA carboxylase complex) and pyc (encoding pyruvate carboxylase) caused efficient overproduction of l-glutamate in C. glutamicum; production was much higher than that of the wild-type strain and ΔdtsR1 strain under glutamate-inducing conditions. In the absence of any inducing conditions, the amount of glutamate produced by the double-deletion strain ΔdtsR1Δpyc was more than that of the mutant ΔdtsR1. The activity of phosphoenolpyruvate carboxylase (PEPC) was found to be higher in the ΔdtsR1Δpyc strain than in the ΔdtsR1 strain and the wild-type strain. Therefore, PEPC appears to be an important anaplerotic enzyme for glutamate synthesis in ΔdtsR1 derivatives. Moreover, this conclusion was confirmed by overexpression of ppc and pyc in the two double-deletion strains (ΔdtsR1Δppc and ΔdtsR1Δpyc), respectively. Based on the data generated in this investigation, we suggest a new method that will improve glutamate production strains and provide a better understanding of the interaction(s) between the anaplerotic pathway and fatty acid synthesis.     

Electrostatic attraction between cytochrome <Emphasis Type="Italic">bc</Emphasis><Subscript>1</Subscript> and cytochrome <Emphasis Type="Italic">c</Emphasis> affects kinetics of cytochrome <Emphasis Type="Italic">c</Emphasis> reduction

The kinetics of the ubiquinol-cytochrome c reductase reaction was examined using membrane fragments and purified bc(1) complexes derived from a wild-type (WT) and a newly constructed mutant (MUT) strains of Paracoccus denitrificans. The cytochrome c(1) of the WT samples possessed an additional stretch of acidic amino acids, which was lacking in the mutant. The reaction was followed with positively charged mitochondrial and negatively charged bacterial cytochromes c, and specific activities, apparent k(cat) values, and first-order rate constant values were compared. These values were distinctly lower for the MUT fractions using mitochondrial cytochrome c but differed only slightly with the bacterial species. The MUT preparations were less sensitive to changes of ionic strength of the reaction media and showed pure first-order kinetics with both samples of cytochrome c. The reaction of the WT enzyme was first order only with bacterial cytochrome c but proceeded with a non-linear profile with mitochondrial cytochrome c. The analysis of the reaction pattern revealed a rapid onset of the reaction with a successively declining rate. Experiments performed in the absence of an electron donor indicated that electrostatic attraction could directly participate in cytochrome c reduction.     

Characterization of the cytochrome <Emphasis Type="Italic">b</Emphasis><Subscript>6</Subscript><Emphasis Type="Italic">f</Emphasis> complex from marine green alga, <Emphasis Type="Italic">Bryopsis corticulans</Emphasis>

A pure, active cytochrome b ₆ f was isolated from the chloroplasts of the marine green alga, Bryopsis corticulans. To investigate and characterize this cytochrome b ₆ f complex, sodium dodecyl sulfate–polyacrylamide gel electrophoresis (SDS–PAGE), absorption spectra measurement and HPLC were employed. It was shown that this purified complex contained four large subunits with apparent molecular masses of 34.8, 24, 18.7 and 16.7 kD. The ratio of Cyt b ₆ to Cytf was 2.01 : 1. The cytochromeb ₆ f was shown to catalyze the transfer of 73 electrons from decylplastoquinol to plastocyanin–ferricyanide per Cyt f per second. α-Carotene, one kind of carotenoid that has not been found to present in cytochrome b ₆ f complex, was discovered in this preparation by reversed phase HPLC. It was different from β-carotene usually found in cytochrome b ₆ f complex. The configuration of the major α-carotene component was assigned to be 9-cis by resonance Raman spectroscopy. Different from the previous reports, the configuration of this α-carotene in dissociated state was determined to be all-trans. Besides this carotene, chlorophyll a was also found in this complex. It was shown that the molecular ratios of chlorophylla, cis and all-trans-α-carotene to Cyt f in this complex were 1.2, 0.7 and 0.2, respectively.     

<Emphasis Type="Italic">Agrobacterium</Emphasis><Emphasis Type="Italic">tumefaciens</Emphasis>-mediated transformation of callus cells of <Emphasis Type="Italic">Crataegus</Emphasis><Emphasis Type="Italic">aronia</Emphasis>

A genetic transformation system has been developed for callus cells of Crataegus aronia using Agrobacterium tumefaciens. Callus culture was established from internodal stem segments incubated on Murashige and Skoog (MS) medium supplemented with 5 mg l⁻¹ Indole-3-butyric acid (IBA) and 0.5 mg l⁻¹ 6-benzyladenine (BA). In order to optimize the callus culture system with respect to callus growth and coloration, different types and concentrations of plant growth regulators were tested. Results indicated that the best average fresh weight of red colored callus was obtained on MS medium supplemented with 2 mg l⁻¹ 2,4-dichlorophenoxyacetic acid (2,4-D) and 1.5 mg l⁻¹ kinetin (Kin) (callus maintenance medium). Callus cells were co-cultivated with Agrobacterium harboring the binary plasmid pCAMBIA1302 carrying the mgfp5 and hygromycin phosphotransferase (hptII) genes conferring green fluorescent protein (GFP) activity and hygromycin resistance, respectively. Putative transgenic calli were obtained 4 weeks after incubation of the co-cultivated explants onto maintenance medium supplemented with 50 mg l⁻¹ hygromycin. Molecular analysis confirmed the integration of the transgenes in transformed callus. To our knowledge, this is the first time to report an Agrobacterium-mediated transformation system in Crataegus aronia.     

Revision of the <Emphasis Type="Italic">Serrulati</Emphasis> species of <Emphasis Type="Italic">Penstemon</Emphasis> section <Emphasis Type="Italic">Saccanthera</Emphasis> subsection <Emphasis Type="Italic">Serrulati</Emphasis>

A revision of Penstemon sect. Saccanthera subsect. Serrulati includes a new species (P. salmonensis), a new variety (P. triphyllus var. infernalis), and the elevation of a subspecies to species (P. curtiflorus), bringing the total number of species to eight, which are keyed and described, complete with nomenclature and type citations.     

The <Emphasis Type="Italic">Caenorhabditis</Emphasis><Emphasis Type="Italic">briggsae</Emphasis> genome contains active <Emphasis Type="Italic">CbmaT1</Emphasis> and <Emphasis Type="Italic">Tcb1</Emphasis> transposons

The maT clade of transposons is a group of transposable elements intermediate in sequence and predicted protein structure to mariner and Tc transposons, with a distribution thus far limited to a few invertebrate species. We present evidence, based on searches of publicly available databases, that the nematode Caenorhabditis briggsae has several maT-like transposons, which we have designated as CbmaT elements, dispersed throughout its genome. We also describe two additional transposon sequences that probably share their evolutionary history with the CbmaT transposons. One resembles a fold back variant of a CbmaT element, with long (380-bp) inverted terminal repeats (ITRs) that show a high degree (71%) of identity to CbmaT1. The other, which shares only the 26-bp ITR sequences with one of the CbmaT variants, is present in eight nearly identical copies, but does not have a transposase gene and may therefore be cross mobilised by a CbmaT transposase. Using PCR-based mobility assays, we show that CbmaT1 transposons are capable of excising from the C. briggsae genome. CbmaT1 excised approximately 500 times less frequently than Tcb1 in the reference strain AF16, but both CbmaT1 and Tcb1 excised at extremely high frequencies in the HK105 strain. The HK105 strain also exhibited a high frequency of spontaneous induction of unc-22 mutants, suggesting that it may be a mutator strain of C. briggsae.     

pH and CO<Subscript>2</Subscript> effects on <Emphasis Type="Italic">Coelastrella</Emphasis> (<Emphasis Type="Italic">Scotiellopsis</Emphasis>) <Emphasis Type="Italic">rubescens</Emphasis> growth and metabolism

We studied effects of рН and СО₂ enrichment on the physiological condition and biochemical composition of a carotenogenic microalga Coelastrella (Scotiellopsis Vinatzer) rubescens Kaufnerová et Eliás (Scenedesmaceae, Sphaeropleales, Chlorophyceae), a promising source of natural astaxanthin. The microalga was grown at a constant pH (5, 6, 7 or 8) maintained by direct СО₂ injection. The air-sparged culture served as the control. Cell division rate and size, dry biomass productivity, the rates of nitrogen and phosphorus uptake as well as photosynthetic pigment and total lipid content and fatty acid composition were followed. С. rubescens possessed a narrow-range рН tolerance (the optimum рН 6–7). Under these conditions, the highest values of the maximum (1.0–1.1 1/day) and average (0.3–0.35 1/day) specific growth rate, chlorophyll а (4.8–4.9%) and total carotenoid dry weight percentages (1.7–1.8%) were recorded. Cell lipid fatty acid unsaturation index (1.851) and polyunsaturated fatty acid percentage (36–39%) and С18:3 ω3/С18:1 ω9 ratio (3.8–4.5) were also the highest under these conditions. A decline of рН to 5 brought about severe stress manifesting itself as a cell division cessation, photosynthetic apparatus reduction, two-fold increase in cell volume, accumulation of dry weight and lipids and a considerable decline in fatty acid unsaturation. Cultivation of С. rubescens without СО₂ enrichment resulted in a rapid alkalization of the medium to рН 9.5–10.5 impairing the physiological condition of the cells. Reasons of the deteriorative effects of suboptimal pH values on the physiological condition of C. rubescens are discussed.     

<Emphasis Type="SmallCaps">l</Emphasis>-Proline nutrition and catabolism in <Emphasis Type="Italic">Staphylococcus saprophyticus</Emphasis>

Staphylococcus saprophyticus strains ATCC 15305, ATCC 35552, and ATCC 49907 were found to require l-proline but not l-arginine for growth in a defined culture medium. All three strains could utilize l-ornithine as a proline source and contained l-ornithine aminotransferase and Δ¹-pyrroline-5-carboxylate reductase activities; strains ATCC 35552 and ATCC 49907 could use l-arginine as a proline source and had l-arginase activity. The proline requirement also could be met by l-prolinamide, l-proline methyl ester, and the dipeptides l-alanyl-l-proline and l-leucyl-l-proline. The bacteria exhibited l-proline degradative activity as measured by the formation of Δ¹-pyrroline-5-carboxylate. The specific activity of proline degradation was not affected by addition of l-proline or NaCl but was highest in strain ATCC 49907 after growth in Mueller–Hinton broth. A membrane fraction from this strain had l-proline dehydrogenase activity as detected both by reaction of Δ¹-pyrroline-5-carboxylate with 2-aminobenzaldehyde (0.79 nmol min⁻¹ mg⁻¹) and by the proline-dependent reduction of p-iodonitrotetrazolium (20.1 nmol min⁻¹ mg⁻¹). A soluble fraction from this strain had Δ¹-pyrroline-5-carboxylate dehydrogenase activity (88.8 nmol min⁻¹ mg⁻¹) as determined by the NAD⁺-dependent oxidation of dl-Δ¹-pyrroline-5-carboxylate. Addition of l-proline to several culture media did not increase the growth rate or final yield of bacteria but did stimulate growth during osmotic stress. When grown with l-ornithine as the proline source, S. saprophyticus was most susceptible to the proline analogues L-azetidine-2-carboylate, 3,4-dehydro-dl-proline, dl-thiazolidine-2-carboxylate, and l-thiazolidine-4-carboxylate. These results indicate that proline uptake and metabolism may be a potential target of antimicrobial therapy for this organism.     

Biotechnological production of <Emphasis Type="SmallCaps">l</Emphasis>-ribose from <Emphasis Type="SmallCaps">l</Emphasis>-arabinose

l-Ribose is a rare and expensive sugar that can be used as a precursor for the production of l-nucleoside analogues, which are used as antiviral drugs. In this work, we describe a novel way of producing l-ribose from the readily available raw material l-arabinose. This was achieved by introducing l-ribose isomerase activity into l-ribulokinase-deficient Escherichia coli UP1110 and Lactobacillus plantarum BPT197 strains. The process for l-ribose production by resting cells was investigated. The initial l-ribose production rates at 39°C and pH 8 were 0.46 ± 0.01 g g⁻¹ h⁻¹ (1.84 ± 0.03 g l⁻¹ h⁻¹) and 0.27 ± 0.01 g g⁻¹ h⁻¹ (1.91 ± 0.1 g l⁻¹ h⁻¹) for E. coli and for L. plantarum, respectively. Conversions were around 20% at their highest in the experiments. Also partially purified protein precipitates having both l-arabinose isomerase and l-ribose isomerase activity were successfully used for converting l-arabinose to l-ribose.     

Stable feature selection based on the ensemble <Emphasis Type="Italic">L</Emphasis><Subscript><Emphasis Type="Italic">1</Emphasis></Subscript>-norm support vector machine for biomarker discovery

Background

Lately, biomarker discovery has become one of the most significant research issues in the biomedical field. Owing to the presence of high-throughput technologies, genomic data, such as microarray data and RNA-seq, have become widely available. Many kinds of feature selection techniques have been applied to retrieve significant biomarkers from these kinds of data. However, they tend to be noisy with high-dimensional features and consist of a small number of samples; thus, conventional feature selection approaches might be problematic in terms of reproducibility.

Results

In this article, we propose a stable feature selection method for high-dimensional datasets. We apply an ensemble L ₁ -norm support vector machine to efficiently reduce irrelevant features, considering the stability of features. We define the stability score for each feature by aggregating the ensemble results, and utilize backward feature elimination on a purified feature set based on this score; therefore, it is possible to acquire an optimal set of features for performance without the need to set a specific threshold. The proposed methodology is evaluated by classifying the binary stage of renal clear cell carcinoma with RNA-seq data.

Conclusion

A comparison with established algorithms, i.e., a fast correlation-based filter, random forest, and an ensemble version of an L ₂ -norm support vector machine-based recursive feature elimination, enabled us to prove the superior performance of our method in terms of classification as well as stability in general. It is also shown that the proposed approach performs moderately on high-dimensional datasets consisting of a very large number of features and a smaller number of samples. The proposed approach is expected to be applicable to many other researches aimed at biomarker discovery.

     

<Emphasis Type="Italic">Pl</Emphasis><Subscript><Emphasis Type="Italic">Arg</Emphasis></Subscript> from <Emphasis Type="Italic">Helianthus argophyllus</Emphasis> is unlinked to other known downy mildew resistance genes in sunflower

The Pl _Arg locus in the sunflower (Helianthus annuus L.) inbred line Arg1575-2 conferring resistance to at least four tested races (300, 700, 730, 770) of downy mildew (Plasmopara halstedii) was localized by the use of simple sequence repeat (SSR) markers. Bulked segregant analysis (BSA) was conducted on 126 individuals of an F₂ progeny from a cross between a downy mildew susceptible line, CmsHA342, and Arg1575-2. Twelve SSR markers linked to the Pl _Arg locus were identified. All markers were located proximal to Pl _Arg on linkage group LG1 based on the map of Yu et al. (2003) in a window of 9.3 cM. Since Pl _Arg was mapped to a linkage group different from all other Pl genes previously mapped with SSRs, it can be concluded that Pl _Arg provides a new source of resistance against P. halstedii in sunflower.     

High-light-induced superoxide anion radical formation in cytochrome <Emphasis Type="Italic">b</Emphasis><Subscript><Emphasis Type="Italic">6</Emphasis></Subscript><Emphasis Type="Italic">f</Emphasis> complex from spinach as detected by EPR spectroscopy

The generation of superoxide anion radical (O₂ ^·−) in the cytochrome b ₆ f complex (Cyt b ₆ f) of spinach under high-light illumination was studied using electron paramagnetic resonance spectroscopy. The generation of O₂ ^·− was lost in the absence of molecular oxygen. It was also suppressed in the presence of NaN₃ and could be scavenged by extraneous antioxidants such as ascorbate, β-carotene, and glutathione. The results also indicate that O₂ ^·−, which is produced under high-light illumination of the Cyt b ₆ f from spinach, might be generated from a reaction involing ¹O₂, and the Rieske Fe-S protein could serve as the electron donor in the O₂ ^·− production. The mechanism of photoprotection of the Cyt b ₆ f complex by antioxidants is discussed.     

Production of <Emphasis Type="SmallCaps">l</Emphasis>-alanine by metabolically engineered <Emphasis Type="Italic">Escherichia coli</Emphasis>

Escherichia coli W was genetically engineered to produce l-alanine as the primary fermentation product from sugars by replacing the native d-lactate dehydrogenase of E. coli SZ194 with alanine dehydrogenase from Geobacillus stearothermophilus. As a result, the heterologous alanine dehydrogenase gene was integrated under the regulation of the native d-lactate dehydrogenase (ldhA) promoter. This homologous promoter is growth-regulated and provides high levels of expression during anaerobic fermentation. Strain XZ111 accumulated alanine as the primary product during glucose fermentation. The methylglyoxal synthase gene (mgsA) was deleted to eliminate low levels of lactate and improve growth, and the catabolic alanine racemase gene (dadX) was deleted to minimize conversion of l-alanine to d-alanine. In these strains, reduced nicotinamide adenine dinucleotide oxidation during alanine biosynthesis is obligately linked to adenosine triphosphate production and cell growth. This linkage provided a basis for metabolic evolution where selection for improvements in growth coselected for increased glycolytic flux and alanine production. The resulting strain, XZ132, produced 1,279 mmol alanine from 120 g l⁻¹ glucose within 48 h during batch fermentation in the mineral salts medium. The alanine yield was 95% on a weight basis (g g⁻¹ glucose) with a chiral purity greater than 99.5% l-alanine. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.