Construction and characterization of a BAC library for functional genomics in Xenopus tropicalis.

Large insert genomic DNA libraries are useful resources for genomic studies. Although the genome of Xenopus tropicalis stands as the amphibian reference genome because it benefitted from large-scale sequencing studies, physical mapping resources such as BAC libraries are lagging behind. Here we present the construction and characterization of a BAC library that covers the whole X. tropicalis genome. We prepared this BAC library from the genomic DNA of X. tropicalis females of the Adiopodoume strain. We characterized BAC clones by screening for specific loci, by chromosomal localization using FISH and by systematic BAC end sequencing. The median insert size is about 110kbp and the library coverage is around six genome equivalents. We obtained a total of 163,787 BAC end sequences with mate pairs for 77,711 BAC clones. We mapped all BAC end sequences to the reference X. tropicalis genome assembly to enable the identification of BAC clones covering specific loci. Overall, this BAC library resource complements the knowledge of the X. tropicalis genome and should further promote its use as a reference genome for developmental biology studies and amphibian comparative genomics.

Exploring the genome of the salt-marsh Spartina maritima (Poaceae, Chloridoideae) through BAC end sequence analysis.

Spartina species play an important ecological role on salt marshes. Spartina maritima is an Old-World species distributed along the European and North-African Atlantic coasts. This hexaploid species (2n = 6x = 60, 2C = 3,700 Mb) hybridized with different Spartina species introduced from the American coasts, which resulted in the formation of new invasive hybrids and allopolyploids. Thus, S. maritima raises evolutionary and ecological interests. However, genomic information is dramatically lacking in this genus. In an effort to develop genomic resources, we analysed 40,641 high-quality bacterial artificial chromosome-end sequences (BESs), representing 26.7 Mb of the S. maritima genome. BESs were searched for sequence homology against known databases. A fraction of 16.91% of the BESs represents known repeats including a majority of long terminal repeat (LTR) retrotransposons (13.67%). Non-LTR retrotransposons represent 0.75%, DNA transposons 0.99%, whereas small RNA, simple repeats and low-complexity sequences account for 1.38% of the analysed BESs. In addition, 4,285 simple sequence repeats were detected. Using the coding sequence database of Sorghum bicolor, 6,809 BESs found homology accounting for 17.1% of all BESs. Comparative genomics with related genera reveals that the microsynteny is better conserved with S. bicolor compared to other sequenced Poaceae, where 37.6% of the paired matching BESs are correctly orientated on the chromosomes. We did not observe large macrosyntenic rearrangements using the mapping strategy employed. However, some regions appeared to have experienced rearrangements when comparing Spartina to Sorghum and to Oryza. This work represents the first overview of S. maritima genome regarding the respective coding and repetitive components. The syntenic relationships with other grass genomes examined here help clarifying evolution in Poaceae, S. maritima being a part of the poorly-known Chloridoideae sub-family.

Transcriptome de novo assembly from next-generation sequencing and comparative analyses in the hexaploid salt marsh species Spartina maritima and Spartina alterniflora (Poaceae).

Spartina species have a critical ecological role in salt marshes and represent an excellent system to investigate recurrent polyploid speciation. Using the 454 GS-FLX pyrosequencer, we assembled and annotated the first reference transcriptome (from roots and leaves) for two related hexaploid Spartina species that hybridize in Western Europe, the East American invasive Spartina alterniflora and the Euro-African S. maritima. The de novo read assembly generated 38 478 consensus sequences and 99% found an annotation using Poaceae databases, representing a total of 16 753 non-redundant genes. Spartina expressed sequence tags were mapped onto the Sorghum bicolor genome, where they were distributed among the subtelomeric arms of the 10 S. bicolor chromosomes, with high gene density correlation. Normalization of the complementary DNA library improved the number of annotated genes. Ecologically relevant genes were identified among GO biological function categories in salt and heavy metal stress response, C4 photosynthesis and in lignin and cellulose metabolism. Expression of some of these genes had been found to be altered by hybridization and genome duplication in a previous microarray-based study in Spartina. As these species are hexaploid, up to three duplicated homoeologs may be expected per locus. When analyzing sequence polymorphism at four different loci in S. maritima and S. alterniflora, we found up to four haplotypes per locus, suggesting the presence of two expressed homoeologous sequences with one or two allelic variants each. This reference transcriptome will allow analysis of specific Spartina genes of ecological or evolutionary interest, estimation of homoeologous gene expression variation using RNA-seq and further gene expression evolution analyses in natural populations.

Estimate of human gene number provided by genome-wide analysis using Tetraodon nigroviridis DNA sequence.

The number of genes in the human genome is unknown, with estimates ranging from 50,000 to 90,000 (refs 1, 2), and to more than 140,000 according to unpublished sources. We have developed 'Exofish', a procedure based on homology searches, to identify human genes quickly and reliably. This method relies on the sequence of another vertebrate, the pufferfish Tetraodon nigroviridis, to detect conserved sequences with a very low background. Similar to Fugu rubripes, a marine pufferfish proposed by Brenner et al. as a model for genomic studies, T. nigroviridis is a more practical alternative with a genome also eight times more compact than that of human. Many comparisons have been made between F. rubripes and human DNA that demonstrate the potential of comparative genomics using the pufferfish genome. Application of Exofish to the December version of the working draft sequence of the human genome and to Unigene showed that the human genome contains 28,000-34,000 genes, and that Unigene contains less than 40% of the protein-coding fraction of the human genome.

Spastin, a new AAA protein, is altered in the most frequent form of autosomal dominant spastic paraplegia.

Autosomal dominant hereditary spastic paraplegia (AD-HSP) is a genetically heterogeneous neurodegenerative disorder characterized by progressive spasticity of the lower limbs. Among the four loci causing AD-HSP identified so far, the SPG4 locus at chromosome 2p2-1p22 has been shown to account for 40-50% of all AD-HSP families. Using a positional cloning strategy based on obtaining sequence of the entire SPG4 interval, we identified a candidate gene encoding a new member of the AAA protein family, which we named spastin. Sequence analysis of this gene in seven SPG4-linked pedigrees revealed several DNA modifications, including missense, nonsense and splice-site mutations. Both SPG4 and its mouse orthologue were shown to be expressed early and ubiquitously in fetal and adult tissues. The sequence homologies and putative subcellular localization of spastin suggest that this ATPase is involved in the assembly or function of nuclear protein complexes.

DNA from soil mirrors plant taxonomic and growth form diversity.

Ecosystems across the globe are threatened by climate change and human activities. New rapid survey approaches for monitoring biodiversity would greatly advance assessment and understanding of these threats. Taking advantage of next-generation DNA sequencing, we tested an approach we call metabarcoding: high-throughput and simultaneous taxa identification based on a very short (usually <100 base pairs) but informative DNA fragment. Short DNA fragments allow the use of degraded DNA from environmental samples. All analyses included amplification using plant-specific versatile primers, sequencing and estimation of taxonomic diversity. We tested in three steps whether degraded DNA from dead material in soil has the potential of efficiently assessing biodiversity in different biomes. First, soil DNA from eight boreal plant communities located in two different vegetation types (meadow and heath) was amplified. Plant diversity detected from boreal soil was highly consistent with plant taxonomic and growth form diversity estimated from conventional above-ground surveys. Second, we assessed DNA persistence using samples from formerly cultivated soils in temperate environments. We found that the number of crop DNA sequences retrieved strongly varied with years since last cultivation, and crop sequences were absent from nearby, uncultivated plots. Third, we assessed the universal applicability of DNA metabarcoding using soil samples from tropical environments: a large proportion of species and families from the study site were efficiently recovered. The results open unprecedented opportunities for large-scale DNA-based biodiversity studies across a range of taxonomic groups using standardized metabarcoding approaches.

Deep RNA sequencing improved the structural annotation of the Tuber melanosporum transcriptome.

• The functional complexity of the Tuber melanosporum transcriptome has not yet been fully elucidated. Here, we applied high-throughput Illumina RNA-sequencing (RNA-Seq) to the transcriptome of T. melanosporum at different major developmental stages, that is free-living mycelium, fruiting body and ectomycorrhiza. • Sequencing of cDNA libraries generated a total of c. 24 million sequence reads representing > 882 Mb of sequence data. To construct a coverage signal profile across the genome, all reads were then aligned to the reference genome assembly of T. melanosporum Mel28. • We were able to identify a substantial number of novel transcripts, antisense transcripts, new exons, untranslated regions (UTRs), alternative upstream initiation codons and upstream open reading frames. • This RNA-Seq analysis allowed us to improve the genome annotation. It also provided us with a genome-wide view of the transcriptional and post-transcriptional mechanisms generating an increased number of transcript isoforms during major developmental transitions in T. melanosporum.

Finding candidate genes under positive selection in Non-model species: examples of genes involved in host specialization in pathogens.

Numerous genes in diverse organisms have been shown to be under positive selection, especially genes involved in reproduction, adaptation to contrasting environments, hybrid inviability, and host-pathogen interactions. Looking for genes under positive selection in pathogens has been a priority in efforts to investigate coevolution dynamics and to develop vaccines or drugs. To elucidate the functions involved in host specialization, here we aimed at identifying candidate sequences that could have evolved under positive selection among closely related pathogens specialized on different hosts. For this goal, we sequenced c. 17,000-32,000 ESTs from each of four Microbotryum species, which are fungal pathogens responsible for anther smut disease on host plants in the Caryophyllaceae. Forty-two of the 372 predicted orthologous genes showed significant signal of positive selection, which represents a good number of candidate genes for further investigation. Sequencing 16 of these genes in 9 additional Microbotryum species confirmed that they have indeed been rapidly evolving in the pathogen species specialized on different hosts. The genes showing significant signals of positive selection were putatively involved in nutrient uptake from the host, secondary metabolite synthesis and secretion, respiration under stressful conditions and stress response, hyphal growth and differentiation, and regulation of expression by other genes. Many of these genes had transmembrane domains and may therefore also be involved in pathogen recognition by the host. Our approach thus revealed fruitful and should be feasible for many non-model organisms for which candidate genes for diversifying selection are needed.

Cure of short- and long-term experimental Chagas' disease using D0870.

Chagas' disease, a protozoan infection by the kinetoplastid Trypanosoma cruzi, constitutes a major public health problem in Latin America. With the use of mouse models of both short- and long-term forms of the disease, the efficacy of D0870, a bis-triazole derivative, was tested. D0870 was able to prevent death and induced parasitological cure in 70 to 90 percent of animals, in both the short- and long-term disease. In contrast, currently used drugs such as nifurtimox or ketoconazole prolonged survival but did not induce significant curing effects. D0870 may be useful in the treatment of human long-term Chagas' disease, a condition that is currently incurable.

Single cell genomics yields a wide diversity of small planktonic protists across major ocean ecosystems.

Marine planktonic protists are critical components of ocean ecosystems and are highly diverse. Molecular sequencing methods are being used to describe this diversity and reveal new associations and metabolisms that are important to how these ecosystems function. We describe here the use of the single cell genomics approach to sample and interrogate the diversity of the smaller (pico- and nano-sized) protists from a range of oceanic samples. We created over 900 single amplified genomes (SAGs) from 8 Tara Ocean samples across the Indian Ocean and the Mediterranean Sea. We show that flow cytometric sorting of single cells effectively distinguishes plastidic and aplastidic cell types that agree with our understanding of protist phylogeny. Yields of genomic DNA with PCR-identifiable 18S rRNA gene sequence from single cells was low (15% of aplastidic cell sorts, and 7% of plastidic sorts) and tests with alternate primers and comparisons to metabarcoding did not reveal phylogenetic bias in the major protist groups. There was little evidence of significant bias against or in favor of any phylogenetic group expected or known to be present. The four open ocean stations in the Indian Ocean had similar communities, despite ranging from 14°N to 20°S latitude, and they differed from the Mediterranean station. Single cell genomics of protists suggests that the taxonomic diversity of the dominant taxa found in only several hundreds of microliters of surface seawater is similar to that found in molecular surveys where liters of sample are filtered.

A direct method for the chromosomal assignment of DNA markers in Leishmania.

A simple method for the chromosomal assignment of any DNA marker would be an important tool for the ongoing project to map the genome of the protozoan parasite Leishmania. The Leishmania chromosomes enter pulsed field gel electrophoresis (PFGE) gels under current electrophoretic conditions, but their direct identification in a given strain is hampered by their stacking in a few chromosomal bands, and by the very frequent size variations of the same chromosome among parasite strains. To overcome these problems. we determined the complete karyotypes of 12 Old World Leishmania cloned strains. This enabled us to select three of these strains that display great chromosome size polymorphisms, such that every chromosome can be individualized by a specific pattern after hybridization onto these three karyotypes. The complete resolution of the genomes of these three strains can be carried out with only three electrophoretic conditions. This makes a series of three blots sufficient for the assignment of any new marker on a particular Leishmania chromosome.

Conserved linkage groups associated with large-scale chromosomal rearrangements between Old World and New World Leishmania genomes.

The genus Leishmania can be taxonomically separated into three main groups: the Old World subgenus L. (Leishmania), the New World subgenus L. (Leishmania) and the New World subgenus L. (Viannia). The haploid genome of Old World Leishmania species has been shown to contain 36 chromosomes defined as physical linkage groups; the latter were found entirely conserved across species. In the present study, we tried to verify whether this conservation of the genome structure extends to the New World species of Leishmania. 300 loci were explored by hybridization on optimized pulsed field gel electrophoresis separations of the chromosomes of polymorphic strains of the six main pathogenic Leishmania species of the New World. When comparing these New World karyotypes with their Old World counterparts, 32 out of 36 linkage groups were found conserved among all species. Four chromosomal rearrangements were found. All species belonging to the L. (Viannia) subgenus were characterized by the presence (i) of a short sequence exchange between chromosomes 26 and 35, and (ii) more importantly, of a fused version of chromosomes 20 and 34 which are separated in all Old World species. 69 additional markers were isolated from a plasmid library specifically constructed from the rearranged chromosomes 20+34 in an attempt to detect mechanisms other than a fusion or breakage: only two markers out of 40 did not belong to the linkage groups 20 and 34. On the other hand, all strains belonging to the New World subgenus L. (Leishmania) were characterized by two different chromosomal rearrangements of the same type (fusion/breakage) as above as compared with Old World species: chromosomes 8+29 and 20+36. Consequently, these two groups of species have 35 and 34 heterologous chromosomes, respectively. Overall, these results show that large-scale chromosomal rearrangements occurred during the evolution of the genus Leishmania, and that the three main groups of pathogenic species are characterized by different chromosome numbers. Nevertheless, translocations seem particularly rare, and the conservation of the major linkage groups should be an essential feature for the compared genetics between species of this parasite.

Genomic exploration of the hemiascomycetous yeasts: 16. Candida tropicalis.

The genome of the diploid hemiascomycetous yeast Candida tropicalis, an opportunistic human pathogen and an important organism for industrial applications, was explored by the analysis of 2541 Random Sequenced Tags (RSTs) covering about 20% of its genome. Comparison of these sequences with Saccharomyces cerevisiae and other species permitted the identification and the analysis of a total of more than 1000 novel genetic elements of C. tropicalis. Moreover, the present study confirms that in C. tropicalis, the rare CUG codon is read as a serine and not a leucine. The sequences have been deposited at EMBL with the accession numbers AL438875-AL441602.

Genomic exploration of the hemiascomycetous yeasts: 5. Saccharomyces bayanus var. uvarum.

Saccharomyces bayanus var. uvarum investigated here is the species closest to Saccharomyces cerevisiae. Random sequence tags (RSTs) allowed us to identify homologues to 2789 open reading frames (ORFs) in S. cerevisiae, ORFs duplicated in S. uvarum but not in S. cerevisiae, centromeres, tRNAs, homologues of Ty1/2 and Ty4 retrotransposons, and a complete rDNA repeat. Only 13 RSTs seem to be homologous to sequences in other organisms but not in S. cerevisiae. As the synteny between the two species is very high, cases in which synteny is lost suggest special mechanisms of genome evolution. The corresponding RSTs revealed that S. uvarum can exist without any S. cerevisiae DNA introgression. Accession numbers are from AL397139 to AL402278 in the EMBL databank.

Genomic exploration of the hemiascomycetous yeasts: 6. Saccharomyces exiguus.

Random sequence tags were obtained from a genomic DNA library of Saccharomyces exiguus. The mitochondrial genome appeared to be at least 25.7 kb in size, with a different organization compared to Saccharomyces cerevisiae. An unusual putative 953 bp long terminal repeated element associated to Ty3 was found. A set of 1451 genes was identified homologous to S. cerevisiae open reading frames. Only five genes were identified outside the S. cerevisiae taxon, confirming that S. exiguus is phylogenetically closely related to S. cerevisiae. Unexpectedly, numerous duplicated genes were found whereas they are unique in S. cerevisiae. The sequences are deposited at EMBL under the accession numbers: AL407377-AL409955.

Genomic exploration of the hemiascomycetous yeasts: 7. Saccharomyces servazzii.

The genome of Saccharomyces servazzii was analyzed with 2570 random sequence tags totalling 2.3 Mb. BLASTX comparisons revealed a minimum of 1420 putative open reading frames with significant homology to Saccharomyces cerevisiae (58% aa identity on average), two with Schizosaccharomyces pombe and one with a human protein, confirming that S. servazzii is closely related to S. cerevisiae. About 25% of the S. servazzii genes were identified, assuming that the gene complement is identical in both yeasts. S. servazzii carries very few transposable elements related to Ty elements in S. cerevisiae. Most of the mitochondrial genes were identified in eight contigs altogether spanning 25 kb for a predicted size of 29 kb. A significant match with the Kluyveromyces lactis linear DNA plasmid pGKL-1 encoded RF4 killer protein suggests that a related plasmid exists in S. servazzii. The sequences have been deposited with EMBL under the accession numbers AL402279-AL404848.

Genomic exploration of the hemiascomycetous yeasts: 9. Saccharomyces kluyveri.

The genome of Saccharomyces kluyveri was explored through 2528 random sequence tags with an average length of 981 bp. The complete nuclear ribosomal DNA unit was found to be 8656 bp in length. Sequences homologous to retroelements of the gypsy and copia types were identified as well as numerous solo long terminal repeats. We identified at least 1406 genes homologous to Saccharomyces cerevisiae open reading frames, with on average 58.1% and 72.4% amino acid identity and similarity, respectively. In addition, by comparison with completely sequenced genomes and the SwissProt database, we found 27 novel S. kluyveri genes. Most of these genes belong to pathways or have functions absent from S. cerevisiae, such as the catabolic pathway of purines or pyrimidines, melibiose fermentation, sorbitol utilization, or degradation of pollutants. The sequences are deposited in EMBL under the accession numbers AL404849-AL407376.

Genomic exploration of the hemiascomycetous yeasts: 10. Kluyveromyces thermotolerans.

A genomic exploration of Kluyveromyces thermotolerans was performed by random sequence tag (RST) analysis. We sequenced 2653 RSTs corresponding to inserts sequenced from both ends. We performed a systematic comparison with a complete set of proteins from Saccharomyces cerevisiae, other completely sequenced genomes and SwissProt. We identified six mitochondrial genes and 1358-1496 nuclear genes by comparison with S. cerevisiae. In addition, 25 genes were identified by comparison with other organisms. This corresponds to about 24% of the estimated gene content of this organism. A lower level of conservation is observed with orthologues to genes of S. cerevisiae previously classified as orphans. Gene order was found to be conserved between S. cerevisiae and K. thermotolerans in 56.5% of studied cases.

Genomic exploration of the hemiascomycetous yeasts: 12. Kluyveromyces marxianus var. marxianus.

As part of the comparative genomics project 'GENOLEVURES', we studied the Kluyveromyces marxianus var. marxianus strain CBS712 using a partial random sequencing strategy. With a 0.2 x genome equivalent coverage, we identified ca. 1300 novel genes encoding proteins, some containing spliceosomal introns with consensus splice sites identical to those of Saccharomyces cerevisiae, 28 tRNA genes, the whole rDNA repeat, and retrotransposons of the Ty1/2 family of S. cerevisiae with diverged Long Terminal Repeats. Functional classification of the K. marxianus genes, as well as the analysis of the paralogous gene families revealed few differences with respect to S. cerevisiae. Only 42 K. marxianus identified genes are without detectable homolog in the baker's yeast. However, we identified several genetic rearrangements between these two yeast species.

Genomic exploration of the hemiascomycetous yeasts: 13. Pichia angusta.

As part of a comparative genomics project on 13 hemiascomycetous yeasts, the Pichia angusta type strain was studied using a partial random sequencing strategy. With coverage of 0.5 genome equivalents, about 2500 novel protein-coding genes were identified, probably corresponding to more than half of the P. angusta protein-coding genes, 6% of which do not have homologs in Saccharomyces cerevisiae. Some of them contain one or two introns, on average three times shorter than those in S. cerevisiae. We also identified 28 tRNA genes, a few retrotransposons similar to Ty5 of S. cerevisiae, solo long terminal repeats, the whole ribosomal DNA cluster, and segments of mitochondrial DNA. The P. angusta sequences were deposited in EMBL under the accession numbers AL430961 to AL436044.