Metabolic Blockade-Based Genome Mining of Sea Anemone-Associated Streptomyces sp. S1502 Identifies Atypical Angucyclines WS-5995 A-E: Isolation, Identification, Biosynthetic Investigation, and Bioactivities

Marine symbiotic and epiphyte microorganisms are sources of bioactive or structurally novel natural products. Metabolic blockade-based genome mining has been proven to be an effective strategy to accelerate the discovery of natural products from both terrestrial and marine microorganisms. Here, the metabolic blockade-based genome mining strategy was applied to the discovery of other metabolites in a sea anemone-associated Streptomyces sp. S1502. We constructed a mutant Streptomyces sp. S1502/Δstp1 that switched to producing the atypical angucyclines WS-5995 A-E, among which WS-5995 E is a new compound. A biosynthetic gene cluster (wsm) of the angucyclines was identified through gene knock-out and heterologous expression studies. The biosynthetic pathways of WS-5995 A-E were proposed, the roles of some tailoring and regulatory genes were investigated, and the biological activities of WS-5995 A-E were evaluated. WS-5995 A has significant anti-Eimeria tenell activity with an IC50 value of 2.21 μM. The production of antibacterial streptopyrroles and anticoccidial WS-5995 A-E may play a protective role in the mutual relationship between Streptomyces sp. S1502 and its host.

A 69 kbp Deletion at the Berry Color Locus Is Responsible for Berry Color Recovery in Vitis vinifera L. Cultivar 'Riesling Rot'

‘Riesling Weiss’ is a white grapevine variety famous worldwide for fruity wines with higher acidity. Hardly known is ‘Riesling Rot’, a red-berried variant of ‘Riesling Weiss’ that disappeared from commercial cultivation but has increased in awareness in the last decades. The question arises of which variant, white or red, is the original and, consequently, which cultivar is the true ancestor. Sequencing the berry color locus of ‘Riesling Rot’ revealed a new VvmybA gene variant in one of the two haplophases called VvmybA3/1RR. The allele displays homologous recombination of VvmybA3 and VvmybA1 with a deletion of about 69 kbp between both genes that restores VvmybA1 transcripts. Furthermore, analysis of ‘Riesling Weiss’, ‘Riesling Rot’, and the ancestor ‘Heunisch Weiss’ along chromosome 2 using SSR (simple sequence repeat) markers elucidated that the haplophase of ‘Riesling Weiss’ was inherited from the white-berried parent variety ‘Heunisch Weiss’. Since no color mutants of ‘Heunisch Weiss’ are described that could have served as allele donors, we concluded that, in contrast to the public opinion, ‘Riesling Rot’ resulted from a mutational event in ‘Riesling Weiss’ and not vice versa.

Identification and Heterologous Expression of the Kendomycin B Biosynthetic Gene Cluster from Verrucosispora sp. SCSIO 07399

Verrucosispora sp. SCSIO 07399, a rare marine-derived actinomycete, produces a set of ansamycin-like polyketides kendomycin B–D (1–3) which possess potent antibacterial activities and moderate tumor cytotoxicity. Structurally, kendomycin B–D contain a unique aliphatic macrocyclic ansa scaffold in which the highly substituted pyran ring is connected to the quinone moiety. In this work, a type I/type III polyketide synthase (PKS) hybrid biosynthetic gene cluster coding for assembly of kendomycin B (kmy), and covering 33 open reading frames, was identified from Verrucosispora sp. SCSIO 07399. The kmy cluster was found to be essential for kendomycin B biosynthesis as verified by gene disruption and heterologous expression. Correspondingly, a biosynthetic pathway was proposed based on bioinformatics, cluster alignments, and previous research. Additionally, the role of type III PKS for generating the precursor unit 3,5-dihydroxybenzoic acid (3,5-DHBA) was demonstrated by chemical complementation, and type I PKS executed the polyketide chain elongation. The kmy cluster was found to contain a positive regulatory gene kmy4 whose regulatory effect was identified using real-time quantitative PCR (RT-qPCR). These advances shed important new insights into kendomycin B biosynthesis and help to set the foundation for further research aimed at understanding and exploiting the carbacylic ansa scaffold.

The Operon Encoding Hydrolytic Dehalogenation of 4-Chlorobenzoate Is Transcriptionally Regulated by the TetR-Type Repressor FcbR and Its Ligand 4-Chlorobenzoyl Coenzyme A

The bacterial hydrolytic dehalogenation of 4-chlorobenzoate (4CBA) is a coenzyme A (CoA)-activation-type catabolic pathway that is usually a common part of the microbial mineralization of chlorinated aromatic compounds. Previous studies have shown that the transport and dehalogenation genes for 4CBA are typically clustered as an fcbBAT1T2T3C operon and inducibly expressed in response to 4CBA. However, the associated molecular mechanism remains unknown. In this study, a gene (fcbR) adjacent to the fcb operon was predicted to encode a TetR-type transcriptional regulator in Comamonas sediminis strain CD-2. The fcbR knockout strain exhibited constitutive expression of the fcb cluster. In the host Escherichia coli, the expression of the P_fcb -fused green fluorescent protein (gfp) reporter was repressed by the introduction of the fcbR gene, and genetic studies combining various catabolic genes suggest that the ligand for FcbR may be an intermediate metabolite. Purified FcbR could bind to the P_fcb DNA probe in vitro, and the metabolite 4-chlorobenzyl-CoA (4CBA-CoA) prevented FcbR binding to the P _fcb DNA probe. Isothermal titration calorimetry (ITC) measurements showed that 4CBA-CoA could bind to FcbR at a 1:1 molar ratio. DNase I footprinting showed that FcbR protected a 42-bp DNA motif (5'-GGAAATCAATAGGTCCATAGAAAATCTATTGACTAATCGAAT-3') that consists of two sequence repeats containing four pseudopalindromic sequences (5'-TCNATNGA-3'). This binding motif overlaps with the -35 box of P_fcb and was proposed to prevent the binding of RNA polymerase. This study characterizes a transcriptional repressor of the fcb operon, together with its ligand, thus identifying halogenated benzoyl-CoA as belonging to the class of ligands of transcriptional regulators.IMPORTANCE The bacterial hydrolytic dehalogenation of 4CBA is a special CoA-activation-type catabolic pathway that plays an important role in the biodegradation of polychlorinated biphenyls and some herbicides. With genetic and biochemical approaches, the present study identified the transcriptional repressor and its cognate effector of a 4CBA hydrolytic dehalogenation operon. This work extends halogenated benzoyl-CoA as a new member of CoA-derived effector compounds that mediate allosteric regulation of transcriptional regulators.

Discovery, Bioactivity Evaluation, Biosynthetic Gene Cluster Identification, and Heterologous Expression of Novel Albofungin Derivatives

The crude extract of Streptomyces chrestomyceticus exhibited strong and broad activities against most “ESKAPE pathogens.” We conducted a comprehensive chemical investigation for secondary metabolites from the S. chrestomyceticus strain and identified two novel albofungin (alb) derivatives, i.e., albofungins A (1) and B (2), along with two known compounds, i.e., albofungin (3) and chloroalbofungin (4). The chemical structures of the novel compounds were elucidated using HRMS, 1D and 2D NMR, and electronic circular dichroism spectroscopy. The draft genome of S. chrestomyceticus was sequenced, and a 72 kb albofungin (alb) gene cluster with 72 open reading frames encoding type II polyketide synthases (PKSs), regulators, and transporters, and tailoring enzymes were identified using bioinformatics analysis. The alb gene cluster was confirmed using the heterologous expression in Streptomyces coelicolor, which successfully produced the compounds 3 and 4. Furthermore, compounds 1–4 displayed remarkable activities against Gram-positive bacteria and antitumor activities toward various cancer cells. Notably, compounds 1 and 3 showed potent activities against Gram-negative pathogenic bacteria. The terminal deoxynucleotidyl transferase (dUTP) nick-end labeling and flow cytometry analysis verified that compound 1 inhibited cancer cell proliferation by inducing cellular apoptosis. These results indicated that albofungins might be potential candidates for the development of antibiotics and antitumor drugs.

Revisiting the Pig IGHC Gene Locus in Different Breeds Uncovers Nine Distinct IGHG Genes

IgG subclass diversification is common in placental mammals. It has been well documented in humans and mice that different IgG subclasses, with diversified functions, synergistically regulate humoral immunity. However, our knowledge on the genomic and functional diversification of IgG subclasses in the pig, a mammalian species with high agricultural and biomedical importance, is incomplete. Using bacterial artificial chromosome sequencing and newly assembled genomes generated by the PacBio sequencing approach, we characterized and mapped the IgH C region gene locus in three indigenous Chinese breeds (Erhualian, Xiang, and Luchuan) and compared them to that of Duroc. Our data revealed that IGHG genes in Chinese pigs differ from the Duroc, whereas the IGHM, IGHD, IGHA, and IGHE genes were all single copy and highly conserved in the pig breeds examined. Most striking were differences in numbers of IGHG genes: there are seven genes in Erhualian pigs, six in the Duroc, but only five in Xiang pigs. Phylogenetic analysis suggested that all reported porcine IGHG genes could be classified into nine subclasses: IGHG1, IGHG2a, IGHG2b, IGHG2c, IGHG3, IGHG4, IGHG5a, IGHG5b, and IGHG5c. Using sequence information, we developed a mouse mAb specific for IgG3. This study offers a starting point to investigate the structure-function relationship of IgG subclasses in pigs.

High-throughput long paired-end sequencing of a Fosmid library by PacBio

BACKGROUND

Large insert paired-end sequencing technologies are important tools for assembling genomes, delineating associated breakpoints and detecting structural rearrangements. To facilitate the comprehensive detection of inter- and intra-chromosomal structural rearrangements or variants (SVs) and complex genome assembly with long repeats and segmental duplications, we developed a new method based on single-molecule real-time synthesis sequencing technology for generating long paired-end sequences of large insert DNA libraries.

RESULTS

A Fosmid vector, pHZAUFOS3, was developed with the following new features: (1) two 18-bp non-palindromic I-SceI sites flank the cloning site, and another two sites are present in the skeleton of the vector, allowing long DNA inserts (and the long paired-ends in this paper) to be recovered as single fragments and the vector (~ 8 kb) to be fragmented into 2-3 kb fragments by I-SceI digestion and therefore was effectively removed from the long paired-ends (5-10 kb); (2) the chloramphenicol (Cm) resistance gene and replicon (oriV), necessary for colony growth, are located near the two sides of the cloning site, helping to increase the proportion of the paired-end fragments to single-end fragments in the paired-end libraries. Paired-end libraries were constructed by ligating the size-selected, mechanically sheared pooled Fosmid DNA fragments to the Ampicillin (Amp) resistance gene fragment and screening the colonies with Cm and Amp. We tested this method on yeast and Setaria italica Yugu1. Fosmid-size paired-ends with an average length longer than 2 kb for each end were generated. The N50 scaffold lengths of the de novo assemblies of the yeast and S. italica Yugu1 genomes were significantly improved. Five large and five small structural rearrangements or assembly errors spanning tens of bp to tens of kb were identified in S. italica Yugu1 including deletions, inversions, duplications and translocations.

CONCLUSIONS

We developed a new method for long paired-end sequencing of large insert libraries, which can efficiently improve the quality of de novo genome assembly and identify large and small structural rearrangements or assembly errors.

Construction of a Genomic Bacterial Artificial Chromosome (BAC) Library for the Prawn Macrobrachium rosenbergii and Initial Analysis of ZW Chromosome-Derived BAC Inserts

Knowledge on sex determination has proven valuable for commercial production of the prawn Macrobrachium rosenbergii due to sex dimorphism of the male and female individuals. Previous studies indicated that prawn sex is determined by a ZW-ZZ chromosomal system, but no genomic information is available for the sex chromosome. Herein, we constructed a genomic bacterial artificial chromosome (BAC) library and identified the ZW-derived BAC clones for initial analysis of the sex chromosomal DNA sequence. The arrayed BAC library contains 200,448 clones with average insert size of 115.4 kb, corresponding to ∼ 4× coverage of the estimated 5.38 Gb genome. Based on a short female-specific marker, a Z- and a W-fragment were retrieved with the genomic walking method. Screening the BAC library using a ZW-specific marker as probe resulted in 12 positive clones. From these, a Z-derived (P331M17) and a W-derived (P122G2) BAC clones were randomly selected and sequenced by PacBio method. We report the construction of a large insert, deep-coverage, and high-quality BAC library for M. rosenbergii that provides a useful resource for positional cloning of target genes, genomic organization, and comparative genomics analysis. Our study not only confirmed the ZW/ZZ system but also discovered sex-linked genes on ZW chromosomes for the first time, contributing to a comprehensive understanding of the genomic structure of sex chromosomes in M. rosenbergii.

Heterologous expression of Avermectins biosynthetic gene cluster by construction of a Bacterial Artificial Chromosome library of the producers

Avermectins, a group of polyketide natural products, are widely used as anthelmintics in agriculture. Metabolic engineering and combinatorial biosynthesis were extensively employed to improve Avermectins production and create novel Avermectin derivatives, including Ivermectin and Doramectin. It is labor intensive and time cost to genetically manipulate Avermectins producer Streptomyces avermitilis in vivo. Cloning and heterologous expression of Avermectins biosynthetic gene cluster will make it possible to tailor the cluster in vitro. We constructed a Bacterial Artificial Chromosome (BAC) library of S. avermitilis ATCC 31267 with inserted DNA fragments ranged from 100 to 130 Kb. Five recombinant BAC clones which carried the Avermectins biosynthetic gene cluster ave (81 Kb in size) were screened out from the library. Then, ave was hetero-expressed in S. lividans. Three Avermectin components, A2a, B1a and A1a were detected from the cell extracts of recombinant strains. It will facilitate the development of Avermectin derivatives by polyketide synthase domain swapping and provide functional element for Avermectins synthetic biology study.

Global genomic diversity of Oryza sativa varieties revealed by comparative physical mapping

Bacterial artificial chromosome (BAC) physical maps embedding a large number of BAC end sequences (BESs) were generated for Oryza sativa ssp. indica varieties Minghui 63 (MH63) and Zhenshan 97 (ZS97) and were compared with the genome sequences of O. sativa spp. japonica cv. Nipponbare and O. sativa ssp. indica cv. 93-11. The comparisons exhibited substantial diversities in terms of large structural variations and small substitutions and indels. Genome-wide BAC-sized and contig-sized structural variations were detected, and the shared variations were analyzed. In the expansion regions of the Nipponbare reference sequence, in comparison to the MH63 and ZS97 physical maps, as well as to the previously constructed 93-11 physical map, the amounts and types of the repeat contents, and the outputs of gene ontology analysis, were significantly different from those of the whole genome. Using the physical maps of four wild Oryza species from OMAP (http://www.omap.org) as a control, we detected many conserved and divergent regions related to the evolution process of O. sativa. Between the BESs of MH63 and ZS97 and the two reference sequences, a total of 1532 polymorphic simple sequence repeats (SSRs), 71,383 SNPs, 1767 multiple nucleotide polymorphisms, 6340 insertions, and 9137 deletions were identified. This study provides independent whole-genome resources for intra- and intersubspecies comparisons and functional genomics studies in O. sativa. Both the comparative physical maps and the GBrowse, which integrated the QTL and molecular markers from GRAMENE (http://www.gramene.org) with our physical maps and analysis results, are open to the public through our Web site (http://gresource.hzau.edu.cn/resource/resource.html).

Genomic resources for gene discovery, functional genome annotation, and evolutionary studies of maize and its close relatives

Maize is one of the most important food crops and a key model for genetics and developmental biology. A genetically anchored and high-quality draft genome sequence of maize inbred B73 has been obtained to serve as a reference sequence. To facilitate evolutionary studies in maize and its close relatives, much like the Oryza Map Alignment Project (OMAP) (www.OMAP.org) bacterial artificial chromosome (BAC) resource did for the rice community, we constructed BAC libraries for maize inbred lines Zheng58, Chang7-2, and Mo17 and maize wild relatives Zea mays ssp. parviglumis and Tripsacum dactyloides. Furthermore, to extend functional genomic studies to maize and sorghum, we also constructed binary BAC (BIBAC) libraries for the maize inbred B73 and the sorghum landrace Nengsi-1. The BAC/BIBAC vectors facilitate transfer of large intact DNA inserts from BAC clones to the BIBAC vector and functional complementation of large DNA fragments. These seven Zea Map Alignment Project (ZMAP) BAC/BIBAC libraries have average insert sizes ranging from 92 to 148 kb, organellar DNA from 0.17 to 2.3%, empty vector rates between 0.35 and 5.56%, and genome equivalents of 4.7- to 8.4-fold. The usefulness of the Parviglumis and Tripsacum BAC libraries was demonstrated by mapping clones to the reference genome. Novel genes and alleles present in these ZMAP libraries can now be used for functional complementation studies and positional or homology-based cloning of genes for translational genomics.

Dynamic intra-japonica subspecies variation and resource application

We constructed a physical map of O. sativa ssp. japonica cv. ZH11 and compared it and its random sample sequences with the Nipponbare RefSeq derived from the same subspecies. This comparison showed that the two japonica genomes were highly syntenic but revealed substantial differences in terms of structural variations, rates of substitutions and indels, and transposable element content. For example, contractions/expansions as large as 450 kb and repeat sequences that were present in high copy numbers only in ZH11 were detected. In tri-alignment regions using the indica variety 93-11 sequence as an outgroup, we found that: (1) the substitution rates of the two japonica-indica inter-subspecies comparison combinations were close but almost a magnitude higher than the substitution rate between the japonica rice varieties ZH11 and Nipponbare; (2) of the substitutions found between ZH11 and Nipponbare, 47.2% occurred in ZH11 and 52.6% in Nipponbare; (3) of the indels found between ZH11 and Nipponbare, the indels that occurred in ZH11 were 15.8 times of those in Nipponbare. Of the indels that occurred in ZH11, 75.67% were insertions and 24.33% deletions. Of the indels that occurred in Nipponbare, 48.23% were insertions and 51.77% were deletions. The ZH11 comparative map covered four Nipponbare physical gaps, detected assembly errors in the Nipponbare sequence, and was integrated with the FSTs of a large ZH11 T-DNA insertion mutant library. ZH11 BAC clones can be browsed, searched, and obtained at our website, http://GResource.hzau.edu.cn.

Advancing Eucalyptus genomics: identification and sequencing of lignin biosynthesis genes from deep-coverage BAC libraries

Background

Eucalyptus species are among the most planted hardwoods in the world because of their rapid growth, adaptability and valuable wood properties. The development and integration of genomic resources into breeding practice will be increasingly important in the decades to come. Bacterial artificial chromosome (BAC) libraries are key genomic tools that enable positional cloning of important traits, synteny evaluation, and the development of genome framework physical maps for genetic linkage and genome sequencing.

Results

We describe the construction and characterization of two deep-coverage BAC libraries EG_Ba and EG_Bb obtained from nuclear DNA fragments of E. grandis (clone BRASUZ1) digested with HindIII and BstYI, respectively. Genome coverages of 17 and 15 haploid genome equivalents were estimated for EG_Ba and EG_Bb, respectively. Both libraries contained large inserts, with average sizes ranging from 135 Kb (Eg_Bb) to 157 Kb (Eg_Ba), very low extra-nuclear genome contamination providing a probability of finding a single copy gene ≥ 99.99%. Libraries were screened for the presence of several genes of interest via hybridizations to high-density BAC filters followed by PCR validation. Five selected BAC clones were sequenced and assembled using the Roche GS FLX technology providing the whole sequence of the E. grandis chloroplast genome, and complete genomic sequences of important lignin biosynthesis genes.

Conclusions

The two E. grandis BAC libraries described in this study represent an important milestone for the advancement of Eucalyptus genomics and forest tree research. These BAC resources have a highly redundant genome coverage (> 15×), contain large average inserts and have a very low percentage of clones with organellar DNA or empty vectors. These publicly available BAC libraries are thus suitable for a broad range of applications in genetic and genomic research in Eucalyptus and possibly in related species of Myrtaceae, including genome sequencing, gene isolation, functional and comparative genomics. Because they have been constructed using the same tree (E. grandis BRASUZ1) whose full genome is being sequenced, they should prove instrumental for assembly and gap filling of the upcoming Eucalyptus reference genome sequence.

Genome-wide BAC-end sequencing of Cucumis melo 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.

A genetic map of melon highly enriched with fruit quality QTLs and EST markers, including sugar and carotenoid metabolism genes

A genetic map of melon enriched for fruit traits was constructed, using a recombinant inbred (RI) population developed from a cross between representatives of the two subspecies of Cucumis melo L.: PI 414723 (subspecies agrestis) and 'Dulce' (subspecies melo). Phenotyping of 99 RI lines was conducted over three seasons in two locations in Israel and the US. The map includes 668 DNA markers (386 SSRs, 76 SNPs, six INDELs and 200 AFLPs), of which 160 were newly developed from fruit ESTs. These ESTs include candidate genes encoding for enzymes of sugar and carotenoid metabolic pathways that were cloned from melon cDNA or identified through mining of the International Cucurbit Genomics Initiative database (http://www.icugi.org/). The map covers 1,222 cM with an average of 2.672 cM between markers. In addition, a skeleton physical map was initiated and 29 melon BACs harboring fruit ESTs were localized to the 12 linkage groups of the map. Altogether, 44 fruit QTLs were identified: 25 confirming QTLs described using other populations and 19 newly described QTLs. The map includes QTLs for fruit sugar content, particularly sucrose, the major sugar affecting sweetness in melon fruit. Six QTLs interacting in an additive manner account for nearly all the difference in sugar content between the two genotypes. Three QTLs for fruit flesh color and carotenoid content were identified. Interestingly, no clear colocalization of QTLs for either sugar or carotenoid content was observed with over 40 genes encoding for enzymes involved in their metabolism. The RI population described here provides a useful resource for further genomics and metabolomics studies in melon, as well as useful markers for breeding for fruit quality.

Isoeugenol monooxygenase and its putative regulatory gene are located in the eugenol metabolic gene cluster in Pseudomonas nitroreducens Jin1

The plant-derived phenylpropanoids eugenol and isoeugenol have been proposed as useful precursors for the production of natural vanillin. Genes involved in the metabolism of eugenol and isoeugenol were clustered in region of about a 30 kb of Pseudomonas nitroreducens Jin1. Two of the 23 ORFs in this region, ORFs 26 (iemR) and 27 (iem), were predicted to be involved in the conversion of isoeugenol to vanillin. The deduced amino acid sequence of isoeugenol monooxygenase (Iem) of strain Jin1 had 81.4% identity to isoeugenol monooxygenase from Pseudomonas putida IE27, which also transforms isoeugenol to vanillin. Iem was expressed in E. coli BL21(DE3) and was found to lead to isoeugenol to vanillin transformation. Deletion and cloning analyses indicated that the gene iemR, located upstream of iem, is required for expression of iem in the presence of isoeugenol, suggesting it to be the iem regulatory gene. Reverse transcription, real-time PCR analyses indicated that the genes involved in the metabolism of eugenol and isoeugenol were differently induced by isoeugenol, eugenol, and vanillin.

Exploiting synteny in Cucumis for mapping of Psm: a unique locus controlling paternal mitochondrial sorting

The three genomes of cucumber show different modes of transmission, nuclear DNA bi-parentally, plastid DNA maternally, and mitochondrial DNA paternally. The mosaic (MSC) phenotype of cucumber is associated with mitochondrial DNA rearrangements and is a valuable tool for studying mitochondrial transmission. A nuclear locus (Psm) has been identified in cucumber that controls sorting of paternally transmitted mitochondrial DNA. Comparative sequencing and mapping of cucumber and melon revealed extensive synteny on the recombinational and sequence levels near Psm and placed this locus on linkage group R of cucumber and G10 of melon. However, the cucumber genomic region near Psm was surprisingly monomorphic with an average of one SNP every 25 kb, requiring that a family from a more diverse cross is produced for fine mapping and eventual cloning of Psm. The cucumber ortholog of Arabidopsis mismatch repair (MSH1) was cloned and it segregated independently of Psm, revealing that this candidate gene is not Psm.

Construction of a hexaploid wheat (Triticum aestivum L.) bacterial artificial chromosome library for cloning genes for stripe rust resistance

A hexaploid wheat (Triticum aestivum L.) bacterial artificial chromosome (BAC) library was constructed for cloning Yr5 and other genes conferring resistance to stripe rust (Puccinia striiformis f. sp. tritici). Intact nuclei from a Yr5 near-isogenic line were used to isolate high molecular weight DNA, which was partially cleaved with HindIII and cloned into pECBAC1 and pIndigoBAC-5 vectors. The wheat BAC library consisted of 422,400 clones arrayed in 1100 micro-titer plates (each plate with 384 wells). Random sampling of 300 BAC clones indicated an average insert size of 140 kb, with a size range from 25 to 365 kb. Ninety percent of the clones in the library had an insert size greater than 100 kb and fewer than 5% of the clones did not contain inserts. Based on an estimated genome size of 15,966 Mb for hexaploid wheat, the BAC library was estimated to have a total coverage of 3.58x wheat genome equivalents, giving approximately 96% probability of identifying a clone representing any given wheat DNA sequence. Twelve BAC clones containing an Yr5 locus-specific marker (Yr5STS7/8) were successfully selected by PCR screening of 3-dimensional BAC pools. The results demonstrated that the T. aestivum BAC library is a valuable genomic resource for positional cloning of Yr5. The library also should be useful in cloning other genes for stripe rust resistance and other traits of interest in hexaploid wheat.

New resources for marine genomics: bacterial artificial chromosome libraries for the Eastern and Pacific oysters (Crassostrea virginica and C. gigas)

Large-insert genomic bacterial artificial chromosome (BAC) libraries of two culturally and economically important oyster species, Crassostrea virginica and C. gigas, have been developed as part of an international effort to develop tools and reagents that will advance our ability to conduct genetic and genomic research. A total of 73,728 C. gigas clones with an average insert size of 152 kb were picked and arrayed representing an 11.8-fold genome coverage. A total of 55,296 clones with an average insert size of 150 kb were picked and arrayed for C. virginica, also representing an 11.8-fold genome coverage. The C. gigas and C. virginica libraries were screened with probes derived from selected oyster genes using high-density BAC colony filter arrays. The probes identified 4 to 25 clones per gene for C. virginica and 5 to 50 clones per gene for C. gigas. We conducted a preliminary analysis of genetic polymorphism represented in the C. gigas library. The results suggest that the degree of divergence among similar sequences is highly variable and concentrated in intronic regions. Evidence supporting allelic polymorphism is reported for two genes and allelic and/or locus specific polymorphism for several others. Classical inheritance studies are needed to confirm the nature of these polymorphisms. The oyster BAC libraries are publicly available to the research community on a cost-recovery basis at (www.genome.clemson.edu).

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.1x) 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.

Construction of a watermelon BAC library and identification of SSRs anchored to melon or Arabidopsis genomes

A bacterial artificial chromosome (BAC) library was constructed for watermelon (Citrullus lanatus (Thunb.) Matsum. & Nakai var. lanatus) with an average insert-size of 106 kb, providing 21 haploid genome equivalents. The library was used to identify BAC clones that are anchored to probes evenly distributed on the genomes of melon or Arabidopsis. Twenty eight probes (representing 66% of the tested probes) from melon and 30 probes (65%) from Arabidopsis identified positive BAC clones. Two methods were implemented to identify SSRs from the positively hybridizing BAC clones. First, analysis of BAC end sequences revealed 37 SSRs. For the second method, pooled DNA of BACs identified by the melon probes was used to develop a shotgun library. The library was then screened with synthetic SSR oligonucleotides by hybridization. Sequence analysis of positively hybridizing shotgun clones revealed 142 different SSRs. Thirty eight SSRs were characterized using three watermelon cultivars, five plant introduction (PI) accessions of C. lanatus var lanatus and four PIs of C. lanatus var citroides. Of these, 36 (95%) were found to be polymorphic with up to six alleles per marker. Polymorphism information content values for polymorphic markers varied between 0.22 and 0.79 with an average of 0.53. The methods described herein will be valuable for the construction of a watermelon linkage map with SSRs evenly distributed on its genome that is anchored to the genomes of melon and Arabidopsis.

Construction of a quinoa (Chenopodium quinoa Willd.) BAC library and its use in identifying genes encoding seed storage proteins

Quinoa (Chenopodium quinoa Willd.) is adapted to the harsh environments of the Andean Altiplano region. Its seeds have a well-balanced amino acid composition and exceptionally high protein content with respect to human nutrition. Quinoa grain is a staple in the diet of some of the most impoverished people in the world. The plant is an allotetraploid displaying disomic inheritance (2n=4x=36) with a di-haploid genome of 967 Mbp (megabase pair), or 2C=2.01 pg. We constructed two quinoa BAC libraries using BamHI (26,880 clones) and EcoRI (48,000 clones) restriction endonucleases. Cloned inserts in the BamHI library average 113 kb (kilobase) with approximately 2% of the clones lacking inserts, whereas cloned inserts in the EcoRI library average 130 kb and approximately 1% lack inserts. Three plastid genes used as probes of high-density arrayed blots of 73,728 BACs identified approximately 2.8% of the clones as containing plastid DNA inserts. We estimate that the combined quinoa libraries represent at least 9.0 di-haploid nuclear genome equivalents. An average of 12.2 positive clones per probe were identified with 13 quinoa single-copy ESTs as probes of the high-density arrayed blots, suggesting that the estimate of 9.0x coverage of the genome is conservative. Utility of the BAC libraries for gene identification was demonstrated by probing the library with a partial sequence of the 11S globulin seed storage protein gene and identifying multiple positive clones. The presence of the 11S globulin gene in four of the clones was verified by direct comparison with quinoa genomic DNA on a Southern blot. Besides serving as a useful tool for gene identification, the quinoa BAC libraries will be an important resource for physical mapping of the quinoa genome.

Physical map of the wheat high-grain protein content gene Gpc-B1 and development of a high-throughput molecular marker

Grain protein content (GPC) is important for human nutrition and has a strong influence on pasta and bread quality. A quantitative trait locus, derived from a Triticum turgidum ssp. dicoccoides accession (DIC), with an average increase in GPC of 14 g kg(-1) was mapped on chromosome 6BS. Using the wheat-rice colinearity, a high-density map of the wheat region was developed and the quantitative trait locus was mapped as a simple Mendelian locus designated Gpc-B1. A physical map of approx. 250 kb of the Gpc-B1 region was developed using a tetraploid wheat bacterial artificial chromosome library. The constructed physical map included the two Gpc-B1 flanking markers and one potential candidate gene from the colinear rice region completely linked to Gpc-B1. The relationship between physical and genetic distances and the feasibility of isolating genes by positional cloning in wheat are discussed. A high-throughput codominant marker, Xuhw89, was developed. A 4-bp deletion present in the DIC allele was absent in a collection of 117 cultivated tetraploid and hexaploid wheat germplasm, suggesting that this marker will be useful to incorporate the high GPC allele from the DIC accession studied here into commercial wheat varieties.

The Oryza bacterial artificial chromosome library resource: construction and analysis of 12 deep-coverage large-insert BAC libraries that represent the 10 genome types of the genus Oryza

Rice (Oryza sativa L.) is the most important food crop in the world and a model system for plant biology. With the completion of a finished genome sequence we must now functionally characterize the rice genome by a variety of methods, including comparative genomic analysis between cereal species and within the genus Oryza. Oryza contains two cultivated and 22 wild species that represent 10 distinct genome types. The wild species contain an essentially untapped reservoir of agriculturally important genes that must be harnessed if we are to maintain a safe and secure food supply for the 21st century. As a first step to functionally characterize the rice genome from a comparative standpoint, we report the construction and analysis of a comprehensive set of 12 BAC libraries that represent the 10 genome types of Oryza. To estimate the number of clones required to generate 10 genome equivalent BAC libraries we determined the genome sizes of nine of the 12 species using flow cytometry. Each library represents a minimum of 10 genome equivalents, has an average insert size range between 123 and 161 kb, an average organellar content of 0.4%-4.1% and nonrecombinant content between 0% and 5%. Genome coverage was estimated mathematically and empirically by hybridization and extensive contig and BAC end sequence analysis. A preliminary analysis of BAC end sequences of clones from these libraries indicated that LTR retrotransposons are the predominant class of repeat elements in Oryza and a roughly linear relationship of these elements with genome size was observed.

Random sheared fosmid library as a new genomic tool to accelerate complete finishing of rice (Oryza sativa spp. Nipponbare) genome sequence: sequencing of gap-specific fosmid clones uncovers new euchromatic portions of the genome

The International Rice Genome Sequencing Project has recently announced the high-quality finished sequence that covers nearly 95% of the japonica rice genome representing 370 Mbp. Nevertheless, the current physical map of japonica rice contains 62 physical gaps corresponding to approximately 5% of the genome, that have not been identified/represented in the comprehensive array of publicly available BAC, PAC and other genomic library resources. Without finishing these gaps, it is impossible to identify the complete complement of genes encoded by rice genome and will also leave us ignorant of some 5% of the genome and its unknown functions. In this article, we report the construction and characterization of a tenfold redundant, 40 kbp insert fosmid library generated by random mechanical shearing. We demonstrated its utility in refining the physical map of rice by identifying and in silico mapping 22 gap-specific fosmid clones with particular emphasis on chromosomes 1, 2, 6, 7, 8, 9 and 10. Further sequencing of 12 of the gap-specific fosmid clones uncovered unique rice genome sequence that was not previously reported in the finished IRGSP sequence and emphasizes the need to complete finishing of the rice genome.

Construction and characterization of a peanut HindIII BAC library

Bacterial artificial chromosome (BAC) libraries have been an essential tool for physical analyses of genomes of many crops. We constructed and characterized the first large-insert DNA library for Arachis hypogaea L. The HindIII BAC library contains 182,784 clones; only 5,484 (3%) had no inserts; and the average insert size is 104.05 kb. Chloroplast DNA contamination was very low, only nine clones, and r-DNA content was 1,208, 0.66% of clones. The depth of coverage is estimated to be 6.5 genome-equivalents, allowing the isolation of virtually any single-copy locus. This rate of coverage was confirmed with the application of 20 overgos, which identified 305 positive clones from the library. The identification of multiple loci by most probes in polyploids complicates anchoring of physical and genetic maps. We explored the practicality of a hybridization-based approach for determination of map locations of BAC clones in peanut by analyzing 94 clones detected by seven different overgos. The banding patterns on Southern blots were good predictors of contig composition; that is, the clones that shared the same size bands and ascribed to the same overgos usually also located in the same contigs. This BAC library has great potential to advance future research about the peanut genome.

Construction of two BAC libraries from cucumber (Cucumis sativus L.) and identification of clones linked to yield component quantitative trait loci

Two bacterial artificial chromosome (BAC) libraries were constructed from an inbred line derived from a cultivar of cucumber (Cucumis sativus L.). Intact nuclei were isolated and embedded in agarose plugs, and high-molecular-weight DNA was subsequently partially digested with BamHI or EcoRI. Ligation of double size-selected DNA fragments with the pECBAC1 vector yielded two libraries containing 23,040 BamHI and 18,432 EcoRI clones. The average BamHI and EcoRI insert sizes were estimated to be 107.0 kb and 100.8 kb, respectively, and BAC clones lacking inserts were 1.3% and 14.5% in the BamHI and EcoRI libraries, respectively. The two libraries together represent approximately 10.8 haploid cucumber genomes. Hybridization with a C(0)t-1 DNA probe revealed that approximately 36% of BAC clones likely carried repetitive sequence-enriched DNA. The frequencies of BAC clones that carry chloroplast or mitochondrial DNA range from 0.20% to 0.47%. Four sequence-characterized amplified region (SCAR), four simple sequence repeat, and an randomly amplified polymorphic DNA marker linked with yield component quantitative trait loci were used either as probes to hybridize high-density colony filters prepared from both libraries or as primers to screen an ordered array of pooled BAC DNA prepared from the BamHI library. Positive BAC clones were identified in predicted numbers, as screening by polymerase chain reaction amplification effectively overcame the problems associated with an overabundance of positives from hybridization with two SCAR markers. The BAC clones identified herein that are linked to the de (determinate habit) and F (gynoecy) locus will be useful for positional cloning of these economically important genes. These BAC libraries will also facilitate physical mapping of the cucumber genome and comparative genome analyses with other plant species.

Comparative mapping of ZYMV resistances in cucumber (Cucumis sativus L.) and melon (Cucumis melo L.)

Zucchini yellow mosaic virus (ZYMV) routinely causes significant losses in cucumber ( Cucumis sativus L.) and melon ( Cucumis melo L.). ZYMV resistances from the cucumber population 'TMG1' and the melon plant introduction (PI) 414723 show different modes of inheritance and their genetic relationships are unknown. We used molecular markers tightly linked to ZYMV resistances from cucumber and melon for comparative mapping. A 5-kb genomic region (YCZ-5) cosegregating with the zym locus of cucumber was cloned and sequenced to reveal single nucleotide polymorphisms and indels distinguishing alleles from ZYMV-resistant (TMG1) and susceptible (Straight 8) cucumbers. A low-copy region of the YCZ-5 clone was hybridized to bacterial artificial chromosome (BAC) clones of melon and a 180-kb contig assembled. One end of this melon contig was mapped in cucumber and cosegregated with ZYMV resistance, demonstrating that physically linked regions in melon show genetic linkage in cucumber. However the YCZ-5 region segregated independently of ZYMV resistance loci in two melon families. These results establish that these sources of ZYMV resistances from cucumber TMG1 and melon PI414723 are likely non-syntenic.

The Fusarium wilt resistance locus Fom-2 of melon contains a single resistance gene with complex features

The soil-borne fungus Fusarium oxysporum f.sp. melonis causes significant losses in the cultivated melon, a key member of the economically important family, the Cucurbitaceae. Here, we report the map-based cloning and characterization of the resistance gene Fom-2 that confers resistance to race 0 and 1 of this plant pathogen. Two recombination events, 75 kb apart, were found to bracket Fom-2 after screening approximately 1324 gametes with PCR-based markers. Sequence analysis of the Fom-2 interval revealed the presence of two candidate genes. One candidate gene showed significant similarity to previously characterized resistance genes. Sequence analysis of this gene revealed clear polymorphisms between resistant and susceptible materials and was therefore designated as Fom-2. Analysis of susceptible breeding lines (BL) presenting a haplotype very similar to the resistant cultivar MR-1 indicated that a gene conversion had occurred in Fom-2, resulting in a significant rearrangement of this gene. The second candidate gene which shared high similarity to an essential gene in Arabidopsis, presented an almost identical sequence in MR-1 and BL, further supporting Fom-2 identity. The gene conversion in Fom-2 produced a truncated R gene, revealing new insights into R gene evolution. Fom-2 was predicted to encode an NBS-LRR type R protein of the non-TIR subfamily. In contrast to most members of this class a coiled-coil structure was predicted within the LRR region rather than in the N-terminal. The Fom-2 physical region contained retroelement-like sequences and truncated genes, suggesting that this locus is complex.

Cloning of disease-resistance homologues in end sequences of BAC clones linked to Fom-2, a gene conferring resistance to Fusarium wilt in melon (Cucumis melo L.)

Disease resistance has not yet been characterized at the molecular level in cucurbits, a group of high-value, nutritious, horticultural plants. Previously, we genetically mapped the Fom-2 gene that confers resistance to Fusarium wilt races 0 and I of melon. In this paper, two cosegregating codominant markers (AM, AFLP marker; FM, Fusarium marker) were used to screen a melon bacterial artificial chromosome (BAC) library. Identified clones were fingerprinted and end sequenced. Fingerprinting analysis showed that clones identified by each marker assembled into two separate contigs at high stringency. GenBank searches produced matches to leucine-rich repeats (LRRs) of resistance genes (R genes); to retroelements and to cellulose synthase in clones identified by FM; and to nucleotide-binding sites (NBSs) of R genes, retroelements, and cytochrome P-450 in clones identified by AM. A 6.5-kb fragment containing both NBS and LRR sequences was found to share high homology to TIR (Toll-interleukin-1 receptor)-NBS-LRR R genes, such as N, with 42% identity and 58% similarity in the TIR-NBS and LRR regions. The sequence information may be useful for identifying NBS-LRR class of R genes in other cucurbits.