Complete genome of the mutualistic, N2-fixing grass endophyte Azoarcus sp. strain BH72

Azoarcus sp. strain BH72, a mutualistic endophyte of rice and other grasses, is of agrobiotechnological interest because it supplies biologically fixed nitrogen to its host and colonizes plants in remarkably high numbers without eliciting disease symptoms. The complete genome sequence is 4,376,040-bp long and contains 3,992 predicted protein-coding sequences. Genome comparison with the Azoarcus-related soil bacterium strain EbN1 revealed a surprisingly low degree of synteny. Coding sequences involved in the synthesis of surface components potentially important for plant-microbe interactions were more closely related to those of plant-associated bacteria. Strain BH72 appears to be 'disarmed' compared to plant pathogens, having only a few enzymes that degrade plant cell walls; it lacks type III and IV secretion systems, related toxins and an N-acyl homoserine lactones-based communication system. The genome contains remarkably few mobile elements, indicating a low rate of recent gene transfer that is presumably due to adaptation to a stable, low-stress microenvironment.     

Transcriptomic analysis of responses to exudates reveal genes required for rhizosphere competence of the endophyte Azoarcus sp. strain BH72

Endophytic colonization is a very complex process which is not yet completely understood. Molecules exuded by the plants may act as signals which influence the ability of the microbe to colonize the host or survive in the rhizosphere. Here we used the whole genome microarray approach to investigate the response of the diazotrophic model endophyte, Azoarcus sp. strain BH72, to exudates of O.?sativa cv. Nipponbare in order to identify differentially regulated genes. On exposure to exudates, an overall expression of 4.4% of the 3992 protein coding genes of Azoarcus sp. strain BH72 was altered, out of which 2.4% was upregulated and 2.0% was downregulated. Genes with modulated expression included a few whose involvement in plant-microbe interaction had already been established, whereas a large fraction comprised of genes encoding proteins with putative or unknown functions. Mutational analysis of several differentially regulated genes like those encoding a minor pilin PilX, signal transduction proteins containing GGDEF domains and a serine-threonine kinase as a putative component of the type IV secretion system (T6SS), revealed their role in host colonization. Our data suggest that strain BH72 may be primed for the endophytic lifestyle by exudates, as the expression of bacterial genes relevant for endophytic colonization of roots is induced by root exudates.     

The Korea brassica genome project: a glimpse of the brassica genome based on comparative genome analysis with Arabidopsis

A complete genome sequence provides unlimited information in the sequenced organism as well as in related taxa. According to the guidance of the Multinational Brassica Genome Project (MBGP), the Korea Brassica Genome Project (KBGP) is sequencing chromosome 1 (cytogenetically oriented chromosome #1) of Brassica rapa. We have selected 48 seed BACs on chromosome 1 using EST genetic markers and FISH analyses. Among them, 30 BAC clones have been sequenced and 18 are on the way. Comparative genome analyses of the EST sequences and sequenced BAC clones from Brassica chromosome 1 revealed their homeologous partner regions on the Arabidopsis genome and a syntenic comparative map between Brassica chromosome 1 and Arabidopsis chromosomes. In silico chromosome walking and clone validation have been successfully applied to extending sequence contigs based on the comparative map and BAC end sequences. In addition, we have defined the (peri)centromeric heterochromatin blocks with centromeric tandem repeats, rDNA and centromeric retrotransposons. In-depth sequence analyses of five homeologous BAC clones and an Arabidopsis chromosomal region reveal overall co-linearity, with 82% sequence similarity. The data indicate that the Brassica genome has undergone triplication and subsequent gene losses after the divergence of Arabidopsis and Brassica. Based on in-depth comparative genome analyses, we propose a comparative genomics approach for conquering the Brassica genome. In 2005 we intend to construct an integrated physical map, including sequence information from 500 BAC clones and integration of fingerprinting data and end sequence data of more than 100 000 BAC clones. The sequences have been submitted to GenBank with accession numbers: 10 204 BAC ends of the KBrH library (CW978640-CW988843); KBrH138P04, AC155338; KBrH117N09, AC155337; KBrH097M21, AC155348; KBrH093K03, AC155347; KBrH081N08, AC155346; KBrH080L24, AC155345; KBrH077A05, AC155343; KBrH020D15, AC155340; KBrH015H17, AC155339; KBrH001H24, AC155335; KBrH080A08, AC155344; KBrH004D11, AC155341; KBrH117M18, AC146875; KBrH052O08, AC155342.     

Complete genome sequence of the denitrifying and N2O-reducing bacterium Azoarcus sp. strain KH32C

We report the finished and annotated genome sequence of a denitrifying and N(2)O-reducing betaproteobacterium, Azoarcus sp. strain KH32C. The genome is composed of one chromosome and one megaplasmid and contains genes for plant-microbe interactions and the gene clusters for aromatic-compound degradations.     

A BAC-based physical map of the channel catfish genome

Catfish is the major aquaculture species in the United States. To enhance its genome studies involving genetic linkage and comparative mapping, a bacterial artificial chromosome (BAC) contig-based physical map of the channel catfish (Ictalurus punctatus) genome was generated using four-color fluorescence-based fingerprints. Fingerprints of 34,580 BAC clones (5.6x genome coverage) were generated for the FPC assembly of the BAC contigs. A total of 3307 contigs were assembled using a cutoff value of 1x10(-20). Each contig contains an average of 9.25 clones with an average size of 292 kb. The combined contig size for all contigs was 0.965 Gb, approximately the genome size of the channel catfish. The reliability of the contig assembly was assessed by both hybridization of gene probes to BAC clones contained in the fingerprinted assembly and validation of randomly selected contigs using overgo probes designed from BAC end sequences. The presented physical map should greatly enhance genome research in the catfish, particularly aiding in the identification of genomic regions containing genes underlying important performance traits.     

Root colonization and systemic spreading of Azoarcus sp. strain BH72 in grasses.

The invasive properties of Azoarcus sp. strain BH72, an endorhizospheric isolate of Kallar grass, on gnotobiotically grown seedlings of Oryza sativa IR36 and Leptochloa fusca (L.) Kunth were studied. Additionally, Azoarcus spp. were localized in roots of field-grown Kallar grass. To facilitate localization and to assure identity of bacteria, genetically engineered microorganisms expressing beta-glucuronidase were also used as inocula. beta-Glucuronidase staining indicated that the apical region of the root behind the meristem was the most intensively colonized. Light and electron microscopy showed that strain BH72 penetrated the rhizoplane preferentially in the zones of elongation and differentiation and colonized the root interior inter- and intracellularly. In addition to the root cortex, stelar tissue was also colonized; bacteria were found in the xylem. No evidence was obtained that Azoarcus spp. could reside in living plant cells; rather, plant cells were apparently destroyed after bacteria had penetrated the cell wall. A common pathogenicity test on tobacco leaves provided no evidence that representative strains of Azoarcus spp. are phytopathogenic. Compared with the control, inoculation with strain BH72 significantly promoted growth of rice seedlings. This effect was reversed when the plant medium was supplemented with malate (0.2 g/liter). N2 fixation was apparently not involved, because the same response was obtained with a nifK mutant of strain BH72, which has a Nif- phenotype. Also, Western blot (immunoblot) analysis of protein extracts from rice seedlings gave no indication that nitrogenase was present. PCR and Western immunoblotting, using primers specific for eubacteria and antibodies recognizing type-specific antigens, respectively, indicated that strain BH72 could colonize rice plants systemically, probably mediated by longitudinal spreading through vessels.     

Twitching motility is essential for endophytic rice colonization by the N2-fixing endophyte Azoarcus sp. strain BH72

Azoarcus sp. strain BH72, as an endophyte of grasses, depends on successful host colonization. Type IV pili are essential for mediating the initial interaction with rice roots. In the genome sequence analysis, the pilT gene was identified, which encodes for a putative type IV pilus retraction protein. PilT of Azoarcus sp. BH72 shares high similarity to PilT of the human pathogen Pseudomonas aeruginosa PAO1 (77% amino acid sequence identity) and contains a predicted nucleotide-binding motif. To gain more insights into the role of the type IV pili in the colonization process of Azoarcus spp., we constructed an insertional mutant of pilT and a deletion mutant of pilA, the major structural component of the pilus structure. The pilT mutant, as the pilin deletion mutant deltapilA, was abolished in twitching motility. Western blot analyses and electron microscopy studies demonstrated an enhanced piliation of the Azoarcus pilT mutant strain compared with the wild type, indicating that, indeed, PilT has a role in pilus retraction. Studies on rice root colonization in gnotobiotic cultures revealed that the establishment of microcolonies on the root surface was strongly reduced in the deltapilA mutant, whereas the surface colonization was reduced by only 50% in the nontwitching pilT mutant. However, endophytic colonization of rice roots was strongly reduced in both mutants. These results demonstrate that the retractile force mediated by PilT is not essential for the bacterial colonization of the plant surface, but that twitching motility is necessary for invasion of and establishment inside the plant. Thus, a novel determinant for endophytic interactions with grasses was identified.     

BioNano genome mapping of individual chromosomes supports physical mapping and sequence assembly in complex plant genomes

The assembly of a reference genome sequence of bread wheat is challenging due to its specific features such as the genome size of 17 Gbp, polyploid nature and prevalence of repetitive sequences. BAC‐by‐BAC sequencing based on chromosomal physical maps, adopted by the International Wheat Genome Sequencing Consortium as the key strategy, reduces problems caused by the genome complexity and polyploidy, but the repeat content still hampers the sequence assembly. Availability of a high‐resolution genomic map to guide sequence scaffolding and validate physical map and sequence assemblies would be highly beneficial to obtaining an accurate and complete genome sequence. Here, we chose the short arm of chromosome 7D (7DS) as a model to demonstrate for the first time that it is possible to couple chromosome flow sorting with genome mapping in nanochannel arrays and create a de novo genome map of a wheat chromosome. We constructed a high‐resolution chromosome map composed of 371 contigs with an N50 of 1.3 Mb. Long DNA molecules achieved by our approach facilitated chromosome‐scale analysis of repetitive sequences and revealed a ~800‐kb array of tandem repeats intractable to current DNA sequencing technologies. Anchoring 7DS sequence assemblies obtained by clone‐by‐clone sequencing to the 7DS genome map provided a valuable tool to improve the BAC‐contig physical map and validate sequence assembly on a chromosome‐arm scale. Our results indicate that creating genome maps for the whole wheat genome in a chromosome‐by‐chromosome manner is feasible and that they will be an affordable tool to support the production of improved pseudomolecules.     

A bacterial artificial chromosome based physical map of the Ustilago maydis genome.

Ustilago maydis, a basidiomycete, is a model organism among phytopathogenic fungi. A physical map of U. maydis strain 521 was developed from bacterial artificial chromosome (BAC) clones. BAC fingerprints used polyacrylamide gel electrophoresis to separate restriction fragments. Fragments were labeled at the HindIII site and co-digested with HaeIII to reduce fragments to 50-750 bp. Contiguous overlapping sets of clones (contigs) were assembled at nine stringencies (from P < or = 1 x 10(-6) to 1 x 10(-24)). Each assembly nucleated contigs with different percentages of bands overlapping between clones (from 20% to 97%). The number of clones per contig decreased linearly from 41 to 12 from P < or = 1 x 10(-7) to 1 x 10 (-12). The number of separate contigs increased from 56 to 150 over the same range. A hybridization-based physical map of the same BAC clones was compared with the fingerprint contigs built at P < or = 1 x 10(-7). The two methods provided consistent physical maps that were largely validated by genome sequence. The combined hybridization and fingerprint physical map provided a minimum tile path composed of 258 BAC clones (18-20 Mbp) distributed among 28 merged contigs. The genome of U. maydis was estimated to be 20.5 Mbp by pulsed-field gel electrophoresis and 24 Mbp by BAC fingerprints. There were 23 separate chromosomes inferred by both pulsed-field gel electrophoresis and fingerprint contigs. Only 11 of the tile path BAC clones contained recognizable centromere, telomere, and subtelomere repeats (high-copy DNA), suggesting that repeats caused some false merges. There were 247 tile path BAC clones that encompassed about 17.5 Mbp of low-copy DNA sequence. BAC clones are available for repeat and unique gene cluster analysis including tDNA-mediated transformation. Program FingerPrint Contigs maps aligned with each chromosome can be viewed at http://www.siu.edu/~meksem/ustilago_maydis/.     

A physical map of the highly heterozygous Populus genome: integration with the genome sequence and genetic map and analysis of haplotype variation

As part of a larger project to sequence the Populus genome and generate genomic resources for this emerging model tree, we constructed a physical map of the Populus genome, representing one of the few such maps of an undomesticated, highly heterozygous plant species. The physical map, consisting of 2802 contigs, was constructed from fingerprinted bacterial artificial chromosome (BAC) clones. The map represents approximately 9.4-fold coverage of the Populus genome, which has been estimated from the genome sequence assembly to be 485 ± 10 Mb in size. BAC ends were sequenced to assist long-range assembly of whole-genome shotgun sequence scaffolds and to anchor the physical map to the genome sequence. Simple sequence repeat-based markers were derived from the end sequences and used to initiate integration of the BAC and genetic maps. A total of 2411 physical map contigs, representing 97% of all clones assigned to contigs, were aligned to the sequence assembly (JGI Populus trichocarpa , version 1.0). These alignments represent a total coverage of 384 Mb (79%) of the entire poplar sequence assembly and 295 Mb (96%) of linkage group sequence assemblies. A striking result of the physical map contig alignments to the sequence assembly was the co-localization of multiple contigs across numerous regions of the 19 linkage groups. Targeted sequencing of BAC clones and genetic analysis in a small number of representative regions showed that these co-aligning contigs represent distinct haplotypes in the heterozygous individual sequenced, and revealed the nature of these haplotype sequence differences.     

An integrated BAC and genome sequence physical map of Phytophthora sojae

Phytophthora spp. are serious pathogens that threaten numerous cultivated crops, trees, and natural vegetation worldwide. The soybean pathogen P. sojae has been developed as a model oomycete. Here, we report a bacterial artificial chromosome (BAC)-based, integrated physical map of the P. sojae genome. We constructed two BAC libraries, digested 8,681 BACs with seven restriction enzymes, end labeled the digested fragments with four dyes, and analyzed them with capillary electrophoresis. Fifteen data sets were constructed from the fingerprints, using individual dyes and all possible combinations, and were evaluated for contig assembly. In all, 257 contigs were assembled from the XhoI data set, collectively spanning approximately 132 Mb in physical length. The BAC contigs were integrated with the draft genome sequence of P. sojae by end sequencing a total of 1,440 BACs that formed a minimal tiling path. This enabled the 257 contigs of the BAC map to be merged with 207 sequence scaffolds to form an integrated map consisting of 79 superscaffolds. The map represents the first genome-wide physical map of a Phytophthora sp. and provides a valuable resource for genomics and molecular biology research in P. sojae and other Phytophthora spp. In one illustration of this value, we have placed the 350 members of a superfamily of putative pathogenicity effector genes onto the map, revealing extensive clustering of these genes.     

An overview of the apple genome through BAC end sequence analysis

The apple, Malus x domestica Borkh., is one of the most important fruit trees grown worldwide. A bacterial artificial chromosome (BAC)-based physical map of the apple genome has been recently constructed. Based on this physical map, a total of approximately 2,100 clones from different contigs (overlapping BAC clones) have been selected and sequenced at both ends, generating 3,744 high-quality BAC end sequences (BESs) including 1,717 BAC end pairs. Approximately 8.5% of BESs contain simple sequence repeats (SSRs), most of which are AT/TA dimer repeats. Potential transposable elements are identified in approximately 21% of BESs, and most of these elements are retrotransposons. About 11% of BESs have homology to the Arabidopsis protein database. The matched proteins cover a broad range of categories. The average GC content of the predicted coding regions of BESs is 42.4%; while, that of the whole BESs is 39%. A small number of BES pairs were mapped to neighboring chromosome regions of A. thaliana and Populus trichocarpa; whereas, no pairs are mapped to the Oryza sativa genome. The apple has a higher degree of synteny with the closely related Populus than with the distantly related Arabidopsis. BAC end sequencing can be used to anchor a small proportion of the apple genome to the Populus and possibly to the Arabidopsis genomes.     

A comparative map of bovine chromosome 19 based on a combination of mapping on a bacterial artificial chromosome scaffold map, a whole genome radiation hybrid panel and the human draft sequence

We have constructed a medium density physical map of bovine chromosome 19 using a combination of mapping loci on both a bovine bacterial artificial chromosome (BAC) scaffold map and a whole genome radiation hybrid (WGRH) panel. The resulting map contains 70 loci spanning the length of bovine chromosome 19. Three contiguous groups of BACs were identified on the basis of multiple loci mapping to individual BAC clones. Bovine chromosome 19 was found in this study to be comprised almost entirely from regions of human chromosome 17, with a small region putatively assigned to human chromosome 10. Fourteen breakpoints between the bovine and human chromosomes were detected, with a possibility of five more based on ordering of the WGRH map.     

Using comparative genomics to reorder the human genome sequence into a virtual sheep genome

Background

Is it possible to construct an accurate and detailed subgene-level map of a genome using bacterial artificial chromosome (BAC) end sequences, a sparse marker map, and the sequences of other genomes?

Results

A sheep BAC library, CHORI-243, was constructed and the BAC end sequences were determined and mapped with high sensitivity and low specificity onto the frameworks of the human, dog, and cow genomes. To maximize genome coverage, the coordinates of all BAC end sequence hits to the cow and dog genomes were also converted to the equivalent human genome coordinates. The 84,624 sheep BACs (about 5.4-fold genome coverage) with paired ends in the correct orientation (tail-to-tail) and spacing, combined with information from sheep BAC comparative genome contigs (CGCs) built separately on the dog and cow genomes, were used to construct 1,172 sheep BAC-CGCs, covering 91.2% of the human genome. Clustered non-tail-to-tail and outsize BACs located close to the ends of many BAC-CGCs linked BAC-CGCs covering about 70% of the genome to at least one other BAC-CGC on the same chromosome. Using the BAC-CGCs, the intrachromosomal and interchromosomal BAC-CGC linkage information, human/cow and vertebrate synteny, and the sheep marker map, a virtual sheep genome was constructed. To identify BACs potentially located in gaps between BAC-CGCs, an additional set of 55,668 sheep BACs were positioned on the sheep genome with lower confidence. A coordinate conversion process allowed us to transfer human genes and other genome features to the virtual sheep genome to display on a sheep genome browser.

Conclusion

We demonstrate that limited sequencing of BACs combined with positioning on a well assembled genome and integrating locations from other less well assembled genomes can yield extensive, detailed subgene-level maps of mammalian genomes, for which genomic resources are currently limited.     

Comparative physical mapping between Oryza sativa (AA genome type) and O. punctata (BB genome type)

A comparative physical map of the AA genome (Oryza sativa) and the BB genome (O. punctata) was constructed by aligning a physical map of O. punctata, deduced from 63,942 BAC end sequences (BESs) and 34,224 fingerprints, onto the O. sativa genome sequence. The level of conservation of each chromosome between the two species was determined by calculating a ratio of BES alignments. The alignment result suggests more divergence of intergenic and repeat regions in comparison to gene-rich regions. Further, this characteristic enabled localization of heterochromatic and euchromatic regions for each chromosome of both species. The alignment identified 16 locations containing expansions, contractions, inversions, and transpositions. By aligning 40% of the punctata BES on the map, 87% of the punctata FPC map covered 98% of the O. sativa genome sequence. The genome size of O. punctata was estimated to be 8% larger than that of O. sativa with individual chromosome differences of 1.5-16.5%. The sum of expansions and contractions observed in regions >500 kb were similar, suggesting that most of the contractions/expansions contributing to the genome size difference between the two species are small, thus preserving the macro-collinearity between these species, which diverged approximately 2 million years ago.     

Integration of physical and genetic maps in apple confirms whole-genome and segmental duplications in the apple genome

A total of 355 simple sequence repeat (SSR) markers were developed, based on expressed sequence tag (EST) and bacterial artificial chromosome (BAC)-end sequence databases, and successfully used to construct an SSR-based genetic linkage map of the apple. The consensus linkage map spanned 1143 cM, with an average density of 2.5 cM per marker. Newly developed SSR markers along with 279 SSR markers previously published by the HiDRAS project were further used to integrate physical and genetic maps of the apple using a PCR-based BAC library screening approach. A total of 470 contigs were unambiguously anchored onto all 17 linkage groups of the apple genome, and 158 contigs contained two or more molecular markers. The genetically mapped contigs spanned ～421 Mb in cumulative physical length, representing 60.0% of the genome. The sizes of anchored contigs ranged from 97 kb to 4.0 Mb, with an average of 995 kb. The average physical length of anchored contigs on each linkage group was ～24.8 Mb, ranging from 17.0 Mb to 37.73 Mb. Using BAC DNA as templates, PCR screening of the BAC library amplified fragments of highly homologous sequences from homoeologous chromosomes. Upon integrating physical and genetic maps of the apple, the presence of not only homoeologous chromosome pairs, but also of multiple locus markers mapped to adjacent sites on the same chromosome was detected. These findings demonstrated the presence of both genome-wide and segmental duplications in the apple genome and provided further insights into the complex polyploid ancestral origin of the apple.     

A comprehensive radiation hybrid map of the bovine genome comprising 5593 loci

A bovine whole genome 7000-rad radiation hybrid (RH) panel, SUNbRH(7000-rad), was constructed to build a high-resolution RH map. The Shirakawa-USDA linkage map served as a scaffold to construct a framework map of 3216 microsatellites on which 2377 ESTs were ordered. The resulting RH map provided essentially complete coverage across the genome, with 1 cR7000 corresponding to 114 kb, and a cattle-human comparative map of 1716 bovine genes and sequences annotated in the human genome, which covered 79 and 72% of the bovine and human genomes, respectively. We then integrated the bovine RH and comparative maps with BAC fingerprint information in to construct a detailed, BAC-based physical map covering a reported 40-cM quantitative trait locus region for intramuscular fat or "marbling" on BTA 4. In summary, the new, high-resolution SUNbRH7000-rad, comparative, Shirakawa-USDA linkage, and BAC fingerprint maps provide a set of genomic tools for fine mapping regions of interest in cattle.