Comparing the whole-genome-shotgun and map-based sequences of the rice genome

A Comprehensive Identification and Expression Analysis of VQ Motif-Containing Proteins in Sugarcane (Saccharum spontaneum L.) under Phytohormone Treatment and Cold Stress

The VQ motif-containing proteins play a vital role in various processes such as growth, resistance to biotic and abiotic stresses and development. However, there is currently no report on the VQ genes in sugarcane (Saccharum spp.). Herein, 78 VQ genes in Saccharum spontaneum were identified and classified into nine subgroups (I-IX) by comparative genomic analyses. Each subgroup had a similar structural and conservative motif. These VQ genes expanded mainly through whole-genome segmental duplication. The cis-regulatory elements (CREs) of the VQ genes were widely involved in stress responses, phytohormone responses and physiological regulation. The RNA-seq data showed that SsVQ gene expression patterns in 10 different samples, including different developmental stages, revealed distinct temporal and spatial patterns. A total of 23 SsVQ genes were expressed in all tissues, whereas 13 SsVQ genes were not expressed in any tissues. Sequence Read Archive (SRA) data showed that the majority of SsVQs responded to cold and drought stress. In addition, quantitative real-time PCR analysis showed that the SsVQs were variously expressed under salicylic acid (SA), jasmonic acid (JA), abscisic acid (ABA) and cold treatment. This study conducted a full-scale analysis of the VQ gene family in sugarcane, which could be beneficial for the functional characterization of sugarcane VQ genes and provide candidate genes for molecular resistance breeding in cultivated sugarcane in the future.

Genome-Wide Identification of the VQ Protein Gene Family of Tobacco (Nicotiana tabacum L.) and Analysis of Its Expression in Response to Phytohormones and Abiotic and Biotic Stresses

VQ motif-containing proteins (VQ proteins) are transcriptional regulators that work independently or in combination with other transcription factors (TFs) to control plant growth and development and responses to biotic and abiotic stresses. VQ proteins contain a conserved FxxhVQxhTG amino acid motif that is the main element of its interaction with WRKY TFs. We identified 59 members of the tobacco (Nicotiana tabacum L.) NtVQ gene family by in silico analysis and examined their differential expression in response to phytohormonal treatments and following exposure to biotic and abiotic stressors. NtVQ proteins clustered into eight groups based upon their amino acid sequence and presence of various conserved domains. Groups II, IV, V, VI, and VIII contained the largest proportion of NtVQ gene family members differentially expressed in response to one or more phytohormone, and NtVQ proteins with similar domain structures had similar patterns of response to different phytohormones. NtVQ genes differentially expressed in response to temperature alterations and mechanical wounding were also identified. Over half of the NtVQ genes were significantly induced in response to Ralstonia solanacearum infection. This first comprehensive characterization of the NtVQ genes in tobacco lays the foundation for further studies of the NtVQ-mediated regulatory network in plant growth, developmental, and stress-related processes.

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.

Comparative BAC-based physical mapping of Oryza sativa ssp. indica var. 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.

Dynamic intra-japonica subspecies variation and resource application

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

A physical map of the heterozygous grapevine 'Cabernet Sauvignon' allows mapping candidate genes for disease resistance

BACKGROUND

Whole-genome physical maps facilitate genome sequencing, sequence assembly, mapping of candidate genes, and the design of targeted genetic markers. An automated protocol was used to construct a Vitis vinifera 'Cabernet Sauvignon' physical map. The quality of the result was addressed with regard to the effect of high heterozygosity on the accuracy of contig assembly. Its usefulness for the genome-wide mapping of genes for disease resistance, which is an important trait for grapevine, was then assessed.

RESULTS

The physical map included 29,727 BAC clones assembled into 1,770 contigs, spanning 715,684 kbp, and corresponding to 1.5-fold the genome size. Map inflation was due to high heterozygosity, which caused either the separation of allelic BACs in two different contigs, or local mis-assembly in contigs containing BACs from the two haplotypes. Genetic markers anchored 395 contigs or 255,476 kbp to chromosomes. The fully automated assembly and anchorage procedures were validated by BAC-by-BAC blast of the end sequences against the grape genome sequence, unveiling 7.3% of chimerical contigs. The distribution across the physical map of candidate genes for non-host and host resistance, and for defence signalling pathways was then studied. NBS-LRR and RLK genes for host resistance were found in 424 contigs, 133 of them (32%) were assigned to chromosomes, on which they are mostly organised in clusters. Non-host and defence signalling genes were found in 99 contigs dispersed without a discernable pattern across the genome.

CONCLUSION

Despite some limitations that interfere with the correct assembly of heterozygous clones into contigs, the 'Cabernet Sauvignon' physical map is a useful and reliable intermediary step between a genetic map and the genome sequence. This tool was successfully exploited for a quick mapping of complex families of genes, and it strengthened previous clues of co-localisation of major NBS-LRR clusters and disease resistance loci in grapevine.

Gene duplication and paleopolyploidy in soybean and the implications for whole genome sequencing

BACKGROUND

Soybean, Glycine max (L.) Merr., is a well documented paleopolyploid. What remains relatively under characterized is the level of sequence identity in retained homeologous regions of the genome. Recently, the Department of Energy Joint Genome Institute and United States Department of Agriculture jointly announced the sequencing of the soybean genome. One of the initial concerns is to what extent sequence identity in homeologous regions would have on whole genome shotgun sequence assembly.

RESULTS

Seventeen BACs representing approximately 2.03 Mb were sequenced as representative potential homeologous regions from the soybean genome. Genetic mapping of each BAC shows that 11 of the 20 chromosomes are represented. Sequence comparisons between homeologous BACs shows that the soybean genome is a mosaic of retained paleopolyploid regions. Some regions appear to be highly conserved while other regions have diverged significantly. Large-scale "batch" reassembly of all 17 BACs combined showed that even the most homeologous BACs with upwards of 95% sequence identity resolve into their respective homeologous sequences. Potential assembly errors were generated by tandemly duplicated pentatricopeptide repeat containing genes and long simple sequence repeats. Analysis of a whole-genome shotgun assembly of 80,000 randomly chosen JGI-DOE sequence traces reveals some new soybean-specific repeat sequences.

CONCLUSION

This analysis investigated both the structure of the paleopolyploid soybean genome and the potential effects retained homeology will have on assembling the whole genome shotgun sequence. Based upon these results, homeologous regions similar to those characterized here will not cause major assembly issues.

Gene duplication and paleopolyploidy in soybean and the implications for whole genome sequencing

BACKGROUND

Soybean, Glycine max (L.) Merr., is a well documented paleopolyploid. What remains relatively under characterized is the level of sequence identity in retained homeologous regions of the genome. Recently, the Department of Energy Joint Genome Institute and United States Department of Agriculture jointly announced the sequencing of the soybean genome. One of the initial concerns is to what extent sequence identity in homeologous regions would have on whole genome shotgun sequence assembly.

RESULTS

Seventeen BACs representing approximately 2.03 Mb were sequenced as representative potential homeologous regions from the soybean genome. Genetic mapping of each BAC shows that 11 of the 20 chromosomes are represented. Sequence comparisons between homeologous BACs shows that the soybean genome is a mosaic of retained paleopolyploid regions. Some regions appear to be highly conserved while other regions have diverged significantly. Large-scale "batch" reassembly of all 17 BACs combined showed that even the most homeologous BACs with upwards of 95% sequence identity resolve into their respective homeologous sequences. Potential assembly errors were generated by tandemly duplicated pentatricopeptide repeat containing genes and long simple sequence repeats. Analysis of a whole-genome shotgun assembly of 80,000 randomly chosen JGI-DOE sequence traces reveals some new soybean-specific repeat sequences.

CONCLUSION

This analysis investigated both the structure of the paleopolyploid soybean genome and the potential effects retained homeology will have on assembling the whole genome shotgun sequence. Based upon these results, homeologous regions similar to those characterized here will not cause major assembly issues.

A genome-wide microsatellite polymorphism database for the indica and japonica rice

Microsatellite (MS) polymorphism is an important source of genetic diversity, providing support for map-based cloning and molecular breeding. We have developed a new database that contains 52 845 polymorphic MS loci between indica and japonica, composed of ample Class II MS markers, and integrated 18 828 MS loci from IRGSP and genetic markers from RGP. Based on genetic marker positions on the rice genome (http://rise.genomics.org.cn/rice2/index.jsp ), we determined the approximate genetic distances of these MS loci and validated 100 randomly selected markers experimentally with 90% success rate. In addition, we recorded polymorphic MS positions in indica cv. 9311 that is the most important paternal parent of the two-line hybrid rice in China. Our database will undoubtedly facilitate the application of MS markers in genetic researches and marker-assisted breeding. The data set is freely available from www.wigs.zju.edu.cn/achievment/polySSR.

A molecular isolation mechanism associated with high intra-specific diversity in rice

High levels of inter-specific diversity are expected due to genetic isolation, the reproductive or geographical barriers, which lead to the accumulation of nucleotide variation. However, high levels of genetic variation are repeatedly observed even within species, notably at loci of the human major histocompatability complex and of plant resistance genes. Are molecular isolations responsible for the high intra-specific variation? To address this issue, we performed a genome-wide survey of the relationship between the possible factors that could cause genetic isolation, and the level of polymorphism, based on two rice genome comparisons. Here, we show that the levels of polymorphism in rice genes are positively correlated with the proportions of non-alignable flanking sequences, and that the correlation is observed even in single-copy genes. The physical locations of the genes were also investigated, and a strong association between the asymmetric architecture of genomes and the levels of polymorphism was revealed. These results suggest that the flank heterogeneity and the asymmetric architecture between genomes serve as isolation mechanisms at the molecular level that result in accumulation of higher genetic variation. This mechanism is of fundamental importance to understand natural genetic variation within species.