Structure of two melon regions reveals high microsynteny with sequenced plant species

In this study, two melon bacterial artificial chromosome (BAC) clones have been sequenced and annotated. BAC 1-21-10 spans 92 kb and contains the nsv locus conferring resistance to the Melon Necrotic Spot Virus (MNSV) in melon linkage group 11. BAC 13J4 spans 98 kb and belongs to a BAC contig containing resistance gene homologues, extending a previous sequenced region of 117 kb in linkage group 4. Both regions have microsyntenic relationships to the model plant species Arabidopsis thaliana, and to Medicago truncatula and Populus trichocarpa. The network of synteny found between melon and each of the sequenced genomes reflects the polyploid structure of Arabidopsis, Populus, and Medicago genomes due to whole genome duplications (WGD). A detailed analysis revealed that both melon regions have a lower relative syntenic quality with Arabidopsis (eurosid II) than when compared to Populus and Medicago (eurosid I). Although phylogenetically Cucurbitales seem to be closer to Fabales than to Malphigiales, synteny was higher between both melon regions and Populus. Presented data imply that the recently completed Populus genome sequence could preferentially be used to obtain positional information in melon, based on microsynteny. Wim Deleu and Víctor González contributed equally to this work. The nucleotide sequences of BACs 1-21-10 and 13J4 are available in the DDBJ/EMBL/GenBank databases under the accession numbers EF188258 and EF657230, respectively.     

Genome-Wide Analysis of Repetitive Elements in Papaya

Papaya (Carica papaya L.) is an important fruit crop cultivated in tropical and subtropical regions worldwide. A first draft of its genome sequence has been recently released. Together with Arabidopsis, rice, poplar, grapevine and other genomes in the pipeline, it represents a good opportunity to gain insight into the organization of plant genomes. Here we report a detailed analysis of repetitive elements in the papaya genome, including transposable elements (TEs), tandemly-arrayed sequences, and high copy number genes. These repetitive sequences account for ～56% of the papaya genome with TEs being the most abundant at 52%, tandem repeats at 1.3% and high copy number genes at 3%. Most common types of TEs are represented in the papaya genome with retrotransposons being the dominant class, accounting for 40% of the genome. The most prevalent retrotransposons are Ty3-gypsy (27.8%) and Ty1-copia (5.5%). Among the tandem repeats, microsatellites are the most abundant in number, but represent only 0.19% of the genome. Minisatellites and satellites are less abundant, but represent 0.68% and 0.43% of the genome, respectively, due to greater repeat length. Despite an overall smaller gene repertoire in papaya than many other angiosperms, a significant fraction of genes (>2%) are present in large gene families with copy number greater than 20. This repeat database clarified a major part of the papaya genome organization and partly explained the lower gene repertoire in papaya than in Arabidopsis.     

Exploratory integration of peanut genetic and physical maps and possible contributions from Arabidopsis

Arachis hypogaea is a widely cultivated crop both as an oilseed and protein source. The genomic analysis of Arachis species hitherto has been limited to the construction of genetic maps; the most comprehensive one contains 370 loci over 2,210 cM in length. However, no attempt has been made to analyze the physical structure of the peanut genome. To investigate the practicality of physical mapping in peanut, we applied a total of 117 oligonucleotide-based probes (overgos) derived from genetically mapped RFLP probes onto peanut BAC filters containing 182,784 peanut large-insert DNA clones in a multiplex experimental design; 91.5% of the overgos identified at least one BAC clone. In order to gain insights into the potential value of Arabidopsis genome sequence for studies in divergent species with complex genomes such as peanut, we employed 576 Arabidopsis-derived overgos selected on the basis of maximum homology to orthologous sequences in other plant taxa to screen the peanut BAC library. A total of 353 (61.3%) overgos detected at least one peanut BAC clone. This experiment represents the first steps toward the creation of a physical map in peanut and illustrates the potential value of leveraging information from distantly related species such as Arabidopsis for both practical applications such as comparative map-based cloning and shedding light on evolutionary relationships. We also evaluated the possible correlation between functional categories of Arabidopsis overgos and their success rates in hybridization to the peanut BAC library.Electronic Supplementary Material Supplementary material is available for this article at     

Significant microsynteny with new evolutionary highlights is detected between <Emphasis Type="Italic">Arabidopsis</Emphasis> and legume model plants despite the lack of macrosynteny

The increased amount of data produced by large genome sequencing projects allows scientists to carry out important syntenic studies to a great extent. Detailed genetic maps and entirely or partially sequenced genomes are compared, and macro- and microsyntenic relations can be determined for different species. In our study, the syntenic relationships between key legume plants and two model plants, Arabidopsis thaliana and Populus trichocarpa were investigated. The comparison of the map position of 172 gene-based Medicago sativa markers to the organization of homologous A. thaliana genes could not identify any sign of macrosynteny between the two genomes. A 276 kb long section of chromosome 5 of the model legume Medicago truncatula was used to investigate potential microsynteny with the other legume Lotus japonicus, as well as with Arabidopsis and Populus. Besides the overall correlation found between the legume plants, the comparison revealed several microsyntenic regions in the two more distant plants with significant resemblance. Despite the large phylogenetic distance, clear microsyntenic regions between Medicago and Arabidopsis or Populus were detected unraveling new intragenomic evolutionary relations in Arabidopsis. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Characterizing the walnut genome through analyses of BAC end sequences

Persian walnut (Juglans regia L.) is an economically important tree for its nut crop and timber. To gain insight into the structure and evolution of the walnut genome, we constructed two bacterial artificial chromosome (BAC) libraries, containing a total of 129,024 clones, from in vitro-grown shoots of J. regia cv. Chandler using the HindIII and MboI cloning sites. A total of 48,218 high-quality BAC end sequences (BESs) were generated, with an accumulated sequence length of 31.2?Mb, representing approximately 5.1% of the walnut genome. Analysis of repeat DNA content in BESs revealed that approximately 15.42% of the genome consists of known repetitive DNA, while walnut-unique repetitive DNA identified in this study constitutes 13.5% of the genome. Among the walnut-unique repetitive DNA, Julia SINE and JrTRIM elements represent the first identified walnut short interspersed element (SINE) and terminal-repeat retrotransposon in miniature (TRIM) element, respectively; both types of elements are abundant in the genome. As in other species, these SINEs and TRIM elements could be exploited for developing repeat DNA-based molecular markers in walnut. Simple sequence repeats (SSR) from BESs were analyzed and found to be more abundant in BESs than in expressed sequence tags. The density of SSR in the walnut genome analyzed was also slightly higher than that in poplar and papaya. Sequence analysis of BESs indicated that approximately 11.5% of the walnut genome represents a coding sequence. This study is an initial characterization of the walnut genome and provides the largest genomic resource currently available; as such, it will be a valuable tool in studies aimed at genetically improving walnut.     

The Complete Mitochondrial Genome of the Entomopathogenic Nematode Steinernema carpocapsae: Insights into Nematode Mitochondrial DNA Evolution and Phylogeny

We determined the complete sequence of the mitochondrial DNA of the entomopathogenic nematode Steinernema carpocapsae and analyzed its structure and composition as well as the secondary structures predicted for its tRNAs and rRNAs. Almost the complete genome has been amplified in one fragment with long PCR and sequenced using a shotgun strategy. The 13,925-bp genome contains genes for 2 rRNAs, 22 tRNAs, and 12 proteins and lacks an ORF encoding ATPase subunit 8. Four initiation codons were inferred, TTT, TTA, ATA, and ATT, most of the genes ended with TAA or TAG, and only two had a T as an incomplete stop codon. All predicted tRNAs showed the nonconventional secondary structure typical of Secernentea. Although we were able to fold the sequences of trnN, trnD, and trnC into more conventional cloverleaf structures after adding adjacent nucleotides, northern blot experiments showed that the nonstandard tRNAs are actually expressed. Phylogenetic and comparative analyses showed that the mitochondrial genome of S. carpocapsae is more closely related to the genomes of A. suum and C. elegans than to that of Strongyloides stercoralis. This finding does not support the phylogeny based on nuclear small subunit ribosomal DNA sequences previously published. This discrepancy may result from differential reproductive strategies and/or differential selective pressure acting on nuclear and mitochondrial genes. The distinctive characteristics observed among mitochondrial genomes of Secernentea may have arisen to counteract the deleterious effects of Muller’s ratchet, which is probably enhanced by the reproductive strategies and selective pressures referred to above. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Rafael Zardoya]     

Detection and correction of assembly errors of rice Nipponbare reference sequence

A complete and high‐quality genome reference sequence of an organism provides a solid foundation for a wide research community and determines the outcomes of relevant genomic, genetic, molecular and evolutionary research. Rice is an important food crop and a model plant for grasses, and therefore was the first chosen crop plant for whole genome sequencing. The genome of the japonica representative rice variety, Nipponbare, was sequenced using a gold standard, map‐based clone‐by‐clone strategy. However, although the Nipponbare reference sequence (RefSeq) has the best quality for existing crop genome sequences, it still contains many assembly errors and gaps. To improve the Nipponbare RefSeq, first a robust method is required to detect the hidden assembly errors. Through alignments between BAC‐end sequences (BESs) embedded in the Nipponbare bacterial artificial chromosome (BAC) physical map and the Nipponbare RefSeq, we detected locations on the Nipponbare RefSeq that were inversely matched with BESs and could therefore be candidates for spurious inversions of assembly. We performed further analysis of five potential locations and confirmed assembly errors at those locations; four of them, two on chr4 and two on chr11 of the Nipponbare RefSeq (IRGSP build 5), were found to be caused by reverse repetitive sequences flanking the locations. Our approach is effective in detecting spurious inversions in the Nipponbare RefSeq and can be applied for improving the sequence qualities of other genomes as well.     

Cloning and characterization of microRNAs from wheat (<Emphasis Type="Italic">Triticum aestivum</Emphasis> L.)

Background  

MicroRNAs (miRNAs) are a class of small, non-coding regulatory RNAs that regulate gene expression by guiding target mRNA cleavage or translational inhibition. So far, identification of miRNAs has been limited to a few model plant species, such as Arabidopsis, rice and Populus, whose genomes have been sequenced. Wheat is one of the most important cereal crops worldwide. To date, only a few conserved miRNAs have been predicted in wheat and the computational identification of wheat miRNAs requires the genome sequence, which is unknown.     

Intron-flanking EST–PCR markers: from genetic marker development to gene structure analysis in Rhododendron

With a long-term goal of constructing a linkage map of Rhododendron enriched with gene-specific markers, we utilized Rhododendron catawbiense ESTs for the development of high-efficiency (in terms of generating polymorphism frequency) PCR-based markers. Using the gene-sequence alignment between Rhododendron ESTs and the genomic sequences of Arabidopsis homologs, we developed ‘intron-flanking‘ EST–PCR-based primers that would anneal in conserved exon regions and amplify across the more highly diverged introns. These primers resulted in increased efficiency (61% vs. 13%; 4.7-fold) of polymorphism-detection compared with conventional EST–PCR methods, supporting the assumption that intron regions are more diverged than exons. Significantly, this study demonstrates that Arabidopsis genome database can be useful in developing gene-specific PCR-based markers for other non-model plant species for which the EST data are available but genomic sequences are not. The comparative analysis of intron sizes between Rhododendron and Arabidopsis (made possible in this study by aligning of Rhododendron ESTs with Arabidopsis genomic sequences and the sequencing of Rhododendron genomic PCR products) provides the first insight into the gene structure of Rhododendron. Electronic Supplementary Material Supplementary material is available for this article at     

BAHD superfamily of acyl-CoA dependent acyltransferases in <Emphasis Type="Italic">Populus</Emphasis> and <Emphasis Type="Italic">Arabidopsis</Emphasis>: bioinformatics and gene expression

Plant acyl-CoA dependent acyltransferases constitute a large specific protein superfamily, named BAHD. Using the conserved sequence motifs of BAHD members, we searched the genome sequences of Populus and Arabidopsis, and identified, respectively, 94- and 61-putative genes. Subsequently, we analyzed the phylogeny, gene structure, and chromosomal distribution of BAHD members of both species; then, we profiled expression patterns of BAHD genes by “in silico” northern- and microarray-analyses based on public databases, and by RT-PCR. While our genomic- and bioinformatic- analyses provided full sets of BAHD superfamily genes, and cleaned up a few existing annotation errors, importantly it led to our recognizing several unique Arabidopsis BAHD genes that inversely overlapped with their neighboring genes on the genome, and disclosing a potential natural anti-sense regulation for gene expressions. Systemic gene-expression profiling of BAHD members revealed distinct tissue-specific/preferential expression patterns, indicating their diverse biological functions. Our study affords a strong knowledge base for understanding BAHD members’ evolutionary relationships and gene functions implicated in plant growth, development and metabolism. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Genome-wide identification of NBS-encoding resistance genes in <Emphasis Type="Italic">Brassica rapa</Emphasis>

Nucleotide-binding site (NBS)-encoding resistance genes are key plant disease-resistance genes and are abundant in plant genomes, comprising up to 2% of all genes. The availability of genome sequences from several plant models enables the identification and cloning of NBS-encoding genes from closely related species based on a comparative genomics approach. In this study, we used the genome sequence of Brassica rapa to identify NBS-encoding genes in the Brassica genome. We identified 92 non-redundant NBS-encoding genes [30 CC-NBS-LRR (CNL) and 62 TIR-NBS-LRR (TNL) genes] in approximately 100 Mbp of B. rapa euchromatic genome sequence. Despite the fact that B. rapa has a significantly larger genome than Arabidopsis thaliana due to a recent whole genome triplication event after speciation, B. rapa contains relatively small number of NBS-encoding genes compared to A. thaliana, presumably because of deletion of redundant genes related to genome diploidization. Phylogenetic and evolutionary analyses suggest that relatively higher relaxation of selective constraints on the TNL group after the old duplication event resulted in greater accumulation of TNLs than CNLs in both Arabidopsis and Brassica genomes. Recent tandem duplication and ectopic deletion are likely to have played a role in the generation of novel Brassica lineage-specific resistance genes.     

Patterns of tandem repetition in plant whole genome assemblies

Tandem repeats often confound large genome assemblies. A survey of tandemly arrayed repetitive sequences was carried out in whole genome sequences of the green alga Chlamydomonas reinhardtii, the moss Physcomitrella patens, the monocots rice and sorghum, and the dicots Arabidopsis thaliana, poplar, grapevine, and papaya, in order to test how these assemblies deal with this fraction of DNA. Our results suggest that plant genome assemblies preferentially include tandem repeats composed of shorter monomeric units (especially dinucleotide and 9–30-bp repeats), while higher repetitive units pose more difficulties to assemble. Nevertheless, notwithstanding that currently available sequencing technologies struggle with higher arrays of repeated DNA, major well-known repetitive elements including centromeric and telomeric repeats as well as high copy-number genes, were found to be reasonably well represented. A database including all tandem repeat sequences characterized here was created to benefit future comparative genomic analyses.     

Genome Sequences of Populus tremula Chloroplast and Mitochondrion: Implications for Holistic Poplar Breeding

Complete Populus genome sequences are available for the nucleus (P. trichocarpa; section Tacamahaca) and for chloroplasts (seven species), but not for mitochondria. Here, we provide the complete genome sequences of the chloroplast and the mitochondrion for the clones P. tremula W52 and P. tremula x P. alba 717-1B4 (section Populus). The organization of the chloroplast genomes of both Populus clones is described. A phylogenetic tree constructed from all available complete chloroplast DNA sequences of Populus was not congruent with the assignment of the related species to different Populus sections. In total, 3,024 variable nucleotide positions were identified among all compared Populus chloroplast DNA sequences. The 5-prime part of the LSC from trnH to atpA showed the highest frequency of variations. The variable positions included 163 positions with SNPs allowing for differentiating the two clones with P. tremula chloroplast genomes (W52, 717-1B4) from the other seven Populus individuals. These potential P. tremula-specific SNPs were displayed as a whole-plastome barcode on the P. tremula W52 chloroplast DNA sequence. Three of these SNPs and one InDel in the trnH-psbA linker were successfully validated by Sanger sequencing in an extended set of Populus individuals. The complete mitochondrial genome sequence of P. tremula is the first in the family of Salicaceae. The mitochondrial genomes of the two clones are 783,442 bp (W52) and 783,513 bp (717-1B4) in size, structurally very similar and organized as single circles. DNA sequence regions with high similarity to the W52 chloroplast sequence account for about 2% of the W52 mitochondrial genome. The mean SNP frequency was found to be nearly six fold higher in the chloroplast than in the mitochondrial genome when comparing 717-1B4 with W52. The availability of the genomic information of all three DNA-containing cell organelles will allow a holistic approach in poplar molecular breeding in the future.     

Nuclear DNA content of some important plant species

Nuclear DNA contents of more than 100 important plant species were measured by flow cytometry of isolated nuclei stained with propidium iodide.Arabidopsis exhibits developmentally regulated multiploidy and has a 2C nuclear DNA content of 0.30 pg (145 Mbp/1C), twice the value usually cited. The 2C value for rice is only about three times that ofArabidopsis. Tomato has a 2C value of about 2.0 pg, larger than commonly cited. This survey identified several horticultural crops in a variety of families with genomes only two or three times as large asArabidopsis; these include several fruit trees (a pricot, cherry, mango, orange, papaya, and peach). The small genome sizes of rice and the horticultural plants should facilitate molecular studies of these crops.     

Alignment of a <Emphasis Type="Italic">Salix</Emphasis> linkage map to the <Emphasis Type="Italic">Populus</Emphasis> genomic sequence reveals macrosynteny between willow and poplar genomes

Poplars (genus Populus) and willows (genus Salix) are members of the Salicaceae, a family of catkin-bearing trees, shrubs and sub-shrubs. Poplar is considered the model system for biological studies in trees and considerable genetic and genomic resources have become available in recent years. The transfer of information to research studies in willow, for which fewer resources are currently available, would be highly beneficial. However, the extent of conservation between poplar and willow genomes has not yet been extensively studied. To address this, we have constructed a linkage map of willow based on a large mapping population derived from a cross between two Salix viminalis × (S. viminalis × S. schwerinii) hybrid sibs, and aligned this to the publicly available poplar genome sequence. A set of genome-wide, expressed poplar sequences was selected and used to design primer sets that efficiently amplified homeologous regions in willow. Direct sequencing of the willow products confirmed homology with the poplar target in the majority of instances and allowed identification of single nucleotide polymorphisms (SNPs) that were used to map these loci. In total, 202 amplified fragment length polymorphisms (AFLPs), 75 microsatellites and 79 SNPs were used to construct a willow consensus map that spanned 1,856.7 cM with an average interval between markers of 6.3 cM. Poplar sequences homologous to those of the mapped willow microsatellite loci were identified and used in addition to the SNP markers to putatively align all but two minor linkage groups to the poplar genome sequence. A high degree of macrosynteny was revealed.     

Processed Pseudogenes, Processed Genes, and Spontaneous Mutations in the Arabidopsis Genome

We identified 411 processed sequences in the Arabidopsis thaliana genome based on the fact that they have lost their intron(s) and have a length that is at least 95% of the length of the gene that gave rise to them. These sequences were generated by 230 different genes and clearly originated from retrotranspositons events because most of them (91%) have a poly(A)-tail. They are composed of 376 sequences with frame shifts and/or premature stop codons (processed pseudogenes) and 35 sequences without disablements (processed genes). Eleven of these processed genes are likely functional retrotransposed genes because they have low Ka/Ks ratios and high Ks values, and their sequences match numerous Arabidopsis ESTs. Processed sequences are mostly randomly distributed in the Arabidopsis genome and their rate of accumulation has steadily been decreasing since it peaked some 50 MYA. In contrast with the situation observed in mammals, the processed sequences found in the Arabidopsis genome originate from genes with high copy numbers and not from highly expressed genes. The patterns of spontaneous mutations in Arabidopsis are slightly different than those of mammals but are similar to those observed in Drosophila. This suggests that methylated cytosine deamination is less frequent in Arabidopsis than in mammals. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Juergen Brosius]     

Comparative genomic mapping between a 754 kb region flanking DREB1A in Arabidopsis thaliana and maize

Comparative mapping between model plant species for which the complete genome sequence is known and crop species has been suggested as a new strategy for the isolation of agronomically valuable genes. In this study, we tested whether comparative mapping between Arabidopsisand maize of a small region (754 kb) surrounding the DREB1A gene in Arabidopsis could lead to the identification of an orthologous region in maize containing the DREB1A homologue. The genomic sequence information available for Arabidopsis allowed for the selection of conserved, low-copy genes that were used for the identification of maize homologues in a large EST database. In total, 17 maize homologues were mapped. A second BLAST comparison of these genes to the recently completed Arabidopsis sequence revealed that 15 homologues are likely to be orthologous as the highest similarity score was obtained either with the original Arabidopsis gene or with a highly similar Arabidopsis gene localized on a duplication of the investigated region on chromosome 5. The map position of these genes showed a significant degree of orthology with the Arabidopsis region. Nevertheless, extensive duplications and rearrangements in the Arabidopsisand maize genomes as well as the evolutionary distance between Arabidopsis and maize make it unlikely that orthology and collinearity between these two species are sufficient to aid gene prediction and cloning in maize.