Inbreeding coefficient estimation with dense SNP data: comparison of strategies and application to HapMap III.

BACKGROUND/AIMS: If the parents of an individual are related, it is possible for the individual to have received at 1 locus 2 identical-by-descent alleles that are copies of a single allele carried by the parents' common ancestor. The inbreeding coefficient measures the probability of this event and increases with increasing relatedness between the parents. It is traditionally computed from the observed inbreeding loops in the genealogies and its accuracy thus depends on the depth and reliability of the genealogies. With the availability of genome-wide genetic data, it has become possible to compute a genome-based inbreeding coefficient f, and different methods have been developed to estimate f and identify inbred individuals in a sample from the observed patterns of homozygosity at markers. METHODS: For this paper, we performed simulations with known genealogies using different SNP panels with different levels of linkage disequilibrium (LD) to compare several estimators of f, including single-point estimates, methods based on the length of runs of homozygosity (ROHs) and different methods that use hidden Markov models (HMMs). We also compared the performances of some of these estimators to identify inbred individuals in a sample using either HMM likelihood ratio tests or an adapted version of the ERSA software. RESULTS: Single-point methods were found to have higher standard deviations than other methods. ROHs gave the best estimates provided the correct length threshold is known. HMMs on sparse data gave equivalent or better results than HMMs modeling LD. Provided LD is correctly accounted for, the inbreeding estimates were very similar using the different SNP panels. The HMM likelihood ratio tests were found to perform better at detecting inbred individuals in a sample than the adapted ERSA. All methods accurately detected inbreeding up to second-cousin offspring. We applied the best method on release 3 of the HapMap phase III project, found up to 4% of inbred individuals, and created HAP1067, an unrelated and outbred dataset of this release. CONCLUSIONS: We recommend using HMMs on multiple sparse maps to estimate and detect inbreeding in large samples. If the sample of individuals is too small to estimate allele frequencies, we advise to estimate them on reference panels or to use 1,500-kb ROHs. Finally, we suggest to investigators using HapMap to be careful with inbred individuals, especially in the GIH (Gujarati Indians from Houston in Texas) population.

Mobyle: a new full web bioinformatics framework.

MOTIVATION: For the biologist, running bioinformatics analyses involves a time-consuming management of data and tools. Users need support to organize their work, retrieve parameters and reproduce their analyses. They also need to be able to combine their analytic tools using a safe data flow software mechanism. Finally, given that scientific tools can be difficult to install, it is particularly helpful for biologists to be able to use these tools through a web user interface. However, providing a web interface for a set of tools raises the problem that a single web portal cannot offer all the existing and possible services: it is the user, again, who has to cope with data copy among a number of different services. A framework enabling portal administrators to build a network of cooperating services would therefore clearly be beneficial. RESULTS: We have designed a system, Mobyle, to provide a flexible and usable Web environment for defining and running bioinformatics analyses. It embeds simple yet powerful data management features that allow the user to reproduce analyses and to combine tools using a hierarchical typing system. Mobyle offers invocation of services distributed over remote Mobyle servers, thus enabling a federated network of curated bioinformatics portals without the user having to learn complex concepts or to install sophisticated software. While being focused on the end user, the Mobyle system also addresses the need, for the bioinfomatician, to automate remote services execution: PlayMOBY is a companion tool that automates the publication of BioMOBY web services, using Mobyle program definitions. AVAILABILITY: The Mobyle system is distributed under the terms of the GNU GPLv2 on the project web site (http://bioweb2.pasteur.fr/projects/mobyle/). It is already deployed on three servers: http://mobyle.pasteur.fr, http://mobyle.rpbs.univ-paris-diderot.fr and http://lipm-bioinfo.toulouse.inra.fr/Mobyle. The PlayMOBY companion is distributed under the terms of the CeCILL license, and is available at http://lipm-bioinfo.toulouse.inra.fr/biomoby/PlayMOBY/.

Narcisse: a mirror view of conserved syntenies.

New methods and tools are needed to exploit the unprecedented source of information made available by the completed and ongoing whole genome sequencing projects. The Narcisse database is dedicated to the study of genome conservation, from sequence similarities to conserved chromosomal segments or conserved syntenies, for a large number of animals, plants and bacterial completely sequenced genomes. The query interface, a comparative genome browser, enables to navigate between genome dotplots, comparative maps and sequence alignments. The Narcisse database can be accessed at http://narcisse.toulouse.inra.fr.

Genome Size Variation and Comparative Genomics Reveal Intraspecific Diversity in Brassica rapa.

Traditionally, reference genomes in crop species rely on the assembly of one accession, thus occulting most of intraspecific diversity. However, rearrangements, gene duplications, and transposable element content may have a large impact on the genomic structure, which could generate new phenotypic traits. Comparing two Brassica rapa genomes recently sequenced and assembled using long-read technology and optical mapping, we investigated structural variants and repetitive content between the two accessions and genome size variation among a core collection. We explored the structural consequences of the presence of large repeated sequences in B. rapa 'Z1' genome vs. the B. rapa 'Chiifu' genome, using comparative genomics and cytogenetic approaches. First, we showed that large genomic variants on chromosomes A05, A06, A09, and A10 are due to large insertions and inversions when comparing B. rapa 'Z1' and B. rapa 'Chiifu' at the origin of important length differences in some chromosomes. For instance, lengths of 'Z1' and 'Chiifu' A06 chromosomes were estimated in silico to be 55 and 29 Mb, respectively. To validate these observations, we compared using fluorescent in situ hybridization (FISH) the two A06 chromosomes present in an F1 hybrid produced by crossing these two varieties. We confirmed a length difference of 17.6% between the A06 chromosomes of 'Z1' compared to 'Chiifu.' Alternatively, using a copy number variation approach, we were able to quantify the presence of a higher number of rDNA and gypsy elements in 'Z1' genome compared to 'Chiifu' on different chromosomes including A06. Using flow cytometry, the total genome size of 12 Brassica accessions corresponding to a B. rapa available core collection was estimated and revealed a genome size variation of up to 16% between these accessions as well as some shared inversions. This study revealed the contribution of long-read sequencing of new accessions belonging to different cultigroups of B. rapa and highlighted the potential impact of differential insertion of repeat elements and inversions of large genomic regions in genome size intraspecific variability.

Associations between nitric oxide synthase genes and exhaled NO-related phenotypes according to asthma status.

BACKGROUND: The nitric oxide (NO) pathway is involved in asthma, and eosinophils participate in the regulation of the NO pool in pulmonary tissues. We investigated associations between single nucleotide polymorphisms (SNPs) of NO synthase genes (NOS) and biological NO-related phenotypes measured in two compartments (exhaled breath condensate and plasma) and blood eosinophil counts. METHODOLOGY: SNPs (N = 121) belonging to NOS1, NOS2 and NOS3 genes were genotyped in 1277 adults from the French Epidemiological study on the Genetics and Environment of Asthma (EGEA). Association analyses were conducted on four quantitative phenotypes: the exhaled fraction of NO (Fe(NO)), plasma and exhaled breath condensate (EBC) nitrite-nitrate levels (NO2-NO3) and blood eosinophils in asthmatics and non-asthmatics separately. Genetic heterogeneity of these phenotypes between asthmatics and non-asthmatics was also investigated. PRINCIPAL FINDINGS: In non-asthmatics, after correction for multiple comparisons, we found significant associations of Fe(NO) levels with three SNPs in NOS3 and NOS2 (P ≤ 0.002), and of EBC NO2-NO3 level with NOS2 (P = 0.002). In asthmatics, a single significant association was detected between Fe(NO) levels and one SNP in NOS3 (P = 0.004). Moreover, there was significant heterogeneity of NOS3 SNP effect on Fe(NO) between asthmatics and non-asthmatics (P = 0.0002 to 0.005). No significant association was found between any SNP and NO2-NO3 plasma levels or blood eosinophil counts. CONCLUSIONS: Variants in NO synthase genes influence Fe(NO) and EBC NO2-NO3 levels in adults. These genetic determinants differ according to asthma status. Significant associations were only detected for exhaled phenotypes, highlighting the critical relevance to have access to specific phenotypes measured in relevant biological fluid.

Transcription reprogramming during root nodule development in Medicago truncatula.

Many genes which are associated with root nodule development and activity in the model legume Medicago truncatula have been described. However information on precise stages of activation of these genes and their corresponding transcriptional regulators is often lacking. Whether these regulators are shared with other plant developmental programs also remains an open question. Here detailed microarray analyses have been used to study the transcriptome of root nodules induced by either wild type or mutant strains of Sinorhizobium meliloti. In this way we have defined eight major activation patterns in nodules and identified associated potential regulatory genes. We have shown that transcription reprogramming during consecutive stages of nodule differentiation occurs in four major phases, respectively associated with (i) early signalling events and/or bacterial infection; plant cell differentiation that is either (ii) independent or (iii) dependent on bacteroid differentiation; (iv) nitrogen fixation. Differential expression of several genes involved in cytokinin biosynthesis was observed in early symbiotic nodule zones, suggesting that cytokinin levels are actively controlled in this region. Taking advantage of databases recently developed for M. truncatula, we identified a small subset of gene expression regulators that were exclusively or predominantly expressed in nodules, whereas most other regulators were also activated under other conditions, and notably in response to abiotic or biotic stresses. We found evidence suggesting the activation of the jasmonate pathway in both wild type and mutant nodules, thus raising questions about the role of jasmonate during nodule development. Finally, quantitative RT-PCR was used to analyse the expression of a series of nodule regulator and marker genes at early symbiotic stages in roots and allowed us to distinguish several early stages of gene expression activation or repression.