Sander A. Peters

Complete genome sequence of Lactobacillus plantarum WCFS1

The 3,308,274-bp sequence of the chromosome of Lactobacillus plantarum strain WCFS1, a single colony isolate of strain NCIMB8826 that was originally isolated from human saliva, has been determined, and contains 3,052 predicted protein-encoding genes. Putative biological functions could be assigned to 2,120 (70%) of the predicted proteins. Consistent with the classification of L. plantarum as a facultative heterofermentative lactic acid bacterium, the genome encodes all enzymes required for the glycolysis and phosphoketolase pathways, all of which appear to belong to the class of potentially highly expressed genes in this organism, as was evident from the codon-adaptation index of individual genes. Moreover, L. plantarum encodes a large pyruvate-dissipating potential, leading to various end-products of fermentation. L. plantarum is a species that is encountered in many different environmental niches, and this flexible and adaptive behavior is reflected by the relatively large number of regulatory and transport functions, including 25 complete PTS sugar transport systems. Moreover, the chromosome encodes >200 extracellular proteins, many of which are predicted to be bound to the cell envelope. A large proportion of the genes encoding sugar transport and utilization, as well as genes encoding extracellular functions, appear to be clustered in a 600-kb region near the origin of replication. Many of these genes display deviation of nucleotide composition, consistent with a foreign origin. These findings suggest that these genes, which provide an important part of the interaction of L. plantarum with its environment, form a lifestyle adaptation region in the chromosome.

Exploring genetic variation in the tomato (Solanum section Lycopersicon) clade by whole-genome sequencing.

We explored genetic variation by sequencing a selection of 84 tomato accessions and related wild species representative of the Lycopersicon, Arcanum, Eriopersicon and Neolycopersicon groups, which has yielded a huge amount of precious data on sequence diversity in the tomato clade. Three new reference genomes were reconstructed to support our comparative genome analyses. Comparative sequence alignment revealed group-, species- and accession-specific polymorphisms, explaining characteristic fruit traits and growth habits in the various cultivars. Using gene models from the annotated Heinz 1706 reference genome, we observed differences in the ratio between non-synonymous and synonymous SNPs (dN/dS) in fruit diversification and plant growth genes compared to a random set of genes, indicating positive selection and differences in selection pressure between crop accessions and wild species. In wild species, the number of single-nucleotide polymorphisms (SNPs) exceeds 10 million, i.e. 20-fold higher than found in most of the crop accessions, indicating dramatic genetic erosion of crop and heirloom tomatoes. In addition, the highest levels of heterozygosity were found for allogamous self-incompatible wild species, while facultative and autogamous self-compatible species display a lower heterozygosity level. Using whole-genome SNP information for maximum-likelihood analysis, we achieved complete tree resolution, whereas maximum-likelihood trees based on SNPs from ten fruit and growth genes show incomplete resolution for the crop accessions, partly due to the effect of heterozygous SNPs. Finally, results suggest that phylogenetic relationships are correlated with habitat, indicating the occurrence of geographical races within these groups, which is of practical importance for Solanum genome evolution studies.

Complete Sequences of Four Plasmids of Lactococcus lactis subsp. cremoris SK11 Reveal Extensive Adaptation to the Dairy Environment

ABSTRACT Lactococcus lactis strains are known to carry plasmids encoding industrially important traits. L. lactis subsp. cremoris SK11 is widely used by the dairy industry in cheese making. Its complete plasmid complement was sequenced and found to contain the plasmids pSK11A (10,372 bp), pSK11B (13,332 bp), pSK11L (47,165 bp), and pSK11P (75,814 bp). Six highly homologous repB-containing replicons were found, all belonging to the family of lactococcal theta-type replicons. Twenty-three complete insertion sequence elements segment the plasmids into numerous modules, many of which can be identified as functional units or containing functionally related genes. Plasmid-encoded functions previously known to reside on L. lactis SK11 plasmids were now mapped in detail, e.g., lactose utilization (lacR-lacABCDFEGX), the proteolytic system (prtM-prtP, pepO, pepF), and the oligopeptide permease system (oppDFBCA). Newly identified plasmid-encoded functions could facilitate the uptake of various cations, while the pabA and pabB genes could be essential for folate biosynthesis. A competitive advantage could be obtained by using the putative flavin adenine dinucleotide-dependent d-lactate dehydrogenase and oxalate:formate antiporter for enhanced ATP synthesis, while the activity of the predicted α-acetolactate decarboxylase may contribute to the formation of an additional electron sink. Various stress response proteins are plasmid encoded, which could enhance strain robustness. A substantial number of these “adaptation” genes have not been described before on L. lactis plasmids. Moreover, several genes were identified for the first time in L. lactis, possibly reflecting horizontal gene transfer.

Cross-Species Bacterial Artificial Chromosome–Fluorescence in Situ Hybridization Painting of the Tomato and Potato Chromosome 6 Reveals Undescribed Chromosomal Rearrangements

Ongoing genomics projects of tomato (Solanum lycopersicum) and potato (S. tuberosum) are providing unique tools for comparative mapping studies in Solanaceae. At the chromosomal level, bacterial artificial chromosomes (BACs) can be positioned on pachytene complements by fluorescence in situ hybridization (FISH) on homeologous chromosomes of related species. Here we present results of such a cross-species multicolor cytogenetic mapping of tomato BACs on potato chromosomes 6 and vice versa. The experiments were performed under low hybridization stringency, while blocking with Cot-100 was essential in suppressing excessive hybridization of repeat signals in both within-species FISH and cross-species FISH of tomato BACs. In the short arm we detected a large paracentric inversion that covers the whole euchromatin part with breakpoints close to the telomeric heterochromatin and at the border of the short arm pericentromere. The long arm BACs revealed no deviation in the colinearity between tomato and potato. Further comparison between tomato cultivars Cherry VFNT and Heinz 1706 revealed colinearity of the tested tomato BACs, whereas one of the six potato clones (RH98-856-18) showed minor putative rearrangements within the inversion. Our results present cross-species multicolor BAC–FISH as a unique tool for comparative genetic studies across Solanum species.

Complete genome sequence and taxonomic position of anguillid herpesvirus 1.

Eel herpesvirus or anguillid herpesvirus 1 (AngHV1) frequently causes disease in freshwater eels. The complete genome sequence of AngHV1 and its taxonomic position within the family Alloherpesviridae were determined. Shotgun sequencing revealed a 249 kbp genome including an 11 kbp terminal direct repeat that contains 7 of the 136 predicted protein-coding open reading frames. Twelve of these genes are conserved among other members of the family Alloherpesviridae and another 28 genes have clear homologues in cyprinid herpesvirus 3. Phylogenetic analyses based on amino acid sequences of five conserved genes, including the ATPase subunit of the terminase, confirm the position of AngHV1 within the family Alloherpesviridae, where it is most closely related to the cyprinid herpesviruses. Our analyses support a recent proposal to subdivide the family Alloherpesviridae into two sister clades, one containing AngHV1 and the cyprinid herpesviruses and the other containing Ictalurid herpesvirus 1 and the ranid herpesviruses.

Structural homology in the Solanaceae: analysis of genomic regions in support of synteny studies in tomato, potato and pepper

We have analysed the structural homology in euchromatin regions of tomato, potato and pepper with special attention for the long arm of chromosome 2 (2L). Molecular organization and colinear junctions were delineated using multi-color BAC FISH analysis and comparative sequence alignment. We found large-scale rearrangements including inversions and segmental translocations that were not reported in previous comparative studies. Some of the structural rearrangements are specific for the tomato clade, and differentiate tomato from potato, pepper and other Solanaceous species. Although local gene vicinity is largely preserved, there are many small-scale synteny perturbations. Gene adjacency in the aligned segments was frequently disrupted for 47% of the ortholog pairs as a result of gene and LTR retrotransposon insertions, and occasionally by single gene inversions and translocations. Our data also suggests that long distance intra-chromosomal rearrangements and local gene rearrangements have evolved frequently during speciation in the Solanum genus, and that small changes are more prevalent than large-scale differences. The occurrence of sonata and harbinger transposable elements and other repeats near or at junction breaks is considered in the light of repeat-mediated rearrangements and a reconstruction scenario for an ancestral 2L topology is discussed.

Solanum lycopersicum cv. Heinz 1706 chromosome 6: distribution and abundance of genes and retrotransposable elements

We studied the physical and genetic organization of chromosome 6 of tomato (Solanum lycopersicum) cv. Heinz 1706 by combining bacterial artificial chromosome (BAC) sequence analysis, high-information-content fingerprinting, genetic analysis, and BAC-fluorescent in situ hybridization (FISH) mapping data. The chromosome positions of 81 anchored seed and extension BACs corresponded in most cases with the linear marker order on the high-density EXPEN 2000 linkage map. We assembled 25 BAC contigs and eight singleton BACs spanning 2.0 Mb of the short-arm euchromatin, 1.8 Mb of the pericentromeric heterochromatin and 6.9 Mb of the long-arm euchromatin. Sequence data were combined with their corresponding genetic and pachytene chromosome positions into an integrated map that covers approximately a third of the chromosome 6 euchromatin and a small part of the pericentromeric heterochromatin. We then compared physical length (Mb), genetic (cM) and chromosome distances (microm) for determining gap sizes between contigs, revealing relative hot and cold spots of recombination. Through sequence annotation we identified several clusters of functionally related genes and an uneven distribution of both gene and repeat sequences between heterochromatin and euchromatin domains. Although a greater number of the non-transposon genes were located in the euchromatin, the highly repetitive (22.4%) pericentromeric heterochromatin displayed an unexpectedly high gene content of one gene per 36.7 kb. Surprisingly, the short-arm euchromatin was relatively rich in repeats as well, with a repeat content of 13.4%, yet the ratio of Ty3/Gypsy and Ty1/Copia retrotransposable elements across the chromosome clearly distinguished euchromatin (2:3) from heterochromatin (3:2).

SNP markers retrieval for a non-model species: a practical approach.

BackgroundSNP (Single Nucleotide Polymorphism) markers are rapidly becoming the markers of choice for applications in breeding because of next generation sequencing technology developments. For SNP development by NGS technologies, correct assembly of the huge amounts of sequence data generated is essential. Little is known about assemblers performance, especially when dealing with highly heterogeneous species that show a high genome complexity and what the possible consequences are of differences in assemblies on SNP retrieval. This study tested two assemblers (CAP3 and CLC) on 454 data from four lily genotypes and compared results with respect to SNP retrieval.ResultsCAP3 assembly resulted in higher numbers of contigs, lower numbers of reads per contig, and shorter average read lengths compared to CLC. Blast comparisons showed that CAP3 contigs were highly redundant. Contrastingly, CLC in rare cases combined paralogs in one contig. Redundant and chimeric contigs may lead to erroneous SNPs. Filtering for redundancy can be done by blasting selected SNP markers to the contigs and discarding all the SNP markers that show more than one blast hit. Results on chimeric contigs showed that only four out of 2,421 SNP markers were selected from chimeric contigs.ConclusionIn practice, CLC performs better in assembling highly heterogeneous genome sequences compared to CAP3, and consequently SNP retrieval is more efficient. Additionally a simple flow scheme is suggested for SNP marker retrieval that can be valid for all non-model species.

Genomic organisation of the Mal d 1 gene cluster on linkage group 16 in apple

European populations exhibit progressive sensitisation to food allergens, and apples are one of the foods for which sensitisation is observed most frequently. Apple cultivars vary greatly in their allergenic characteristics, and a better understanding of the genetic basis of low allergenicity may therefore allow allergic individuals to increase their fruit intake. Mal d 1 is considered to be a major apple allergen, and this protein is encoded by the most complex allergen gene family. Not all Mal d 1 members are likely to be involved in allergenicity. Therefore, additional knowledge about the existence and characteristics of the different Mal d 1 genes is required. In the present study, we investigated the genomic organisation of the Mal d 1 gene cluster in linkage group 16 of apple through the sequencing of two bacterial artificial chromosome clones. The results provided new information on the composition of this family with respect to the number and orientation of functional and pseudogenes and their physical distances. The results were compared with the apple and peach genome sequences that have recently been made available. A broad analysis of the whole apple genome revealed the presence of new genes in this family, and a complete list of the observed Mal d 1 genes is supplied. Thus, this study provides an important contribution towards a better understanding of the genetics of the Mal d 1 family and establishes the basis for further research on allelic diversity among cultivars in relation to variation in allergenicity.

TOPAAS, a Tomato and Potato Assembly Assistance System for Selection and Finishing of Bacterial Artificial Chromosomes

We have developed the software package Tomato and Potato Assembly Assistance System (TOPAAS), which automates the assembly and scaffolding of contig sequences for low-coverage sequencing projects. The order of contigs predicted by TOPAAS is based on read pair information; alignments between genomic, expressed sequence tags, and bacterial artificial chromosome (BAC) end sequences; and annotated genes. The contig scaffold is used by TOPAAS for automated design of nonredundant sequence gap-flanking PCR primers. We show that TOPAAS builds reliable scaffolds for tomato (Solanum lycopersicum) and potato (Solanum tuberosum) BAC contigs that were assembled from shotgun sequences covering the target at 6- to 8-fold coverage. More than 90% of the gaps are closed by sequence PCR, based on the predicted ordering information. TOPAAS also assists the selection of large genomic insert clones from BAC libraries for walking. For this, tomato BACs are screened by automated BLAST analysis and in parallel, high-density nonselective amplified fragment length polymorphism fingerprinting is used for constructing a high-resolution BAC physical map. BLAST and amplified fragment length polymorphism analysis are then used together to determine the precise overlap. Assembly onto the seed BAC consensus confirms the BACs are properly selected for having an extremely short overlap and largest extending insert. This method will be particularly applicable where related or syntenic genomes are sequenced, as shown here for the Solanaceae, and potentially useful for the monocots Brassicaceae and Leguminosea.