Characterization of a deep-coverage carrot (<Emphasis Type="Italic">Daucus carota</Emphasis> L.) BAC library and initial analysis of BAC-end sequences

Carrot is the most economically important member of the Apiaceae family and a major source of provitamin A carotenoids in the human diet. However, carrot molecular resources are relatively underdeveloped, hampering a number of genetic studies. Here, we report on the synthesis and characterization of a bacterial artificial chromosome (BAC) library of carrot. The library is 17.3-fold redundant and consists of 92,160 clones with an average insert size of 121 kb. To provide an overview of the composition and organization of the carrot nuclear genome we generated and analyzed 2,696 BAC-end sequences (BES) from nearly 2,000 BACs, totaling 1.74 Mb of BES. This analysis revealed that 14% of the BES consists of known repetitive elements, with transposable elements representing more than 80% of this fraction. Eleven novel carrot repetitive elements were identified, covering 8.5% of the BES. Analysis of microsatellites showed a comparably low frequency for these elements in the carrot BES. Comparisons of the translated BES with protein databases indicated that approximately 10% of the carrot genome represents coding sequences. Moreover, among eight dicot species used for comparison purposes, carrot BES had highest homology to protein-coding sequences from tomato. This deep-coverage library will aid carrot breeding and genetics. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. Nucleotide sequence data reported are available in the DDBJ/EMBL/GenBank databases under the accession numbers FJ147695–FJ150390.     

Intraspecific sequence comparisons reveal similar rates of non-collinear gene insertion in the B and D genomes of bread wheat

ABSTRACT: BACKGROUND: Polyploidization is considered one of the main mechanisms of plant genome evolution. The presence of multiple copies of the same gene reduces selection pressure and permits sub-functionalization and neo-functionalization leading to plant diversification, adaptation and speciation. In bread wheat, polyploidization and the prevalence of transposable elements resulted in massive gene duplication and movement. As a result, the number of genes which are non-collinear to genomes of related species seems markedly increased in wheat. RESULTS: We used new-generation sequencing (NGS) to generate sequence of a Mb-sized region from wheat chromosome arm 3DS. Sequence assembly of 24 BAC clones resulted in two scaffolds of 1,264,820 and 333,768 bases. The sequence was annotated and compared to the homoeologous region on wheat chromosome 3B and orthologous loci of Brachypodium distachyon and rice. Among 39 coding sequences in the 3DS scaffolds, 32 have a homoeolog on chromosome 3B. In contrast, only fifteen and fourteen orthologs were identified in the corresponding regions in rice and Brachypodium, respectively. Interestingly, five pseudogenes were identified among the non-collinear coding sequences at the 3B locus, while none was found at the 3DS locus. CONCLUSION: Direct comparison of two Mb-sized regions of the B and D genomes of bread wheat revealed similar rates of non-collinear gene insertion in both genomes with a majority of gene duplications occurring before their divergence. Relatively low proportion of pseudogenes was identified among non-collinear coding sequences. Our data suggest that the pseudogenes did not originate from insertion of non-functional copies, but were formed later during the evolution of hexaploid wheat. Some evidence was found for gene erosion along the B genome locus.     

Genome-wide BAC-end sequencing of <Emphasis Type="Italic">Cucumis melo</Emphasis> using two BAC libraries

Background

Although melon (Cucumis melo L.) is an economically important fruit crop, no genome-wide sequence information is openly available at the current time. We therefore sequenced BAC-ends representing a total of 33,024 clones, half of them from a previously described melon BAC library generated with restriction endonucleases and the remainder from a new random-shear BAC library.

Results

We generated a total of 47,140 high-quality BAC-end sequences (BES), 91.7% of which were paired-BES. Both libraries were assembled independently and then cross-assembled to obtain a final set of 33,372 non-redundant, high-quality sequences. These were grouped into 6,411 contigs (4.5 Mb) and 26,961 non-assembled BES (14.4 Mb), representing ~4.2% of the melon genome. The sequences were used to screen genomic databases, identifying 7,198 simple sequence repeats (corresponding to one microsatellite every 2.6 kb) and 2,484 additional repeats of which 95.9% represented transposable elements. The sequences were also used to screen expressed sequence tag (EST) databases, revealing 11,372 BES that were homologous to ESTs. This suggests that ~30% of the melon genome consists of coding DNA. We observed regions of microsynteny between melon paired-BES and six other dicotyledonous plant genomes.

Conclusion

The analysis of nearly 50,000 BES from two complementary genomic libraries covered ~4.2% of the melon genome, providing insight into properties such as microsatellite and transposable element distribution, and the percentage of coding DNA. The observed synteny between melon paired-BES and six other plant genomes showed that useful comparative genomic data can be derived through large scale BAC-end sequencing by anchoring a small proportion of the melon genome to other sequenced genomes.

     

Next-generation sequencing and syntenic integration of flow-sorted arms of wheat chromosome 4A exposes the chromosome structure and gene content

Wheat is the third most important crop for human nutrition in the world. The availability of high-resolution genetic and physical maps and ultimately a complete genome sequence holds great promise for breeding improved varieties to cope with increasing food demand under the conditions of changing global climate. However, the large size of the bread wheat (Triticum aestivum) genome (approximately 17 Gb/1C) and the triplication of genic sequence resulting from its hexaploid status have impeded genome sequencing of this important crop species. Here we describe the use of mitotic chromosome flow sorting to separately purify and then shotgun-sequence a pair of telocentric chromosomes that together form chromosome 4A (856 Mb/1C) of wheat. The isolation of this much reduced template and the consequent avoidance of the problem of sequence duplication, in conjunction with synteny-based comparisons with other grass genomes, have facilitated construction of an ordered gene map of chromosome 4A, embracing ≥85% of its total gene content, and have enabled precise localization of the various translocation and inversion breakpoints on chromosome 4A that differentiate it from its progenitor chromosome in the A genome diploid donor. The gene map of chromosome 4A, together with the emerging sequences of homoeologous wheat chromosome groups 4, 5 and 7, represent unique resources that will allow us to obtain new insights into the evolutionary dynamics between homoeologous chromosomes and syntenic chromosomal regions.     

Lr34 multi-pathogen resistance ABC transporter: molecular analysis of homoeologous and orthologous genes in hexaploid wheat and other grass species

The Triticum aestivum (bread wheat) disease resistance gene Lr34 confers durable, race non-specific protection against three fungal pathogens, and has been a highly relevant gene for wheat breeding since the green revolution. Lr34, located on chromosome 7D, encodes an ATP-binding cassette (ABC) transporter. Both wheat cultivars with and without Lr34-based resistance encode a putatively functional protein that differ by only two amino acid polymorphisms. In this study, we focused on the identification and characterization of homoeologous and orthologous Lr34 genes in hexaploid wheat and other grasses. In hexaploid wheat we found an expressed and putatively functional Lr34 homoeolog located on chromosome 4A, designated Lr34-B. Another homoeologous Lr34 copy, located on chromosome 7A, was disrupted by the insertion of repetitive elements. Protein sequences of LR34-B and LR34 were 97% identical. Orthologous Lr34 genes were detected in the genomes of Oryza sativa (rice) and Sorghum bicolor (sorghum). Zea mays (maize), Brachypodium distachyon and Hordeum vulgare (barley) lacked Lr34 orthologs, indicating independent deletion of this particular ABC transporter. Lr34 was part of a gene-rich island on the wheat D genome. We found gene colinearity on the homoeologous A and B genomes of hexaploid wheat, but little microcolinearity in other grasses. The homoeologous LR34-B protein and the orthologs from rice and sorghum have the susceptible haplotype for the two critical polymorphisms distinguishing the LR34 proteins from susceptible and resistant wheat cultivars. We conclude that the particular Lr34-haplotype found in resistant wheat cultivars is unique. It probably resulted from functional gene diversification that occurred after the polyploidization event that was at the origin of cultivated bread wheat.     

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.     

Comparative analysis of Stowaway-like miniature inverted repeat transposable elements in wheat group 7 chromosomes: Abundance,composition, and evolution

Miniature inverted repeat transposable elements (MITEs) are the most ubiquitous transposable elements in eukaryotic genomes; they play a prominent role in sequence divergence and genome evolution. There are many well-characterized Stowaway-like MITE families in wheat, but their distribution, abundance, and composition at the chromosome level are still not well understood. In this study, we systematically investigated the Stowaway-like MITEs in wheat group 7 chromosomes based on the survey sequences of isolated wheat chromosomes, to compare them at the chromosome level and to reveal their evolutionary role on wheat polyploidization. In summary, 2026 MITEs were identified, of which 587, 714, and 725 were distributed on 7A, 7B, and 7D chromosomes, respectively. There are more MITEs present on 7D, compared to 7A and 7B, suggesting A and B subgenomes eliminated some repetitive elements during two hybridization processes. Furthermore, some chromosome/arm-specific MITEs were also identified, providing information on the function and evolution of MITEs in wheat genomes. The sequence diversity of the MITE insertions was also investigated. This study for the first time investigated the abundance and composition of MITEs at the chromosome level, which will be beneficial to improve our understanding of the distribution of wheat MITEs and their evolutionary role in polyploidization.     

Updating of transposable element annotations from large wheat genomic sequences reveals diverse activities and gene associations

Triticeae species (including wheat, barley and rye) have huge and complex genomes due to polyploidization and a high content of transposable elements (TEs). TEs are known to play a major role in the structure and evolutionary dynamics of Triticeae genomes. During the last 5 years, substantial stretches of contiguous genomic sequence from various species of Triticeae have been generated, making it necessary to update and standardize TE annotations and nomenclature. In this study we propose standard procedures for these tasks, based on structure, nucleic acid and protein sequence homologies. We report statistical analyses of TE composition and distribution in large blocks of genomic sequences from wheat and barley. Altogether, 3.8 Mb of wheat sequence available in the databases was analyzed or re-analyzed, and compared with 1.3 Mb of re-annotated genomic sequences from barley. The wheat sequences were relatively gene-rich (one gene per 23.9 kb), although wheat gene-derived sequences represented only 7.8% (159 elements) of the total, while the remainder mainly comprised coding sequences found in TEs (54.7%, 751 elements). Class I elements [mainly long terminal repeat (LTR) retrotransposons] accounted for the major proportion of TEs, in terms of sequence length as well as element number (83.6% and 498, respectively). In addition, we show that the gene-rich sequences of wheat genome A seem to have a higher TE content than those of genomes B and D, or of barley gene-rich sequences. Moreover, among the various TE groups, MITEs were most often associated with genes: 43.1% of MITEs fell into this category. Finally, the TRIM and copia elements were shown to be the most active TEs in the wheat genome. The implications of these results for the evolution of diploid and polyploid wheat species are discussed. Electronic Supplementary Material Supplementary material is available for this article at     

A genome-wide survey of switchgrass genome structure and organization

The perennial grass, switchgrass (Panicum virgatum L.), is a promising bioenergy crop and the target of whole genome sequencing. We constructed two bacterial artificial chromosome (BAC) libraries from the AP13 clone of switchgrass to gain insight into the genome structure and organization, initiate functional and comparative genomic studies, and assist with genome assembly. Together representing 16 haploid genome equivalents of switchgrass, each library comprises 101,376 clones with average insert sizes of 144 (HindIII-generated) and 110 kb (BstYI-generated). A total of 330,297 high quality BAC-end sequences (BES) were generated, accounting for 263.2 Mbp (16.4%) of the switchgrass genome. Analysis of the BES identified 279,099 known repetitive elements, >50,000 SSRs, and 2,528 novel repeat elements, named switchgrass repetitive elements (SREs). Comparative mapping of 47 full-length BAC sequences and 330K BES revealed high levels of synteny with the grass genomes sorghum, rice, maize, and Brachypodium. Our data indicate that the sorghum genome has retained larger microsyntenous regions with switchgrass besides high gene order conservation with rice. The resources generated in this effort will be useful for a broad range of applications.     

Sequence analysis of the long arm of rice chromosome 11 for rice–wheat synteny

The DNA sequence of 106 BAC/PAC clones in the minimum tiling path (MTP) of the long arm of rice chromosome 11, between map positions 57.3 and 116.2 cM, has been assembled to phase 2 or PLN level. This region has been sequenced to 10× redundancy by the Indian Initiative for Rice Genome Sequencing (IIRGS) and is now publicly available in GenBank. The region, excluding overlaps, has been predicted to contain 2,932 genes using different software. A gene-by-gene BLASTN search of the NCBI wheat EST database of over 420,000 cDNA sequences revealed that 1,143 of the predicted rice genes (38.9%) have significant homology to wheat ESTs (bit score 100). Further BLASTN search of these 1,143 rice genes with the GrainGenes database of sequence contigs containing bin-mapped wheat ESTs allowed 113 of the genes to be placed in bins located on wheat chromosomes of different homoeologous groups. The largest number of genes, about one-third, mapped to the homoeologous group 4 chromosomes of wheat, suggesting a common evolutionary origin. The remaining genes were located on wheat chromosomes of different groups with significantly higher numbers for groups 3 and 5. Location of bin-mapped wheat contigs to chromosomes of all the seven homoeologous groups can be ascribed to movement of genes (transpositions) or chromosome segments (translocations) within rice or the hexaploid wheat genomes. Alternatively, it could be due to ancient duplications in the common ancestral genome of wheat and rice followed by selective elimination of genes in the wheat and rice genomes. While there exists definite conservation of gene sequences and the ancestral chromosomal identity between rice and wheat, there is no obvious conservation of the gene order at this level of resolution. Lack of extensive colinearity between rice and wheat genomes suggests that there have been many insertions, deletions, duplications and translocations that make the synteny comparisons much more complicated than earlier thought. However, enhanced resolution of comparative sequence analysis may reveal smaller conserved regions of colinearity, which will facilitate selection of markers for saturation mapping and sequencing of the gene-rich regions of the wheat genome.     

Chromosomal location of a K/Na discrimination character in the D genome of wheat

Summary K/Na ratios have been determined in the leaves of salt-treated plants of 14 disomic substitution lines in which each of the D-genome chromosomes replaces the homoeologous A- or B-genome chromosome in the tetraploid wheat variety Langdon (AABB genome). Aneuploid lines of hexaploid bread wheat (cv Chinese Spring) having a reduced or an enhanced complement of chromosome 4D have also been examined. These investigations show that the gene(s) determining K/Na ratios in the leaves of wheat plants grown in the presence of salt is located on the long arm of chromosome 4D.     

BAC-FISH in wheat identifies chromosome landmarks consisting of different types of transposable elements

Fluorescence in situ hybridization (FISH) has been widely used in the physical mapping of genes and chromosome landmarks in plants and animals. Bacterial artificial chromosomes (BACs) contain large inserts making them amenable for FISH mapping. We used BAC-FISH to study genome organization and evolution in hexaploid wheat and its relatives. We selected 56 restriction fragment length polymorphism (RFLP) locus-specific BAC clones from libraries of Aegilops tauschii (the D-genome donor of hexaploid wheat) and A-genome diploid Triticum monococcum. Different types of repetitive sequences were identified using BAC-FISH. Two BAC clones gave FISH patterns similar to the repetitive DNA family pSc119; one BAC clone gave a FISH pattern similar to the repetitive DNA family pAs1. In addition, we identified several novel classes of repetitive sequences: one BAC clone hybridized to the centromeric regions of wheat and other cereal species, except rice; one BAC clone hybridized to all subtelomeric chromosome regions in wheat, rye, barley and oat; one BAC clone contained a localized tandem repeat and hybridized to five D-genome chromosome pairs in wheat; and four BAC clones hybridized only to a proximal region in the long arm of chromosome 4A of hexaploid wheat. These repeats are valuable markers for defined chromosome regions and can also be used for chromosome identification. Sequencing results revealed that all these repeats are transposable elements (TEs), indicating the important role of TEs, especially retrotransposons, in genome evolution of wheat.Communicated by P.B. Moens     

Molecular characterization and chromosome-specific TRAP-marker development for Langdon durum D-genome disomic substitution lines.

The aneuploid stocks of durum wheat (Triticum turgidum L. subsp. durum (Desf.) Husnot) and common wheat (T. aestivum L.) have been developed mainly in 'Langdon' (LDN) and 'Chinese Spring' (CS) cultivars, respectively. The LDN-CS D-genome chromosome disomic substitution (LDN-DS) lines, where a pair of CS D-genome chromosomes substitute for a corresponding homoeologous A- or B-genome chromosome pair of LDN, have been widely used to determine the chromosomal locations of genes in tetraploid wheat. The LDN-DS lines were originally developed by crossing CS nulli-tetrasomics with LDN, followed by 6 backcrosses with LDN. They have subsequently been improved with 5 additional backcrosses with LDN. The objectives of this study were to characterize a set of the 14 most recent LDN-DS lines and to develop chromosome-specific markers, using the newly developed TRAP (target region amplification polymorphism)-marker technique. A total of 307 polymorphic DNA fragments were amplified from LDN and CS, and 302 of them were assigned to individual chromosomes. Most of the markers (95.5%) were present on a single chromosome as chromosome-specific markers, but 4.5% of the markers mapped to 2 or more chromosomes. The number of markers per chromosome varied, from a low of 10 (chromosomes 1A and 6D) to a high of 24 (chromosome 3A). There was an average of 16.6, 16.6, and 15.9 markers per chromosome assigned to the A-, B-, and D-genome chromosomes, respectively, suggesting that TRAP markers were detected at a nearly equal frequency on the 3 genomes. A comparison of the source of the expressed sequence tags (ESTs), used to derive the fixed primers, with the chromosomal location of markers revealed that 15.5% of the TRAP markers were located on the same chromosomes as the ESTs used to generate the fixed primers. A fixed primer designed from an EST mapped on a chromosome or a homoeologous group amplified at least 1 fragment specific to that chromosome or group, suggesting that the fixed primers might generate markers from target regions. TRAP-marker analysis verified the retention of at least 13 pairs of A- or B-genome chromosomes from LDN and 1 pair of D-genome chromosomes from CS in each of the LDN-DS lines. The chromosome-specific markers developed in this study provide an identity for each of the chromosomes, and they will facilitate molecular and genetic characterization of the individual chromosomes, including genetic mapping and gene identification.     

Genetic mapping of the wheat dimeric α-amylase inhibitor multi-gene family using allele-specific primers based on intergenomic SNPs

Single nucleotide polymorphisms (SNPs) identified in EST sequences can be used to map expressed genes. Though SNPs are useful markers for genetic mapping, SNP mapping of genes in common wheat is challenging because the genetic complement of wheat consists of three closely related genomes (designated A, B, and D), and most genes are present in triplicate sets. Mapping multi-gene family members is further complicated by the fact that it is difficult to distinguish SNP differences between the various paralogs from those between the different genomes. We have developed a PCR-based method for assigning wheat EST sequences to their proper genetic loci by first identifying and mapping SNPs that distinguish the three genomes. To develop this method, we focused on EST sequences encoding the dimeric α-amylase inhibitors (WDAI), The WDAI coding regions of hexaploid wheat were aligned. The sequences were classified into three groups based on nucleotide variations. Twenty-two SNPs were identified that distinguish the three groups. Group-specific primers based on these SNPs were designed to permit selective amplification of each group. The chromosomal location of each group was then determined using Group 3 ditelosomic lines of Chinese Spring. Groups 1 and 2 were assigned to chromosome locations 3DS and 3BS, respectively, whereas no sequence could be assigned to 3AS. A remarkable feature of this method is the ability to discriminate the location of homoeologous multigenes in the three genomes of wheat. This strategy can be useful for assigning unknown wheat EST sequences to specific chromosomes.     

A whole-genome snapshot of 454 sequences exposes the composition of the barley genome and provides evidence for parallel evolution of genome size in wheat and barley

The genomes of barley and wheat, two of the world's most important crops, are very large and complex due to their high content of repetitive DNA. In order to obtain a whole-genome sequence sample, we performed two runs of 454 (GS20) sequencing on genomic DNA of barley cv. Morex, which yielded approximately 1% of a haploid genome equivalent. Almost 60% of the sequences comprised known transposable element (TE) families, and another 9% represented novel repetitive sequences. We also discovered high amounts of low-complexity DNA and non-genic low-copy DNA. We identified almost 2300 protein coding gene sequences and more than 660 putative conserved non-coding sequences. Comparison of the 454 reads with previously published genomic sequences suggested that TE families are distributed unequally along chromosomes. This was confirmed by in situ hybridizations of selected TEs. A comparison of these data for the barley genome with a large sample of publicly available wheat sequences showed that several TE families that are highly abundant in wheat are absent from the barley genome. This finding implies that the TE composition of their genomes differs dramatically, despite their very similar genome size and their close phylogenetic relationship.     

The characterization and comparative analysis of high-molecular-weight glutenin genes from genomes A and B of a hexaploid bread wheat

Summary Two high-molecular-weight subunit (HMWS) glutenin genes from the A and B genomes of the hexaploid bread wheat Triticum aestivum L. cv Cheyenne have been isolated and sequenced. Both of these genes are of the high M_r class (x-type) of HMW glutenins, and have not been previously reported. The entire set of six HMW genes from cultivar Cheyenne have now been isolated and characterized. An analysis of the Ax and Bx sequences shows that the Ax sequence is similar to the homoeologous gene from the D genome, while the Bx repeat structure is significantly different. The repetitive region of these proteins can be modelled as a series of interspersed copies of repeat modifs of 6, 9, and 15 amino acid residues. The evolution of these genes includes single-base substitutions over the entire coding region, plus insertion/deletions of single or blocks of repeats in the central repetitive domain.     

Origin of a Giant Sex Chromosome

Chromosome size and morphology vary within and among species, but little is known about the proximate or ultimate causes of these differences. Cichlid fish species in the tribe Oreochromini share an unusual giant chromosome that is ∼3 times longer than the other chromosomes. This giant chromosome functions as a sex chromosome in some of these species. We test two hypotheses of how this giant sex chromosome may have evolved. The first hypothesis proposes that it evolved by accumulating repetitive elements as recombination was reduced around a dominant sex determination locus, as suggested by canonical models of sex chromosome evolution. An alternative hypothesis is that the giant sex chromosome originated via the fusion of an autosome with a highly repetitive B chromosome, one of which carried a sex determination locus. We test these hypotheses using comparative analysis of chromosome-scale cichlid and teleost genomes. We find that the giant sex chromosome consists of three distinct regions based on patterns of recombination, gene and transposable element content, and synteny to the ancestral autosome. The WZ sex determination locus encompasses the last ∼105 Mb of the 134-Mb giant chromosome. The last 47 Mb of the giant chromosome shares no obvious homology to any ancestral chromosome. Comparisons across 69 teleost genomes reveal that the giant sex chromosome contains unparalleled amounts of endogenous retroviral elements, immunoglobulin genes, and long noncoding RNAs. The results favor the B chromosome fusion hypothesis for the origin of the giant chromosome.