Conservation phylogeography: does historical diversity contribute to regional vulnerability in European tree frogs (Hyla arborea)?

Documenting and preserving the genetic diversity of populations, which conditions their long-term survival, have become a major issue in conservation biology. The loss of diversity often documented in declining populations is usually assumed to result from human disturbances; however, historical biogeographic events, otherwise known to strongly impact diversity, are rarely considered in this context. We apply a multilocus phylogeographic study to investigate the late-Quaternary history of a tree frog (Hyla arborea) with declining populations in the northern and western part of its distribution range. Mitochondrial and nuclear polymorphisms reveal high genetic diversity in the Balkan Peninsula, with a spatial structure moulded by the last glaciations. While two of the main refugial lineages remained limited to the Balkans (Adriatic coast, southern Balkans), a third one expanded to recolonize Northern and Western Europe, loosing much of its diversity in the process. Our findings show that mobile and a priori homogeneous taxa may also display substructure within glacial refugia ('refugia within refugia') and emphasize the importance of the Balkans as a major European biodiversity centre. Moreover, the distribution of diversity roughly coincides with regional conservation situations, consistent with the idea that historically impoverished genetic diversity may interact with anthropogenic disturbances, and increase the vulnerability of populations. Phylogeographic models seem important to fully appreciate the risks of local declines and inform conservation strategies.

Evolutionary history of almond tree domestication in the Mediterranean basin.

Genetic diversity of contemporary domesticated species is shaped by both natural and human-driven processes. However, until now, little is known about how domestication has imprinted the variation of fruit tree species. In this study, we reconstruct the recent evolutionary history of the domesticated almond tree, Prunus dulcis, around the Mediterranean basin, using a combination of nuclear and chloroplast microsatellites [i.e. simple sequence repeat (SSRs)] to investigate patterns of genetic diversity. Whereas conservative chloroplast SSRs show a widespread haplotype and rare locally distributed variants, nuclear SSRs show a pattern of isolation by distance with clines of diversity from the East to the West of the Mediterranean basin, while Bayesian genetic clustering reveals a substantial longitudinal genetic structure. Both kinds of markers thus support a single domestication event, in the eastern side of the Mediterranean basin. In addition, model-based estimation of the timing of genetic divergence among those clusters is estimated sometime during the Holocene, a result that is compatible with human-mediated dispersal of almond tree out of its centre of origin. Still, the detection of region-specific alleles suggests that gene flow from relictual wild preglacial populations (in North Africa) or from wild counterparts (in the Near East) could account for a fraction of the diversity observed.

Hierarchical analyses of genetic differentiation in a hybrid zone of Sorex araneus (Insectivora : Soricidae)

Microsatellites are used to unravel the fine-scale genetic structure of a hybrid zone between chromosome races Valais and Cordon of the common shrew (Sorex araneus) located in the French Alps. A total of 269 individuals collected between 1992 and 1995 was typed for seven microsatellite loci. A modified version of the classical multiple correspondence analysis is carried out. This analysis clearly shows the dichotomy between the two races. Several approaches are used to study genetic structuring. Gene flow is clearly reduced between these chromosome races and is estimated at one migrant every two generations using X-statistics and one migrant per generation using F-statistics. Hierarchical F- and R-statistics are compared and their efficiency to detect inter- and intraracial patterns of divergence is discussed. Within-race genetic structuring is significant, but remains weak. F-ST displays similar values on both sides of the hybrid zone, although no environmental barriers are found on the Cordon side, whereas the Valais side is divided by several mountain rivers. We introduce the exact G-test to microsatellite data which proved to be a powerful test to detect genetic differentiation within as well as among races. The genetic background of karyotypic hybrids was compared with the genetic background of pure parental forms using a CRT-MCA. Our results indicate that, without knowledge of the karyotypes, we would not have been able to distinguish these hybrids from karyotypically pure samples.

Use of the frc gene as a molecular marker to characterize oxalate-oxidizing bacterial abundance and diversity structure in soil.

Oxalate catabolism, which can have both medical and environmental implications, is performed by phylogenetically diverse bacteria. The formyl-CoA-transferase gene was chosen as a molecular marker of the oxalotrophic function. Degenerated primers were deduced from an alignment of frc gene sequences available in databases. The specificity of primers was tested on a variety of frc-containing and frc-lacking bacteria. The frc-primers were then used to develop PCR-DGGE and real-time SybrGreen PCR assays in soils containing various amounts of oxalate. Some PCR products from pure cultures and from soil samples were cloned and sequenced. Data were used to generate a phylogenetic tree showing that environmental PCR products belonged to the target physiological group. The extent of diversity visualised on DGGE pattern was higher for soil samples containing carbonate resulting from oxalate catabolism. Moreover, the amount of frc gene copies in the investigated soils was detected in the range of 1.64x10(7) to 1.75x10(8)/g of dry soil under oxalogenic tree (representing 0.5 to 1.2% of total 16S rRNA gene copies), whereas the number of frc gene copies in the reference soil was 6.4x10(6) (or 0.2% of 16S rRNA gene copies). This indicates that oxalotrophic bacteria are numerous and widespread in soils and that a relationship exists between the presence of the oxalogenic trees Milicia excelsa and Afzelia africana and the relative abundance of oxalotrophic guilds in the total bacterial communities. This is obviously related to the accomplishment of the oxalate-carbonate pathway, which explains the alkalinization and calcium carbonate accumulation occurring below these trees in an otherwise acidic soil. The molecular tools developed in this study will allow in-depth understanding of the functional implication of these bacteria on carbonate accumulation as a way of atmospheric CO(2) sequestration.

Shifts in species and phylogenetic diversity between sapling and tree communities indicate negative density dependence in a lowland rain forest

1. As trees in a given cohort progress through ontogeny, many individuals die. This risk of mortality is unevenly distributed across species because of many processes such as habitat filtering, interspecific competition and negative density dependence. Here, we predict and test the patterns that such ecological processes should inscribe on both species and phylogenetic diversity as plants recruit from saplings to the canopy. 2. We compared species and phylogenetic diversity of sapling and tree communities at two sites in French Guiana. We surveyed 2084 adult trees in four 1-ha tree plots and 943 saplings in sixteen 16-m2 subplots nested within the tree plots. Species diversity was measured using Fisher's alpha (species richness) and Simpson's index (species evenness). Phylogenetic diversity was measured using Faith's phylogenetic diversity (phylogenetic richness) and Rao's quadratic entropy index (phylogenetic evenness). The phylogenetic diversity indices were inferred using four phylogenetic hypotheses: two based on rbcLa plastid DNA sequences obtained from the inventoried individuals with different branch lengths, a global phylogeny available from the Angiosperm Phylogeny Group, and a combination of both. 3. Taxonomic identification of the saplings was performed by combining morphological and DNA barcoding techniques using three plant DNA barcodes (psbA-trnH, rpoC1 and rbcLa). DNA barcoding enabled us to increase species assignment and to assign unidentified saplings to molecular operational taxonomic units. 4. Species richness was similar between saplings and trees, but in about half of our comparisons, species evenness was higher in trees than in saplings. This suggests that negative density dependence plays an important role during the sapling-to-tree transition. 5. Phylogenetic richness increased between saplings and trees in about half of the comparisons. Phylogenetic evenness increased significantly between saplings and trees in a few cases (4 out of 16) and only with the most resolved phylogeny. These results suggest that negative density dependence operates largely independently of the phylogenetic structure of communities. 6. Synthesis. By contrasting species richness and evenness across size classes, we suggest that negative density dependence drives shifts in composition during the sapling-to-tree transition. In addition, we found little evidence for a change in phylogenetic diversity across age classes, suggesting that the observed patterns are not phylogenetically constrained.

A simple permutation test for clusteredness

Hierarchical clustering is a popular method for finding structure in multivariate data,resulting in a binary tree constructed on the particular objects of the study, usually samplingunits. The user faces the decision where to cut the binary tree in order to determine the numberof clusters to interpret and there are various ad hoc rules for arriving at a decision. A simplepermutation test is presented that diagnoses whether non-random levels of clustering are presentin the set of objects and, if so, indicates the specific level at which the tree can be cut. The test isvalidated against random matrices to verify the type I error probability and a power study isperformed on data sets with known clusteredness to study the type II error.

Terrestrial arthropods from tree canopies in the Pantanal of Mato Grosso, Brazil

Terrestrial arthropods from tree canopies in the Pantanal of Mato Grosso, Brazil. This study represents a contribution to the knowledge of the diversity of arthropods associated to the canopy of Vochysia divergens Pohl (Vochysiaceae). Three trees individuals were sampled during two seasonal periods in this region: a) by spraying one tree canopy during high water (February); b) by fogging two tree canopies during low water (September/October). The 15,744 arthropods (183.2±38.9 individuals/m²) obtained from all three trees (86 m²) represented 20 taxonomic orders, 87.1% were Insecta, and 12.9% Arachnida. The dominant groups were Hymenoptera (48.5%; 88.9 individuals/m²), mostly Formicidae (44.5%; 81.4 individuals/m²), followed by Coleoptera (14.0%; 25.5 individuals/m²) and Araneae (10.2%; 19.5 individuals/m²), together representing 62.5% of the total catch. Fourteen (70%) of all orders occurred on three trees. Dermaptera, Isoptera, Neuroptera, Odonata, Plecoptera and Trichoptera were collected from only one tree. Of the total, 2,197 adult Coleoptera collected (25.5±11.3 individuals/m²), 99% were assigned to 32 families and 256 morphospecies. Nitidulidae (17.9% of the total catch; 4.6 individuals/m²), Anobiidae (16.7%; 4.3 individuals/m²), Curculionidae (13.2%; 3.4 individuals/m²) and Meloidae (11.4%; 2.9 individuals/m²) dominated. The communitiy of adult Coleoptera on V. divergens indicated a dominance of herbivores (37.8% of the total catch, 127 spp.) and predators (35.2%, 82 spp.), followed by saprophages (16.2%, 32 spp.) and fungivores (10.8%, 15 spp.). The influence of the flood pulse on the community of arboreal arthropods in V. divergens is indicated by the seasonal variation in evaluated groups, causing changes in their structure and composition.

Further insights on the French WISC-IV factor structure through Bayesian structural equation modeling

The interpretation of the Wechsler Intelligence Scale for Children-Fourth Edition (WISC-IV) is based on a 4-factor model, which is only partially compatible with the mainstream Cattell-Horn-Carroll (CHC) model of intelligence measurement. The structure of cognitive batteries is frequently analyzed via exploratory factor analysis and/or confirmatory factor analysis. With classical confirmatory factor analysis, almost all crossloadings between latent variables and measures are fixed to zero in order to allow the model to be identified. However, inappropriate zero cross-loadings can contribute to poor model fit, distorted factors, and biased factor correlations; most important, they do not necessarily faithfully reflect theory. To deal with these methodological and theoretical limitations, we used a new statistical approach, Bayesian structural equation modeling (BSEM), among a sample of 249 French-speaking Swiss children (8-12 years). With BSEM, zero-fixed cross-loadings between latent variables and measures are replaced by approximate zeros, based on informative, small-variance priors. Results indicated that a direct hierarchical CHC-based model with 5 factors plus a general intelligence factor better represented the structure of the WISC-IV than did the 4-factor structure and the higher order models. Because a direct hierarchical CHC model was more adequate, it was concluded that the general factor should be considered as a breadth rather than a superordinate factor. Because it was possible for us to estimate the influence of each of the latent variables on the 15 subtest scores, BSEM allowed improvement of the understanding of the structure of intelligence tests and the clinical interpretation of the subtest scores.

HIERFSTAT, a package for R to compute and test hierarchical F-statistics

The package HIERFSTAT for the statistical software R, created by the R Development Core Team, allows the estimate of hierarchical F-statistics from a hierarchy with any numbers of levels. In addition, it allows testing the statistical significance of population differentiation for these different levels, using a generalized likelihood-ratio test. The package HIERFSTAT is available at http://www.unil.ch/popgen/softwares/hierfstat.htm.

Structure-function characterization and optimization of a plant-derived antibacterial peptide.

Crushed seeds of the Moringa oleifera tree have been used traditionally as natural flocculants to clarify drinking water. We previously showed that one of the seed peptides mediates both the sedimentation of suspended particles such as bacterial cells and a direct bactericidal activity, raising the possibility that the two activities might be related. In this study, the conformational modeling of the peptide was coupled to a functional analysis of synthetic derivatives. This indicated that partly overlapping structural determinants mediate the sedimentation and antibacterial activities. Sedimentation requires a positively charged, glutamine-rich portion of the peptide that aggregates bacterial cells. The bactericidal activity was localized to a sequence prone to form a helix-loop-helix structural motif. Amino acid substitution showed that the bactericidal activity requires hydrophobic proline residues within the protruding loop. Vital dye staining indicated that treatment with peptides containing this motif results in bacterial membrane damage. Assembly of multiple copies of this structural motif into a branched peptide enhanced antibacterial activity, since low concentrations effectively kill bacteria such as Pseudomonas aeruginosa and Streptococcus pyogenes without displaying a toxic effect on human red blood cells. This study thus identifies a synthetic peptide with potent antibacterial activity against specific human pathogens. It also suggests partly distinct molecular mechanisms for each activity. Sedimentation may result from coupled flocculation and coagulation effects, while the bactericidal activity would require bacterial membrane destabilization by a hydrophobic loop.

Genetic differentiation for size at first reproduction through male versus female functions in the widespread Mediterranean tree Pinus pinaster.

BACKGROUND AND AIMS: The study of local adaptation in plant reproductive traits has received substantial attention in short-lived species, but studies conducted on forest trees are scarce. This lack of research on long-lived species represents an important gap in our knowledge, because inferences about selection on the reproduction and life history of short-lived species cannot necessarily be extrapolated to trees. This study considers whether the size for first reproduction is locally adapted across a broad geographical range of the Mediterranean conifer species Pinus pinaster. In particular, the study investigates whether this monoecious species varies genetically among populations in terms of whether individuals start to reproduce through their male function, their female function or both sexual functions simultaneously. Whether differences among populations could be attributed to local adaptation across a climatic gradient is then considered. METHODS: Male and female reproduction and growth were measured during early stages of sexual maturity of a P. pinaster common garden comprising 23 populations sampled across the species range. Generalized linear mixed models were used to assess genetic variability of early reproductive life-history traits. Environmental correlations with reproductive life-history traits were tested after controlling for neutral genetic structure provided by 12 nuclear simple sequence repeat markers. KEY RESULTS: Trees tended to reproduce first through their male function, at a size (height) that varied little among source populations. The transition to female reproduction was slower, showed higher levels of variability and was negatively correlated with vegetative growth traits. Several female reproductive traits were correlated with a gradient of growth conditions, even after accounting for neutral genetic structure, with populations from more unfavourable sites tending to commence female reproduction at a lower individual size. CONCLUSIONS: The study represents the first report of genetic variability among populations for differences in the threshold size for first reproduction between male and female sexual functions in a tree species. The relatively uniform size at which individuals begin reproducing through their male function probably represents the fact that pollen dispersal is also relatively invariant among sites. However, the genetic variability in the timing of female reproduction probably reflects environment-dependent costs of cone production. The results also suggest that early sex allocation in this species might evolve under constraints that do not apply to other conifers.

BIO-INSPIRED COMPUTATIONAL TECHNIQUES APPLIED TO THE CLUSTERING AND VISUALIZATION OF SPATIO-TEMPORAL GEOSPATIAL DATA

The coverage and volume of geo-referenced datasets are extensive and incessantly¦growing. The systematic capture of geo-referenced information generates large volumes¦of spatio-temporal data to be analyzed. Clustering and visualization play a key¦role in the exploratory data analysis and the extraction of knowledge embedded in¦these data. However, new challenges in visualization and clustering are posed when¦dealing with the special characteristics of this data. For instance, its complex structures,¦large quantity of samples, variables involved in a temporal context, high dimensionality¦and large variability in cluster shapes.¦The central aim of my thesis is to propose new algorithms and methodologies for¦clustering and visualization, in order to assist the knowledge extraction from spatiotemporal¦geo-referenced data, thus improving making decision processes.¦I present two original algorithms, one for clustering: the Fuzzy Growing Hierarchical¦Self-Organizing Networks (FGHSON), and the second for exploratory visual data analysis:¦the Tree-structured Self-organizing Maps Component Planes. In addition, I present¦methodologies that combined with FGHSON and the Tree-structured SOM Component¦Planes allow the integration of space and time seamlessly and simultaneously in¦order to extract knowledge embedded in a temporal context.¦The originality of the FGHSON lies in its capability to reflect the underlying structure¦of a dataset in a hierarchical fuzzy way. A hierarchical fuzzy representation of¦clusters is crucial when data include complex structures with large variability of cluster¦shapes, variances, densities and number of clusters. The most important characteristics¦of the FGHSON include: (1) It does not require an a-priori setup of the number¦of clusters. (2) The algorithm executes several self-organizing processes in parallel.¦Hence, when dealing with large datasets the processes can be distributed reducing the¦computational cost. (3) Only three parameters are necessary to set up the algorithm.¦In the case of the Tree-structured SOM Component Planes, the novelty of this algorithm¦lies in its ability to create a structure that allows the visual exploratory data analysis¦of large high-dimensional datasets. This algorithm creates a hierarchical structure¦of Self-Organizing Map Component Planes, arranging similar variables' projections in¦the same branches of the tree. Hence, similarities on variables' behavior can be easily¦detected (e.g. local correlations, maximal and minimal values and outliers).¦Both FGHSON and the Tree-structured SOM Component Planes were applied in¦several agroecological problems proving to be very efficient in the exploratory analysis¦and clustering of spatio-temporal datasets.¦In this thesis I also tested three soft competitive learning algorithms. Two of them¦well-known non supervised soft competitive algorithms, namely the Self-Organizing¦Maps (SOMs) and the Growing Hierarchical Self-Organizing Maps (GHSOMs); and the¦third was our original contribution, the FGHSON. Although the algorithms presented¦here have been used in several areas, to my knowledge there is not any work applying¦and comparing the performance of those techniques when dealing with spatiotemporal¦geospatial data, as it is presented in this thesis.¦I propose original methodologies to explore spatio-temporal geo-referenced datasets¦through time. Our approach uses time windows to capture temporal similarities and¦variations by using the FGHSON clustering algorithm. The developed methodologies¦are used in two case studies. In the first, the objective was to find similar agroecozones¦through time and in the second one it was to find similar environmental patterns¦shifted in time.¦Several results presented in this thesis have led to new contributions to agroecological¦knowledge, for instance, in sugar cane, and blackberry production.¦Finally, in the framework of this thesis we developed several software tools: (1)¦a Matlab toolbox that implements the FGHSON algorithm, and (2) a program called¦BIS (Bio-inspired Identification of Similar agroecozones) an interactive graphical user¦interface tool which integrates the FGHSON algorithm with Google Earth in order to¦show zones with similar agroecological characteristics.

Impact of community structure on information transfer

The observation that real complex networks have internal structure has important implication for dynamic processes occurring on such topologies. Here we investigate the impact of community structure on a model of information transfer able to deal with both search and congestion simultaneously. We show that networks with fuzzy community structure are more efficient in terms of packet delivery than those with pronounced community structure. We also propose an alternative packet routing algorithm which takes advantage of the knowledge of communities to improve information transfer and show that in the context of the model an intermediate level of community structure is optimal. Finally, we show that in a hierarchical network setting, providing knowledge of communities at the level of highest modularity will improve network capacity by the largest amount.

Inferring epidemic contact structure from phylogenetic trees.

Contact structure is believed to have a large impact on epidemic spreading and consequently using networks to model such contact structure continues to gain interest in epidemiology. However, detailed knowledge of the exact contact structure underlying real epidemics is limited. Here we address the question whether the structure of the contact network leaves a detectable genetic fingerprint in the pathogen population. To this end we compare phylogenies generated by disease outbreaks in simulated populations with different types of contact networks. We find that the shape of these phylogenies strongly depends on contact structure. In particular, measures of tree imbalance allow us to quantify to what extent the contact structure underlying an epidemic deviates from a null model contact network and illustrate this in the case of random mixing. Using a phylogeny from the Swiss HIV epidemic, we show that this epidemic has a significantly more unbalanced tree than would be expected from random mixing.