Population-specific metabolic phenotypes of Arabidopsis lyrata ssp. petraea

Plant populations growing at the margin of their range may exhibit traits that indicate genetic differentiation and adaptation to their local abiotic environment. Here, it was investigated whether geographically separated marginal populations of Arabidopsis lyrata ssp. petraea have distinct metabolic phenotypes within the plant foliage. Seeds of A. petraea were obtained from populations along a latitudinal gradient (49-64 N), namely Germany, Wales, Sweden and Iceland and grown in a controlled cabinet environment. Targeted metabolic profiles and fingerprints were obtained at the same initial developmental stage. The free amino acid compositions were population specific, with fold differences in arginine, aspartic acid, asparagines, glycine, phenylalanine, alanine, threonine, histidine, serine and gamma-aminobutyric acid (GABA) concentrations. Sucrose, mannose and fructose concentrations were also different between populations but polyhydric alcohol concentrations were not. Principal component analysis (PCA) of metabolite fingerprints revealed metabolic phenotypes for each population. It is suggested that glucosinolates were responsible for discriminating populations within the PCA. Metabolite fingerprinting and profiling has proved to be sufficiently sensitive to identify metabolic differences between plant populations. These findings show that there is significant natural variation in metabolism among populations of A. petraea.     

Expression of a novel <Emphasis Type="Italic">OSPGYRP</Emphasis> (rice proline-, glycine- and tyrosine-rich protein) gene,which is involved in vesicle trafficking,enhanced cold tolerance in <Emphasis Type="Italic">E. coli</Emphasis>

A novel OSPGYRP gene encoding a rice proline-, glycine- and tyrosine-rich protein was isolated from cold-stress treated rice seedlings using suppression subtractive hybridization. Both amino acid sequence analysis and subcellular localization confirm that OsPGYRP is a novel protein involved in vesicle trafficking. The expression of the OSPGYRP gene was induced by cold, salt, and osmotic stress. In addition, expression of the OSPGYRP gene in E. coli increased the resistance to cold stress. These results show that OsPGYRP is a novel protein involved in vesicle trafficking and plays an important role in plant adaptation to stress. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

贵州特有药用植物毕节小檗的形态变异

为探讨贵州特有药用植物毕节小檗(Berberis guizhouensis)形态多样性的变异式样及其适应机制,对其野外居群营养器官和花部形态进行研究。对毕节小檗一直空缺的花部性状进行了描述。结果表明,毕节小檗的性状在同一居群或同一性状在不同居群中都存在较丰富的变异,变异系数为0.03~0.57。从总体来看,营养器官性状的变异程度要高于花部性状。单因素方差分析表明居群间部分花部性状存在显著差异(P0.05)。聚类分析表明小白岩居群形成1个独立的分支,这种分化很可能是环境压力选择的结果。主成分分析表明花部性状与营养性状在维持形态结构的稳定性具有协调一致的生态适应。     

Utility of the methylation-sensitive amplified polymorphism (MSAP) marker for detection of DNA methylation polymorphism and epigenetic population structure in a wild barley species (<Emphasis Type="Italic">Hordeum brevisubulatum</Emphasis>)

We report here that by using a modified scoring criterion, the methylation-sensitive amplified polymorphism or MSAP marker can be used effectively to detect polymorphism in DNA methylation patterns within and among populations of a perennial wild barley species, Hordeum brevisubulatum. Twenty-four selected individual genotypes representing four natural populations of H. brevisubulatum distributed in the Songnen Prairie in northeastern China were studied. The utility of MSAP was evidenced by its detection of high levels of polymorphism in DNA methylation patterns between individuals within a given population, and the clear inter-population differentiation in methylation patterns (methylation-based epigenetic population structure) revealed among the four populations. The resolving power of MSAP to detect DNA methylation polymorphism was found to be comparable with that of a retrotransposon-based sequence-specific amplified polymorphism marker, or SSAP, to detect genetic polymorphism in the same set of plants, suggesting that MSAP with a modified scoring criterion can be used efficiently to detect DNA methylation polymorphism and assess epigenetic population structure in natural plant populations. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Flooding effects on starch partitioning during early growth of two oak species

Pedunculate (Quercus robur L.) and sessile (Q. petraea [Matt.] Liebl.) oaks are the most common oak species in Western Europe. They are known to display different ecological requirements, particularly relative to root hypoxia induced by flooding: In a glasshouse study of seedlings, we quantified the effects of flooding on starch mobilization from cotyledons and starch partitioning. Growth and distribution of lateral roots were also measured. The above-ground growth of Q. robur was less affected by flooding than that of Q. petraea which failed to develop a second flush. Root growth was also severely inhibited, particularly in Q. petraea. In Q. robur, lateral root initiation as well as elongation was restricted to the soil surface layer. Flooding markedly reduced total growth and concentrations of in all components except stems. Starch mobilization from cotyledons was delayed by flooding, especially in Q. robur seedlings. Under flooding, the decrease of cotyledons dry mass and starch content in Q. robur was lower than in Q. petraea, whereas Q. robur displayed larger growth than Q. petraea. The features of carbohydrate management may be crucial in the observed differences in flooding tolerance of these species.     

Overexpression of a cyanobacterial phospho<Emphasis Type="Italic">enol</Emphasis>pyruvate carboxylase with diminished sensitivity to feedback inhibition in<Emphasis Type="Italic"> Arabidopsis</Emphasis> changes amino acid metabolism

Phosphoenolpyruvate carboxylase (EC 4.1.1.31) from Synechococcus vulcanus (SvPEPC) is a unique enzyme, being almost insensitive to feedback inhibition at neutral pH. In order to assess its usefulness in metabolic engineering of plants, SvPEPC was expressed in Arabidopsis thaliana (L.) Heynh. under the control of the cauliflower mosaic virus 35S promoter. About one-third of the transformants of the T1 generation showed severe visible phenotypes such as leaf bleaching and were infertile when grown on soil. However, no such phenotype was observed with Arabidopsis transformed with Zea mays L. PEPC (ZmPEPC) for C4 photosynthesis, which is normally sensitive to a feedback inhibitor, l-malate. For the SvPEPC transformants of the T2 generation, which had been derived from fertile T1 transformants, three kinds of phenotype were observed when plants were grown on an agar medium containing sucrose: Type-I plants showed poor growth and a block in true leaf development; Type-II plants produced a few true leaves, which were partially bleached; Type-III plants were apparently normal. In Type-I plants, total PEPC activity was increased about 2-fold over the control plant but there was no such increase in Type-III plants. The phenotypes of Type-I plants were rescued when the sucrose-containing agar medium was supplemented with aromatic amino acids. Measurement of the free amino acid content in whole seedlings of Type-I transformants revealed that the levels of the aromatic amino acids Phe and Tyr were lowered significantly as compared with the control plants. In contrast, the levels of several amino acids of the aspartic and glutamic families, such as Asn, Gln and Arg, were markedly enhanced (4- to 8-fold per plant fresh weight). However, when the medium was supplemented with aromatic amino acids, the levels of Asn, Gln, and Arg decreased to levels slightly higher than those of control plants, accompanied by growth recovery. Taken together, it can be envisaged that SvPEPC is capable of efficiently exerting its activity in the plant cell environment so as to cause imbalance between aromatic and non-aromatic amino acid syntheses. The growth inhibition of Type-I plants was presumed to be primarily due to a decreased availability of phosphoenolpyruvate, one of the precursors for the shikimate pathway for the synthesis of aromatic amino acids and phenylpropanoids. The possible usefulness of SvPEPC as one of the key components for installing the C₄-like pathway is proposed.Abbreviations CaMV Cauliflower mosaic virus - GUS -Glucuronidase - Kan Kanamycin - 2-ME 2-Mercaptoethanol - MS/G medium 1/2 Murashige–Skoog and 1/2 Gamborg mixed medium - PEP Phosphoenolpyruvate - PEPC Phosphoenolpyruvate carboxylase - Sv Synechococcus vulcanus - ZmPEPC Maize PEPC involved in C₄ photosynthesis     

冬虫夏草菌菌丝体中的氨基酸及其衍生物的代谢组学分析

为探究冬虫夏草菌（Ophiocordyceps sinensis）菌丝体氨基酸及其衍生物的差异性,采用超高效液相色谱串联质谱（Ultra performance liquid Chromatography, tandem mass spectrometry, UPLC-MS/MS）技术对3株冬虫夏草菌（玉树菌株XSH、达日菌株DR、贵德菌株LJ）菌丝体的氨基酸及其衍生物进行靶向定量检测,利用主成分分析（Principal component analysis, PCA）和正交偏最小二乘判别分析（Orthogonal signal correction and partial least squares-discriminant analysis, OPLS-DA）考察样品分类情况、筛选差异代谢物。结果表明,3株菌株菌丝体中均检测到72种氨基酸及其衍生物,发现除组氨酸外可合成蛋白质的19种氨基酸,菌株XSH的总氨基酸及必需氨基酸绝对含量均高于其他菌株。菌株XSH相比菌株DR和菌株LJ分别含有29种和28种差异氨基酸及其衍生物,菌株DR相比菌株LJ含有10种差异氨基酸及其衍生物,此外,3株菌株菌丝体含有3个共有的差异氨基酸及其衍生物,分别为高精氨酸、N-乙酰-L-谷氨酰胺、肌氨酸。京都基因与基因组百科全书（Kyoto encyclopedia of genes and genomes,KEGG）通路富集分析显示,差异氨基酸及其衍生物在精氨酸和脯氨酸代谢通路中更活跃,不同产区冬虫夏草菌株的氨基酸含量可能受精氨酸和脯氨酸代谢通路的影响。基于研究结果发现不同产区的冬虫夏草菌株的氨基酸及其衍生物存在较大的差异。     

Mitochondrial DNA sequence variability in red-legged partridge, <Emphasis Type="Italic">Alectoris rufa</Emphasis>, Spanish populations and the origins of genetic contamination from <Emphasis Type="Italic">A. chukar</Emphasis>

The analysis of 135 mitochondrial D-loop sequences of the Iberian autochthonous red-legged partridge (Alectoris rufa) from wild population hunting bags from various locations and fowl runs in Spain yielded 37 different haplotypes. Among these, three haplotypes correspond to chukar partridges (Alectoris chukar), indicating genetic introgression from birds illegally introduced for restocking: three individuals carrying such haplotypes where found in natural populations, one appeared among those sampled on a mass reproduction farm and the remaining 10 in another fowl-run. The geographical origin of the contaminating chukar haplotypes could be assigned to the most easterly area of the chukar partridge geographical distribution in China. Molecular diversity parameters in the A. rufa samples indicate a considerable amount of genetic variation. Φ_ST showed significant differences among populations that are not explained by geographical distance alone. Particularly, one northern population (Palencia) shows a certain degree of genetic differentiation that could reflect a previously suggested subspecies division. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Proteomic insights into adaptive responses of Saccharomyces cerevisiae to the repeated vacuum fermentation

The responses and adaptation mechanisms of the industrial Saccharomyces cerevisiae to vacuum fermentation were explored using proteomic approach. After qualitative and quantitative analyses, a total of 106 spots corresponding to 68 different proteins were identified by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. The differentially expressed proteins were involved in amino acid and carbohydrate metabolisms, various signal pathways (Ras/MAPK, Ras–cyclic adenosine monophosphate, and HOG pathway), and heat shock and oxidative responses. Among them, alternations in levels of 17 proteins associated with carbohydrate metabolisms, in particular, the upregulations of proteins involved in glycolysis, trehalose biosynthesis, and the pentose phosphate pathway, suggested vacuum-induced redistribution of the metabolic fluxes. The upregulation of 17 heat stress and oxidative response proteins indicated that multifactors contributed to oxidative stresses by affecting cell redox homeostasis. Taken together with upregulation in 14-3-3 proteins levels, 22 proteins were detected in multispots, respectively, indicating that vacuum might have promoted posttranslational modifications of some proteins in S. cerevisiae. Further investigation revealed that the elevations of the differentially expressed proteins were mainly derived from vacuum stress rather than the absence of oxygen. These findings provide new molecular mechanisms for understanding of adaptation and tolerance of yeast to vacuum fermentation. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Comparative morphology of populations of Monstera Adans. (Araceae) from natural forest fragments in Northeast Brazil using elliptic Fourier Analysis of leaf outlines

Summary  A comparative study of the leaf outline morphometrics of Monstera adansonii var. klotzschiana, M. adansonii var. laniata and M. praetermissa was carried out. The study focused on populations in isolated montane humid (brejo) forests of Ceará state in Northeast Brazil and compared them with populations from Amazonia and the Brazilian Atlantic Forest. Digitised outlines were prepared from a total of 1,695 field-collected leaf images from 20 populations, and elliptic Fourier analysis was used to generate matrices of coefficients, from which six shape variables (principal components) were extracted using Principal Components Analysis. Intra-population variability and inter-population differences were analysed with multivariate distance methods. Separate analyses were carried out for each of three leaf size classes (juvenile, submature, mature) because of the strong heteroblasty typical of this genus. Juvenile leaves were the least variable size class within populations of M. adansonii var. klotzschiana. The shape variables expressed very similar types of variation in all three size classes. The Ceará brejo populations of M. adansonii var. klotzschiana showed significant differences between mature leaf outlines in all pairwise comparisons; the Pacatuba population was the most distinct. The Ceará populations did not cluster together exclusively. In all three size classes, populations clustered together into their taxonomic groups, most clearly so in mature leaves. No correlation between morphological and geographic distance matrices was found, nor between morphological and molecular distance. The study showed that leaf outline shape is a practicable and useful quantitative trait for studying morphological variability at species, varietal and population levels. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Mitochondrial <Emphasis Type="Italic">atpA</Emphasis> gene is altered in a new <Emphasis Type="Italic">orf220</Emphasis>-type cytoplasmic male-sterile line of stem mustard (<Emphasis Type="Italic">Brassica juncea</Emphasis>)

The purpose of this research is to identify the probable mitochondrial factor associated with cytoplasmic male sterility (cms) by comparative analysis of cms and its isogenic maintainer lines in stem mustards. Dramatic variations in the morphology of floral organs were observed in cms stem mustard. Mitochondrial atpA gene was shown to be altered in cms compared with that in its maintainer line, of which mitochondrial atpA gene from its maintainer line was sequenced to encode 507 amino acids. It was indicative of high homology with mitochondrial atpA genes from other species, even as high as 94% in similarity with Oryza sativa in terms of amino acid constituents. However, only 429 amino acids were deduced in cms showing 83% similarity with atpA gene from its maintainer line. Two copies were observed in its maintainer line, but only one was found in cms. Such numerous differences of mitochondrial atpA gene between cms and its maintainer lines may not be the results of evolutionary divergence but the rearrangements of mitochondria. Expression of mitochondrial atpA gene was shown to be down-regulated in cms by using Northern blot. Consequently, mitochondrial ATP synthesis was severely decreased more than one fold in cms stem mustard indicating deficiency in mitochondrial ATP synthesis in this type of cms. Therefore, we deduced that mitochondrial atpA gene altered in cms could be associated with male-sterility in this type of cms. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. Jing-Hua Yang and Yan Huai contributed equally to this work.     

Comparative metabolomic analysis on industrial continuous and batch ethanol fermentation processes by GC-TOF-MS

The intracellular metabolic profile characterization of Saccharomyces cerevisiae throughout industrial ethanol fermentation was investigated using gas chromatography coupled to time-of-flight mass spectrometry. A total of 143 and 128 intracellular metabolites in S. cerevisiae were detected and quantified in continuous and batch fermentations, respectively. The two fermentation processes were both clearly distinguished into three main phases by principal components analysis. Furthermore, the levels of some metabolites involved in central carbon metabolism varied significantly throughout both processes. Glycerol and phosphoric acid were principally responsible for discriminating seed, main and final phases of continuous fermentation, while lactic acid and glycerol contributed mostly to telling different phases of batch fermentation. In addition, the levels of some amino acids such as glycine varied significantly during both processes. These findings provide new insights into the metabolomic characteristics during industrial ethanol fermentation processes. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Population structure and genetic diversity distribution in wild and cultivated populations of the traditional Chinese medicinal plant Magnolia officinalis subsp. biloba (Magnoliaceae)

Magnolia officinalis subsp. biloba, a traditional Chinese medicinal plant, experienced severe declines in the number of populations and the number of individuals in the late 20th century due to the widespread harvest of the subspecies. A large-scale cultivation program was initiated and cultivated populations rapidly recovered the loss in individual plant numbers, but wild populations remained small as a consequence of cutting. In this study, the levels of genetic variation and genetic structure of seven wild populations and five domestic populations of M. officinalis subsp. biloba were estimated employing an AFLP methodology. The plant exhibited a relatively high level of intra-population genetic diversity (h = 0.208 and H _j = 0.268). The cultivated populations maintained approximately 95% of the variation exhibited in wild populations, indicating a slight genetic bottleneck in the cultivated populations. The analysis of genetic differentiation revealed that most of the AFLP diversity resided within populations both for the wild group (78.22%) and the cultivated group (85.92%). Genetic differentiation among populations in the wild group was significant (F _ST = 0.1092, P < 0.005), suggesting wild population level genetic structure. Principal coordinates analysis (PCO) did not discern among wild and cultivated populations, indicating that alleles from the wild population were maintained in the cultivated gene pool. Results from the present study provide important baseline data for effectively conserving the genetic resources of this medicinal subspecies.     

Comprehensive metabolite profiling of <Emphasis Type="Italic">Sinorhizobium meliloti</Emphasis> using gas chromatography-mass spectrometry

A metabolite analysis of the soil bacterium Sinorhizobium meliloti was established as a first step towards a better understanding of the symbiosis with its host plant Medicago truncatula. A crucial step was the development of fast harvesting and extraction methods for the bacterial metabolites because of rapid changes in their composition. S. meliloti 1021 cell cultures grown in minimal medium were harvested by centrifugation, filtration or immediate freezing in liquid nitrogen followed by a lyophilisation step. Bacteria were lysed mechanically in methanol and hydrophilic compounds were analysed after methoxymation and silylisation via GC-MS. The different compounds were identified by comparison with the NIST 98 database and available standards. From about 200 peaks in each chromatogram 65 compounds have been identified so far. A comparison of the different extraction methods giving the metabolite composition revealed clear changes in several amino acids and amino acid precursor pools. A principal component analysis (PCA) was able to distinguish S. meliloti cells grown on different carbon sources based on their metabolite profile. A comparison of the metabolite composition of a S. meliloti leucine auxotrophic mutant with the wild type revealed a marked accumulation of 2-isopropylmalate in the mutant. Interestingly, the accumulated metabolite is not the direct substrate of the mutated enzyme, 3-isopropylmalate dehydrogenase, but the substrate of isopropylmalate isomerase, which acts one step further upstream in the biosynthetic pathway of leucine. This finding further emphasises the importance of integrating metabolic data into post-genomic research.     

GA-20 oxidase as a candidate for the semidwarf gene sdw1/denso in barley

The barley sdw1/denso gene not only controls plant height but also yield and quality. The sdw1/denso gene was mapped to the long arm of chromosome 3H. Comparative genomic analysis revealed that the sdw1/denso gene was located in the syntenic region of the rice semidwarf gene sd1 on chromosome 1. The sd1 gene encodes a gibberellic acid (GA)-20 oxidase enzyme. The gene ortholog of rice sd1 was isolated from barley using polymerase chain reaction. The barley and rice genes showed a similar gene structure consisting of three exons and two introns. Both genes share 88.3% genomic sequence similarity and 89% amino acid sequence identity. A single nucleotide polymorphism was identified in intron 2 between barley varieties Baudin and AC Metcalfe with Baudin known to contain the denso semidwarf gene. The single nucleotide polymorphism (SNP) marker was mapped to chromosome 3H in a doubled haploid population of Baudin × AC Metcalfe with 178 DH lines. Quantitative trait locus analysis revealed that plant height cosegregated with the SNP. The sdw1/denso gene in barley is the most likely ortholog of the sd1 in rice. The result will facilitate understanding of the molecular mechanism controlling semidwarf phenotype and provide a diagnostic marker for selection of semidwarf gene in barley. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Population cytogenetics of <Emphasis Type="Italic">Chironomus circumdatus</Emphasis> Kieffer, 1921 (Diptera,Chironomidae) from Thailand

The objectives of this study were to explore cytogenetic variation and the role of chromosomal change on local adaptation and genetic differentiation of Chironomus circumdatus Kieffer from Thailand. A total of 1,505 larvae from 24 populations were examined cytologically. Twelve chromosomal inversions were found and most of these (9 of 12) were rare inversions. All populations were in Hardy–Weinberg equilibrium. Significant association (P < 0.001) between the A2 and B5 inversions was detected in one population. Population genetic structure analysis indicated significant genetic differentiation between populations (F _ST = 0.037, P < 0.001). Geographic distance was the principal factor limiting gene flow between populations. Nei’s modified genetic distance (D _A) between populations ranged from 0.001 to 0.011 with an average of 0.003. An UPGMA population phenogram depicting relationship between populations based on D _A values revealed three groups of populations, group I, II and III each characterized by different inversions/inversion frequencies. Significant correlation of inversion C3 and water temperature suggested that this inversion might have a role to play on adaptation to high temperature habitat. However, if detection of significant population subdivision and relationship between genetic and geographic distance are taken into account, relationship between C3 and water temperature will also be due to the effect of migration/drift alone without the effect of selection.     

Two stable QTL involved in adult plant resistance to powdery mildew in the winter wheat line RE714 are expressed at different times along the growing season

Powdery mildew (Blumeria graminis f. sp. tritici) is one of the major diseases of wheat (Triticum aestivum). Adult plant resistance (APR) to powdery mildew is considered more durable than resistance conferred by major race-specific resistance genes. The objective of the present study was a better understanding of the genetic basis of APR in RE714 by means of QTL analysis of several resistance scores along the growing season. A population of 160 recombinant inbred lines obtained from the cross between RE714 and Hardi (susceptible) was assessed for APR under natural infection conditions during 3 years and a genetic map with whole genome coverage was developed with microsatellite and AFLP markers in this population. Two major QTL on chromosomes 5D and 6A were detected each year, and 6 minor QTL were detected only in 1 or 2 years. The QTL on chromosome 5D was detected during all the growing season each year and its R ² value varied between 8.5 and 56.3%, whereas the QTL on chromosome 6A was detected at 1–4 scoring dates in the 3 years, and its R ² value varied between 6.1 and 20.5%. The two QTL explained between 24.4 and 52.1% of the phenotypic variance for AUDPC, depending on the year. The models including QTL and cofactors in the composite interval mapping explained between 29 and 72% of the variance. The molecular markers linked to the two major QTL could be used in marker-assisted selection for adult plant resistance to powdery mildew. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Microsatellite analysis of the natterjack toad (Bufo calamita) in Denmark: populations are islands in a fragmented landscape

The European natterjack toad (Bufo calamita) has declined rapidly in recent years, primarily due to loss of habitat, and in Denmark it is estimated that 50% of the isolated populations are lost each decade. To efficiently manage and conserve this species and its genetic diversity, knowledge of the genetic structure is crucial. Based on nine polymorphic microsatellite loci, the genetic diversity, genetic structure and gene flow were investigated at 12 sites representing 5–10% of the natterjack toad localities presently known in Denmark. The expected heterozygosity (H _E) within each locality was generally low (range: 0.18–0.43). Further analyses failed to significantly correlate genetic diversity with population size, degree of isolation and increasing northern latitude, indicating a more complex combination of factors in determining the present genetic profile. Genetic differentiation was high (overall θ = 0.29) and analyses based on a Bayesian clustering method revealed that the dataset constituted 11 genetic clusters, defining nearly all sampling sites as distinct populations. Contemporary gene flow among populations was undetectable in nearly all cases, and the failure to detect a pattern of isolation by distance within major regions supported this apparent lack of a gene flow continuum. Indications of a genetic bottleneck were found in three populations. The analyses suggest that the remaining Bufo calamita populations in Denmark are genetically isolated, and represent independent units in a highly fragmented gene pool. Future conservation management of this species is discussed in light of these results. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.