The situation of open access institutional repositories in Spain: 2009 report

Introduction. The DRIVER I project drew up a detailed report of European repositories based on data gathered in a survey in which Spain's participation was very low. This created a highly distorted image of the implementation of repositories in Spain. This study aims to analyse the current state of Spanish open-access institutional repositories and to describe their characteristics. Method. The data were gathered through a Web survey. The questionnaire was based on that used by DRIVER I: coverage; technical infrastructure and technical issues; institutional policies; services created; and stimulators and inhibitors for establishing, filling and maintaining their digital institutional repositories. Analysis. Data were tabulated and analysed systematically according responses obtained from the questionnaire and grouped by coverage. Results. Responses were obtained from 38 of the 104 institutions contacted, which had 29 institutional repositories. This represents 78.3% of the Spanish repositories according to the BuscaRepositorios directory. Spanish repositories contained mainly full-text materials (journal articles and doctoral theses) together with metadata. The software most used was DSpace, followed by EPrints. The metadata standard most used was Dublin Core. Spanish repositories offered more usage statistics and fewer author-oriented services than the European average. The priorities for the future development of the repositories are the need for clear policies on access to scientific production based on public funding and the need for quality control indicators. Conclusions.This is the first detailed study of Spanish institutional repositories. The key stimulants for establishing, filling and maintaining were, in order of importance, the increase of visibility and citation, the interest of decision-makers, simplicity of use and search services. On the other hand the main inhibitors identified were the absence of policies, the lack of integration with other national and international systems and the lack of awareness efforts among academia.

Metallogenic model of the Trepca Pb-Zn-Ag skarn deposit, Kosovo: Evidence from fluid inclusions, Rare Earth Elements, and stable isotope data

The Trepca Pb-Zn-Ag skarn deposit (29 Mt of ore at 3.45% Pb, 2.30% Zn, and 80 g/t Ag) is located in the Kopaonik block of the western Vardar zone, Kosovo. The mineralization, hosted by recrystallized limestone of Upper Triassic age, was structurally and lithologically controlled. Ore deposition is spatially and temporally related with the postcollisional magmatism of Oligocene age (23-26 Ma). The deposit was formed during two distinct mineralization stages: an early prograde closed-system and a later retrograde open-system stage. The prograde mineralization consisting mainly of pyroxenes (Hd(54-100)Jo(0-45)Di(0-45)) resulted from the interaction of magmatic fluids associated with Oligocene (23-26 Ma) postcollisional magmatism. Whereas there is no direct contact between magmatic rocks and the mineralization, the deposit is classified as a distal Pb-Zn-Ag skarn. Abundant pyroxene reflects low oxygen fugacity (<10(-31) bar) and anhydrous environment. Fluid inclusion data and mineral assemblage limit the prograde stage within a temperature range between 390 degrees and 475 degrees C. Formation pressure is estimated below 900 bars. Isotopic composition of aqueous fluid, inclusions hosted by hedenbergite (delta D = -108 to -130 parts per thousand; delta O-18 = 7.5-8.0 parts per thousand), Mn-enriched mineralogy and high REE content of the host carbonates at the contact with the skarn mineralization suggest that a magmatic fluid was modified during its infiltration through the country rocks. The retrograde mineral assemblage comprises ilvaite, magnetite, arsenopyrite, pyrrhotite, marcasite, pyrite, quartz, and various carbonates. Increases in oxygen and sulfur fugacities, as well as a hydrous character of mineralization, require an open-system model. The opening of the system is related to phreatomagmatic explosion and formation of the breccia. Arsenopyrite geothermometer limits the retrograde stage within the temperature range between 350 degrees and 380 degrees C and sulfur fugacity between 10(-8.8) and 10(-7.2) bars. The principal ore minerals, galena, sphalerite, pyrite, and minor chalcopyrite, were deposited from a moderately saline Ca-Na chloride fluid at around 350 degrees C. According to the isotopic composition of fluid inclusions hosted by sphalerite (delta D = -55 to -74 parts per thousand; delta O-18 = -9.6 to -13.6 parts per thousand), the fluid responsible for ore deposition was dominantly meteoric in origin. The delta S-31 values of the sulfides spanning between -5.5 and +10 parts per thousand point to a magmatic origin of sulfur. Ore deposition appears to have been largely contemporaneous with the retrograde stage of the skarn development. Postore stage accompanied the precipitation of significant amount of carbonates including the travertine deposits at the deposit surface. Mineralogical composition of travertine varies from calcite to siderite and all carbonates contain significant amounts of Mn. Decreased formation temperature and depletion in the REE content point to an influence of pH-neutralized cold ground water and dying magmatic system.

Cilengitide combined with standard treatment for patients with newly diagnosed glioblastoma with methylated MGMT promoter (CENTRIC EORTC 26071-22072 study): a multicentre, randomised, open-label, phase 3 trial.

BACKGROUND: Cilengitide is a selective αvβ3 and αvβ5 integrin inhibitor. Data from phase 2 trials suggest that it has antitumour activity as a single agent in recurrent glioblastoma and in combination with standard temozolomide chemoradiotherapy in newly diagnosed glioblastoma (particularly in tumours with methylated MGMT promoter). We aimed to assess cilengitide combined with temozolomide chemoradiotherapy in patients with newly diagnosed glioblastoma with methylated MGMT promoter. METHODS: In this multicentre, open-label, phase 3 study, we investigated the efficacy of cilengitide in patients from 146 study sites in 25 countries. Eligible patients (newly diagnosed, histologically proven supratentorial glioblastoma, methylated MGMT promoter, and age ≥18 years) were stratified for prognostic Radiation Therapy Oncology Group recursive partitioning analysis class and geographic region and centrally randomised in a 1:1 ratio with interactive voice response system to receive temozolomide chemoradiotherapy with cilengitide 2000 mg intravenously twice weekly (cilengitide group) or temozolomide chemoradiotherapy alone (control group). Patients and investigators were unmasked to treatment allocation. Maintenance temozolomide was given for up to six cycles, and cilengitide was given for up to 18 months or until disease progression or unacceptable toxic effects. The primary endpoint was overall survival. We analysed survival outcomes by intention to treat. This study is registered with ClinicalTrials.gov, number NCT00689221. FINDINGS: Overall, 3471 patients were screened. Of these patients, 3060 had tumour MGMT status tested; 926 patients had a methylated MGMT promoter, and 545 were randomly assigned to the cilengitide (n=272) or control groups (n=273) between Oct 31, 2008, and May 12, 2011. Median overall survival was 26·3 months (95% CI 23·8-28·8) in the cilengitide group and 26·3 months (23·9-34·7) in the control group (hazard ratio 1·02, 95% CI 0·81-1·29, p=0·86). None of the predefined clinical subgroups showed a benefit from cilengitide. We noted no overall additional toxic effects with cilengitide treatment. The most commonly reported adverse events of grade 3 or worse in the safety population were lymphopenia (31 [12%] in the cilengitide group vs 26 [10%] in the control group), thrombocytopenia (28 [11%] vs 46 [18%]), neutropenia (19 [7%] vs 24 [9%]), leucopenia (18 [7%] vs 20 [8%]), and convulsion (14 [5%] vs 15 [6%]). INTERPRETATION: The addition of cilengitide to temozolomide chemoradiotherapy did not improve outcomes; cilengitide will not be further developed as an anticancer drug. Nevertheless, integrins remain a potential treatment target for glioblastoma. FUNDING: Merck KGaA, Darmstadt, Germany.

Endovascular aortic repair versus open surgical repair for descending thoracic aortic disease a systematic review and meta-analysis of comparative studies.

OBJECTIVES: The purpose of this study was to determine whether thoracic endovascular aortic repair (TEVAR) reduces death and morbidity compared with open surgical repair for descending thoracic aortic disease. BACKGROUND: The role of TEVAR versus open surgery remains unclear. Metaregression can be used to maximally inform adoption of new technologies by utilizing evidence from existing trials. METHODS: Data from comparative studies of TEVAR versus open repair of the descending aorta were combined through meta-analysis. Metaregression was performed to account for baseline risk factor imbalances, study design, and thoracic pathology. Due to significant heterogeneity, registry data were analyzed separately from comparative studies. RESULTS: Forty-two nonrandomized studies involving 5,888 patients were included (38 comparative studies, 4 registries). Patient characteristics were balanced except for age, as TEVAR patients were usually older than open surgery patients (p = 0.001). Registry data suggested overall perioperative complications were reduced. In comparative studies, all-cause mortality at 30 days (odds ratio [OR]: 0.44, 95% confidence interval [CI]: 0.33 to 0.59) and paraplegia (OR: 0.42, 95% CI: 0.28 to 0.63) were reduced for TEVAR versus open surgery. In addition, cardiac complications, transfusions, reoperation for bleeding, renal dysfunction, pneumonia, and length of stay were reduced. There was no significant difference in stroke, myocardial infarction, aortic reintervention, and mortality beyond 1 year. Metaregression to adjust for age imbalance, study design, and pathology did not materially change the results. CONCLUSIONS: Current data from nonrandomized studies suggest that TEVAR may reduce early death, paraplegia, renal insufficiency, transfusions, reoperation for bleeding, cardiac complications, pneumonia, and length of stay compared with open surgery. Sustained benefits on survival have not been proven.

MADAP, a flexible clustering tool for the interpretation of one-dimensional genome annotation data.

A recurring task in the analysis of mass genome annotation data from high-throughput technologies is the identification of peaks or clusters in a noisy signal profile. Examples of such applications are the definition of promoters on the basis of transcription start site profiles, the mapping of transcription factor binding sites based on ChIP-chip data and the identification of quantitative trait loci (QTL) from whole genome SNP profiles. Input to such an analysis is a set of genome coordinates associated with counts or intensities. The output consists of a discrete number of peaks with respective volumes, extensions and center positions. We have developed for this purpose a flexible one-dimensional clustering tool, called MADAP, which we make available as a web server and as standalone program. A set of parameters enables the user to customize the procedure to a specific problem. The web server, which returns results in textual and graphical form, is useful for small to medium-scale applications, as well as for evaluation and parameter tuning in view of large-scale applications, requiring a local installation. The program written in C++ can be freely downloaded from ftp://ftp.epd.unil.ch/pub/software/unix/madap. The MADAP web server can be accessed at http://www.isrec.isb-sib.ch/madap/.

Smed454 dataset: unravelling the transcriptome of Schmidtea mediterranea

Background: Freshwater planarians are an attractive model for regeneration and stem cell research and have become a promising tool in the field of regenerative medicine. With the availability of a sequenced planarian genome, the recent application of modern genetic and high-throughput tools has resulted in revitalized interest in these animals, long known for their amazing regenerative capabilities, which enable them to regrow even a new head after decapitation. However, a detailed description of the planarian transcriptome is essential for future investigation into regenerative processes using planarians as a model system. Results: In order to complement and improve existing gene annotations, we used a 454 pyrosequencing approach to analyze the transcriptome of the planarian species Schmidtea mediterranea Altogether, 598,435 454-sequencing reads, with an average length of 327 bp, were assembled together with the ~10,000 sequences of the S. mediterranea UniGene set using different similarity cutoffs. The assembly was then mapped onto the current genome data. Remarkably, our Smed454 dataset contains more than 3 million novel transcribed nucleotides sequenced for the first time. A descriptive analysis of planarian splice sites was conducted on those Smed454 contigs that mapped univocally to the current genome assembly. Sequence analysis allowed us to identify genes encoding putative proteins with defined structural properties, such as transmembrane domains. Moreover, we annotated the Smed454 dataset using Gene Ontology, and identified putative homologues of several gene families that may play a key role during regeneration, such as neurotransmitter and hormone receptors, homeobox-containing genes, and genes related to eye function. Conclusions: We report the first planarian transcript dataset, Smed454, as an open resource tool that can be accessed via a web interface. Smed454 contains significant novel sequence information about most expressed genes of S. mediterranea. Analysis of the annotated data promises to contribute to identification of gene families poorly characterized at a functional level. The Smed454 transcriptome data will assist in the molecular characterization of S. mediterranea as a model organism, which will be useful to a broad scientific community.

Smed454 dataset: unravelling the transcriptome of Schmidtea mediterranea

Background Freshwater planarians are an attractive model for regeneration and stem cell research and have become a promising tool in the field of regenerative medicine. With the availability of a sequenced planarian genome, the recent application of modern genetic and high-throughput tools has resulted in revitalized interest in these animals, long known for their amazing regenerative capabilities, which enable them to regrow even a new head after decapitation. However, a detailed description of the planarian transcriptome is essential for future investigation into regenerative processes using planarians as a model system. Results In order to complement and improve existing gene annotations, we used a 454 pyrosequencing approach to analyze the transcriptome of the planarian species Schmidtea mediterranea Altogether, 598,435 454-sequencing reads, with an average length of 327 bp, were assembled together with the ~10,000 sequences of the S. mediterranea UniGene set using different similarity cutoffs. The assembly was then mapped onto the current genome data. Remarkably, our Smed454 dataset contains more than 3 million novel transcribed nucleotides sequenced for the first time. A descriptive analysis of planarian splice sites was conducted on those Smed454 contigs that mapped univocally to the current genome assembly. Sequence analysis allowed us to identify genes encoding putative proteins with defined structural properties, such as transmembrane domains. Moreover, we annotated the Smed454 dataset using Gene Ontology, and identified putative homologues of several gene families that may play a key role during regeneration, such as neurotransmitter and hormone receptors, homeobox-containing genes, and genes related to eye function. Conclusions We report the first planarian transcript dataset, Smed454, as an open resource tool that can be accessed via a web interface. Smed454 contains significant novel sequence information about most expressed genes of S. mediterranea. Analysis of the annotated data promises to contribute to identification of gene families poorly characterized at a functional level. The Smed454 transcriptome data will assist in the molecular characterization of S. mediterranea as a model organism, which will be useful to a broad scientific community.

First day of an oil spill on the open sea: early mass transfers of hydrocarbons to air and water.

During the first hours after release of petroleum at sea, crude oil hydrocarbons partition rapidly into air and water. However, limited information is available about very early evaporation and dissolution processes. We report on the composition of the oil slick during the first day after a permitted, unrestrained 4.3 m(3) oil release conducted on the North Sea. Rapid mass transfers of volatile and soluble hydrocarbons were observed, with >50% of ≤C17 hydrocarbons disappearing within 25 h from this oil slick of <10 km(2) area and <10 μm thickness. For oil sheen, >50% losses of ≤C16 hydrocarbons were observed after 1 h. We developed a mass transfer model to describe the evolution of oil slick chemical composition and water column hydrocarbon concentrations. The model was parametrized based on environmental conditions and hydrocarbon partitioning properties estimated from comprehensive two-dimensional gas chromatography (GC×GC) retention data. The model correctly predicted the observed fractionation of petroleum hydrocarbons in the oil slick resulting from evaporation and dissolution. This is the first report on the broad-spectrum compositional changes in oil during the first day of a spill at the sea surface. Expected outcomes under other environmental conditions are discussed, as well as comparisons to other models.

Development of a Database for Drilled Shaft Foundation Testing (DSHAFT), June 2012

Drilled shafts have been used in the US for more than 100 years in bridges and buildings as a deep foundation alternative. For many of these applications, the drilled shafts were designed using the Working Stress Design (WSD) approach. Even though WSD has been used successfully in the past, a move toward Load Resistance Factor Design (LRFD) for foundation applications began when the Federal Highway Administration (FHWA) issued a policy memorandum on June 28, 2000.The policy memorandum requires all new bridges initiated after October 1, 2007, to be designed according to the LRFD approach. This ensures compatibility between the superstructure and substructure designs, and provides a means of consistently incorporating sources of uncertainty into each load and resistance component. Regionally-calibrated LRFD resistance factors are permitted by the American Association of State Highway and Transportation Officials (AASHTO) to improve the economy and competitiveness of drilled shafts. To achieve this goal, a database for Drilled SHAft Foundation Testing (DSHAFT) has been developed. DSHAFT is aimed at assimilating high quality drilled shaft test data from Iowa and the surrounding regions, and identifying the need for further tests in suitable soil profiles. This report introduces DSHAFT and demonstrates its features and capabilities, such as an easy-to-use storage and sharing tool for providing access to key information (e.g., soil classification details and cross-hole sonic logging reports). DSHAFT embodies a model for effective, regional LRFD calibration procedures consistent with PIle LOad Test (PILOT) database, which contains driven pile load tests accumulated from the state of Iowa. PILOT is now available for broader use at the project website: http://srg.cce.iastate.edu/lrfd/. DSHAFT, available in electronic form at http://srg.cce.iastate.edu/dshaft/, is currently comprised of 32 separate load tests provided by Illinois, Iowa, Minnesota, Missouri and Nebraska state departments of transportation and/or department of roads. In addition to serving as a manual for DSHAFT and providing a summary of the available data, this report provides a preliminary analysis of the load test data from Iowa, and will open up opportunities for others to share their data through this quality–assured process, thereby providing a platform to improve LRFD approach to drilled shafts, especially in the Midwest region.

The Brazilian electronic theses and dissertations digital library: providing open access for scholarly information

This paper describes a project led by the Instituto Brasileiro de Informações em Ciência e Tecnologia (Ibict), a government institution, to build a national digital library for electronic theses and dissertations - Bibliteca Digital de Teses e Dissertações (BDTD). The project has been a collaborative effort among Ibict, universities and other research centers in Brazil. The developers adopted a system architecture based on the Open Archives Initiative (OAI) in which universities and research centers act as data providers and Ibict as a service provider. A Brazilian metadata standard for electronic theses and dissertations was developed for the digital library. A toolkit including open source package was also developed by Ibict to be distributed to potential data providers. BDTD has been integrated with the international initiative: the Networked Digital Library of Thesis and Dissertation (NDLTD). Discussions in the paper address various issues related to project design, development and management as well as the role played by Ibict. Conclusions highlight some important lessons learned to date and challenges for the future in expanding the BDTD project.

The "open" chimney graft technique for juxtarenal aortic aneurysms with discrepant renal arteries.

OBJECTIVES: A straightforward original Chimney Graft (CG) protocol has been developed at our institution in selected cases of juxtarenal aortic aneurysm (JRAA). The aim of this study was to present our clinical experience of consecutive series with use of uncovered self-expanding stent (SES) as "Open Chimney" (OCh) in the endovascular repair (EVAR) of JRAA. METHODS: A standard endograft with suprarenal fixation struts is delivered with its proximal covered edge just below the highest RA in JRAA presenting the ostium of the two renal arteries at a different aortic level and the distance between the highest renal artery and the beginning of the aneurysm (improved landing zone) ≥10 mm. The low-lying renal artery is maintained patent by the OCh graft (standard SES) delivered from left brachial access (6 Fr). All clinical, anatomical, and operative data were prospectively collected and retrieved for the study analysis. RESULTS: From July 2010 to November 2012, OCh EVAR was offered to 22 consecutive patients considered unfit for JRAA open repair. All procedures were technically successful with aneurysm exclusion and patent OCh graft. One small perioperative type Ia endoleak spontaneously disappeared at the 3-month CT control. One patient died because of acute decompensated heart failure. One patient presented a left hemispheric stroke. The median follow-up of 18 months (range 7-35) showed aneurysm exclusion in all patients without type I and III endoleaks, SES stenosis, and/or renal impairment. CONCLUSIONS: OCh-EVAR is a straightforward technique that can be employed in selected cases of JRAA, avoiding the more complex and expensive fenestrated EVAR.

The connectome viewer toolkit: an open source framework to manage, analyze, and visualize connectomes.

Advanced neuroinformatics tools are required for methods of connectome mapping, analysis, and visualization. The inherent multi-modality of connectome datasets poses new challenges for data organization, integration, and sharing. We have designed and implemented the Connectome Viewer Toolkit - a set of free and extensible open source neuroimaging tools written in Python. The key components of the toolkit are as follows: (1) The Connectome File Format is an XML-based container format to standardize multi-modal data integration and structured metadata annotation. (2) The Connectome File Format Library enables management and sharing of connectome files. (3) The Connectome Viewer is an integrated research and development environment for visualization and analysis of multi-modal connectome data. The Connectome Viewer's plugin architecture supports extensions with network analysis packages and an interactive scripting shell, to enable easy development and community contributions. Integration with tools from the scientific Python community allows the leveraging of numerous existing libraries for powerful connectome data mining, exploration, and comparison. We demonstrate the applicability of the Connectome Viewer Toolkit using Diffusion MRI datasets processed by the Connectome Mapper. The Connectome Viewer Toolkit is available from http://www.cmtk.org/

ConnectomeViewer - Multi-Modal Multi-Level Network Visualization and Analysis.

Introduction: The field of Connectomic research is growing rapidly, resulting from methodological advances in structural neuroimaging on many spatial scales. Especially progress in Diffusion MRI data acquisition and processing made available macroscopic structural connectivity maps in vivo through Connectome Mapping Pipelines (Hagmann et al, 2008) into so-called Connectomes (Hagmann 2005, Sporns et al, 2005). They exhibit both spatial and topological information that constrain functional imaging studies and are relevant in their interpretation. The need for a special-purpose software tool for both clinical researchers and neuroscientists to support investigations of such connectome data has grown. Methods: We developed the ConnectomeViewer, a powerful, extensible software tool for visualization and analysis in connectomic research. It uses the novel defined container-like Connectome File Format, specifying networks (GraphML), surfaces (Gifti), volumes (Nifti), track data (TrackVis) and metadata. Usage of Python as programming language allows it to by cross-platform and have access to a multitude of scientific libraries. Results: Using a flexible plugin architecture, it is possible to enhance functionality for specific purposes easily. Following features are already implemented: * Ready usage of libraries, e.g. for complex network analysis (NetworkX) and data plotting (Matplotlib). More brain connectivity measures will be implemented in a future release (Rubinov et al, 2009). * 3D View of networks with node positioning based on corresponding ROI surface patch. Other layouts possible. * Picking functionality to select nodes, select edges, get more node information (ConnectomeWiki), toggle surface representations * Interactive thresholding and modality selection of edge properties using filters * Arbitrary metadata can be stored for networks, thereby allowing e.g. group-based analysis or meta-analysis. * Python Shell for scripting. Application data is exposed and can be modified or used for further post-processing. * Visualization pipelines using filters and modules can be composed with Mayavi (Ramachandran et al, 2008). * Interface to TrackVis to visualize track data. Selected nodes are converted to ROIs for fiber filtering The Connectome Mapping Pipeline (Hagmann et al, 2008) processed 20 healthy subjects into an average Connectome dataset. The Figures show the ConnectomeViewer user interface using this dataset. Connections are shown that occur in all 20 subjects. The dataset is freely available from the homepage (connectomeviewer.org). Conclusions: The ConnectomeViewer is a cross-platform, open-source software tool that provides extensive visualization and analysis capabilities for connectomic research. It has a modular architecture, integrates relevant datatypes and is completely scriptable. Visit www.connectomics.org to get involved as user or developer.

GSVA: gene set variation analysis for microarray and RNA-seq data

Gene set enrichment (GSE) analysis is a popular framework for condensing information from gene expression profiles into a pathway or signature summary. The strengths of this approach over single gene analysis include noise and dimension reduction, as well as greater biological interpretability. As molecular profiling experiments move beyond simple case-control studies, robust and flexible GSE methodologies are needed that can model pathway activity within highly heterogeneous data sets. To address this challenge, we introduce Gene Set Variation Analysis (GSVA), a GSE method that estimates variation of pathway activity over a sample population in an unsupervised manner. We demonstrate the robustness of GSVA in a comparison with current state of the art sample-wise enrichment methods. Further, we provide examples of its utility in differential pathway activity and survival analysis. Lastly, we show how GSVA works analogously with data from both microarray and RNA-seq experiments. GSVA provides increased power to detect subtle pathway activity changes over a sample population in comparison to corresponding methods. While GSE methods are generally regarded as end points of a bioinformatic analysis, GSVA constitutes a starting point to build pathway-centric models of biology. Moreover, GSVA contributes to the current need of GSE methods for RNA-seq data. GSVA is an open source software package for R which forms part of the Bioconductor project and can be downloaded at http://www.bioconductor.org.

Reliable selfing rate estimates from imperfect population genetic data.

Genotypic frequencies at codominant marker loci in population samples convey information on mating systems. A classical way to extract this information is to measure heterozygote deficiencies (FIS) and obtain the selfing rate s from FIS = s/(2 - s), assuming inbreeding equilibrium. A major drawback is that heterozygote deficiencies are often present without selfing, owing largely to technical artefacts such as null alleles or partial dominance. We show here that, in the absence of gametic disequilibrium, the multilocus structure can be used to derive estimates of s independent of FIS and free of technical biases. Their statistical power and precision are comparable to those of FIS, although they are sensitive to certain types of gametic disequilibria, a bias shared with progeny-array methods but not FIS. We analyse four real data sets spanning a range of mating systems. In two examples, we obtain s = 0 despite positive FIS, strongly suggesting that the latter are artefactual. In the remaining examples, all estimates are consistent. All the computations have been implemented in a open-access and user-friendly software called rmes (robust multilocus estimate of selfing) available at http://ftp.cefe.cnrs.fr, and can be used on any multilocus data. Being able to extract the reliable information from imperfect data, our method opens the way to make use of the ever-growing number of published population genetic studies, in addition to the more demanding progeny-array approaches, to investigate selfing rates.