Effect of training-sample size and classification difficulty on the accuracy of genomic predictors.

Introduction: As part of the MicroArray Quality Control (MAQC)-II project, this analysis examines how the choice of univariate feature-selection methods and classification algorithms may influence the performance of genomic predictors under varying degrees of prediction difficulty represented by three clinically relevant endpoints. Methods: We used gene-expression data from 230 breast cancers (grouped into training and independent validation sets), and we examined 40 predictors (five univariate feature-selection methods combined with eight different classifiers) for each of the three endpoints. Their classification performance was estimated on the training set by using two different resampling methods and compared with the accuracy observed in the independent validation set. Results: A ranking of the three classification problems was obtained, and the performance of 120 models was estimated and assessed on an independent validation set. The bootstrapping estimates were closer to the validation performance than were the cross-validation estimates. The required sample size for each endpoint was estimated, and both gene-level and pathway-level analyses were performed on the obtained models. Conclusions: We showed that genomic predictor accuracy is determined largely by an interplay between sample size and classification difficulty. Variations on univariate feature-selection methods and choice of classification algorithm have only a modest impact on predictor performance, and several statistically equally good predictors can be developed for any given classification problem.

Alineamiento múltiple de secuencias con T-Coffee: una aproximación paralela

Las aplicaciones de alineamiento múltiple de secuencias son prototipos de aplicaciones que requieren elevada potencia de cómputo y memoria. Se destacan por la relevancia científica que tienen los resultados que brindan a investigaciones científicas en el campo de la biomedicina, genética y farmacología. Las aplicaciones de alineamiento múltiple tienen la limitante de que no son capaces de procesar miles de secuencias, por lo que se hace necesario crear un modelo para resolver la problemática. Analizando el volumen de datos que se manipulan en el área de las ciencias biológica y la complejidad de los algoritmos de alineamiento de secuencias, la única vía de solución del problema es a través de la utilización de entornos de cómputo paralelos y la computación de altas prestaciones. La investigación realizada por nosotros tiene como objetivo la creación de un modelo paralelo que le permita a los algoritmos de alineamiento múltiple aumentar el número de secuencias a procesar, tratando de mantener la calidad en los resultados para garantizar la precisión científica. El modelo que proponemos emplea como base la clusterización de las secuencias de entrada utilizando criterios biológicos que permiten mantener la calidad de los resultados. Además, el modelo se enfoca en la disminución del tiempo de cómputo y consumo de memoria. Para presentar y validar el modelo utilizamos T-Coffee, como plataforma de desarrollo e investigación. El modelo propuesto pudiera ser aplicado a cualquier otro algoritmo de alineamiento múltiple de secuencias.

Versatile gene-specific sequence tags for Arabidopsis functional genomics: transcript profiling and reverse genetics applications.

Microarray transcript profiling and RNA interference are two new technologies crucial for large-scale gene function studies in multicellular eukaryotes. Both rely on sequence-specific hybridization between complementary nucleic acid strands, inciting us to create a collection of gene-specific sequence tags (GSTs) representing at least 21,500 Arabidopsis genes and which are compatible with both approaches. The GSTs were carefully selected to ensure that each of them shared no significant similarity with any other region in the Arabidopsis genome. They were synthesized by PCR amplification from genomic DNA. Spotted microarrays fabricated from the GSTs show good dynamic range, specificity, and sensitivity in transcript profiling experiments. The GSTs have also been transferred to bacterial plasmid vectors via recombinational cloning protocols. These cloned GSTs constitute the ideal starting point for a variety of functional approaches, including reverse genetics. We have subcloned GSTs on a large scale into vectors designed for gene silencing in plant cells. We show that in planta expression of GST hairpin RNA results in the expected phenotypes in silenced Arabidopsis lines. These versatile GST resources provide novel and powerful tools for functional genomics.

The effect of host social system on parasite population genetic structure: comparative population genetics of two ectoparasitic mites and their bat hosts.

BACKGROUND: The population genetic structure of a parasite, and consequently its ability to adapt to a given host, is strongly linked to its own life history as well as the life history of its host. While the effects of parasite life history on their population genetic structure have received some attention, the effect of host social system has remained largely unstudied. In this study, we investigated the population genetic structure of two closely related parasitic mite species (Spinturnix myoti and Spinturnix bechsteini) with very similar life histories. Their respective hosts, the greater mouse-eared bat (Myotis myotis) and the Bechstein's bat (Myotis bechsteinii) have social systems that differ in several substantial features, such as group size, mating system and dispersal patterns. RESULTS: We found that the two mite species have strongly differing population genetic structures. In S. myoti we found high levels of genetic diversity and very little pairwise differentiation, whereas in S. bechsteini we observed much less diversity, strongly differentiated populations and strong temporal turnover. These differences are likely to be the result of the differences in genetic drift and dispersal opportunities afforded to the two parasites by the different social systems of their hosts. CONCLUSIONS: Our results suggest that host social system can strongly influence parasite population structure. As a result, the evolutionary potential of these two parasites with very similar life histories also differs, thereby affecting the risk and evolutionary pressure exerted by each parasite on its host.

A Survey of Active Learning Algorithms for Supervised Remote Sensing Image Classification

Defining an efficient training set is one of the most delicate phases for the success of remote sensing image classification routines. The complexity of the problem, the limited temporal and financial resources, as well as the high intraclass variance can make an algorithm fail if it is trained with a suboptimal dataset. Active learning aims at building efficient training sets by iteratively improving the model performance through sampling. A user-defined heuristic ranks the unlabeled pixels according to a function of the uncertainty of their class membership and then the user is asked to provide labels for the most uncertain pixels. This paper reviews and tests the main families of active learning algorithms: committee, large margin, and posterior probability-based. For each of them, the most recent advances in the remote sensing community are discussed and some heuristics are detailed and tested. Several challenging remote sensing scenarios are considered, including very high spatial resolution and hyperspectral image classification. Finally, guidelines for choosing the good architecture are provided for new and/or unexperienced user.

Quantitative genetics of sleep in inbred mice.

The timing and the organization of sleep architecture are mainly controlled by the circadian system, while sleep need and intensity are regulated by a homeostatic process. How independent these two systems are in regulating sleep is not well understood. In contrast to the impressive progress in the molecular genetics of circadian rhythms, little is known about the molecular basis of sleep. Nevertheless, as summarized here, phenotypic dissection of sleep into its most basic aspects can be used to identify both the single major genes and small effect quantitative trait loci involved. Although experimental models such as the mouse are more readily amenable to genetic analysis of sleep, similar approaches can be applied to humans.

SYSGENET: a meeting report from a new European network for systems genetics.

The first scientific meeting of the newly established European SYSGENET network took place at the Helmholtz Centre for Infection Research (HZI) in Braunschweig, April 7-9, 2010. About 50 researchers working in the field of systems genetics using mouse genetic reference populations (GRP) participated in the meeting and exchanged their results, phenotyping approaches, and data analysis tools for studying systems genetics. In addition, the future of GRP resources and phenotyping in Europe was discussed.

GWAS of human bitter taste perception identifies new loci and reveals additional complexity of bitter taste genetics.

Human perception of bitterness displays pronounced interindividual variation. This phenotypic variation is mirrored by equally pronounced genetic variation in the family of bitter taste receptor genes. To better understand the effects of common genetic variations on human bitter taste perception, we conducted a genome-wide association study on a discovery panel of 504 subjects and a validation panel of 104 subjects from the general population of São Paulo in Brazil. Correction for general taste-sensitivity allowed us to identify a SNP in the cluster of bitter taste receptors on chr12 (10.88- 11.24 Mb, build 36.1) significantly associated (best SNP: rs2708377, P = 5.31 × 10(-13), r(2) = 8.9%, β = -0.12, s.e. = 0.016) with the perceived bitterness of caffeine. This association overlaps with-but is statistically distinct from-the previously identified SNP rs10772420 influencing the perception of quinine bitterness that falls in the same bitter taste cluster. We replicated this association to quinine perception (P = 4.97 × 10(-37), r(2) = 23.2%, β = 0.25, s.e. = 0.020) and additionally found the effect of this genetic locus to be concentration specific with a strong impact on the perception of low, but no impact on the perception of high concentrations of quinine. Our study, thus, furthers our understanding of the complex genetic architecture of bitter taste perception.

Landscape genetics of the Alpine newt (Mesotriton alpestris) inferred from a strip-based approach

Habitat destruction and fragmentation are known to strongly affect dispersal by altering the quality of the environment between populations. As a consequence, lower landscape connectivity is expected to enhance extinction risks through a decrease in gene flow and the resulting negative effects of genetic drift, accumulation of deleterious mutations and inbreeding depression. Such phenomena are particularly harmful for amphibian species, characterized by disjunct breeding habitats. The dispersal behaviour of amphibians being poorly understood, it is crucial to develop new tools, allowing us to determine the influence of landscape connectivity on the persistence of populations. In this study, we developed a new landscape genetics approach that aims at identifying land-uses affecting genetic differentiation, without a priori assumptions about associated ecological costs. We surveyed genetic variation at seven microsatellite loci for 19 Alpine newt (Mesotriton alpestris) populations in western Switzerland. Using strips of varying widths that define a dispersal corridor between pairs of populations, we were able to identify land-uses that act as dispersal barriers (i.e. urban areas) and corridors (i.e. forests). Our results suggest that habitat destruction and landscape fragmentation might in the near future affect common species such as M. alpestris. In addition, by identifying relevant landscape variables influencing population structure without unrealistic assumptions about dispersal, our method offers a simple and flexible tool of investigation as an alternative to least-cost models and other approaches.

Structural prediction of peptides bound to MHC class I.

An ab initio structure prediction approach adapted to the peptide-major histocompatibility complex (MHC) class I system is presented. Based on structure comparisons of a large set of peptide-MHC class I complexes, a molecular dynamics protocol is proposed using simulated annealing (SA) cycles to sample the conformational space of the peptide in its fixed MHC environment. A set of 14 peptide-human leukocyte antigen (HLA) A0201 and 27 peptide-non-HLA A0201 complexes for which X-ray structures are available is used to test the accuracy of the prediction method. For each complex, 1000 peptide conformers are obtained from the SA sampling. A graph theory clustering algorithm based on heavy atom root-mean-square deviation (RMSD) values is applied to the sampled conformers. The clusters are ranked using cluster size, mean effective or conformational free energies, with solvation free energies computed using Generalized Born MV 2 (GB-MV2) and Poisson-Boltzmann (PB) continuum models. The final conformation is chosen as the center of the best-ranked cluster. With conformational free energies, the overall prediction success is 83% using a 1.00 Angstroms crystal RMSD criterion for main-chain atoms, and 76% using a 1.50 Angstroms RMSD criterion for heavy atoms. The prediction success is even higher for the set of 14 peptide-HLA A0201 complexes: 100% of the peptides have main-chain RMSD values < or =1.00 Angstroms and 93% of the peptides have heavy atom RMSD values < or =1.50 Angstroms. This structure prediction method can be applied to complexes of natural or modified antigenic peptides in their MHC environment with the aim to perform rational structure-based optimizations of tumor vaccines.

Conservation genetics and population monitoring of the Alpine Salamander (Salamandra Atra) :

SUMMARYIn the context of the biodiversity crisis, amphibians are experiencing the most severe worldwide decline of all vertebrates and are in urgent need of better management. Efficient conservation strategies rely on sound knowledge of the species biology and of the genetic and demographic processes that might impair their welfare. Nonetheless, these processes are poorly understood in amphibians. Delineating population boundaries remains consequently problematic for these species, while it is of critical importance to define adequate management units for conservation. In this study, our attention focused on the alpine salamander (Salamandra atra), a species that deserves much interest in terms of both conservation biology and evolution. This endemic alpine species shows peculiar life-history traits (viviparity, reduced activity period, slow maturation) and has a slow population turnover, which might be problematic for its persistence in a changing environment. Due to its elusive behaviour (individuals spend most of their time underground and are unavailable for sampling), dynamic processes of gene and individuals were poorly understood for that species. Consequently, its conservation status could hardly be reliably assessed. Similarly the fire salamander (Salamandra salamandra) also poses special challenges for conservation, as no clear demarcation of geographical populations exists and dispersal patterns are poorly known. Through a phylogeographic analysis, we first studied the evolutionary history of the alpine salamander to better document the distribution of the genetic diversity along its geographical range. This study highlighted the presence of multiple divergent lineages in Italy together with a clear genetic divergence between populations from Northern and Dinaric Alps. These signs of cryptic genetic differentiation, which are not accounted for by the current taxonomy of the species, should not be neglected for further definition of conservation units. In addition, our data supported glacial survival of the species in northern peripheral glacial réfugia and nunataks, a pattern rarely documented for long-lived species. Then, we evaluated the level of gene flow between populations at the local scale and tested for asymmetries in male versus female dispersal using both field-based (mark-recapture) and genetic approaches. This study revealed high level of gene flow between populations, which stems mainly from male dispersal. This corroborated the idea that salamanders are much better dispersers than hitherto thought and provided a well- supported example of male-biased dispersal in amphibians. In a third step, based on a mark- recapture survey, we addressed the problem of sampling unavailability in alpine salamanders and evaluated its impact on two monitoring methods. We showed that about three quarters of individuals were unavailable for sampling during sampling sessions, a proportion that can vary with climatic conditions. If not taken into account, these complexities would result in false assumptions on population trends and misdirect conservation efforts. Finally, regarding the daunting task of delineating management units, our attention was drawn on the fire salamander. We conducted a local population genetic study that revealed high levels of gene flow among sampling sites. Management units for this species should consequently be large. Interestingly, despite the presence of several landscape features often reported to act as barriers, genetic breaks occurred at unexpected places. This suggests that landscape features may rather have idiosyncratic effects on population structure. In conclusion, this work brought new insights on both genetic and demographic processes occurring in salamanders. The results suggest that some biological paradigms should be taken with caution when particular species are in focus. Species- specific studies remain thus fundamental for a better understanding of species evolution and conservation, particularly in the context of current global changes.RESUMEDans le contexte de la crise de la biodiversité actuelle, les amphibiens subissent le déclin le plus important de tous les vertébrés et ont urgemment besoin d'une meilleure protection. L'établissement de stratégies de conservation efficaces repose sur des connaissances solides de la biologie des espèces et des processus génétiques et démographiques pouvant menacer leur survie. Ces processus sont néanmoins encore peu étudiés chez les amphibiens.Dans cette étude, notre attention s'est portée sur la salamandre noire (Salamandra atra), une espèce endémique des Alpes dont les traits d'histoire de vie atypiques (viviparité, phase d'activité réduite, lent turnover des populations) pourraient la rendre très vulnérable face aux changements environnementaux. Par ailleurs, en raison de son comportement cryptique (les individus passent la plupart de leur temps sous terre) la dynamique des gènes et des individus est mal comprise chez cette espèce. Il est donc difficile d'évaluer son statut de conservation de manière fiable. La salamandre tachetée {Salamandra salamandra), pour qui il n'existe aucune démarcation géographique apparente des populations, pose également des problèmes en termes de gestion. Dans un premier temps, nous avons étudié l'histoire évolutive de la salamandre noire afin de mieux décrire la distribution de sa diversité génétique au sein de son aire géographique. Cela a permis de mettre en évidence la présence de multiples lignées en Italie, ainsi qu'une nette divergence entre les populations du nord des Alpes et des Alpes dinariques. Ces résultats seront à prendre en compte lorsqu'il s'agira de définir des unités de conservation pour cette espèce. D'autre part, nos données soutiennent l'hypothèse d'une survie glaciaire dans des refuges nordiques périglaciaires ou dans des nunataks, fait rarement documenté pour une espèce longévive. Nous avons ensuite évalué la différentiation génétique des populations à l'échelle locale, ce qui a révélé d'important flux de gènes, ainsi qu'une asymétrie de dispersion en faveur des mâles. Ces résultats corroborent l'idée que les amphibiens dispersent mieux que ce que l'on pensait, et fournissent un exemple robuste de dispersion biaisée en faveur des mâles chez les amphibiens. Nous avons ensuite abordé le problème de Γ inaccessibilité des individus à la capture. Nous avons montré qu'environ trois quarts des individus sont inaccessibles lors des échantillonnages, une proportion qui peut varier en fonction des conditions climatiques. Ignoré, ce processus pourrait entraîner une mauvaise interprétation des fluctuations de populations ainsi qu'une mauvaise allocation des efforts de conservation. Concernant la définition d'unités de gestion pour la salamandre tachetée, nous avons pu mettre en évidence un flux de gènes important entre les sites échantillonnés. Les unités de gestion pour cette espèce devraient donc être étendues. Etonnamment, malgré la présence de nombreuses barrières potentielles au flux de gènes, les démarcations génétiques sont apparues à des endroits inattendus. En conclusion, ce travail a apporté une meilleure compréhension des processus génétiques et démographiques en action chez les salamandres. Les résultats suggèrent que certains paradigmes biologiques devraient être considérés avec précaution quand il s'agit de les appliquer à des espèces particulières. Les études spécifiques demeurent donc fondamentales pour une meilleure compréhension de l'évolution des espèces et leur conservation, tout particulièrement dans le contexte des changements globaux actuels.