Experimental assessment of the accuracy of genomic selection in sugarcane.

Sugarcane cultivars are interspecific hybrids with an aneuploid, highly heterozygous polyploid genome. The complexity of the sugarcane genome is the main obstacle to the use of marker-assisted selection in sugarcane breeding. Given the promising results of recent studies of plant genomic selection, we explored the feasibility of genomic selection in this complex polyploid crop. Genetic values were predicted in two independent panels, each composed of 167 accessions representing sugarcane genetic diversity worldwide. Accessions were genotyped with 1,499 DArT markers. One panel was phenotyped in Reunion Island and the other in Guadeloupe. Ten traits concerning sugar and bagasse contents, digestibility and composition of the bagasse, plant morphology, and disease resistance were used. We used four statistical predictive models: bayesian LASSO, ridge regression, reproducing kernel Hilbert space, and partial least square regression. The accuracy of the predictions was assessed through the correlation between observed and predicted genetic values by cross validation within each panel and between the two panels. We observed equivalent accuracy among the four predictive models for a given trait, and marked differences were observed among traits. Depending on the trait concerned, within-panel cross validation yielded median correlations ranging from 0.29 to 0.62 in the Reunion Island panel and from 0.11 to 0.5 in the Guadeloupe panel. Cross validation between panels yielded correlations ranging from 0.13 for smut resistance to 0.55 for brix. This level of correlations is promising for future implementations. Our results provide the first validation of genomic selection in sugarcane.

A branch-heterogeneous model of protein evolution for efficient inference of ancestral sequences.

Most models of nucleotide or amino acid substitution used in phylogenetic studies assume that the evolutionary process has been homogeneous across lineages and that composition of nucleotides or amino acids has remained the same throughout the tree. These oversimplified assumptions are refuted by the observation that compositional variability characterizes extant biological sequences. Branch-heterogeneous models of protein evolution that account for compositional variability have been developed, but are not yet in common use because of the large number of parameters required, leading to high computational costs and potential overparameterization. Here, we present a new branch-nonhomogeneous and nonstationary model of protein evolution that captures more accurately the high complexity of sequence evolution. This model, henceforth called Correspondence and likelihood analysis (COaLA), makes use of a correspondence analysis to reduce the number of parameters to be optimized through maximum likelihood, focusing on most of the compositional variation observed in the data. The model was thoroughly tested on both simulated and biological data sets to show its high performance in terms of data fitting and CPU time. COaLA efficiently estimates ancestral amino acid frequencies and sequences, making it relevant for studies aiming at reconstructing and resurrecting ancestral amino acid sequences. Finally, we applied COaLA on a concatenate of universal amino acid sequences to confirm previous results obtained with a nonhomogeneous Bayesian model regarding the early pattern of adaptation to optimal growth temperature, supporting the mesophilic nature of the Last Universal Common Ancestor.

Genetic mapping of a major dominant gene for resistance to Ralstonia solanacearum in eggplant.

Resistance of eggplant against Ralstonia solanacearum phylotype I strains was assessed in a F(6) population of recombinant inbred lines (RILs) derived from a intra-specific cross between S. melongena MM738 (susceptible) and AG91-25 (resistant). Resistance traits were determined as disease score, percentage of wilted plants, and stem-based bacterial colonization index, as assessed in greenhouse experiments conducted in Réunion Island, France. The AG91-25 resistance was highly efficient toward strains CMR134, PSS366 and GMI1000, but only partial toward the highly virulent strain PSS4. The partial resistance found against PSS4 was overcome under high inoculation pressure, with heritability estimates from 0.28 to 0.53, depending on the traits and season. A genetic map was built with 119 AFLP, SSR and SRAP markers positioned on 18 linkage groups (LG), for a total length of 884 cM, and used for quantitative trait loci (QTL) analysis. A major dominant gene, named ERs1, controlled the resistance to strains CMR134, PSS366, and GMI1000. Against strain PSS4, this gene was not detected, but a significant QTL involved in delay of disease progress was detected on another LG. The possible use of the major resistance gene ERs1 in marker-assisted selection and the prospects offered for academic studies of a possible gene for gene system controlling resistance to bacterial wilt in solanaceous plants are discussed.

Species identification of staphylococci by amplification and sequencing of the tuf gene compared to the gap gene and by matrix-assisted laser desorption ionization time-of-flight mass spectrometry.

Staphylococcal species, notably, coagulase-negative staphylococci (CoNS), are frequently misidentified using phenotypic methods. The partial nucleotide sequences of the tuf and gap genes were determined in 47 reference strains to assess their suitability, practicability, and discriminatory power as target molecules for staphylococcal identification. The partial tuf gene sequence was selected and further assessed with a collection of 186 strains, including 35 species and subspecies. Then, to evaluate the efficacy of this genotyping method versus the technology of matrix-assisted laser desorption ionization time-of-flight mass spectrometry (MALDI-TOF-MS), the 186 strains were identified using MALDI-TOF-MS (Axima® Shimadzu) coupled to the SARAMIS® database (AnagnosTec). The French National Reference Center for Staphylococci identification method was used as a reference. One hundred and eighty-four strains (98.9%) were correctly identified by tuf gene sequencing. Only one strain was misidentified and one was unidentified. MALDI-TOF-MS identified correctly 138 isolates (74.2%). Four strains were misidentified, 39 were unidentified, five were identified at the group (hominis/warneri) level, and one strain was identified at the genus level. These results confirm the value of MALDI-TOF-MS identification for common species in clinical laboratory practice and the value of the partial tuf gene sequence for the identification of all staphylococcal species as required in a reference laboratory.

Relationship between morphological taxonomy and molecular divergence within Crustacea: proposal of a molecular threshold to help species delimitation.

With today's technology for production of molecular sequences, DNA taxonomy and barcoding arose as a new tool for evolutionary biology and ecology. However, their validities still need to be empirically evaluated. Of most importance is the strength of the correlation between morphological taxonomy and molecular divergence and the possibility to define some molecular thresholds. Here, we report measurements of this correlation for two mitochondrial genes (COI and 16S rRNA) within the sub-phylum Crustacea. Perl scripts were developed to ensure objectivity, reproducibility, and exhaustiveness of our tests. Our analysis reveals a general correlation between molecular divergence and taxonomy. This correlation is particularly high for shallow taxonomic levels allowing us to propose a COI universal crustacean threshold to help species delimitation. At higher taxonomic levels this correlation decreases, particularly when comparing different families. Those results plead for DNA use in taxonomy and suggest an operational method to help crustacean species delimitation that is linked to the phylogenetic species definition. This pragmatic tool is expected to fine tune the present classification, and not, as some would have believed, to tear it apart.

Phylogeography of a subterranean amphipod reveals cryptic diversity and dynamic evolution in extreme environments.

Extreme conditions in subsurface are suspected to be responsible for morphological convergences, and so to bias biodiversity assessment. Subterranean organisms are also considered as having poor dispersal abilities that in turn generate a large number of endemic species when habitat is fragmented. Here we test these general hypotheses using the subterranean amphipod Niphargus virei. All our phylogenetic analyses (Bayesian, maximum likelihood and distance), based on two independent genes (28S and COI), revealed the same tripartite structure. N. virei populations from Benelux, Jura region and the rest of France appeared as independent evolutionary units. Molecular rates estimated via global or Bayesian relaxed clock suggest that this split is at least 13 million years old and accredit the cryptic diversity hypothesis. Moreover, the geographical distribution of these lineages showed some evidence of recent dispersal through apparent vicariant barrier. In consequence, we argue that future analyses of evolution and biogeography in subsurface, or more generally in extreme environments, should consider dispersal ability as an evolving trait and morphology as a potentially biased marker.

Synthetic approaches to a mononucleotide prodrug of cytarabine.

Synthetic pathways to a mononucleotide prodrug of cytarabine (Ara-C) bearing S-pivaloyl-2-thioethyl (tBuSATE) groups, as biolabile phosphate protections, are reported. Using a common phosphoramidite approach, two different kinds of nucleoside protecting groups have been investigated. During this study, we observed an intermolecular migration of the Boc protecting group in the course of the tert-butyldimethylsilyl ether cleavage using tetrabutyl ammonium fluoride.

Genome sequence and gene compaction of the eukaryote parasite Encephalitozoon cuniculi.

Microsporidia are obligate intracellular parasites infesting many animal groups. Lacking mitochondria and peroxysomes, these unicellular eukaryotes were first considered a deeply branching protist lineage that diverged before the endosymbiotic event that led to mitochondria. The discovery of a gene for a mitochondrial-type chaperone combined with molecular phylogenetic data later implied that microsporidia are atypical fungi that lost mitochondria during evolution. Here we report the DNA sequences of the 11 chromosomes of the approximately 2.9-megabase (Mb) genome of Encephalitozoon cuniculi (1,997 potential protein-coding genes). Genome compaction is reflected by reduced intergenic spacers and by the shortness of most putative proteins relative to their eukaryote orthologues. The strong host dependence is illustrated by the lack of genes for some biosynthetic pathways and for the tricarboxylic acid cycle. Phylogenetic analysis lends substantial credit to the fungal affiliation of microsporidia. Because the E. cuniculi genome contains genes related to some mitochondrial functions (for example, Fe-S cluster assembly), we hypothesize that microsporidia have retained a mitochondrion-derived organelle.

Nrh, a human homologue of Nr-13 associates with Bcl-Xs and is an inhibitor of apoptosis.

In search of human homologues of the anti-apoptotic protein Nr-13, we have characterized a human EST clone that potentially encodes a protein, which is the closest homologue of Nr-13 among the Bcl-2 family members, to date known, in humans. Phylogenetic analyses suggest Human nrh, Mouse diva/boo and Quail nr-13 to be orthologous genes. The nrh gene has the same overall organization as nr-13 and diva/boo with one single intron interrupting the ORF at the level of the Bcl-2-homology domain BH2. RT-PCR-based analysis of nrh expression indicated that this gene is preferentially expressed in the lungs, the liver and the kidneys. Interestingly, two in frame ATG codons can lead potentially to the synthesis of two products, one of them lacking 10 aminoacids at the N-terminal end. Sequence alignment with Nr-13 and Diva/Boo in addition to secondary structure prediction of the nrh transcript suggested that the shortest protein will be preferentially synthetized. Immunohistochemical analyses have revealed that Nrh is associated with mitochondria and the nuclear envelope. Moreover, Nrh preferentially associates with the apoptosis accelerator Bcl-Xs and behaves as an inhibitor of apoptosis both in yeast and vertebrate cells.

Molecular evidence for the existence of two species of Marteilia in Europe.

Marteilia refringens is one of the most significant pathogens of bivalve molluscs. Previous sequencing of the small subunit ribosomal RNA gene of M. refringens isolates derived from the infected mussels (Mytilus edulis and Mytilus galloprovinciallis) and the oyster (Ostrea edulis) in Europe did not reveal genetic polymorphisms despite indications from epizootiological data that distinct types may exist. We investigated the existence of polymorphisms in the internal transcribed spacer region of the ribosomal RNA genes. The sequences of this region proved to be clearly dimorphic among Marteilia from five sampling sites. The distribution of the two genetic types, named "O" and "M", appeared to be linked to the host species, oysters and mussels, respectively. We therefore support the recognition of two species of Marteilia in Europe and propose that the "O" type corresponds to M. refringens and the "M" type to M. maurini.

Lipid phosphate phosphatases in Arabidopsis. Regulation of the AtLPP1 gene in response to stress.

An Arabidopsis thaliana gene (AtLPP1) was isolated on the basis that it was transiently induced by ionizing radiation. The putative AtLPP1 gene product showed homology to the yeast and mammalian lipid phosphate phosphatase enzymes and possessed a phosphatase signature sequence motif. Heterologous expression and biochemical characterization of the AtLPP1 gene in yeast showed that it encoded an enzyme (AtLpp1p) that exhibited both diacylglycerol pyrophosphate phosphatase and phosphatidate phosphatase activities. Kinetic analysis indicated that diacylglycerol pyrophosphate was the preferred substrate for AtLpp1p in vitro. A second Arabidopsis gene (AtLPP2) was identified based on sequence homology to AtLPP1 that was also heterologously expressed in yeast. The AtLpp2p enzyme also utilized diacylglycerol pyrophosphate and phosphatidate but with no preference for either substrate. The AtLpp1p and AtLpp2p enzymes showed differences in their apparent affinities for diacylglycerol pyrophosphate and phosphatidate as well as other enzymological properties. Northern blot analyses showed that the AtLPP1 gene was preferentially expressed in leaves and roots, whereas the AtLPP2 gene was expressed in all tissues examined. AtLPP1, but not AtLPP2, was regulated in response to various stress conditions. The AtLPP1 gene was transiently induced by genotoxic stress (gamma ray or UV-B) and elicitor treatments with mastoparan and harpin. The regulation of the AtLPP1 gene in response to stress was consistent with the hypothesis that its encoded lipid phosphate phosphatase enzyme may attenuate the signaling functions of phosphatidate and/or diacylglycerol pyrophosphate that form in response to stress in plants.

Bacterial molecular phylogeny using supertree approach.

It has been claimed that complete genome sequences would clarify phylogenetic relationships between organisms but, up to now, no satisfying approach has been proposed to use efficiently these data. For instance, if the coding of presence or absence of genes in complete genomes gives interesting results, it does not take into account the phylogenetic information contained in sequences and ignores hidden paralogy by using a similarity-based definition of orthology. Also, concatenation of sequences of different genes takes hardly in consideration the specific evolutionary rate of each gene. At last, building a consensus tree is strongly limited by the low number of genes shared among all organisms. Here, we use a new method based on supertree construction, which permits to cumulate in one supertree the information and statistical support of hundreds of trees from orthologous gene families and to build the phylogeny of 33 prokaryotes and four eukaryotes with completely sequenced genomes. This approach gives a robust supertree, which demonstrates that a phylogeny of prokaryotic species is conceivable and challenges the hypothesis of a thermophilic origin of bacteria and present-day life. The results are compatible with the hypothesis of a core of genes for which lateral transfers are rare but they raise doubts on the widely admitted "complexity hypothesis" which predicts that this core is mainly implicated in informational processes.

Phylogenetic analysis of the small subunit ribosomal RNA of Marteilia refringens validates the existence of phylum Paramyxea (Desportes and Perkins, 1990).

Marteilia refringens is recognized as one of the most significant pathogens of bivalve molluscs. The nucleotide sequence of the small subunit ribosomal RNA gene of Marteilia refringens is used to elucidate the phylogenetic position of the phylum Paramyxea. Genomic DNA was extracted from sporangia of Marteilia, purified from infected blue mussels, Mytilus edulis, and flat oysters, Ostrea edulis. The sequences obtained from Marteilia species purified from both oysters and mussels were identical. The sequence identity was confirmed by in situ hybridization using a DNA probe targeted to a variable region of the ribosomal DNA. The small subunit ribosomal RNA gene sequence of M. refringens is very different from all known sequences of eukaryotic organisms, including those of myxosporeans and haplosporeans. Therefore, the phylum Paramyxea should continue to be recognized as an independent eukaryotic phylum.

HOBACGEN: database system for comparative genomics in bacteria.

We present here HOBACGEN, a database system devoted to comparative genomics in bacteria. HOBACGEN contains all available protein genes from bacteria, archaea, and yeast, taken from SWISS-PROT/TrEMBL and classified into families. It also includes multiple alignments and phylogenetic trees built from these families. The database is organized under a client/server architecture with a client written in Java, which may run on any platform. This client integrates a graphical interface allowing users to select families according to various criteria and notably to select homologs common to a given set of taxa. This interface also allows users to visualize multiple alignments and trees associated to families. In tree displays, protein gene names are colored according to the taxonomy of the corresponding organisms. Users may access all information associated to sequences and multiple alignments by clicking on genes. This graphic tool thus gives a rapid and simple access to all data required to interpret homology relationships between genes and distinguish orthologs from paralogs. Instructions for installation of the client or the server are available at http://pbil.univ-lyon1. fr/databases/hobacgen.html.

Accounting for evolutionary rate variation among sequence sites consistently changes universal phylogenies deduced from rRNA and protein-coding genes.

Phylogenetic analyses of gene and protein sequences have led to two major competing views of the universal phylogeny, the evolutionary tree relating the three kinds of living organisms, Bacteria, Archaea, and Eukarya. In the first scheme, called "the archaebacterial tree, " organisms of the same type are clustered together. In the second scenario, called "the eocyte tree," the archaeal phylum of Crenarchaeota is more closely related to eukaryotes than are other Archaea. A major property of the evolution of functional ribosomal and protein-encoding genes is that the rate of nucleotide and amino acid substitution varies across sequence sites. Here, using distance-based and maximum-likelihood methods, we show that universal phylogenies of ribosomal RNAs and RNA polymerases built by ignoring this variation are biased toward the archaebacterial tree because of attraction between long branches. In contrast, taking among-site rate variability into account gives support for the eocyte tree.

A nonhyperthermophilic common ancestor to extant life forms.

The G+C nucleotide content of ribosomal RNA (rRNA) sequences is strongly correlated with the optimal growth temperature of prokaryotes. This property allows inference of the environmental temperature of the common ancestor to all life forms from knowledge of the G+C content of its rRNA sequences. A model of sequence evolution, assuming varying G+C content among lineages and unequal substitution rates among sites, was devised to estimate ancestral base compositions. This method was applied to rRNA sequences of various species representing the major lineages of life. The inferred G+C content of the common ancestor to extant life forms appears incompatible with survival at high temperature. This finding challenges a widely accepted hypothesis about the origin of life.

Sensitivity of the relative-rate test to taxonomic sampling.

Relative-rate tests may be used to compare substitution rates between more than two sequences, which yields two main questions: What influence does the number of sequences have on relative-rate tests and what is the influence of the sampling strategy as characterized by the phylogenetic relationships between sequences? Using both simulations and analysis of real data from murids (APRT and LCAT nuclear genes), we show that comparing large numbers of species significantly improves the power of the test. This effect is stronger if species are more distantly related. On the other hand, it appears to be less rewarding to increase outgroup sampling than to use the single nearest outgroup sequence. Rates may be compared between paraphyletic ingroups and using paraphyletic outgroups, but unbalanced taxonomic sampling can bias the test. We present a simple phylogenetic weighting scheme which takes taxonomic sampling into account and significantly improves the relative-rate test in cases of unbalanced sampling. The answers are thus: (1) large taxonomic sampling of compared groups improves relative-rate tests, (2) sampling many outgroups does not bring significant improvement, (3) the only constraint on sampling strategy is that the outgroup be valid, and (4) results are more accurate when phylogenetic relationships between the investigated sequences are taken into account. Given current limitations of the maximum-likelihood and nonparametric approaches, the relative-rate test generalized to any number of species with phylogenetic weighting appears to be the most general test available to compare rates between lineages.

Inferring pattern and process: maximum-likelihood implementation of a nonhomogeneous model of DNA sequence evolution for phylogenetic analysis.

A nonhomogeneous, nonstationary stochastic model of DNA sequence evolution allowing varying equilibrium G + C contents among lineages is devised in order to deal with sequences of unequal base compositions. A maximum-likelihood implementation of this model for phylogenetic analyses allows handling of a reasonable number of sequences. The relevance of the model and the accuracy of parameter estimates are theoretically and empirically assessed, using real or simulated data sets. Overall, a significant amount of information about past evolutionary modes can be extracted from DNA sequences, suggesting that process (rates of distinct kinds of nucleotide substitutions) and pattern (the evolutionary tree) can be simultaneously inferred. G + C contents at ancestral nodes are quite accurately estimated. The new method appears to be useful for phylogenetic reconstruction when base composition varies among compared sequences. It may also be suitable for molecular evolution studies.

Microsporidian Encephalitozoon cuniculi, a unicellular eukaryote with an unusual chromosomal dispersion of ribosomal genes and a LSU rRNA reduced to the universal core.

Microsporidia are eukaryotic parasites lacking mitochondria, the ribosomes of which present prokaryote-like features. In order to better understand the structural evolution of rRNA molecules in microsporidia, the 5S and rDNA genes were investigated in Encephalitozoon cuniculi . The genes are not in close proximity. Non-tandemly arranged rDNA units are on every one of the 11 chromosomes. Such a dispersion is also shown in two other Encephalitozoon species. Sequencing of the 5S rRNA coding region reveals a 120 nt long RNA which folds according to the eukaryotic consensus structural shape. In contrast, the LSU rRNA molecule is greatly reduced in length (2487 nt). This dramatic shortening is essentially due to truncation of divergent domains, most of them being removed. Most variable stems of the conserved core are also deleted, reducing the LSU rRNA to only those structural features preserved in all living cells. This suggests that the E.cuniculi LSU rRNA performs only the basic mechanisms of translation. LSU rRNA phylogenetic analysis with the BASEML program favours a relatively recent origin of the fast evolving microsporidian lineage. Therefore, the prokaryote-like ribosomal features, such as the absence of ITS2, may be derived rather than primitive characters.