A bacterial artificial chromosome (BAC) library for sunflower, and identification of clones containing genes for putative transmembrane receptors

Sunflower (Helianthus annuus L.) is an economically important oil seed crop with an estimated genome size of 3000 Mb. We have constructed a bacterial artificial chromosome (BAC) library for sunflower, which represents an estimated 4- to 5-fold coverage of the genome. Nuclei isolated from young leaves were used as a source of high-molecular-weight DNA and a partial restriction endonuclease digestion protocol was used to cleave the DNA. A random sample of 60 clones indicated an average insert size of 80 kb, implying a 95% probability of recovering any specific sequence of interest. The library was screened with chloroplast DNA probes. Only 0.1% of the clones were identified to be of chloroplast origin, indicating that contamination with organellar DNAs is very low. The utility of the library was evaluated by screening for the presence of genes for putative transmembrane receptors sharing epidermal growth factor (EGF) and integrin-like domains. First, a homologous sunflower EST (HaELP1) was obtained by degenerate RT-PCR cloning, using Arabidopsis thaliana genes (AtELP) as a source of consensus sequences. Three different BACs yielded positive hybridization signals when HaELP1 was used as a probe. BAC subcloning and sequencing demonstrated the presence of two different loci putatively homologous to genes for transmembrane proteins with EGF- and integrin-like domains from sunflower. This work demonstrates the suitability of the library for homology map-based cloning of sunflower genes and physical mapping of the sunflower genome.     

Construction of a BAC library of pearl millet, Pennisetum glaucum

A bacterial artificial chromosome (BAC) library was constructed using nuclear DNA from pearl millet (Pennisetum glaucum), and used as a resource for the isolation of microsatellite sequences. The library contains a total of 159,100 clones with an average insert size of 90 kb, and corresponds to 5.8 haploid genome equivalents. The BAC library was pooled for screening by the polymerase chain reaction (PCR) as well as robotically gridded on high-density filters. PCR-based screening of a subset of the library (4.7 haploid genome equivalents) using five sequence-tagged site (STS) and six microsatellite markers identified between 2 and 11 positives superpools (5.4 on average). The frequency of BAC clones carrying inserts of chloroplast DNA was estimated to be less than 1% by hybridisation with a rice chloroplast probe. Received: 30 January 2000 / Accepted: 16 October 2000     

Fosmid Library Construction and Initial Analysis of End Sequences in Female Half-smooth Tongue Sole (Cynoglossus semilaevis)

Half-smooth tongue sole (Cynoglossus semilaevis: Pleuronectiformes) is a commercially important cultured marine flatfish in China and forms an important fishery resource, but the research of its genome is underdeveloped. In this study, we constructed a female C. semilaevis fosmid library and analyzed the fosmid end sequences to provide a preliminary assessment of the genome. The library consists of 49,920 clones with an average insert size of about 39 kb, amounting to 3.23 genome equivalents. Fosmid stability assays indicate that female C. semilaevis DNA was stable during propagation in the fosmid system. Library screening with eight microsatellite markers yielded between two and five positive clones, and none of those tested was absent from the library. End-sequencing of both 5′ and 3′ ends of 1,152 individual clones generated 2,247 sequences after trimming, with an average sequence length of 855 bp. BLASTN searches of the nr and EST databases of GenBank and BLASTX searches of the nr database resulted in 259 (11.53%) and 287 (12.77%) significant hits (E < e ⁻⁵), respectively. Repetitive sequences analysis resulted in 5.23% of base pairs masked using both the Fugu and Danio databases, repetitive elements were composed of retroelements, DNA transposons, satellites, simple repeats, and low-complexity sequences. The fosmid library, in conjunction with the fosmid end sequences, will serve as a useful resource for large-scale genome sequencing, physical mapping, and positional cloning, and provide a better understanding of female C. semilaevis genome.     

Fosmid Library Construction and Initial Analysis of End Sequences in Zhikong Scallop (<Emphasis Type="Italic">Chlamys farreri</Emphasis>)

Zhikong scallop (Chlamys farreri Jones et Preston, 1904) is one of the most commercially important bivalves in China, but research on its genome is underdeveloped. In this study, we constructed the first Zhikong scallop fosmid library, and analyzed the fosmid end sequences to provide a preliminary assessment of the genome. The library consists of 133,851 clones with an average insert size of about 40 kb, amounting to 4.3 genome equivalents. Fosmid stability assays indicate that Zhikong scallop DNA was stable during propagation in the fosmid system. Library screening with two genes and seven microsatellite markers yielded between two and eight positive clones, and none of those tested was absent from the library. End-sequencing of 480 individual clones generated 828 sequences after trimming, with an average sequence length of 624 bp. BLASTN searches of the nr and EST databases of GenBank and BLASTX searches of the nr database resulted in 213 (25.72%) and 44 (5.31%) significant hits (E < e⁻⁵), respectively. Repetitive sequences analysis resulted in 375 repeats, accounting for 15.84% of total length, which were composed of interspersed repetitive sequences, tandem repeats, and low-complexity sequences. The fosmid library, in conjunction with the fosmid end sequences, will serve as a useful resource for physical mapping and positional cloning, and provide a better understanding of the Zhikong scallop genome.     

Fosmid library construction and end sequences analysis of the Pacific oyster,Crassostrea gigas

The Pacific oyster (Crassostrea gigas) is globally distributed and is one of the most commercially and ecologically important marine organisms. However, little is known about the genome of this species. In this study, a C. gigas fosmid library was constructed that contains 459,936 clones with an average insert size of approximately 40 kb, representing 22.34-fold haploid genome equivalents. End sequencing generated 90,240 fosmid end sequences (FESs) with an average length of 384.27 base pairs (bp), covering approximately 2.58% of the Pacific oyster genome. The FESs were subsequently assembled and annotated, resulting in 6332 sequences with predicted open reading frames≥300 and 1,189,100 bp repeats. Furthermore, a total of 3200 microsatellite repeats were identified, and dinucleotide repeats were found to occur most abundantly, with AG and AAT being the most abundant repeat class of dinucleotides and trinucleotides. We also found that the repeat number was generally negatively proportional to the repeat element length. Microsatellites composition between the transcribed sequences and genomic sequences was shown to be different. Point mutations of microsatellite were non-random and underwent strong selection stress. Overall, a comprehensive sequence resource for the Pacific oyster was created, including annotated transposable elements, tandem repeats, protein coding sequences and microsatellites. These initial findings will serve as resources for further in-depth studies of physical mapping, gene discovery, microsatellite marker developing and evolution studies.     

Construction and characterization of two BAC libraries from Brachypodium distachyon, a new model for grass genomics.

Brachypodium is well suited as a model system for temperate grasses because of its compact genome and a range of biological features. In an effort to develop resources for genome research in this emerging model species, we constructed 2 bacterial artificial chromosome (BAC) libraries from an inbred diploid Brachypodium distachyon line, Bd21, using restriction enzymes HindIII and BamHI. A total of 73,728 clones (36,864 per BAC library) were picked and arrayed in 192,384-well plates. The average insert size for the BamHI and HindIII libraries is estimated to be 100 and 105 kb, respectively, and inserts of chloroplast origin account for 4.4% and 2.4%, respectively. The libraries individually represent 9.4- and 9.9-fold haploid genome equivalents with combined 19.3-fold genome coverage, based on a genome size of 355 Mb reported for the diploid Brachypodium, implying a 99.99% probability that any given specific sequence will be present in each library. Hybridization of the libraries with 8 starch biosynthesis genes was used to empirically evaluate this theoretical genome coverage; the frequency at which these genes were present in the library clones gave an estimated coverage of 11.6- and 19.6-fold genome equivalents. To obtain a first view of the sequence composition of the Brachypodium genome, 2185 BAC end sequences (BES) representing 1.3 Mb of random genomic sequence were compared with the NCBI GenBank database and the GIRI repeat database. Using a cutoff expectation value of E<10-10, only 3.3% of the BESs showed similarity to repetitive sequences in the existing database, whereas 40.0% had matches to the sequences in the EST database, suggesting that a considerable portion of the Brachypodium genome is likely transcribed. When the BESs were compared with individual EST databases, more matches hit wheat than maize, although their EST collections are of a similar size, further supporting the close relationship between Brachypodium and the Triticeae. Moreover, 122 BESs have significant matches to wheat ESTs mapped to individual chromosome bin positions. These BACs represent colinear regions containing the mapped wheat ESTs and would be useful in identifying additional markers for specific wheat chromosome regions.     

Genome-wide analysis of conservation and divergence of microsatellites in rice

Studies on microsatellite distribution and divergence in related genomes contribute towards understanding of genome evolution in eukaryotes. Despite the availability of whole genome sequences of four rice genomes, occurrence and significance of microsatellites in the rice genome has remained a relatively unexplored area of research. We have aligned genomes of two rice subspecies i.e. indica and japonica to understand the trends of microsatellite conservation and divergence in the rice genome. Nearly 62% of the indica microsatellites were also found in the japonica genome. Occurrence of microsatellites showed a negative association with that of retrotransposons. Microsatellites repeat unit length and sequence showed direct influence on the microsatellite locus length. Further, microsatellite allele length was also influenced by the sequence characteristics of the neighbouring regions. CCG repeats were most conserved microsatellite sequences across the different syntenic regions in the two rice genomes and often showed association with CpG islands. Our study suggested that microsatellite distribution is not only governed by a balance between replication slippage and point mutations as proposed earlier, but also by the microsatellite motif sequence and characteristics of microsatellite neighbouring regions in the genome. Thus, this study is likely to prove an important reference for understanding the process of microsatellite evolution and dynamics in the two rice subspecies.     

Characterisation and analysis of microsatellite loci in a mangrove species, Avicennia marina (Forsk.) Vierh. (Avicenniaceae)

An enriched microsatellite library of the mangrove species Avicennia marina was constructed, in which 85.8% of the clones contained microsatellite sequences. Of the microsatellite repeat sequences isolated, 55.0% were di-nucleotides, 34.2% were tri-nucleotides, 50.0% were perfect, 24.2% were imperfect, and 15.0% were compound. Four different di-nucleotide repeats were isolated with repeat lengths ranging from 5 to 33; ten different tri-nucleotide repeats were isolated with repeat lengths ranging from 3 to 25. The most common di-nucleotide was the AC/TG repeat; the most common tri-nucleotide was the CCG/GGC repeat. Sixteen microsatellite sequences were selected for primer design, and 6 primers were selected to investigate the polymorphism detected among 15 individuals of A. marina from three natural populations in Australia. A total of 40 alleles were detected at 6 microsatellite loci. The number of alleles per microsatellite locus ranged from 5 to 13. On average, 7 alleles were detected per locus. All microsatellite loci showed high levels of gene diversity (heterozygosity), with values ranging from 0.53 to 0.88; the mean value of gene diversity was 0.70. Microsatellite loci were also tested for conservation across Avicennia species. There was a decline in amplification success with increasing divergence between Avicennia species. The results indicate that microsatellites are abundant in the Avicennia genome and can be valuable genetic markers for assessing the effects of deforestation and forest fragmentation in mangrove communities, which is an important issue for mangrove conservation and afforestation schemes. Received: 8 June 1999 / Accepted: 21 September 1999     

Construction of BAC and BIBAC libraries from sunflower and identification of linkage group-specific clones by overgo hybridization

Complementary BAC and BIBAC libraries were constructed from nuclear DNA of sunflower cultivar HA 89. The BAC library, constructed with BamHI in the pECBAC1 vector, contains 107,136 clones and has an average insert size of 140 kb. The BIBAC library was constructed with HindIII in the plant-transformation-competent binary vector pCLD04541 and contains 84,864 clones, with an average insert size of 137 kb. The two libraries combined contain 192,000 clones and are equivalent to approximately 8.9 haploid genomes of sunflower (3,000 Mb/1C), and provide a greater than 99% probability of obtaining a clone of interest. The frequencies of BAC and BIBAC clones carrying chloroplast or mitochondrial DNA sequences were estimated to be 2.35 and 0.04%, respectively, and insert-empty clones were less than 0.5%. To facilitate chromosome engineering and anchor the sunflower genetic map to its chromosomes, one to three single- or low-copy RFLP markers from each linkage group of sunflower were used to design pairs of overlapping oligonucleotides (overgos). Thirty-six overgos were designed and pooled as probes to screen a subset (5.1×) of the BAC and BIBAC libraries. Of the 36 overgos, 33 (92%) gave at least one positive clone and 3 (8%) failed to hit any clone. As a result, 195 BAC and BIBAC clones representing 19 linkage groups were identified, including 76 BAC clones and 119 BIBAC clones, further verifying the genome coverage and utility of the libraries. These BAC and BIBAC libraries and linkage group-specific clones provide resources essential for comprehensive research of the sunflower genome.     

Development of Polymorphic EST Markers Suitable for Genetic Linkage Mapping of Catfish

Expressed sequence tag (EST) markers are important for gene mapping and for marker-assisted selection (MAS). To develop EST markers for use in catfish gene mapping, 100 randomly picked complementary DNAs from the channel catfish (Ictalurus punctatus) pituitary library were sequenced. The EST sequences were used to design primers to amplify channel catfish and blue catfish (I. furcatus) genomic DNAs. Polymerase chain reaction products of the ESTs were analyzed to determine length polymorphism between the channel catfish and blue catfish. Eleven polymorphic EST markers were identified. Five of the 11 EST markers were from known genes and the other six were from unidentified ESTs. Seven ESTs were found to be associated with microsatellite sequences. Analysis of channel catfish gene sequences indicated highly biased codon usage, with 16 codons being preferably used. These codons were more preferably used in highly expressed ribosomal protein genes and in highly expressed pituitary hormone genes. G/C-rich codons are less used in channel catfish than those in other vertebrates suggesting AT-richness of the channel catfish genome. Received June 29, 1998; accepted March 29, 1999.     

The first insight into the <Emphasis Type="Italic">Taxus</Emphasis> genome via fosmid library construction and end sequencing

Taxus mairei is a critically endangered and commercially important cultured medicinal gymnosperm in China and forms an important medicinal resource, but the research of its genome is absent. In this study, we constructed a T. mairei fosmid library and analyzed the fosmid end sequences to provide a preliminary assessment of the genome. The library consists of one million clones with an average insert size of about 39 kb, amounting to 3.9 genome equivalents. Fosmid stability assays indicate that T. mairei DNA was stable during propagation in the fosmid system. End sequencing of both 5′ and 3′ ends of 968 individual clones generated 1,923 sequences after trimming, with an average sequence length of 839 bp. BLASTN searches of the nr and EST databases of GenBank and BLASTX searches of the nr database resulted in 560 (29.1%) significant hits (E < e⁻⁵). Repetitive sequences analysis revealed that 20.8% of end sequences are repetitive elements, which were composed of retroelements, DNA transposons, satellites, simple repeats, and low complexity sequences. The distribution pattern of various repeat types was found to be more similar to the gymnosperm Pinus and Picea than to the monocot and dicot. The satellites of T. mairei were significantly longer than those of P. taeda and P. glauca. The tetra-nucleotide repeats of T. mairei were much longer than those of P. glauca and P. taeda. The fosmid library and the fosmid end sequences, for the first time, will serve as a useful resource for large-scale genome sequencing, physical mapping, SSR marker development and positional cloning, and provide a better understanding of the Taxus genome.     

藏波罗花叶绿体基因组分析

藏波罗花（Incarvillea younghusbandii Sprague）是一种传统的补益类中药。其根作草药使用,用于滋补强壮,治产后少乳、久病虚弱、头晕、贫血等症。但目前关于藏波罗花分子遗传信息的研究很少。本研究基于高通量测序技术对藏波罗花叶绿体基因组进行测序、组装和注释,并对其序列特征、密码子偏好性、重复序列、系统发育和分化时间进行分析。结果表明,藏波罗花叶绿体基因组全长为159 323 bp,包含1个大单拷贝区（80 197 bp）、1个小单拷贝区（9 030 bp）和2个反向重复区（35 048 bp）;共注释出120个基因,包括77个蛋白编码基因、8个rRNA基因和35个tRNA基因;密码子偏好性分析显示,AAA是藏波罗花叶绿体基因组中使用最频繁的密码子;从藏波罗花叶绿体基因组中共检测到42个简单重复序列（simple sequence repeats,SSR）;系统发育分析表明,藏波罗花与密生波罗花（Incarvillea compacta）的亲缘关系最近,且在大概466万年前产生分化。本研究对藏波罗花相关资源的科学保护和开发具有重要的现实意义,也可以为后续角蒿属（Incarvillea）的物种鉴定、紫葳科（Bignoniaceae）的种群遗传多样性研究提供基本的遗传资源。     

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.     

The determination of multiple microsatellite genotypes and DNA sequences from a single pollen grain

Pollen plays important roles in the reproduction and gene flow of flowering plants, and its haploid DNA sequence provides useful information for studies of plant evolution and genealogy. We describe a new method for multiple microsatellite genotyping and DNA sequencing from a single pollen grain. The haploid DNA was extracted from a single pollen grain by using a simple DNA extraction method, and multiple microsatellite genotypes and DNA sequences of multiple chloroplast loci were determined. Using nine pairs of microsatellite primers, more than 90% of genotypes were successfully determined, and 71% and 100% of DNA sequences were determined at two chloroplast DNA loci, the trnL intron region and the trnL/trnF intergenic spacer region, respectively. This simple method of genetic analysis for a single pollen grain will facilitate detailed study of pollination, evolution and genealogy.     

Genealogical use of chloroplast DNA variation for intraspecific studies of Aegilops tauschii Coss.

Intraspecific patterns of chloroplast DNA variation was studied in Aegilops tauschii Coss., the D-genome progenitor of bread wheat. Nucleotide sequences of ten chloroplast microsatellite loci were analyzed for 63 accessions that cover the central part of the species distribution. As is often the case with nuclear microsatellites, those of chloroplasts of Ae. tauschii bear complex mutations. Several types of mutations other than change in the microsatellite repeat number were found, including base substitutions and length mutations in flanking regions. In total, eight mutations were present in the flanking regions of four loci. Most mutations in the flanking regions of microsatellite repeats are associated with biallelic polymorphisms. Phylogeographic analyses showed that such biallelic polymorphisms are useful to investigate intraspecific patterns of monophyletic lineage divergence. In contrast, most microsatellite repeat sites are multiallelic, variable within intraspecific lineages, and useful to compare degrees of genetic diversity between lineages. These findings show that the chloroplast genome harbors evolutionary variations informative for intraspecific studies of Ae. tauschii and can be analyzed by genealogical approaches.Electronic Supplementary Material Supplementary material is available for this article at     

石刁柏雄性偏向核质体DNA的克隆与分析

该研究以雌雄异株植物石刁柏为材料,利用基因组消减杂交技术对石刁柏雌雄核基因组中的性别差异核质体DNA（nuclear plastid DNA,NUPTs）进行了分离和分析。结果表明：（1）通过构建消减杂交文库共获得了52个雄性偏向序列,序列长度分布在63~297 bp之间,其中有19个差异序列属于叶绿体来源序列（命名为Ao1~Ao19）,且这些序列与石刁柏叶绿体基因组的相似性均大于84%,Ao19与石刁柏叶绿体基因组相似性为100%。（2）利用基因组半定量PCR对19个NUPTs序列的性别差异分析表明,有4条序列为稳定的雄性偏向NUPTs序列,分别为Ao1、Ao3、Ao10和Ao18。（3）序列比对表明,转移到核基因组的NUPTs主要来源于叶绿体基因组的反向重复区（包含IRa和IRb区）,说明石刁柏叶绿体基因组重复区序列更容易向核基因组进行转移形成雄性偏向的NUPTs序列。     

Sequence Analysis of the Genome of an Oil-Bearing Tree, Jatropha curcas L.

The whole genome of Jatropha curcas was sequenced, using a combination of the conventional Sanger method and new-generation multiplex sequencing methods. Total length of the non-redundant sequences thus obtained was 285 858 490 bp consisting of 120 586 contigs and 29 831 singlets. They accounted for ∼95% of the gene-containing regions with the average G + C content was 34.3%. A total of 40 929 complete and partial structures of protein encoding genes have been deduced. Comparison with genes of other plant species indicated that 1529 (4%) of the putative protein-encoding genes are specific to the Euphorbiaceae family. A high degree of microsynteny was observed with the genome of castor bean and, to a lesser extent, with those of soybean and Arabidopsis thaliana. In parallel with genome sequencing, cDNAs derived from leaf and callus tissues were subjected to pyrosequencing, and a total of 21 225 unigene data have been generated. Polymorphism analysis using microsatellite markers developed from the genomic sequence data obtained was performed with 12 J. curcas lines collected from various parts of the world to estimate their genetic diversity. The genomic sequence and accompanying information presented here are expected to serve as valuable resources for the acceleration of fundamental and applied research with J. curcas, especially in the fields of environment-related research such as biofuel production. Further information on the genomic sequences and DNA markers is available at http://www.kazusa.or.jp/jatropha/.     

Direct sequencing of tobacco chloroplast genome by the polymerase chain reaction

We have developed a polymerase chain reaction (PCR) method for sequencing of tobacco chloroplast genome. In a mixture containing chloroplast DNA, 5-end-labeled oligonucleotide primer, Taq DNA polymerase and reaction buffer, we were able to sequence a segment of chloroplast 16S rRNA gene. The results showed that the 750 bp of DNA sequenced were identical to the sequence reported, indicating that direct sequencing method that we have developed is useful for the sequencing of chloroplast genome. To analyze the chloroplast genome more rapidly in those in vitro grown plantlets, we also developed a simple method which is applicable for the amplifications and sequencing of chloroplast 16S rRNA fragment from either 0.15 g of tobacco leaf or stem tissue. The readable sequences obtained from the presented methods were consistent with the published sequence.     

Expressed sequence tag analysis of the diapausing queen of the bumblebee Bombus ignitus

We constructed a full‐length cDNA library from diapausing queens of the bumblebee Bombus ignitus. A total of 480 randomly selected clones was sequenced by single‐run 5′‐end sequencing. Of these, there were 437 high quality clones, 23 poor quality clones and 20 read‐fail clones. Each high quality clone sequence was searched against a public protein database. The most frequently found matching genes were ribosomal proteins (12.5%), p10 (3.58%), cytochrome P450 monooxygenase (3.13%) and sensory appendage protein (2.9%). Sequence similarity analysis between bumblebees and other insect species showed that 72 out of 437 (16.5%) bumblebee expressed sequence tags (EST) matched sequences of Apis mellifera, with matches to Drosophila melanogaster (6.6%), Caenorhabditis briggsae (6.2%), Lysiphlebus testaceipes (4.8%), Periplaneta americana (3.7%) and Anopheles gambiae (3.4%) following, suggesting that sequence similarity of bumblebee EST is closest to that of A. mellifera. Functional classification of EST based on Gene Ontology showed that most genes found by sequencing are associated with physiological processes in the bumblebee. The results of sequencing and analysis of our 437 cDNA demonstrated that high‐throughput EST sequencing and data analysis are powerful means for identifying novel genes and for expression profiling. Our bumblebee EST collection could be a useful platform for further studies of gene expression in diapausing bumblebees.