An integrated BAC and genome sequence physical map of Phytophthora sojae

Karyotype variation, spontaneous genome rearrangements affecting chemical insensitivity, and expression level polymorphisms in the plant pathogen Phytophthora infestans revealed using its first chromosome-scale assembly

Natural isolates of the potato and tomato pathogen Phytophthora infestans exhibit substantial variation in virulence, chemical sensitivity, ploidy, and other traits. A chromosome-scale assembly was developed to expand genomic resources for this oomyceteous microbe, and used to explore the basis of variation. Using PacBio and Illumina data, a long-range linking library, and an optical map, an assembly was created and coalesced into 15 pseudochromosomes spanning 219 Mb using SNP-based genetic linkage data. De novo gene prediction combined with transcript evidence identified 19,981 protein-coding genes, plus about eight thousand tRNA genes. The chromosomes were comprised of a mosaic of gene-rich and gene-sparse regions plus very long centromeres. Genes exhibited a biased distribution across chromosomes, especially members of families encoding RXLR and CRN effectors which clustered on certain chromosomes. Strikingly, half of F1 progeny of diploid parents were polyploid or aneuploid. Substantial expression level polymorphisms between strains were identified, much of which could be attributed to differences in chromosome dosage, transposable element insertions, and adjacency to repetitive DNA. QTL analysis identified a locus on the right arm of chromosome 3 governing sensitivity to the crop protection chemical metalaxyl. Strains heterozygous for resistance often experienced megabase-sized deletions of that part of the chromosome when cultured on metalaxyl, increasing resistance due to loss of the sensitive allele. This study sheds light on diverse phenomena affecting variation in P. infestans and relatives, helps explain the prevalence of polyploidy in natural populations, and provides a new foundation for biologic and genetic investigations.

Identification of Chromosomes and Chromosome Rearrangements in Crop Brassicas and Raphanus sativus: A Cytogenetic Toolkit Using Synthesized Massive Oligonucleotide Libraries

Crop brassicas include three diploid [Brassica rapa (AA; 2n = 2x = 16), B. nigra (BB; 2n = 2x = 18), and B. oleracea (CC; 2n = 2x = 20)] and three derived allotetraploid species. It is difficult to distinguish Brassica chromosomes as they are small and morphologically similar. We aimed to develop a genome-sequence based cytogenetic toolkit for reproducible identification of Brassica chromosomes and their structural variations. A bioinformatic pipeline was used to extract repeat-free sequences from the whole genome assembly of B. rapa. Identified sequences were subsequently used to develop four c. 47-mer oligonucleotide libraries comprising 27,100, 11,084, 9,291, and 16,312 oligonucleotides. We selected these oligonucleotides after removing repeats from 18 identified sites (500-1,000 kb) with 1,997-5,420 oligonucleotides localized at each site in B. rapa. For one set of probes, a new method for amplification or immortalization of the library is described. oligonucleotide probes produced specific and reproducible in situ hybridization patterns for all chromosomes belonging to A, B, C, and R (Raphanus sativus) genomes. The probes were able to identify structural changes between the genomes, including translocations, fusions, and deletions. Furthermore, the probes were able to identify a structural translocation between a pak choi and turnip cultivar of B. rapa. Overall, the comparative chromosomal mapping helps understand the role of chromosome structural changes during genome evolution and speciation in the family Brassicaceae. The probes can also be used to identify chromosomes in aneuploids such as addition lines used for gene mapping, and to track transfer of chromosomes in hybridization and breeding programs.

Pathotype and Genetic Shifts in a Population of Phytophthora sojae Under Soybean Cultivar Rotation

Changes in pathotype structure of Phytophthora sojae populations have been attributed to deployment of race-specific resistant Rps genes in soybean that have been incorporated into commercial cultivars to reduce losses due to Phytophthora root and stem rot. To test this hypothesis, a cultivar rotation study was established from 2007 through 2010 in microplots at a site in Iowa with no history of soybean cultivation. All microplots were inoculated with P. sojae isolate PR1, race 1 (vir 7) prior to planting in year 1. Six rotations were tested: (i) continuous planting of a P. sojae-susceptible cultivar, (ii) continuous planting of a cultivar with high partial resistance to the pathogen, (iii) continuous planting of a cultivar with the Rps 1k gene, (iv) annual rotation of a susceptible with a resistant cultivar, (v) annual rotation of a partially resistant cultivar with a cultivar with the Rps 1k gene, and (vi) 4-year rotation of cultivars with Rps 1k, 1c, 3a, and 1k genes in year one, two, three, and four, respectively. The diversity of 121 isolates of P. sojae that were recovered by baiting from soil samples collected from the experiment were assessed using pathotyping and eight microsatellite markers, and compared with PR1. Changes in pathotype and multilocus genotypes (MLGs) were recorded at the second sampling date, indicating that P. sojae has the ability to evolve quickly. In total, 14 pathotypes and 21 MLGs were recovered over the 4-year experiment, and only 49 and 22% of the isolates had the same pathotype and MLG, respectively, as PR1. The number of isolates of P. sojae recovered varied among rotations, with more isolates recovered from rotations that included a cultivar with partial resistance. Gain of virulence was detected on Rps 1a, 1b, 1c, 1d, and 3a and was not dependent on rotation. Using simple-sequence repeat analysis, 10 alleles that were different from those of PR1 were detected throughout the 4-year period. Cultivar rotation affected the genetic structure of the P. sojae population. Recovery of isolates with different MLGs, genotypic diversity (G = 4.7), and gene diversity (UHe = 0.45) were greater under continuous rotation with partial resistance. Phytophthora root and stem rot causes economic losses in the north-central region of the United States annually. An improved understanding of the effect of Rps gene deployment on P. sojae diversity would lead to improved management practices and reduced losses.

Genome physical mapping of polyploids: a BIBAC physical map of cultivated tetraploid cotton, Gossypium hirsutum L

Polyploids account for approximately 70% of flowering plants, including many field, horticulture and forage crops. Cottons are a world-leading fiber and important oilseed crop, and a model species for study of plant polyploidization, cellulose biosynthesis and cell wall biogenesis. This study has addressed the concerns of physical mapping of polyploids with BACs and/or BIBACs by constructing a physical map of the tetraploid cotton, Gossypium hirsutum L. The physical map consists of 3,450 BIBAC contigs with an N50 contig size of 863 kb, collectively spanning 2,244 Mb. We sorted the map contigs according to their origin of subgenome, showing that we assembled physical maps for the A- and D-subgenomes of the tetraploid cotton, separately. We also identified the BIBACs in the map minimal tilling path, which consists of 15,277 clones. Moreover, we have marked the physical map with nearly 10,000 BIBAC ends (BESs), making one BES in approximately 250 kb. This physical map provides a line of evidence and a strategy for physical mapping of polyploids, and a platform for advanced research of the tetraploid cotton genome, particularly fine mapping and cloning the cotton agronomic genes and QTLs, and sequencing and assembling the cotton genome using the modern next-generation sequencing technology.

Construction of BIBAC and BAC libraries from a variety of organisms for advanced genomics research

Large-insert BAC (bacterial artificial chromosome) and BIBAC (binary BAC) libraries are essential for modern genomics research for all organisms. We helped pioneer the BAC and BIBAC technologies, and by using them we have constructed hundreds of BAC and BIBAC libraries for different species of plants, animals, marine animals, insects, algae and microbes. These libraries have been used globally for different aspects of genomics research. Here we describe the procedure with the latest improvements that we have made and used for construction of BIBAC libraries. The procedure includes the preparation of BIBAC vectors, the preparation of clonable fragments of the desired size from the source DNA, the construction and transformation of BIBACs and, finally, the characterization and assembly of BIBAC libraries. We also specify the modifications necessary for construction of BAC libraries using the protocol. The entire protocol takes ∼7 d.

Preparation of megabase-sized DNA from a variety of organisms using the nuclei method for advanced genomics research

Megabase-sized DNA is crucial to modern genomics research of all organisms. Among the preparation methods developed, the nuclei method is the simplest and most widely used for preparing high-quality megabase-sized DNA from divergent organisms. In this method, nuclei are first isolated by physically grinding the source tissues. The nontarget cytoplast organellar genomes and metabolites are removed by centrifugation and washing, thus maximizing the utility of the method and substantially improving the digestibility and clonability of the resultant DNA. The nuclei are then embedded in an agarose matrix containing numerous pores, allowing the access of restriction enzymes while preventing the DNA from physical shearing. DNA is extracted from the nuclei, purified and subsequently manipulated in the agarose matrix. Here we describe the nuclei method that we have successfully used to prepare high-quality megabase-sized DNA from hundreds of plant, animal, fish, insect, algal and microbial species. The entire protocol takes ∼3 d.

A BAC-based physical map of Zhikong scallop (Chlamys farreri Jones et Preston)

Zhikong scallop (Chlamys farreri) is one of the most economically important aquaculture species in China. Physical maps are crucial tools for genome sequencing, gene mapping and cloning, genetic improvement and selective breeding. In this study, we have developed a genome-wide, BAC-based physical map for the species. A total of 81,408 clones from two BAC libraries of the scallop were fingerprinted using an ABI 3130xl Genetic Analyzer and a fingerprinting kit developed in our laboratory. After data processing, 63,641 (∼5.8× genome coverage) fingerprints were validated and used in the physical map assembly. A total of 3,696 contigs were assembled for the physical map. Each contig contained an average of 10.0 clones, with an average physical size of 490 kb. The combined total physical size of all contigs was 1.81 Gb, equivalent to approximately 1.5 fold of the scallop haploid genome. A total of 10,587 BAC end sequences (BESs) and 167 markers were integrated into the physical map. We evaluated the physical map by overgo hybridization, BAC-FISH (fluorescence in situ hybridization), contig BAC pool screening and source BAC library screening. The results have provided evidence of the high reliability of the contig physical map. This is the first physical map in mollusc; therefore, it provides an important platform for advanced research of genomics and genetics, and mapping of genes and QTL of economical importance, thus facilitating the genetic improvement and selective breeding of the scallop and other marine molluscs.

A major invasion of transposable elements accounts for the large size of the Blumeria graminis f.sp. tritici genome

Powdery mildew of wheat (Triticum aestivum L.) is caused by the ascomycete fungus Blumeria graminis f.sp. tritici. Genomic approaches open new ways to study the biology of this obligate biotrophic pathogen. We started the analysis of the Bg tritici genome with the low-pass sequencing of its genome using the 454 technology and the construction of the first genomic bacterial artificial chromosome (BAC) library for this fungus. High-coverage contigs were assembled with the 454 reads. They allowed the characterization of 56 transposable elements and the establishment of the Blumeria repeat database. The BAC library contains 12,288 clones with an average insert size of 115 kb, which represents a maximum of 7.5-fold genome coverage. Sequencing of the BAC ends generated 12.6 Mb of random sequence representative of the genome. Analysis of BAC-end sequences revealed a massive invasion of transposable elements accounting for at least 85% of the genome. This explains the unusually large size of this genome which we estimate to be at least 174 Mb, based on a large-scale physical map constructed through the fingerprinting of the BAC library. Our study represents a crucial step in the perspective of the determination and study of the whole Bg tritici genome sequence.

A BAC/BIBAC-based physical map of chickpea, Cicer arietinum L

Background

Chickpea (Cicer arietinum L.) is the third most important pulse crop worldwide. Despite its importance, relatively little is known about its genome. The availability of a genome-wide physical map allows rapid fine mapping of QTL, development of high-density genome maps, and sequencing of the entire genome. However, no such a physical map has been developed in chickpea.

Results

We present a genome-wide, BAC/BIBAC-based physical map of chickpea developed by fingerprint analysis. Four chickpea BAC and BIBAC libraries, two of which were constructed in this study, were used. A total of 67,584 clones were fingerprinted, and 64,211 (~11.7 ×) of the fingerprints validated and used in the physical map assembly. The physical map consists of 1,945 BAC/BIBAC contigs, with each containing an average of 28.3 clones and having an average physical length of 559 kb. The contigs collectively span approximately 1,088 Mb. By using the physical map, we identified the BAC/BIBAC contigs containing or closely linked to QTL4.1 for resistance to Didymella rabiei (RDR) and QTL8 for days to first flower (DTF), thus further verifying the physical map and confirming its utility in fine mapping and cloning of QTL.

Conclusion

The physical map represents the first genome-wide, BAC/BIBAC-based physical map of chickpea. This map, along with other genomic resources previously developed in the species and the genome sequences of related species (soybean, Medicago and Lotus), will provide a foundation necessary for many areas of advanced genomics research in chickpea and other legume species. The inclusion of transformation-ready BIBACs in the map greatly facilitates its utility in functional analysis of the legume genomes.

Identification and occurrence of the LTR-Copia-like retrotransposon, PSCR and other Copia-like elements in the genome of Phytophthora sojae

Sequence analysis of the genomic region of Phytophthora sojae close to the Avr4/6 locus specifying virulence on soybean identified a Ty1/Copia-like retrotransposon that we have named Phytophthora sojae Copia-like retrotransposon (PSCR). Twelve near-complete homologs of PSCR were found in the published P. sojae genome sequence, none of which encoded a full-length polyprotein characteristic of Copia-like retrotransposons, or appears to exhibit transcriptional activity or show evidence of recent movement, suggesting they are non-functional and unlikely to have caused pathogenic variability. However, reconstructed consensus PSCR sequence encoding a full-length polyprotein resembles a functional, ancestral retroelement within P. sojae. Homologs were also found in sequence databases of other Phytophthora species. Database searches found other families of Copia-like elements in genomes of P. sojae, P. ramorum and P. infestans that were different from members of the PSCR family and from Copia-like elements reported in other organisms. It is possible that the various families of Copia-like retroelements identified in this study represent introgressions into the genome of ancient ancestor(s) of current Phytophthora species, where they have evolved and diverged considerably during the speciation. Some Copia-like families are transcriptionally active with the potential to transpose and contribute to pathogenic variation in current populations of P. sojae.

Copy number variation and transcriptional polymorphisms of Phytophthora sojae RXLR effector genes Avr1a and Avr3a

The importance of segmental duplications and copy number variants as a source of genetic and phenotypic variation is gaining greater appreciation, in a variety of organisms. Now, we have identified the Phytophthora sojae avirulence genes Avr1a and Avr3a and demonstrate how each of these Avr genes display copy number variation in different strains of P. sojae. The Avr1a locus is a tandem array of four near-identical copies of a 5.2 kb DNA segment. Two copies encoding Avr1a are deleted in some P. sojae strains, causing changes in virulence. In other P. sojae strains, differences in transcription of Avr1a result in gain of virulence. For Avr3a, there are four copies or one copy of this gene, depending on the P. sojae strain. In P. sojae strains with multiple copies of Avr3a, this gene occurs within a 10.8 kb segmental duplication that includes four other genes. Transcriptional differences of the Avr3a gene among P. sojae strains cause changes in virulence. To determine the extent of duplication within the superfamily of secreted proteins that includes Avr1a and Avr3a, predicted RXLR effector genes from the P. sojae and the P. ramorum genomes were compared by counting trace file matches from whole genome shotgun sequences. The results indicate that multiple, near-identical copies of RXLR effector genes are prevalent in oomycete genomes. We propose that multiple copies of particular RXLR effectors may contribute to pathogen fitness. However, recognition of these effectors by plant immune systems results in selection for pathogen strains with deleted or transcriptionally silenced gene copies.

The first generation of a BAC-based physical map of Brassica rapa

Background

The genus Brassica includes the most extensively cultivated vegetable crops worldwide. Investigation of the Brassica genome presents excellent challenges to study plant genome evolution and divergence of gene function associated with polyploidy and genome hybridization. A physical map of the B. rapa genome is a fundamental tool for analysis of Brassica "A" genome structure. Integration of a physical map with an existing genetic map by linking genetic markers and BAC clones in the sequencing pipeline provides a crucial resource for the ongoing genome sequencing effort and assembly of whole genome sequences.

Results

A genome-wide physical map of the B. rapa genome was constructed by the capillary electrophoresis-based fingerprinting of 67,468 Bacterial Artificial Chromosome (BAC) clones using the five restriction enzyme SNaPshot technique. The clones were assembled into contigs by means of FPC v8.5.3. After contig validation and manual editing, the resulting contig assembly consists of 1,428 contigs and is estimated to span 717 Mb in physical length. This map provides 242 anchored contigs on 10 linkage groups to be served as seed points from which to continue bidirectional chromosome extension for genome sequencing.

Conclusion

The map reported here is the first physical map for Brassica "A" genome based on the High Information Content Fingerprinting (HICF) technique. This physical map will serve as a fundamental genomic resource for accelerating genome sequencing, assembly of BAC sequences, and comparative genomics between Brassica genomes. The current build of the B. rapa physical map is available at the B. rapa Genome Project website for the user community.

RXLR effector reservoir in two Phytophthora species is dominated by a single rapidly evolving superfamily with more than 700 members

Pathogens secrete effector molecules that facilitate the infection of their hosts. A number of effectors identified in plant pathogenic Phytophthora species possess N-terminal motifs (RXLR-dEER) required for targeting these effectors into host cells. Here, we bioinformatically identify >370 candidate effector genes in each of the genomes of P. sojae and P. ramorum. A single superfamily, termed avirulence homolog (Avh) genes, accounts for most of the effectors. The Avh proteins show extensive sequence divergence but are all related and likely evolved from a common ancestor by rapid duplication and divergence. More than half of the Avh proteins contain conserved C-terminal motifs (termed W, Y, and L) that are usually arranged as a module that can be repeated up to eight times. The Avh genes belong to the most rapidly evolving part of the genome, and they are nearly always located at synteny breakpoints. The superfamily includes all experimentally identified oomycete effector and avirulence genes, and its rapid pace of evolution is consistent with a role for Avh proteins in interaction with plant hosts.

Bacterial artificial chromosome-based physical map of Gibberella zeae (Fusarium graminearum)

Fusarium graminearum is the primary causal pathogen of Fusarium head blight of wheat and barley. To accelerate genomic analysis of F. graminearum, we developed a bacterial artificial chromosome (BAC)-based physical map and integrated it with the genome sequence and genetic map. One BAC library, developed in the HindIII restriction enzyme site, consists of 4608 clones with an insert size of approximately 107 kb and covers about 13.5 genome equivalents. The other library, developed in the BamHI restriction enzyme site, consists of 3072 clones with an insert size of approximately 95 kb and covers about 8.0 genome equivalents. We fingerprinted 2688 clones from the HindIII library and 1536 clones from the BamHI library and developed a physical map of F. graminearum consisting of 26 contigs covering 39.2 Mb. Comparison of our map with the F. graminearum genome sequence showed that the size of our physical map is equivalent to the 36.1 Mb of the genome sequence. We used 31 sequence-based genetic markers, randomly spaced throughout the genome, to integrate the physical map with the genetic map. We also end-sequenced 17 BamHI BAC clones and identified 4 clones that spanned gaps in the genome sequence. Our new integrated map is highly reliable and useful for a variety of genomics studies.

Oomycete genomics: new insights and future directions

The oomycetes form a distinct phylogenetic lineage of fungus-like eukaryotic microorganisms that are relatively closely related to photosynthetic algae such as brown algae and diatoms. Plant pathogenic species, notably those of the genus Phytophthora, are the best-studied oomycetes. The genomes of four Phytophthora and one downy mildew species were recently sequenced resulting in novel insights on the evolution and pathogenesis of oomycetes. This review highlights key findings that emerged from these studies and discusses the future challenges for oomycete research.

Phytophthora sojae: root rot pathogen of soybean and model oomycete

SUMMARY Phytophthora sojae is an oomycete pathogen of soybean, classified in the kingdom Stramenopiles. It causes 'damping off' of seedlings and root rot of older plants, with an annual cost worldwide of $1-2 billion. Owing to its economic importance, this species, along with P. infestans, has been developed as a model species for the study of oomycete plant pathogens. It is readily transformed with DNA enabling over-expression and silencing of selected genes, genetic maps have been constructed and large expressed sequence tag sequence libraries have been developed. A draft genome sequence has recently been completed. This review briefly summarizes current information about the pathogenicity, evolution, molecular biology and genomics of P. sojae.

TAXONOMY

Phytophthora sojae (Kaufman & Gerdman): superkingdom Eukaryota; kingdom Stramenopila; phylum Oomycota; class Peronosporomycetidae; order Pythiales; family Pythiaceae; genus Phytophthora.

HOST RANGE

Soybean is the only economically important host. Several species of lupins have also been reported as hosts. Disease symptoms and signs: All parts of the soybean plant are susceptible to infection by P. sojae, from germinating seedlings to mature plants. In the field, P. sojae causes damping off of soybean seedlings and a root and stem rot of established plants. Leaves can be infected in the field as a result of rain splash or by deliberate inoculation in the laboratory. Damping off can affect germinating seeds or emerged seedlings and is most severe when the spring is very wet and warm (25-30 degrees C). Established plants can become infected when the soil is wet for extended periods, especially if the soil is poorly drained. Both the cortex and the vascular tissue are colonized by P. sojae, and the infection can spread rapidly along the vascular tissues in susceptible cultivars.

USEFUL WEBSITES

http://pmgn.vbi.vt.edu, http://phytophthora.vbi.vt.edu, http://www.jgi.doe.gov/Psojae, http://www.jgi.doe.gov/Pramorum, http://www.pfgd.org, http://pamgo.vbi.vt.edu, http://soy.vbi.vt.edu, https://www.vbi.vt.edu/article/articleview/78, http://plantpath.osu.edu/faculty/dorrance.php.

Phytophthora genomics: the plant destroyers' genome decoded

The year 2004 was an exciting one for the Phytophthora research community. The United States Department of Energy Joint Genome Institute (JGI) completed the draft genome sequence of two Phytophthora species, Phytophthora sojae and Phytophthora ramorum. In August of that year over 50 people gathered at JGI in Walnut Creek, California, for an annotation jamboree and searched for the secrets and surprises that the two genomes have in petto. This culminated in a paper in Science in September of this year describing the highlights of the sequencing project and emphasizing the power of having the genome sequences of two closely related organisms. This MPMI Focus issue on Phytophthora genomics contains a number of more specialized manuscripts centered on gene annotation and genome organization, and complemented with manuscripts that rely on genomics resources.