Genome-wide BAC-end sequencing of Musa acuminata DH Pahang reveals further insights into the genome organization of banana

Banana and plantain (Musa spp.) are grown in more than 120 countries in tropical and subtropical regions and constitute an important staple food for millions of people. A Musa acuminata ssp. malaccencis DH Pahang bacterial artificial chromosome (BAC) library (MAMB) was submitted for BAC-end sequencing. MAMB consists of 23,040 clones, with a 140-kbp average insert size, accounting for a five times coverage of the banana genome. A total of 46,080 reads were generated, and 42,750 (92.8%) high-quality sequences were obtained after trimming for vector and quality. Analysis of these data shows a GC content of 41.39%, whereas interspersed repeats comprise 32.3%. The most common repeated sequences found show homology to ribosomal RNA genes, particularly 18S rRNA, while the Ty3/gypsy type monkey retrotransposon is the most common retro element. The sequence data were used to generate a banana-specific repeat library containing 54 new repetitive elements which accounted for 11.86% of the total nucleotides. Simple sequence repeats represent 0.7% of the sequence data and allowed the identification of 2,455 potentially useful marker sites. Functional annotation identified 2,705 sequences that could code for proteins of known function. Microsynteny analysis shows a higher number of co-linear matches to Oryza sativa, in contrast to Arabidopsis thaliana. This database of BAC-end sequences is useful for the assembly of the complete banana genome sequence and is important for identification in functional genomics experiments.     

Evolutionary history of mitogen-activated protein kinase (MAPK) genes in Lotus,Medicago, and Phaseolus

Mitogen-Activated Protein Kinase (MAPK) genes encode proteins that mediate various signaling pathways associated with biotic and abiotic stress responses in eukaryotes. The MAPK genes form a 3-tier signal transduction cascade between cellular stimuli and physiological responses. Recent identification of soybean MAPKs and availability of genome sequences from other legume species allowed us to identify their MAPK genes. The main objectives of this study were to identify MAPKs in 3 legume species, Lotus japonicus, Medicago truncatula, and Phaseolus vulgaris, and to assess their phylogenetic relationships. We used approaches in comparative genomics for MAPK gene identification and named the newly identified genes following Arabidopsis MAPK nomenclature model. We identified 19, 18, and 15 MAPKs and 7, 4, and 9 MAPKKs in the genome of Lotus japonicus, Medicago truncatula, and Phaseolus vulgaris, respectively. Within clade placement of MAPKs and MAPKKs in the 3 legume species were consistent with those in soybean and Arabidopsis. Among 5 clades of MAPKs, 4 founder clades were consistent to MAPKs of other plant species and orthologs of MAPK genes in the fifth clade-"Clade E" were consistent with those in soybean. Our results also indicated that some gene duplication events might have occurred prior to eudicot-monocot divergence. Highly diversified MAPKs in soybean relative to those in 3 other legume species are attributable to the polyploidization events in soybean. The identification of the MAPK genes in the legume species is important for the legume crop improvement; and evolutionary relationships and functional divergence of these gene members provide insights into plant genome evolution.     

Construction of a BAC library of Korean ginseng and initial analysis of BAC-end sequences

We estimated the genome size of Korean ginseng ( Panax ginseng C.A. Meyer), a medicinal herb, constructed a Hin dIII BAC library, and analyzed BAC-end sequences to provide an initial characterization of the library. The 1C nuclear DNA content of Korean ginseng was estimated to be 3.33 pg (3.12×10³ Mb). The BAC library consists of 106,368 clones with an average size of 98.61 kb, amounting to 3.34 genome equivalents. Sequencing of 2167 BAC clones generated 2492 BAC-end sequences with an average length of 400 bp. Analysis using BLAST and motif searches revealed that 10.2%, 20.9% and 3.8% of the BAC-end sequences contained protein-coding regions, transposable elements and microsatellites, respectively. A comparison of the functional categories represented by the protein-coding regions found in BAC-end sequences with those of Arabidopsis revealed that proteins pertaining to energy metabolism, subcellular localization, cofactor requirement and transport facilitation were more highly represented in the P. ginseng sample. In addition, a sequence encoding a glucosyltransferase-like protein implicated in the ginsenoside biosynthesis pathway was also found. The majority of the transposable element sequences found belonged to the gypsy type (67.6%), followed by copia (11.7%) and LINE (8.0%) retrotransposons, whereas DNA transposons accounted for only 2.1% of the total in our sequence sample. Higher levels of transposable elements than protein-coding regions suggest that mobile elements have played an important role in the evolution of the genome of Korean ginseng, and contributed significantly to its complexity. We also identified 103 microsatellites with 3–38 repeats in their motifs. The BAC library and BAC-end sequences will serve as a useful resource for physical mapping, positional cloning and genome sequencing of P. ginseng.Electronic Supplementary Material Supplementary material is available in the online version of this article at Communicated by M.-A. Grandbastien     

Construction, characterization, and preliminary BAC-end sequencing analysis of a bacterial artificial chromosome library of white clover (Trifolium repens L.).

White clover (Trifolium repens L.) is a forage legume widely used in combination with grass in pastures because of its ability to fix nitrogen. We have constructed a bacterial artificial chromosome (BAC) library of an advanced breeding line of white clover. The library contains 37 248 clones with an average insert size of approximately 85 kb, representing an approximate 3-fold coverage of the white clover genome based on an estimated genome size of 960 Mb. The BAC library was pooled and screened by polymerase chain reaction (PCR) amplification using both white clover microsatellites and PCR-based markers derived from Medicago truncatula, resulting in an average of 6 hits per marker; this supports the estimated 3-fold genome coverage in this allotetraploid species. PCR-based screening of 766 clones with a multiplex set of chloroplast primers showed that only 0.5% of BAC clones contained chloroplast-derived inserts. The library was further evaluated by sequencing both ends of 724 of the clover BACs. These were analysed with respect to their sequence content and their homology to the contents of a range of plant gene, expressed sequence tag, and repeat element databases. Forty-three microsatellites were discovered in the BAC-end sequences (BESs) and investigated as potential genetic markers in white clover. The BESs were also compared with the partially sequenced genome of the model legume M. truncatula with the specific intention of identifying putative comparative-tile BACs, which represent potential regions of microsynteny between the 2 species; 14 such BACs were discovered. The results suggest that a large-scale BAC-end sequencing strategy has the potential to anchor a significant proportion of the genome of white clover onto the gene-space sequence of M. truncatula.     

Sequence variation in dispersed repetitive sequences in Saccharomyces cerevisiae

Two new dispersed repetitive DNA sequences related to the transposable element Tyl have been isolated from the genome of Saccharomyces cerevisiae. One sequence, designated Tyl-17, is present at about six copies per haploid genome, and one copy is located approximately 1000 base-pairs from the LEU2 locus on chromosome III. Tyl-17 is about the same size as Tyl (Cameron et al., 1979) and is flanked by δ sequences, but differs from Tyl by the presence of two large substitutions representing about 50% of the sequence. Tyl and Tyl-17 are found in a ‘head-to-head’ array in at least one cloned region of the yeast genome. Another sequence, designated Tyl-161, is situated about 9000 base-pairs from the PGK locus of chromosome III, and is structurally identical to Tyl except for the presence of a 1200 base-pair insertion near one end of the sequence element.     

Begin at the beginning: A BAC-end view of the passion fruit (Passiflora) genome

Background

The passion fruit (Passiflora edulis) is a tropical crop of economic importance both for juice production and consumption as fresh fruit. The juice is also used in concentrate blends that are consumed worldwide. However, very little is known about the genome of the species. Therefore, improving our understanding of passion fruit genomics is essential and to some degree a pre-requisite if its genetic resources are to be used more efficiently. In this study, we have constructed a large-insert BAC library and provided the first view on the structure and content of the passion fruit genome, using BAC-end sequence (BES) data as a major resource.

Results

The library consisted of 82,944 clones and its levels of organellar DNA were very low. The library represents six haploid genome equivalents, and the average insert size was 108 kb. To check its utility for gene isolation, successful macroarray screening experiments were carried out with probes complementary to eight Passiflora gene sequences available in public databases. BACs harbouring those genes were used in fluorescent in situ hybridizations and unique signals were detected for four BACs in three chromosomes (n = 9). Then, we explored 10,000 BES and we identified reads likely to contain repetitive mobile elements (19.6% of all BES), simple sequence repeats and putative proteins, and to estimate the GC content (~42%) of the reads. Around 9.6% of all BES were found to have high levels of similarity to plant genes and ontological terms were assigned to more than half of the sequences analysed (940). The vast majority of the top-hits made by our sequences were to Populus trichocarpa (24.8% of the total occurrences), Theobroma cacao (21.6%), Ricinus communis (14.3%), Vitis vinifera (6.5%) and Prunus persica (3.8%).

Conclusions

We generated the first large-insert library for a member of Passifloraceae. This BAC library provides a new resource for genetic and genomic studies, as well as it represents a valuable tool for future whole genome study. Remarkably, a number of BAC-end pair sequences could be mapped to intervals of the sequenced Arabidopsis thaliana, V. vinifera and P. trichocarpa chromosomes, and putative collinear microsyntenic regions were identified.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-816) contains supplementary material, which is available to authorized users.     

Analysis of a <Emphasis Type="Italic">c</Emphasis><Subscript>0</Subscript><Emphasis Type="Italic">t-1</Emphasis> library enables the targeted identification of minisatellite and satellite families in <Emphasis Type="Italic">Beta vulgaris</Emphasis>

Background  

Repetitive DNA is a major fraction of eukaryotic genomes and occurs particularly often in plants. Currently, the sequencing of the sugar beet (Beta vulgaris) genome is under way and knowledge of repetitive DNA sequences is critical for the genome annotation. We generated a c ₀ t-1 library, representing highly to moderately repetitive sequences, for the characterization of the major B. vulgaris repeat families. While highly abundant satellites are well-described, minisatellites are only poorly investigated in plants. Therefore, we focused on the identification and characterization of these tandemly repeated sequences.     

The genome of Zea mays,its organization and homology to related grasses

The pattern of genome organization of Zea mays has been analyzed, and the relationship of maize to possible progenitor species assessed by DNADNA hybridization. Reassociation of 470 and 1,350 bp fragments of maize DNA to various C₀t values demonstrates that the genome is composed of 3 major kinetic classes: highly repetitive, mid-repetitive, and unique. Mini-C₀t curves of the repetitive sequences at short fragment length indicate that the highly repetitive sequence class is 20% of the genome and is present at an average reiteration frequency of 800,000 copies; the mid-repetitive sequence class is 40% of the genome and is present at an average reiteration frequency of 1,000 copies. Thermal denaturation studies show that the highly repetitive sequences are 12% divergent and mid-repetitive sequences are 6% divergent. Most of the genome is organized in two interspersion patterns. One, approximately one-third of the genome, is composed of unique sequences of average length 2,100 bp interspersed with mid-repetitive sequences; the other, also one-third of the genome, is mid-repetitive sequences interspersed with highly repetitive sequences. The repetitive sequences are 500 to 1,000 bp by electron microscopic measurement. The remaining third of the genome is unique sequences farther than 5,000 bp from a palindromic or repetitive sequence. Hybridization of maize DNA from Midwestern Dent to popcorn and related grasses indicates that both the unique and repetitive sequence elements have diverged. Teosinte and popcorn are approximately equally divergent from Midwestern Dent whereas Tripsacum is much more divergent. The divergence times calculated from the depression of T_m in heterologous duplexes indicate that the divergence within Zea mays and between maize and near relatives is at least an order of magnitude greater than expected. This high degree of divergence may reflect the pressures of domestication of maize.     

Origin of the main class of repetitive DNA within selected Pennisetum species

In an attempt to identify relationships among genomes of the allotetraploid Pennisetum purpureum Schumach and closely related Pennisetum species with which it can be successfully hybridized, repetitive DNA sequences were examined. Digestion with KpnI revealed two highly repetitive fragments of 140 by and 160 bp. The possibility that these sequences could be used as genome markers was investigated. Average sequences were determined for the 140 by and 160 by KpnI families from P. purpureum and P. squamulatum Fresen. Average sequences (based upon four or five repeats) were determined for the P. glaucum (L.) R. Br. 140 by KpnI family and the diploid P. hohenackeri Hochst. ex Steud. 160 bp KpnI family. The average sequences of the 160 by KpnI families in P. purpureum and P. squamulatum differ by only nine bases. The 140 by KpnI families of the three related species, P. purpureum, P. squamulantum, and P. glaucum are nearly identical, and thus likely represent a recent divergence from a common progenitor or a common genome. Each repetitive sequence may contain internal duplications, which probably diverged following amplification of the original sequence. The 140 by KpnI repeat probably evolved from the 160 by KpnI repeat since the missing 18 by segment is part of the internal duplication that is otherwise conserved in the subrepeats. Tandemly arrayed repetitive sequences in plants are likely to be composed of subrepeats which have been duplicated and amplified.     

Development and analysis of a germline BAC resource for the sea lamprey, a vertebrate that undergoes substantial chromatin diminution

Over the last several years, the sea lamprey (Petromyzon marinus) has grown substantially as a model for understanding the evolutionary fundaments and capacity of vertebrate developmental and genome biology. Recent work on the lamprey genome has resulted in a preliminary assembly of the lamprey genome and led to the realization that nearly all somatic cell lineages undergo extensive programmed rearrangements. Here we describe the development of a bacterial artificial chromosome (BAC) resource for lamprey germline DNA and use sequence information from this resource to probe the subchromosomal structure of the lamprey genome. The arrayed germline BAC library represents ∼10× coverage of the lamprey genome. Analyses of BAC-end sequences reveal that the lamprey genome possesses a high content of repetitive sequences (relative to human), which show strong clustering at the subchromosomal level. This pattern is not unexpected given that the sea lamprey genome is dispersed across a large number of chromosomes (n ∼ 99) and suggests a low-copy DNA targeting strategy for efficiently generating informative paired-BAC-end linkages from highly repetitive genomes. This library therefore represents a new and biologically informed resource for understanding the structure of the lamprey genome and the biology of programmed genome rearrangement.     

Characterizing the composition and evolution of homoeologous genomes in hexaploid wheat through BAC-end sequencing on chromosome 3B

Bread wheat (Triticum aestivum) is one of the most important crops worldwide. However, because of its large, hexaploid, highly repetitive genome it is a challenge to develop efficient means for molecular analysis and genetic improvement in wheat. To better understand the composition and molecular evolution of the hexaploid wheat homoeologous genomes and to evaluate the potential of BAC-end sequences (BES) for marker development, we have followed a chromosome-specific strategy and generated 11 Mb of random BES from chromosome 3B, the largest chromosome of bread wheat. The sequence consisted of about 86% of repetitive elements, 1.2% of coding regions, and 13% remained unknown. With 1.2% of the sequence length corresponding to coding sequences, 6000 genes were estimated for chromosome 3B. New repetitive sequences were identified, including a Triticineae-specific tandem repeat (Fat) that represents 0.6% of the B-genome and has been differentially amplified in the homoeologous genomes before polyploidization. About 10% of the BES contained junctions between nested transposable elements that were used to develop chromosome-specific markers for physical and genetic mapping. Finally, sequence comparison with 2.9 Mb of random sequences from the D-genome of Aegilops tauschii suggested that the larger size of the B-genome is due to a higher content in repetitive elements. It also indicated which families of transposable elements are mostly responsible for differential expansion of the homoeologous wheat genomes during evolution. Our data demonstrate that BAC-end sequencing from flow-sorted chromosomes is a powerful tool for analysing the structure and evolution of polyploid and highly repetitive genomes.     

Structural genes adjacent to interspersed repetitive DNA sequences

The observation that repetitive and single copy sequences are interspersed in animal DNAs has suggested that repetitive sequences are adjacent to single copy structural gene sequences. To test this concept, single copy DNA sequences contiguous to interspersed repetitive sequences were prepared from sea urchin DNA by hydroxyapatite fractionation (repeat-contiguous DNA fraction). These single copy sequences included about one third of the total nonrepetitive sequence in the genome as determined by the amounts recovered during the hydroxyapatite fractionation and by reassociation kinetics. ³H-labeled mRNA from sea urchin gastrula was prepared by puromycin release from polysomes and used in DNA-driven hybridization reactions. The kinetics of mRNA hybridization reactions with excess whole DNA were carefully measured, and the rate of hybridization was found to be 3–5 times slower than the corresponding single copy DNA driver reassociation rate. The mRNA hybridized with excess repeat-contiguous DNA with similar kinetics relative to the driver DNA. At completion 80% of that mRNA hybridizable with whole DNA (approximately 65%) had reacted with the repeat-contiguous DNA fraction (50%). This result shows that 80–100% of the mRNA molecules present in sea urchin embryos are transcribed from single copy DNA sequences adjacent to interspersed repetitive sequences in the genome.     

Sequence organization in Xenopus DNA studied by the electron microscope.

Xenopus laevis DNA was extracted from red blood cells and sheared to a mean length of 2780 nucleotides. The DNA was stripped of foldback-containing fragments and incubated to C₀t 10 (mol · s · l⁻¹), allowing most repetitive sequences to form duplex structures. Duplex-containing fragments were eluted from an hydroxylapatite column and visualized for electron microscopy by spreading from 57% formamide according to the modified Kleinschmidt technique of Davis et al. (1971). The mean length of the fragments observed was 2445 nucleotides. A total of 1700 DNA strands were photographed and studied. Less than 5% of the total strand length was in uninterpretable structures. Every molecule falling within the confines of the plates was included in the sample. Over 50% of the total strand length in the sample was found in structures bearing at least one interspersed repetitive sequence duplex terminated by four single-strand regions. The fraction of DNA present in duplex regions was almost exactly that predicted if the duplex regions represent all the interspersed middle repetitive sequence in the Xenopus genome. Direct measurement of visualized duplexes shows that the mean length of interspersed repetitive sequence elements in this genome is 345 nucleotides. Duplex length was shown to be independent of the length of the strands bearing the duplexes. These observations provide direct confirmation of the length of approximately 300 nucleotides indicated for interspersed repetitive sequences by earlier physical-chemical studies 011 Xenopus DNA. In strands carrying two duplexes terminated by single-strand regions the interduplex, or single-copy sequence element length could be measured. Sequence interspersion curves generated from these data are roughly consistent with those derived earlier from measurements of hydroxylapatite binding as a function of fragment length.     

The nuclear genome of Brachypodium distachyon: analysis of BAC end sequences

Due in part to its small genome (~350 Mb), Brachypodium distachyon is emerging as a model system for temperate grasses, including important crops like wheat and barley. We present the analysis of 10.9% of the Brachypodium genome based on 64,696 bacterial artificial chromosome (BAC) end sequences (BES). Analysis of repeat DNA content in BES revealed that approximately 11.0% of the genome consists of known repetitive DNA. The vast majority of the Brachypodium repetitive elements are LTR retrotransposons. While Bare-1 retrotransposons are common to wheat and barley, Brachypodium repetitive element sequence-1 (BRES-1), closely related to Bare-1, is also abundant in Brachypodium. Moreover, unique Brachypodium repetitive element sequences identified constitute approximately 7.4% of its genome. Simple sequence repeats from BES were analyzed, and flanking primer sequences for SSR detection potentially useful for genetic mapping are available at . Sequence analyses of BES indicated that approximately 21.2% of the Brachypodium genome represents coding sequence. Furthermore, Brachypodium BES have more significant matches to ESTs from wheat than rice or maize, although these species have similar sizes of EST collections. A phylogenetic analysis based on 335 sequences shared among seven grass species further revealed a closer relationship between Brachypodium and Triticeae than Brachypodium and rice or maize. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. N. Huo and G.R. Lazo contributed equally to this work.     

In Depth Characterization of Repetitive DNA in 23 Plant Genomes Reveals Sources of Genome Size Variation in the Legume Tribe Fabeae

The differential accumulation and elimination of repetitive DNA are key drivers of genome size variation in flowering plants, yet there have been few studies which have analysed how different types of repeats in related species contribute to genome size evolution within a phylogenetic context. This question is addressed here by conducting large-scale comparative analysis of repeats in 23 species from four genera of the monophyletic legume tribe Fabeae, representing a 7.6-fold variation in genome size. Phylogenetic analysis and genome size reconstruction revealed that this diversity arose from genome size expansions and contractions in different lineages during the evolution of Fabeae. Employing a combination of low-pass genome sequencing with novel bioinformatic approaches resulted in identification and quantification of repeats making up 55–83% of the investigated genomes. In turn, this enabled an analysis of how each major repeat type contributed to the genome size variation encountered. Differential accumulation of repetitive DNA was found to account for 85% of the genome size differences between the species, and most (57%) of this variation was found to be driven by a single lineage of Ty3/gypsy LTR-retrotransposons, the Ogre elements. Although the amounts of several other lineages of LTR-retrotransposons and the total amount of satellite DNA were also positively correlated with genome size, their contributions to genome size variation were much smaller (up to 6%). Repeat analysis within a phylogenetic framework also revealed profound differences in the extent of sequence conservation between different repeat types across Fabeae. In addition to these findings, the study has provided a proof of concept for the approach combining recent developments in sequencing and bioinformatics to perform comparative analyses of repetitive DNAs in a large number of non-model species without the need to assemble their genomes.