Identification of Fusarium head blight sources of resistance and associated QTLs in historical and modern Canadian spring wheat

Fusarium head blight (FHB) is one the most globally destructive fungal diseases in wheat and other small grains, causing a reduction in grain yield by 10-70%. The present study was conducted in a panel of historical and modern Canadian spring wheat (Triticum aestivum L.) varieties and lines to identify new sources of FHB resistance and map associated quantitative trait loci (QTLs). We evaluated 249 varieties and lines for reaction to disease incidence, severity, and visual rating index (VRI) in seven environments by artificially spraying a mixture of four Fusarium graminearum isolates. A subset of 198 them were genotyped with the Wheat 90K iSelect single nucleotide polymorphisms (SNPs) array. Genome-wide association mapping performed on the overall best linear unbiased estimators (BLUE) computed from all seven environments and the International Wheat Genome Sequencing Consortium (IWGSC) RefSeq v2.0 physical map of 26,449 polymorphic SNPs out of the 90K identified sixteen FHB resistance QTLs that individually accounted for 5.7-10.2% of the phenotypic variance. The positions of two of the FHB resistance QTLs overlapped with plant height and flowering time QTLs. Four of the QTLs (QFhb.dms-3B.1, QFhb.dms-5A.5, QFhb.dms-5A.7, and QFhb.dms-6A.4) were simultaneously associated with disease incidence, severity, and VRI, which accounted for 27.0-33.2% of the total phenotypic variance in the combined environments. Three of the QTLs (QFhb.dms-2A.2, QFhb.dms-2D.2, and QFhb.dms-5B.8) were associated with both incidence and VRI and accounted for 20.5-22.1% of the total phenotypic variance. In comparison with the VRI of the checks, we identified four highly resistant and thirty-three moderately resistant lines and varieties. The new FHB sources of resistance and the physical map of the associated QTLs would provide wheat breeders valuable information towards their efforts in developing improved varieties in western Canada.

A physical map of traits of agronomic importance based on potato and tomato genome sequences

Potato, tomato, pepper, and eggplant are worldwide important crop and vegetable species of the Solanaceae family. Molecular linkage maps of these plants have been constructed and used to map qualitative and quantitative traits of agronomic importance. This research has been undertaken with the vision to identify the molecular basis of agronomic characters on the one hand, and on the other hand, to assist the selection of improved varieties in breeding programs by providing DNA-based markers that are diagnostic for specific agronomic characters. Since 2011, whole genome sequences of tomato and potato became available in public databases. They were used to combine the results of several hundred mapping and map-based cloning studies of phenotypic characters between 1988 and 2022 in physical maps of the twelve tomato and potato chromosomes. The traits evaluated were qualitative and quantitative resistance to pathogenic oomycetes, fungi, bacteria, viruses, nematodes, and insects. Furthermore, quantitative trait loci for yield and sugar content of tomato fruits and potato tubers and maturity or earliness were physically mapped. Cloned genes for pathogen resistance, a few genes underlying quantitative trait loci for yield, sugar content, and maturity, and several hundred candidate genes for these traits were included in the physical maps. The comparison between the physical chromosome maps revealed, in addition to known intrachromosomal inversions, several additional inversions and translocations between the otherwise highly collinear tomato and potato genomes. The integration of the positional information from independent mapping studies revealed the colocalization of qualitative and quantitative loci for resistance to different types of pathogens, called resistance hotspots, suggesting a similar molecular basis. Synteny between potato and tomato with respect to genomic positions of quantitative trait loci was frequently observed, indicating eventual similarity between the underlying genes.

Five Species of Taxus Karyotype Based on Oligo-FISH for 5S rDNA and (AG3T3)3

As a relict plant, Taxus is used in a variety of medicinal ingredients, for instance to treat a variety of cancers. Taxus plants are difficult to distinguish from one another due to their similar morphology; indeed, some species of Taxus cytogenetic data still are unclear. Oligo-FISH can rapidly and efficiently provide insight into the genetic composition and karyotype. This is important for understanding the organization and evolution of chromosomes in Taxus species. We analysed five Taxus species using two oligonucleotide probes. (AG₃T₃)₃ signals were distributed at the chromosome ends and the centromere of five species of Taxus. The 5S rDNA signal was displayed on two chromosomes of five species of Taxus. In addition to Taxus wallichiana var. mairei, 5S rDNA signals were found proximal in the remaining four species, which signals a difference in its location. The karyotype formula of Taxus wallichiana was 2n = 2x = 24m, its karyotype asymmetry index was 55.56%, and its arm ratio was 3.0087. Taxus × media's karyotype formula was 2n = 2x = 24m, its karyotype asymmetry index was 55.09%, and its arm ratio was 3.4198. The karyotype formula of Taxus yunnanensis was 2n = 2x = 24m, its karyotype asymmetry index was 55.56%, and its arm ratio was 2.6402. The karyotype formula of Taxus cuspidate was 2n = 2x = 24m, its karyotype asymmetry index was 54.67%, its arm ratio was 3.0135, and two chromosomes exhibited the 5S rDNA signal. The karyotype formula of T. wallichiana var. mairei was 2n= 2x = 22m + 2sm, its karyotype asymmetry index was 54.33%, and its arm ratio was 2.8716. Our results provide the karyotype analysis and physical genetic map of five species of Taxus, which contributes to providing molecular cytogenetics data for Taxus.

Earliness

In previous studies, we reported quantitative trait loci (QTL) associated with the heading, flowering, and maturity time in four hard red spring wheat recombinant inbred line (RIL) populations but the results are scattered in population-specific genetic maps, which is challenging to exploit efficiently in breeding. Here, we mapped and characterized QTL associated with these three earliness traits using the International Wheat Genome Sequencing Consortium (IWGSC) RefSeq v2.0 physical map. Our data consisted of (i) 6526 single nucleotide polymorphisms (SNPs) and two traits evaluated at five conventionally managed environments in the 'Cutler' × 'AC Barrie' population; (ii) 3158 SNPs and two traits evaluated across three organic and seven conventional managements in the 'Attila' × 'CDC Go' population; (iii) 5731 SilicoDArT and SNP markers and the three traits evaluated at four conventional and organic management systems in the 'Peace' × 'Carberry' population; and (iv) 1058 SNPs and two traits evaluated across two conventionally and organically managed environments in the 'Peace' × 'CDC Stanley' population. Using composite interval mapping, the phenotypic data across all environments, and the IWGSC RefSeq v2.0 physical maps, we identified a total of 44 QTL associated with days to heading (11), flowering (10), and maturity (23). Fifteen of the 44 QTL were common to both conventional and organic management systems, and the remaining QTL were specific to either the conventional (21) or organic (8) management systems. Some QTL harbor known genes, including the Vrn-A1, Vrn-B1, Rht-A1, and Rht-B1 that regulate photoperiodism, flowering time, and plant height in wheat, which lays a solid basis for cloning and further characterization.

A Comprehensive Integrated Genetic Map of the Complete Karyotype of Solea senegalensis (Kaup 1858)

Solea senegalensis aquaculture production has experienced a great increase in the last decade and, consequently, the genome knowledge of the species is gaining attention. In this sense, obtaining a high-density genome mapping of the species could offer clues to the aquaculture improvement in those aspects not resolved so far. In the present article, a review and new processed data have allowed to obtain a high-density BAC-based cytogenetic map of S. senegalensis beside the analysis of the sequences of such BAC clones to achieve integrative data. A total of 93 BAC clones were used to localize the chromosome complement of the species and 588 genes were annotated, thus almost reaching the 2.5% of the S. senegalensis genome sequences. As a result, important data about its genome organization and evolution were obtained, such as the lesser gene density of the large metacentric pair compared with the other metacentric chromosomes, which supports the theory of a sex proto-chromosome pair. In addition, chromosomes with a high number of linked genes that are conserved, even in distant species, were detected. This kind of result widens the knowledge of this species’ chromosome dynamics and evolution.

Ding Chromosomes, Reproduces its Metaphase Dynamics and Distinguishes Its Chromosomes

Chimonanthus campanulatus R.H. Chang & C.S. Ding is a good horticultural tree because of its beautiful yellow flowers and evergreen leaves. In this study, fluorescence in situ hybridization (FISH) was used to analyse mitotic metaphase chromosomes of Ch. campanulatus with 5S rDNA and (AG₃T₃)₃ oligonucleotides. Twenty-two small chromosomes were observed. Weak 5S rDNA signals were observed only in proximal regions of two chromosomes, which were adjacent to the (AG₃T₃)₃ proximal signals. Weak (AG₃T₃)₃ signals were observed on both chromosome ends, which enabled accurate chromosome counts. A pair of satellite bodies was observed. (AG₃T₃)₃ signals displayed quite high diversity, changing in intensity from weak to very strong as follows: far away from the chromosome ends (satellites), ends, subtelomeric regions, and proximal regions. Ten high-quality spreads revealed metaphase dynamics from the beginning to the end and the transition to anaphase. Chromosomes gradually grew larger and thicker into linked chromatids, which grew more significantly in width than in length. Based on the combination of 5S rDNA and (AG₃T₃)₃ signal patterns, ten chromosomes were exclusively distinguished, and the remaining twelve chromosomes were divided into two distinct groups. Our physical map, which can reproduce dynamic metaphase progression and distinguish chromosomes, will powerfully guide cytogenetic research on Chimonanthus and other trees.

Chromosome-scale pseudomolecules refined by optical, physical and genetic maps in flax

Genomes of varying sizes have been sequenced with next-generation sequencing platforms. However, most reference sequences include draft unordered scaffolds containing chimeras caused by mis-scaffolding. A BioNano genome (BNG) optical map was constructed to improve the previously sequenced flax genome (Linum usitatissimum L., 2n = 30, about 373 Mb), which consisted of 3852 scaffolds larger than 1 kb and totalling 300.6 Mb. The high-resolution BNG map of cv. CDC Bethune totalled 317 Mb and consisted of 251 BNG contigs with an N50 of 2.15 Mb. A total of 622 scaffolds (286.6 Mb, 94.9%) aligned to 211 BNG contigs (298.6 Mb, 94.2%). Of those, 99 scaffolds, diagnosed to contain assembly errors, were refined into 225 new scaffolds. Using the newly refined scaffold sequences and the validated bacterial artificial chromosome-based physical map of CDC Bethune, the 211 BNG contigs were scaffolded into 94 super-BNG contigs (N50 of 6.64 Mb) that were further assigned to the 15 flax chromosomes using the genetic map. The pseudomolecules total about 316 Mb, with individual chromosomes of 15.6 to 29.4 Mb, and cover 97% of the annotated genes. Evidence from the chromosome-scale pseudomolecules suggests that flax has undergone palaeopolyploidization and mesopolyploidization events, followed by rearrangements and deletions or fusion of chromosome arms from an ancient progenitor with a haploid chromosome number of eight.

"Targeted Sequencing by Gene Synteny," a New Strategy for Polyploid Species: Sequencing and Physical Structure of a Complex Sugarcane Region

Sugarcane exhibits a complex genome mainly due to its aneuploid nature and high ploidy level, and sequencing of its genome poses a great challenge. Closely related species with well-assembled and annotated genomes can be used to help assemble complex genomes. Here, a stable quantitative trait locus (QTL) related to sugar accumulation in sorghum was successfully transferred to the sugarcane genome. Gene sequences related to this QTL were identified in silico from sugarcane transcriptome data, and molecular markers based on these sequences were developed to select bacterial artificial chromosome (BAC) clones from the sugarcane variety SP80-3280. Sixty-eight BAC clones containing at least two gene sequences associated with the sorghum QTL were sequenced using Pacific Biosciences (PacBio) technology. Twenty BAC sequences were found to be related to the syntenic region, of which nine were sufficient to represent this region. The strategy we propose is called "targeted sequencing by gene synteny," which is a simpler approach to understanding the genome structure of complex genomic regions associated with traits of interest.

Fine Mapping of a Novel Heading Date Gene, TaHdm605, in Hexaploid Wheat

The heading date is critical in determining the adaptability of plants to specific natural environments. Molecular characterization of the wheat genes that regulate heading not only enhances our understanding of the mechanisms underlying wheat heading regulation but also benefits wheat breeding programs by improving heading phenotypes. In this study, we characterized a late heading date mutant, m605, obtained by ethyl methanesulfonate (EMS) mutation. Compared with its wild-type parent, YZ4110, m605 was at least 7 days late in heading when sown in autumn. This late heading trait was controlled by a single recessive gene named TaHdm605. Genetic mapping located the TaHdm605 locus between the molecular markers cfd152 and barc42 on chromosome 3DL using publicly available markers and then further mapped this locus to a 1.86 Mb physical genomic region containing 26 predicted genes. This fine genetic and physical mapping will be helpful for the future map-based cloning of TaHdm605 and for breeders seeking to engineer changes in the wheat heading date trait.

Physical Location of New PCR-Based Markers and Powdery Mildew Resistance Gene(s) on Rye (Secale cereale L.) Chromosome 4 Using 4R Dissection Lines

Rye (Secale cereale L.) 4R chromosome contains elite genes that are applicable for wheat (Triticum aestivum L.) cultivar improvement. PCR-based 4R-specific markers can benefit the detection of elite genes on 4R in wheat backgrounds. In this study, a new fluorescence in situ hybridization (FISH) map of the 4R^Ku chromosome of rye Kustro has been constructed. A set of 4R^Ku dissection lines was obtained and 301 new 4R^Ku-specific markers were developed using specific length amplified fragment sequencing (SLAF-seq) technology. These markers were combined with the 99 4R^Ku-specific markers previously developed, and were physically mapped to 4R^Ku chromosome using the new FISH map and the 4R^Ku dissection lines. A total of 338 of the 400 markers have been successfully mapped to six regions of 4R^Ku chromosome. Additionally, the powdery mildew resistance gene(s) on the 4RL^Ku arm was located to the segment between L.4 and L.8, the same region where 115 4RL^Ku-specific markers were mapped. The markers developed in this study can be used to identify a specific segment of 4R chromatin in wheat backgrounds, help construct a high-density physical map of 4R chromosome, and facilitate the utilization of elite genes on 4R chromosome in wheat breeding programs.

Genome Mapping in Plant Comparative Genomics

Genome mapping produces fingerprints of DNA sequences to construct a physical map of the whole genome. It provides contiguous, long-range information that complements and, in some cases, replaces sequencing data. Recent advances in genome-mapping technology will better allow researchers to detect large (>1kbp) structural variations between plant genomes. Some molecular and informatics complications need to be overcome for this novel technology to achieve its full utility. This technology will be useful for understanding phenotype responses due to DNA rearrangements and will yield insights into genome evolution, particularly in polyploids. In this review, we outline recent advances in genome-mapping technology, including the processes required for data collection and analysis, and applications in plant comparative genomics.

Cytological and molecular analysis of wheat - Agropyron cristatum translocation lines with 6P chromosome fragments conferring superior agronomic traits in common wheat

Agropyron cristatum (2n = 4x = 28, PPPP) is a wild relative of common wheat and confers several desirable agronomic traits to wheat, such as high grain number per spike and enhanced resistance to certain diseases. Development of wheat - A. cristatum 6P translocation lines facilitates its utilization in wheat improvement. In this study, 26 wheat - A. cristatum 6P translocation lines were characterized by in situ hybridization (ISH) and 6P-specific sequence-tagged-site (STS) markers. These translocation lines carried different 6P chromosomal segments, which covered the whole 6P chromosome. FISH results showed that 15, 5, and 6 lines were translocated onto wheat A, B, and D genomes, respectively. Compared with the previous reports, a fine physical map of 6P chromosome was constructed, consisting of 31 chromosomal bins with 255 STS markers. Twelve translocation lines containing 6PS13∼14 chromosomal bins were highly resistant to leaf rust. Two lines showed high grain number per spike, and three lines displayed both enhanced grain number per spike and thousand-grain weight. Development of wheat - A. cristatum 6P translocation lines will not only provide novel wheat germplasm for wheat breeding but also be helpful to broaden the genetic basis of common wheat.

Multiplex sequencing of bacterial artificial chromosomes for assembling complex plant genomes

Hierarchical shotgun sequencing remains the method of choice for assembling high-quality reference sequences of complex plant genomes. The efficient exploitation of current high-throughput technologies and powerful computational facilities for large-insert clone sequencing necessitates the sequencing and assembly of a large number of clones in parallel. We developed a multiplexed pipeline for shotgun sequencing and assembling individual bacterial artificial chromosomes (BACs) using the Illumina sequencing platform. We illustrate our approach by sequencing 668 barley BACs (Hordeum vulgare L.) in a single Illumina HiSeq 2000 lane. Using a newly designed parallelized computational pipeline, we obtained sequence assemblies of individual BACs that consist, on average, of eight sequence scaffolds and represent >98% of the genomic inserts. Our BAC assemblies are clearly superior to a whole-genome shotgun assembly regarding contiguity, completeness and the representation of the gene space. Our methods may be employed to rapidly obtain high-quality assemblies of a large number of clones to assemble map-based reference sequences of plant and animal species with complex genomes by sequencing along a minimum tiling path.

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.

Cytogenetic mapping with centromeric bacterial artificial chromosomes contigs shows that this recombination-poor region comprises more than half of barley chromosome 3H

Genetic maps are based on the frequency of recombination and often show different positions of molecular markers in comparison to physical maps, particularly in the centromere that is generally poor in meiotic recombinations. To decipher the position and order of DNA sequences genetically mapped to the centromere of barley (Hordeum vulgare) chromosome 3H, fluorescence in situ hybridization with mitotic metaphase and meiotic pachytene chromosomes was performed with 70 genomic single-copy probes derived from 65 fingerprinted bacterial artificial chromosomes (BAC) contigs genetically assigned to this recombination cold spot. The total physical distribution of the centromeric 5.5 cM bin of 3H comprises 58% of the mitotic metaphase chromosome length. Mitotic and meiotic chromatin of this recombination-poor region is preferentially marked by a heterochromatin-typical histone mark (H3K9me2), while recombination enriched subterminal chromosome regions are enriched in euchromatin-typical histone marks (H3K4me2, H3K4me3, H3K27me3) suggesting that the meiotic recombination rate could be influenced by the chromatin landscape.

Closing a gap in the physical map of the ovine major histocompatibility complex

A 184 kb gap in an ovine MHC physical map was successfully closed by identification of two overlapping clones (304C7 and 222G18) from a Chinese fine wool merino sheep BAC library. The location and tiling path of the two clones were confirmed by BAC-end sequencing and PCR amplification of loci in overlapping regions. Full-length sequencing of the clones identified 13 novel ovine genes in the gap between loci Notch4 and Btnl2, and eight of them belonging to the Butyrophilin-like (Btn-like or Btnl) gene family. The scattered distribution of the Btnl gene cluster at the gap provided a clue to explain the difficulties previously experienced in closing the gap. Completed BAC contigs of the ovine MHC will facilitate sequencing of the entire ovine leukocyte antigen (OLA) region, providing detailed information for comparative studies of MHC evolution.

High-throughput physical map anchoring via BAC-pool sequencing

BACKGROUND

Physical maps created from large insert DNA libraries, typically cloned in BAC vector, are valuable resources for map-based cloning and de novo genome sequencing. The maps are most useful if contigs of overlapping DNA clones are anchored to chromosome(s), and ordered along them using molecular markers. Here we present a novel approach for anchoring physical maps, based on sequencing three-dimensional pools of BAC clones from minimum tilling path.

RESULTS

We used physical map of wheat chromosome arm 3DS to validate the method with two different DNA sequence datasets. The first comprised 567 genes ordered along the chromosome arm based on syntenic relationship of wheat with the sequenced genomes of Brachypodium, rice and sorghum. The second dataset consisted of 7,136 SNP-containing sequences, which were mapped genetically in Aegilops tauschii, the donor of the wheat D genome. Mapping of sequence reads from individual BAC pools to the first and the second datasets enabled unambiguous anchoring 447 and 311 3DS-specific sequences, respectively, or 758 in total.

CONCLUSIONS

We demonstrate the utility of the novel approach for BAC contig anchoring based on mass parallel sequencing of three-dimensional pools prepared from minimum tilling path of physical map. The existing genetic markers as well as any other DNA sequence could be mapped to BAC clones in a single in silico experiment. The approach reduces significantly the cost and time needed for anchoring and is applicable to any genomic project involving the construction of anchored physical map.

Construction of a high-density, high-resolution genetic map and its integration with BAC-based physical map in channel catfish

Construction of genetic linkage map is essential for genetic and genomic studies. Recent advances in sequencing and genotyping technologies made it possible to generate high-density and high-resolution genetic linkage maps, especially for the organisms lacking extensive genomic resources. In the present work, we constructed a high-density and high-resolution genetic map for channel catfish with three large resource families genotyped using the catfish 250K single-nucleotide polymorphism (SNP) array. A total of 54,342 SNPs were placed on the linkage map, which to our knowledge had the highest marker density among aquaculture species. The estimated genetic size was 3,505.4 cM with a resolution of 0.22 cM for sex-averaged genetic map. The sex-specific linkage maps spanned a total of 4,495.1 cM in females and 2,593.7 cM in males, presenting a ratio of 1.7 : 1 between female and male in recombination fraction. After integration with the previously established physical map, over 87% of physical map contigs were anchored to the linkage groups that covered a physical length of 867 Mb, accounting for ∼90% of the catfish genome. The integrated map provides a valuable tool for validating and improving the catfish whole-genome assembly and facilitates fine-scale QTL mapping and positional cloning of genes responsible for economically important traits.

Selection of a nucleopolyhedrovirus isolate from Helicoverpa armigera as the basis for a biological insecticide

BACKGROUND

The cotton bollworm, Helicoverpa armigera, is an insect that causes damage in a wide range of crops in Spain. Seven isolates of H. armigera single nucleopolyhedrovirus (HearSNPV) from the Iberian Peninsula were subjected to molecular and biological characterization and compared with a Chinese genotype (HearSNPV-G4).

RESULTS

The estimated sizes of the Iberian genomes varied between 116.2 and 132.4 kb, compared to 131.4 kb of the HearSNPV-G4 reference genome. Phylogenetic analysis based on the lef-8, lef-9 and polh genes revealed that the Iberian strains were more closely related to one another than to other HearSNPV isolates. Occlusion body (OB) concentration-mortality responses (LC50 values) did not differ significantly among Iberian isolates when tested against a Helicoverpa armigera colony from Oxford (UK). Despite being the fastest killing isolate, HearSNPV-SP1 was as productive as isolates with lower virulence, with an average yield of 3.1 × 10(9) OBs larva(-1) . OBs of HearSNPV-SP1 and HearSNPV-G4 were similarly pathogenic against a recently established colony from southern Spain, although HearSNPV-SP1 was faster killing than HearSNPV-G4 against a range of instars.

CONCLUSION

The insecticidal properties of HearSNPV-SP1 mean that this strain is likely to prove useful as the basis for a biological insecticide for control of Helicoverpa armigera in Spain.

Detection and correction of assembly errors of rice Nipponbare reference sequence

A complete and high-quality genome reference sequence of an organism provides a solid foundation for a wide research community and determines the outcomes of relevant genomic, genetic, molecular and evolutionary research. Rice is an important food crop and a model plant for grasses, and therefore was the first chosen crop plant for whole genome sequencing. The genome of the japonica representative rice variety, Nipponbare, was sequenced using a gold standard, map-based clone-by-clone strategy. However, although the Nipponbare reference sequence (RefSeq) has the best quality for existing crop genome sequences, it still contains many assembly errors and gaps. To improve the Nipponbare RefSeq, first a robust method is required to detect the hidden assembly errors. Through alignments between BAC-end sequences (BESs) embedded in the Nipponbare bacterial artificial chromosome (BAC) physical map and the Nipponbare RefSeq, we detected locations on the Nipponbare RefSeq that were inversely matched with BESs and could therefore be candidates for spurious inversions of assembly. We performed further analysis of five potential locations and confirmed assembly errors at those locations; four of them, two on chr4 and two on chr11 of the Nipponbare RefSeq (IRGSP build 5), were found to be caused by reverse repetitive sequences flanking the locations. Our approach is effective in detecting spurious inversions in the Nipponbare RefSeq and can be applied for improving the sequence qualities of other genomes as well.

93-11 and evaluation of the two rice reference sequence assemblies

Reference sequences are sequences that are used for public consultation, and therefore must be of high quality. Using the whole-genome shotgun/next-generation sequencing approach, many genome sequences of complex higher plants have been generated in recent years, and are generally considered reference sequences. However, none of these sequences has been experimentally evaluated at the whole-genome sequence assembly level. Rice has a relatively simple plant genome, and the genome sequences for its two sub-species obtained using different sequencing approaches were published approximately 10 years ago. This provides a unique system for a case study to evaluate the qualities and utilities of published plant genome sequences. We constructed a robust BAC physical map embedding a large number of BAC end sequences forrice variety 93-11. Through BAC end sequence alignments and tri-assembly comparisons of the 93-11 physical map and the two reference sequences, we found that the Nipponbare reference sequence generated using the clone-by-clone approach has a high quality but still contains small artifact inversions and missing sequences. In contrast, the 93-11 reference sequence generated using the whole-genome shotgun approach contains many large and varied assembly errors, such as inversions, duplications and translocations, as well as missing sequences. The 93-11 physical map provides an invaluable resource for evaluation and improvements toward completion of both Nipponbare and 93-11 reference sequences.

Integrated physical, genetic and genome map of chickpea (Cicer arietinum L.)

Physical map of chickpea was developed for the reference chickpea genotype (ICC 4958) using bacterial artificial chromosome (BAC) libraries targeting 71,094 clones (~12× coverage). High information content fingerprinting (HICF) of these clones gave high-quality fingerprinting data for 67,483 clones, and 1,174 contigs comprising 46,112 clones and 3,256 singletons were defined. In brief, 574 Mb genome size was assembled in 1,174 contigs with an average of 0.49 Mb per contig and 3,256 singletons represent 407 Mb genome. The physical map was linked with two genetic maps with the help of 245 BAC-end sequence (BES)-derived simple sequence repeat (SSR) markers. This allowed locating some of the BACs in the vicinity of some important quantitative trait loci (QTLs) for drought tolerance and reistance to Fusarium wilt and Ascochyta blight. In addition, fingerprinted contig (FPC) assembly was also integrated with the draft genome sequence of chickpea. As a result, ~965 BACs including 163 minimum tilling path (MTP) clones could be mapped on eight pseudo-molecules of chickpea forming 491 hypothetical contigs representing 54,013,992 bp (~54 Mb) of the draft genome. Comprehensive analysis of markers in abiotic and biotic stress tolerance QTL regions led to identification of 654, 306 and 23 genes in drought tolerance "QTL-hotspot" region, Ascochyta blight resistance QTL region and Fusarium wilt resistance QTL region, respectively. Integrated physical, genetic and genome map should provide a foundation for cloning and isolation of QTLs/genes for molecular dissection of traits as well as markers for molecular breeding for chickpea improvement.

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

Construction of reference chromosome-scale pseudomolecules for potato: integrating the potato genome with genetic and physical maps

The genome of potato, a major global food crop, was recently sequenced. The work presented here details the integration of the potato reference genome (DM) with a new sequence-tagged site marker−based linkage map and other physical and genetic maps of potato and the closely related species tomato. Primary anchoring of the DM genome assembly was accomplished by the use of a diploid segregating population, which was genotyped with several types of molecular genetic markers to construct a new ~936 cM linkage map comprising 2469 marker loci. In silico anchoring approaches used genetic and physical maps from the diploid potato genotype RH89-039-16 (RH) and tomato. This combined approach has allowed 951 superscaffolds to be ordered into pseudomolecules corresponding to the 12 potato chromosomes. These pseudomolecules represent 674 Mb (~93%) of the 723 Mb genome assembly and 37,482 (~96%) of the 39,031 predicted genes. The superscaffold order and orientation within the pseudomolecules are closely collinear with independently constructed high density linkage maps. Comparisons between marker distribution and physical location reveal regions of greater and lesser recombination, as well as regions exhibiting significant segregation distortion. The work presented here has led to a greatly improved ordering of the potato reference genome superscaffolds into chromosomal “pseudomolecules”.

Whole genome profiling physical map and ancestral annotation of tobacco Hicks Broadleaf

Genomics-based breeding of economically important crops such as banana, coffee, cotton, potato, tobacco and wheat is often hampered by genome size, polyploidy and high repeat content. We adapted sequence-based whole-genome profiling (WGP™) technology to obtain insight into the polyploidy of the model plant Nicotiana tabacum (tobacco). N. tabacum is assumed to originate from a hybridization event between ancestors of Nicotiana sylvestris and Nicotiana tomentosiformis approximately 200,000 years ago. This resulted in tobacco having a haploid genome size of 4500 million base pairs, approximately four times larger than the related tomato (Solanum lycopersicum) and potato (Solanum tuberosum) genomes. In this study, a physical map containing 9750 contigs of bacterial artificial chromosomes (BACs) was constructed. The mean contig size was 462 kbp, and the calculated genome coverage equaled the estimated tobacco genome size. We used a method for determination of the ancestral origin of the genome by annotation of WGP sequence tags. This assignment agreed with the ancestral annotation available from the tobacco genetic map, and may be used to investigate the evolution of homoeologous genome segments after polyploidization. The map generated is an essential scaffold for the tobacco genome. We propose the combination of WGP physical mapping technology and tag profiling of ancestral lines as a generally applicable method to elucidate the ancestral origin of genome segments of polyploid species. The physical mapping of genes and their origins will enable application of biotechnology to polyploid plants aimed at accelerating and increasing the precision of breeding for abiotic and biotic stress resistance.

A high-resolution radiation hybrid map of the river buffalo major histocompatibility complex and comparison with BoLA

The major histocompatibility complex (MHC) in mammals codes for antigen-presenting proteins. For this reason, the MHC is of great importance for immune function and animal health. Previous studies revealed this gene-dense and polymorphic region in river buffalo to be on the short arm of chromosome 2, which is homologous to cattle chromosome 23. Using cattle-derived STS markers and a river buffalo radiation hybrid (RH) panel (BBURH5000 ), we generated a high-resolution RH map of the river buffalo MHC region. The buffalo MHC RH map (cR5000 ) was aligned with the cattle MHC RH map (cR12000 ) to compare gene order. The buffalo MHC had similar organization to the cattle MHC, with class II genes distributed in two segments, class IIa and class IIb. Class IIa was closely associated with the class I and class III regions, and class IIb was a separate cluster. A total of 53 markers were distributed into two linkage groups based on a two-point LOD score threshold of ≥8. The first linkage group included 32 markers from class IIa, class I and class III. The second linkage group included 21 markers from class IIb. Bacterial artificial chromosome clones for seven loci were mapped by fluorescence in situ hybridization on metaphase chromosomes using single- and double-color hybridizations. The order of cytogenetically mapped markers in the region corroborated the physical order of markers obtained from the RH map and served as anchor points to align and orient the linkage groups.

Second-generation genetic linkage map of catfish and its integration with the BAC-based physical map

Construction of high-density genetic linkage maps is crucially important for quantitative trait loci (QTL) studies, and they are more useful when integrated with physical maps. Such integrated maps are valuable genome resources for fine mapping of QTL, comparative genomics, and accurate and efficient whole-genome assembly. Previously, we established both linkage maps and a physical map for channel catfish, Ictalurus punctatus, the dominant aquaculture species in the United States. Here we added 2030 BAC end sequence (BES)-derived microsatellites from 1481 physical map contigs, as well as markers from singleton BES, ESTs, anonymous microsatellites, and SNPs, to construct a second-generation linkage map. Average marker density across the 29 linkage groups reached 1.4 cM/marker. The increased marker density highlighted variations in recombination rates within and among catfish chromosomes. This work effectively anchored 44.8% of the catfish BAC physical map contigs, covering ∼52.8% of the genome. The genome size was estimated to be 2546 cM on the linkage map, and the calculated physical distance per centimorgan was 393 Kb. This integrated map should enable comparative studies with teleost model species as well as provide a framework for ordering and assembling whole-genome scaffolds.

Leaf morphology in Cowpea [Vigna unguiculata (L.) Walp]: QTL analysis, physical mapping and identifying a candidate gene using synteny with model legume species

BACKGROUND

Cowpea [Vigna unguiculata (L.) Walp] exhibits a considerable variation in leaf shape. Although cowpea is mostly utilized as a dry grain and animal fodder crop, cowpea leaves are also used as a high-protein pot herb in many countries of Africa.

RESULTS

Leaf morphology was studied in the cowpea RIL population, Sanzi (sub-globose leaf shape) x Vita 7 (hastate leaf shape). A QTL for leaf shape, Hls (hastate leaf shape), was identified on the Sanzi x Vita 7 genetic map spanning from 56.54 cM to 67.54 cM distance on linkage group 15. SNP marker 1_0910 was the most significant over the two experiments, accounting for 74.7% phenotypic variance (LOD 33.82) in a greenhouse experiment and 71.5% phenotypic variance (LOD 30.89) in a field experiment. The corresponding Hls locus was positioned on the cowpea consensus genetic map on linkage group 4, spanning from 25.57 to 35.96 cM. A marker-trait association of the Hls region identified SNP marker 1_0349 alleles co-segregating with either the hastate or sub-globose leaf phenotype. High co-linearity was observed for the syntenic Hls region in Medicago truncatula and Glycine max. One syntenic locus for Hls was identified on Medicago chromosome 7 while syntenic regions for Hls were identified on two soybean chromosomes, 3 and 19. In all three syntenic loci, an ortholog for the EZA1/SWINGER (AT4G02020.1) gene was observed and is the candidate gene for the Hls locus. The Hls locus was identified on the cowpea physical map via SNP markers 1_0910, 1_1013 and 1_0992 which were identified in three BAC contigs; contig926, contig821 and contig25.

CONCLUSIONS

This study has demonstrated how integrated genomic resources can be utilized for a candidate gene approach. Identification of genes which control leaf morphology may be utilized to improve the quality of cowpea leaves for vegetable and or forage markets as well as contribute to more fundamental research understanding the control of leaf shape in legumes.

DaizuBase, an integrated soybean genome database including BAC-based physical maps

Soybean [Glycine max (L) Merrill] is one of the most important leguminous crops and ranks fourth after to rice, wheat and maize in terms of world crop production. Soybean contains abundant protein and oil, which makes it a major source of nutritious food, livestock feed and industrial products. In Japan, soybean is also an important source of traditional staples such as tofu, natto, miso and soy sauce. The soybean genome was determined in 2010. With its enormous size, physical mapping and genome sequencing are the most effective approaches towards understanding the structure and function of the soybean genome. We constructed bacterial artificial chromosome (BAC) libraries from the Japanese soybean cultivar, Enrei. The end-sequences of approximately 100,000 BAC clones were analyzed and used for construction of a BAC-based physical map of the genome. BLAST analysis between Enrei BAC-end sequences and the Williams82 genome was carried out to increase the saturation of the map. This physical map will be used to characterize the genome structure of Japanese soybean cultivars, to develop methods for the isolation of agronomically important genes and to facilitate comparative soybean genome research. The current status of physical mapping of the soybean genome and construction of database are presented.

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.

High-resolution cytogenetic map for the African malaria vector Anopheles gambiae

Cytogenetic and physical maps are indispensible for precise assembly of genome sequences, functional characterization of chromosomal regions, and population genetic and taxonomic studies. We have created a new cytogenetic map for Anopheles gambiae by using a high-pressure squash technique that increases overall band clarity. To link chromosomal regions to the genome sequence, we attached genome coordinates, based on 302 markers of bacterial artificial chromosome, cDNA clones, and PCR-amplified gene fragments, to the chromosomal bands and interbands at approximately a 0.5-1 Mb interval. In addition, we placed the breakpoints of seven common polymorphic inversions on the map and described the chromosomal landmarks for the arm and inversion identification. The map's increased resolution can be used to further enhance physical mapping, improve genome assembly, and stimulate epigenomic studies of malaria vectors.

The rice blast pathosystem as a case study for the development of new tools and raw materials for genome analysis of fungal plant pathogens

Fungi have an astounding and diverse impact on this planet. While they are agents of human diseases and the cause of allergic reactions, factories for the conversion of carbon in environmental and industrially adapted systems, and potential biological weapons, their importance as plant pathogens is unparalleled. In plants alone, fungi cause tens of thousands of different diseases and are responsible for massive losses of food, fiber and forestry at an estimated annual cost of hundreds of billions of dollars. These losses are not only realized in the incomes of individual farmers and state economies, but contribute significantly to world hunger problems and issues relating to safeguarding a global food supply. Our collective understanding of how fungi, particularly plant pathogens, grow, reproduce, identify a host and cause disease is still at a formative stage. There is an equal lack of detailed knowledge about how a plant recognizes that it is being attacked and then mounts an adequate defense response. The advent of genomic technologies has given researchers an unprecedented opportunity to address these mysteries in a powerful and more holistic manner. Where the genetic revolution of only a few years ago allowed for the characterization of single genes, today's genomic technologies are facilitating the evaluation of the entire complement of genes in an organism and the discovery of the suites of genes that act during any one time or particular condition. This review will describe the recent development of tools for whole or partial genome analysis and multigenome comparisons. Th discussion focuses on the rice blast pathosystem as a case study.

The complete nucleotide sequence of the human immunoglobulin heavy chain variable region locus

The complete nucleotide sequence of the 957-kb DNA of the human immunoglobulin heavy chain variable (VH) region locus was determined and 43 novel VH segments were identified. The region contains 123 VH segments classifiable into seven different families, of which 79 are pseudogenes. Of the 44 VH segments with an open reading frame, 39 are expressed as heavy chain proteins and 1 as mRNA, while the remaining 4 are not found in immunoglobulin cDNAs. Combinatorial diversity of VH region was calculated to be approximately 6,000. Conservation of the promoter and recombination signal sequences was observed to be higher in functional VH segments than in pseudogenes. Phylogenetic analysis of 114 VH segments clearly showed clustering of the VH segments of each family. However, an independent branch in the tree contained a single VH, V4-44.1P, sharing similar levels of homology to human VH families and to those of other vertebrates. Comparison between different copies of homologous units that appear repeatedly across the locus clearly demonstrates that dynamic DNA reorganization of the locus took place at least eight times between 133 and 10 million years ago. One nonimmunoglobulin gene of unknown function was identified in the intergenic region.