Correction of the misclassification of species in the Portuguese collection of Cucurbita pepo L. using DNA markers

In this study, the genetic variability among 130 accessions of the Portuguese germplasm collection of Cucurbita pepo L. maintained at the Banco Portugues de Germoplasma Vegetal was assessed using AFLP (amplified fragment length polymorphism) and RAPD (random amplified polymorphic DNA) techniques for the identification of a genetically diverse core group of accessions for field phenotypic analysis. The surprisingly completely different molecular patterns exhibited by multiple accessions was later confirmed in the distribution of the putative C. pepo plants into two clusters drastically separated at a very low level of genetic similarity (DICE coefficient = 0.37). Additional analyses with RAPD and ISSR (inter single sequence repeat) markers and the introduction of standard genotypes of C. maxima L. and C. moschata L. into the analyses allowed the identification of multiple accessions of the last species wrongly included in the C. pepo collection. This study is a good example of the usefulness of DNA markers in the establishment and management of plant germplasm collections.

Occurrence of Mycobacterium avium subspecies paratuberculosis across host species and European countries with evidence for transmission between wildlife and domestic ruminants

BACKGROUND Mycobacterium avium subspecies paratuberculosis (Map) causes an infectious chronic enteritis (paratuberculosis or Johne's disease) principally of ruminants. The epidemiology of Map is poorly understood, particularly with respect to the role of wildlife reservoirs and the controversial issue of zoonotic potential (Crohn's disease). Genotypic discrimination of Map isolates is pivotal to descriptive epidemiology and resolving these issues. This study was undertaken to determine the genetic diversity of Map, enhance our understanding of the host range and distribution and assess the potential for interspecies transmission. RESULTS 164 Map isolates from seven European countries representing 19 different host species were genotyped by standardized IS900--restriction fragment length polymorphism (IS900-RFLP), pulsed-field gel electrophoresis (PFGE), amplified fragment length polymorphisms (AFLP) and mycobacterial interspersed repeat unit-variable number tandem repeat (MIRU-VNTR) analyses. Six PstI and 17 BstEII IS900-RFLP, 31 multiplex [SnaBI-SpeI] PFGE profiles and 23 MIRU-VNTR profiles were detected. AFLP gave insufficient discrimination of isolates for meaningful genetic analysis. Point estimates for Simpson's index of diversity calculated for the individual typing techniques were in the range of 0.636 to 0.664 but a combination of all three methods increased the discriminating power to 0.879, sufficient for investigating transmission dynamics. Two predominant strain types were detected across Europe with all three typing techniques. Evidence for interspecies transmission between wildlife and domestic ruminants on the same property was demonstrated in four cases, between wildlife species on the same property in two cases and between different species of domestic livestock on one property. CONCLUSION The results of this study showed that it is necessary to use multiple genotyping techniques targeting different sources of genetic variation to obtain the level of discrimination necessary to investigate transmission dynamics and trace the source of Map infections. Furthermore, the combination of genotyping techniques may depend on the geographical location of the population to be tested. Identical genotypes were obtained from Map isolated from different host species co-habiting on the same property strongly suggesting that interspecies transmission occurs. Interspecies transmission of Map between wildlife species and domestic livestock on the same property provides further evidence to support a role for wildlife reservoirs of infection.

Functional genomic analysis of heat stress in Vitis vinifera

Grapevine (Vitis vinifera) is one of most agro-economically important fruit crops worldwide, with a special relevance in Portugal where over 300 varieties are used for wine production. Due to global warming, temperature stress is currently a serious issue affecting crop production especially in temperate climates. Mobile genetic elements such as retrotransposons have been shown to be involved in environmental stress induced genetic and epigenetic modifications. In this study, sequences related to Grapevine Retrotransposon 1 (Gret1) were utilized to determine heat induced genomic and transcriptomic modifications in Touriga Nacional, a traditional Portuguese grapevine variety. For this purpose, growing canes were treated to 42 oC for four hours and leaf genomic DNA and RNA was utilized for various techniques to observe possible genomic alterations and variation in transcription levels of coding and non-coding sequences between non-treated plants and treated plants immediately after heat stress (HS-0 h) or after a 24 hour recovery period (HS-24 h). Heat stress was found to induce a significant decrease in Gret1 related sequences in HS-24 h leaves, indicating an effect of heat stress on genomic structure. In order to identify putative heat induced DNA modifications, genome wide approaches such as Amplified Fragment Length Polymorphism were utilized. This resulted in the identification of a polymorphic DNA fragment in HS-0 h and HS-24 h leaves whose sequence mapped to a genomic region flanking a house keeping gene (NADH) that is represented in multiple copies in the Vitis vinifera genome. Heat stress was also found to affect the transcript levels of various non-coding and gene coding sequences. Accordingly, quantitative real time PCR results established that Gret1 related sequences are up regulated immediately after heat stress whereas the level of transcript of genes involved in identification and repair of double strand breaks are significantly down regulated in HS-0 h plants. Taken together, the results of this work demonstrated heat stress affects both genomic integrity and transcription levels.

Genetic diversity and differentiation processes in the ploidy series of Olea europaea L.: a multiscale approach from subspecies to insular populations.

Geographical isolation and polyploidization are central concepts in plant evolution. The hierarchical organization of archipelagos in this study provides a framework for testing the evolutionary consequences for polyploid taxa and populations occurring in isolation. Using amplified fragment length polymorphism and simple sequence repeat markers, we determined the genetic diversity and differentiation patterns at three levels of geographical isolation in Olea europaea: mainland-archipelagos, islands within an archipelago, and populations within an island. At the subspecies scale, the hexaploid ssp. maroccana (southwest Morocco) exhibited higher genetic diversity than the insular counterparts. In contrast, the tetraploid ssp. cerasiformis (Madeira) displayed values similar to those obtained for the diploid ssp. guanchica (Canary Islands). Geographical isolation was associated with a high genetic differentiation at this scale. In the Canarian archipelago, the stepping-stone model of differentiation suggested in a previous study was partially supported. Within the western lineage, an east-to-west differentiation pattern was confirmed. Conversely, the easternmost populations were more related to the mainland ssp. europaea than to the western guanchica lineage. Genetic diversity across the Canarian archipelago was significantly correlated with the date of the last volcanic activity in the area/island where each population occurs. At the island scale, this pattern was not confirmed in older islands (Tenerife and Madeira), where populations were genetically homogeneous. In contrast, founder effects resulted in low genetic diversity and marked genetic differentiation among populations of the youngest island, La Palma.

Decoupled post-glacial history in mutualistic plant-insect interactions: insights from the yellow loosestrife (Lysimachia vulgaris) and its associated oil-collecting bees (Macropis europaea and M-fulvipes)

AimWe take a comparative phylogeographical approach to assess whether three species involved in a specialized oil-rewarding pollination system (i.e. Lysimachia vulgaris and two oil-collecting bees within the genus Macropis) show congruent phylogeographical trajectories during post-glacial colonization processes. Our working hypothesis is that within specialized mutualistic interactions, where each species relies on the co-occurrence of the other for survival and/or reproduction, partners are expected to show congruent evolutionary trajectories, because they are likely to have followed parallel migration routes and to have shared glacial refugia. LocationWestern Palaearctic. MethodsOur analysis relies on the extensive sampling of 104 Western Palaearctic populations (totalling 434, 159 and 74 specimens of Lysimachiavulgaris, Macropiseuropaea and Macropisfulvipes, respectively), genotyped with amplified fragment length polymorphism. Based on this, we evaluated the regional genetic diversity (Shannon diversity and allele rarity index) and genetic structure (assessed using structure, population networks, isolation-by-distance and spatial autocorrelation metrics) of each species. Finally, we compared the general phylogeographical patterns obtained. ResultsContrary to our expectations, the analyses revealed phylogeographical signals suggesting that the investigated organisms demonstrate independent post-glacial trajectories as well as distinct contemporaneous demographic parameters, despite their mutualistic interaction. Main conclusionsThe mutualistic partners investigated here are likely to be experiencing distinct and independent evolutionary dynamics because of their contrasting life-history traits (e.g. dispersal abilities), as well as distinct hubs and migration routes. Such conditions would prevent and/or erase any signature of co-structuring of lineages in space and time. As a result, the lack of phylogeographical congruence driven by differences in life-history traits might have arisen irrespective of the three species having shared similar Pleistocene glacial refugia.

An MLSA-based online scheme for the rapid identification of Stenotrophomonas isolates

An online scheme to assign Stenotrophomonas isolates to genomic groups was developed using the multilocus sequence analysis (MLSA), which is based on the DNA sequencing of selected fragments of the housekeeping genes ATP synthase alpha subunit (atpA), the recombination repair protein (recA), the RNA polymerase alpha subunit (rpoA) and the excision repair beta subunit (uvrB). This MLSA-based scheme was validated using eight of the 10 Stenotrophomonas species that have been previously described. The environmental and nosocomial Stenotrophomonas strains were characterised using MLSA, 16S rRNA sequencing and DNA-DNA hybridisation (DDH) analyses. Strains of the same species were found to have greater than 95% concatenated sequence similarity and specific strains formed cohesive readily recognisable phylogenetic groups. Therefore, MLSA appeared to be an effective alternative methodology to amplified fragment length polymorphism fingerprint and DDH techniques. Strains of Stenotrophomonas can be readily assigned through the open database resource that was developed in the current study (www.steno.lncc.br/).

Characterization of microsatellite loci and reliable genotyping in a polyploid plant, Mercurialis perennis (Euphorbiaceae).

For many applications in population genetics, codominant simple sequence repeats (SSRs) may have substantial advantages over dominant anonymous markers such as amplified fragment length polymorphisms (AFLPs). In high polyploids, however, allele dosage of SSRs cannot easily be determined and alleles are not easily attributable to potentially diploidized loci. Here, we argue that SSRs may nonetheless be better than AFLPs for polyploid taxa if they are analyzed as effectively dominant markers because they are more reliable and more precise. We describe the transfer of SSRs developed for diploid Mercurialis huetii to the clonal dioecious M. perennis. Primers were tested on a set of 54 male and female plants from natural decaploid populations. Eight of 65 tested loci produced polymorphic fragments. Binary profiles from 4 different scoring routines were used to define multilocus lineages (MLLs). Allowing for fragment differences within 1 MLL, all analyses revealed the same 14 MLLs without conflicting with merigenet, sex, or plot assignment. For semiautomatic scoring, a combination of as few as 2 of the 4 most polymorphic loci resulted in unambiguous discrimination of clones. Our study demonstrates that microsatellite fingerprinting of polyploid plants is a cost efficient and reliable alternative to AFLPs, not least because fewer loci are required than for diploids.

Fine-scale genetic structure and marginal processes in an expanding population of Biscutella laevigata L. (Brassicaceae).

Evolutionary processes acting at the expanding margins of a species' range are still poorly understood. Genetic drift is considered prevalent in marginal populations, and the maintenance of genetic diversity during recolonization might seem puzzling. To investigate such processes, a fine-scale investigation of 219 individuals was performed within a population of Biscutella laevigata (Brassicaceae), located at the leading edge of its range. The survey used amplified fragment length polymorphisms (AFLPs). As commonly reported across the whole species distribution range, individual density and genetic diversity decreased along the local axis of recolonization of this expanding population, highlighting the enduring effect of the historical colonization on present-day diversity. The self-incompatibility system of the plant may have prevented local inbreeding in newly found patches and sustained genetic diversity by ensuring gene flow from established populations. Within the more continuously populated region, spatial analysis of genetic structure revealed restricted gene flow among individuals. The distribution of genotypes formed a mosaic of relatively homogenous patches within the continuous population. This pattern could be explained by a history of expansion by long-distance dispersal followed by fine-scale diffusion (that is, a stratified dispersal combination). The secondary contact among expanding patches apparently led to admixture among differentiated genotypes where they met (that is, a reshuffling effect). This type of dynamics could explain the maintenance of genetic diversity during recolonization.

Predicting present and future intra-specific genetic structure through niche hindcasting across 24 millennia.

Paleoclimatic reconstructions coupled with species distribution models and identification of extant spatial genetic structure have the potential to provide insights into the demographic events that shape the distribution of intra-specific genetic variation across time. Using the globeflower Trollius europaeus as a case-study, we combined (1) Amplified Fragment Length Polymorphisms, (2) suites of 1000-years stepwise hindcasted species distributions and (3) a model of diffusion through time over the last 24,000 years, to trace the spatial dynamics that most likely fits the species' current genetic structure. We show that the globeflower comprises four gene pools in Europe which, from the dry period preceding the Last Glacial Maximum, dispersed while tracking the conditions fitting its climatic niche. Among these four gene pools, two are predicted to experience drastic range retraction in the near future. Our interdisciplinary approach, applicable to virtually any taxon, is an advance in inferring how climate change impacts species' genetic structures.

Ecology and life history affect different aspects of the population structure of 27 high-alpine plants.

A plant species' genetic population structure is the result of a complex combination of its life history, ecological preferences, position in the ecosystem and historical factors. As a result, many different statistical methods exist that measure different aspects of species' genetic structure. However, little is known about how these methods are interrelated and how they are related to a species' ecology and life history. In this study, we used the IntraBioDiv amplified fragment length polymorphisms data set from 27 high-alpine species to calculate eight genetic summary statistics that we jointly correlate to a set of six ecological and life-history traits. We found that there is a large amount of redundancy among the calculated summary statistics and that there is a significant association with the matrix of species traits. In a multivariate analysis, two main aspects of population structure were visible among the 27 species. The first aspect is related to the species' dispersal capacities and the second is most likely related to the species' postglacial recolonization of the Alps. Furthermore, we found that some summary statistics, most importantly Mantel's r and Jost's D, show different behaviour than expected based on theory. We therefore advise caution in drawing too strong conclusions from these statistics.

Infectivity of Lactobacillus rhamnosus and Lactobacillus paracasei isolates in a rat model of experimental endocarditis.

The potential pathogenicity of selected (potentially) probiotic and clinical isolates of Lactobacillus rhamnosus and Lactobacillus paracasei was investigated in a rat model of experimental endocarditis. In addition, adhesion properties of the lactobacilli for fibrinogen, fibronectin, collagen and laminin, as well as the killing activity of the platelet-microbicidal proteins fibrinopeptide A (FP-A) and connective tissue activating peptide 3 (CTAP-3), were assessed. The 90 % infective dose (ID(90)) of the L. rhamnosus endocarditis isolates varied between 10(6) and 10(7) c.f.u., whereas four of the six (potentially) probiotic L. rhamnosus isolates showed an ID(90) that was at least 10-fold higher (10(8) c.f.u.) (P<0.001). In contrast, the two other probiotic L. rhamnosus isolates exhibited an ID(90) (10(6) and 10(7) c.f.u.) comparable to the ID(90) of the clinical isolates of this species investigated (P>0.05). Importantly, these two probiotic isolates shared the same fluorescent amplified fragment length polymorphism cluster type as the clinical isolate showing the lowest ID(90) (10(6) c.f.u.). L. paracasei tended to have a lower infectivity than L. rhamnosus (ID(90) of 10(7) to > or =10(8) c.f.u.). All isolates had comparable bacterial counts in cardiac vegetations (P>0.05). Except for one L. paracasei strain adhering to all substrates, all tested lactobacilli adhered only weakly or not at all. The platelet peptide FP-A did not show any microbicidal activity against the tested lactobacilli, whereas CTAP-3 killed the majority of the isolates. In general, these results indicate that probiotic lactobacilli display a lower infectivity in experimental endocarditis compared with true endocarditis pathogens. However, the difference in infectivity between L. rhamnosus endocarditis and (potentially) probiotic isolates could not be explained by differences in adherence or platelet microbicidal protein susceptibility. Other disease-promoting factors may exist in these organisms and warrant further investigation.

Wheat alleles introgress into selfing wild relatives: empirical estimates from approximate Bayesian computation in Aegilops triuncialis.

Extensive gene flow between wheat (Triticum sp.) and several wild relatives of the genus Aegilops has recently been detected despite notoriously high levels of selfing in these species. Here, we assess and model the spread of wheat alleles into natural populations of the barbed goatgrass (Aegilops triuncialis), a wild wheat relative prevailing in the Mediterranean flora. Our sampling, based on an extensive survey of 31 Ae. triuncialis populations collected along a 60 km × 20 km area in southern Spain (Grazalema Mountain chain, Andalousia, totalling 458 specimens), is completed with 33 wheat cultivars representative of the European domesticated pool. All specimens were genotyped with amplified fragment length polymorphism with the aim of estimating wheat admixture levels in Ae. triuncialis populations. This survey first confirmed extensive hybridization and backcrossing of wheat into the wild species. We then used explicit modelling of populations and approximate Bayesian computation to estimate the selfing rate of Ae. triuncialis along with the magnitude, the tempo and the geographical distance over which wheat alleles introgress into Ae. triuncialis populations. These simulations confirmed that extensive introgression of wheat alleles (2.7 × 10(-4) wheat immigrants for each Ae. triuncialis resident, at each generation) into Ae. triuncialis occurs despite a high selfing rate (Fis ≈ 1 and selfing rate = 97%). These results are discussed in the light of risks associated with the release of genetically modified wheat cultivars in Mediterranean agrosystems.

Combining niche modelling and landscape genetics to study local adaptation: A novel approach illustrated using alpine plants

Understanding the factors that shape adaptive genetic variation across species niches has become of paramount importance in evolutionary ecology, especially to understand how adaptation to changing climate affects the geographic range of species. The distribution of adaptive alleles in the ecological niche is determined by the emergence of novel mutations, their fitness consequences and gene flow that connects populations across species niches. Striking demographical differences and source sink dynamics of populations between the centre and the margin of the niche can play a major role in the emergence and spread of adaptive alleles. Although some theoretical predictions have long been proposed, the origin and distribution of adaptive alleles within species niches remain untested. In this paper, we propose and discuss a novel empirical approach that combines landscape genetics with species niche modelling, to test whether alleles that confer local adaptation are more likely to occur in either marginal or central populations of species niches. We illustrate this new approach by using a published data set of 21 alpine plant species genotyped with a total of 2483 amplified fragment length polymorphisms (AFLP), distributed over more than 1733 sampling sites across the Alps. Based on the assumption that alleles that were statistically associated with environmental variables were adaptive, we found that adaptive alleles in the margin of a species niche were also present in the niche centre, which suggests that adaptation originates in the niche centre. These findings corroborate models of species range evolution, in which the centre of the niche contributes to the emergence of novel adaptive alleles, which diffuse towards niche margins and facilitate niche and range expansion through subsequent local adaptation. Although these results need to be confirmed via fitness measurements in natural populations and functionally characterised genetic sequences, this study provides a first step towards understanding how adaptive genetic variation emerges and shapes species niches and geographic ranges along environmental gradients.

Divergence with gene flow and fine-scale phylogeographical structure in the wedge-billed woodcreeper, Glyphorynchus spirurus, a Neotropical rainforest bird.

Determining the relative roles of vicariance and selection in restricting gene flow between populations is of central importance to the evolutionary process of population divergence and speciation. Here we use molecular and morphological data to contrast the effect of isolation (by mountains and geographical distance) with that of ecological factors (altitudinal gradients) in promoting differentiation in the wedge-billed woodcreeper, Glyphorynchus spirurus, a tropical forest bird, in Ecuador. Tarsus length and beak size increased relative to body size with altitude on both sides of the Andes, and were correlated with the amount of moss on tree trunks, suggesting the role of selection in driving adaptive divergence. In contrast, molecular data revealed a considerable degree of admixture along these altitudinal gradients, suggesting that adaptive divergence in morphological traits has occurred in the presence of gene flow. As suggested by mitochondrial DNA sequence data, the Andes act as a barrier to gene flow between ancient subspecific lineages. Genome-wide amplified fragment length polymorphism markers reflected more recent patterns of gene flow and revealed fine-scale patterns of population differentiation that were not detectable with mitochondrial DNA, including the differentiation of isolated coastal populations west of the Andes. Our results support the predominant role of geographical isolation in driving genetic differentiation in G. spirurus, yet suggest the role of selection in driving parallel morphological divergence along ecological gradients.

Reversible phase variation in the phnE gene, which is required for phosphonate metabolism in Escherichia coli K-12

It is known that Escherichia coli K-12 is cryptic (Phn(-)) for utilization of methyl phosphonate (MePn) and that Phn(+) variants can be selected for growth on MePn as the sole P source. Variants arise from deletion via a possible slip strand mechanism of one of three direct 8-bp repeat sequences in phnE, which restores function to a component of a putative ABC type transporter. Here we show that Phn(+) variants are present at the surprisingly high frequency of >10(-2) in K-12 strains. Amplified-fragment length polymorphism analysis was used to monitor instability in phnE in various strains growing under different conditions. This revealed that, once selection for growth on MePn is removed, Phn(+) revertants reappear and accumulate at high levels through reinsertion of the 8-bp repeat element sequence. It appears that, in K-12, phnE contains a high-frequency reversible gene switch, producing phase variation which either allows ("on" form) or blocks ("off" form) MePn utilization. The switch can also block usage of other metabolizable alkyl phosphonates, including the naturally occurring 2-aminoethylphosphonate. All K-12 strains, obtained from collections, appear in the "off" form even when bearing mutations in mutS, mutD, or dnaQ which are known to enhance slip strand events between repetitive sequences. The ability to inactivate the phnE gene appears to be unique to K-12 strains since the B strain is naturally Phn(+) and lacks the inactivating 8-bp insertion in phnE, as do important pathogenic strains for which genome sequences are known and also strains isolated recently from environmental sources.