Genomic and transcriptomic comparison of nucleotide variations for insights into bruchid resistance of mungbean (Vigna radiata [L.] R. Wilczek)

BackgroundMungbean (Vigna radiata [L.] R. Wilczek) is an important legume crop with high nutritional value in South and Southeast Asia. The crop plant is susceptible to a storage pest caused by bruchids (Callosobruchus spp.). Some wild and cultivated mungbean accessions show resistance to bruchids. Genomic and transcriptomic comparison of bruchid-resistant and -susceptible mungbean could reveal bruchid-resistant genes (Br) for this pest and give insights into the bruchid resistance of mungbean.ResultsFlow cytometry showed that the genome size varied by 61 Mb (mega base pairs) among the tested mungbean accessions. Next generation sequencing followed by de novo assembly of the genome of the bruchid-resistant recombinant inbred line 59 (RIL59) revealed more than 42,000 genes. Transcriptomic comparison of bruchid-resistant and -susceptible parental lines and their offspring identified 91 differentially expressed genes (DEGs) classified into 17 major and 74 minor bruchid-resistance–associated genes. We found 408 nucleotide variations (NVs) between bruchid-resistant and -susceptible lines in regions spanning 2 kb (kilo base pairs) of the promoters of 68 DEGs. Furthermore, 282 NVs were identified on exons of 148 sequence-changed-protein genes (SCPs). DEGs and SCPs comprised genes involved in resistant-related, transposable elements (TEs) and conserved metabolic pathways. A large number of these genes were mapped to a region on chromosome 5. Molecular markers designed for variants of putative bruchid-resistance–associated genes were highly diagnostic for the bruchid-resistant genotype.ConclusionsIn addition to identifying bruchid-resistance-associated genes, we found that conserved metabolism and TEs may be modifier factors for bruchid resistance of mungbean. The genome sequence of a bruchid-resistant inbred line, candidate genes and sequence variations in promoter regions and exons putatively conditioning resistance as well as markers detecting these variants could be used for development of bruchid-resistant mungbean varieties.

[1]  R. Horres,et al.  Genome size variation among accessions of Arabidopsis thaliana. , 2004, Annals of botany.

[2]  N. S. Talekar,et al.  RFLP mapping of a major bruchid resistance gene in mungbean (Vigna radiata, L. Wilczek) , 1992, Theoretical and Applied Genetics.

[3]  J. Doležel,et al.  Estimation of nuclear DNA content in plants using flow cytometry , 2007, Nature Protocols.

[4]  M. F. Grossi-de-Sá,et al.  Do legume storage proteins play a role in defending seeds against bruchids? , 2000, Plant physiology.

[5]  P. Perrino,et al.  High antinutrient levels and bruchid resistance of cowpea (Vigna unguiculata) seeds , 2004, Euphytica.

[6]  S. Jackson,et al.  Doubling genome size without polyploidization: dynamics of retrotransposition-driven genomic expansions in Oryza australiensis, a wild relative of rice. , 2006, Genome research.

[7]  啓介 喜多村,et al.  Inheritance of resistance to infestation with azuki bean weevil in Vigna sublobata and successful incorporation to V. radiata. , 1988 .

[8]  R. Flavell.  The Molecular Characterization and Organization of Plant Chromosomal DNA Sequences , 1980 .

[9]  C. Robin Buell,et al.  The TIGR Plant Repeat Databases: a collective resource for the identification of repetitive sequences in plants , 2004, Nucleic Acids Res..

[10]  P. Shewry,et al.  The resistance of cowpea seeds to bruchid beetles is not related to levels of cysteine proteinase inhibitors , 1993, Plant Molecular Biology.

[11]  Rajeev K. Varshney,et al.  Genome sequence of mungbean and insights into evolution within Vigna species , 2014, Nature Communications.

[12]  L Nardi,et al.  Plant Genome Size Estimation by Flow Cytometry: Inter-laboratory Comparison , 1998 .

[13]  E. Kondorosi,et al.  Plant cell-size control: growing by ploidy? , 2000, Current opinion in plant biology.

[14]  Triparna Lahari,et al.  Molecular markers linked with bruchid resistance in Vigna radiata var. Sublobata and their validation , 2011, Journal of Plant Biochemistry and Biotechnology.

[15]  Thomas D. Schmittgen,et al.  Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. , 2001, Methods.

[16]  Huei-Mei Chen,et al.  Development of a molecular marker for a bruchid (Callosobruchus chinensis L.) resistance gene in mungbean , 2007, Euphytica.

[17]  S. Komatsu,et al.  Proteome analysis of the relationship between bruchid-resistant and -susceptible mungbean genotypes , 2003, Plant Genetic Resources.

[18]  Paul Medvedev,et al.  Informed and automated k-mer size selection for genome assembly , 2013, Bioinform..

[19]  H. A. Sober,et al.  Handbook of Biochemistry: Selected Data for Molecular Biology , 1971 .

[20]  David R. Kelley,et al.  Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks , 2012, Nature Protocols.

[21]  M. Ishimoto,et al.  Inheritance of Resistance to Infestation with Azuki Bean Weevil in Vigna sublobata and Successful Incorporation to V.radiata , 1988 .

[22]  E. Meyerowitz,et al.  The DNA of Arabidopsis thaliana , 1984, Molecular and General Genetics MGG.

[23]  Congfen He,et al.  A review of phytochemistry, metabolite changes, and medicinal uses of the common food mung bean and its sprouts (Vigna radiata) , 2014, Chemistry Central Journal.

[24]  M. Ishimoto,et al.  Patterns of resistance to bean weevils (Bruchidae) in Vigna radiata-mungo-sublobata complex inform the breeding of new resistant varieties. , 1989 .

[25]  E. D. Earle,et al.  Nuclear DNA content of some important plant species , 1991, Plant Molecular Biology Reporter.

[26]  A. Meister,et al.  Endopolyploidy in seed plants is differently correlated to systematics, organ, life strategy and genome size , 2003 .

[27]  N. S. Talekar,et al.  Genetics and Breeding for Bruchid Resistance in Asiatic Vigna Species , 1990 .

[28]  J. Bennetzen,et al.  Transposable element contributions to plant gene and genome evolution , 2004, Plant Molecular Biology.

[29]  Thomas R. Gingeras,et al.  STAR: ultrafast universal RNA-seq aligner , 2013, Bioinform..

[30]  Ewan Birney,et al.  Automated generation of heuristics for biological sequence comparison , 2005, BMC Bioinformatics.

[31]  Richard Durbin,et al.  Sequence analysis Fast and accurate short read alignment with Burrows – Wheeler transform , 2009 .

[32]  W. J. Kent,et al.  BLAT--the BLAST-like alignment tool. , 2002, Genome research.

[33]  A. Pingoud,et al.  Real-time PCR-based method for the estimation of genome sizes. , 2003, Nucleic acids research.

[34]  A. Furini.  Retrotransposons and the Eternal Leaves , 2012 .

[35]  Yul-Ho Kim,et al.  Inheritance and Quantitative Trait Loci Analysis of Resistance Genes to Bruchid and Bean Bug in Mungbean (Vigna radiata L. Wilczek) , 2015 .

[36]  M. Ishimoto,et al.  Bruchid resistance of transgenic azuki bean expressing seed α‐amylase inhibitor of common bean , 1996 .

[37]  K. Bains,et al.  High-iron mungbean recipes for north India. , 2003 .

[38]  Michael Ashburner,et al.  Recurrent insertion and duplication generate networks of transposable element sequences in the Drosophila melanogaster genome , 2006, Genome Biology.

[39]  H. Emden,et al.  The importance of the Bruchidae as pests of grain legumes, their distribution and control. , 1978 .

[40]  C. Nusbaum,et al.  ALLPATHS: de novo assembly of whole-genome shotgun microreads. , 2008, Genome research.

[41]  M. Ishimoto,et al.  Inhibitory Effects of Adzuki Bean Weevil-Resistant Mungbean Seeds on Growth of the Bean Bug , 1993 .

[42]  J. Cairney,et al.  A simple and efficient method for isolating RNA from pine trees , 1993, Plant Molecular Biology Reporter.

[43]  David Haussler,et al.  Using native and syntenically mapped cDNA alignments to improve de novo gene finding , 2008, Bioinform..

[44]  Shaoli Wang,et al.  Flow cytometry and K-mer analysis estimates of the genome sizes of Bemisia tabaci B and Q (Hemiptera: Aleyrodidae) , 2015, Front. Physiol..

[45]  L. Pachter,et al.  Streaming fragment assignment for real-time analysis of sequencing experiments , 2012, Nature Methods.

[46]  G. Daglish,et al.  Efficacy of protectants against Callosobruchus phaseoli (Gyll.) and C. maculatus (F.) (Coleoptera: Bruchidae) in mungbeans , 1993 .

[47]  Katinka Weinberger,et al.  Assessment of the Nutritional Impact of Agricultural Research: The Case of Mungbean in Pakistan , 2005, Food and nutrition bulletin.

[48]  G. Holloway,et al.  Legume type and temperature effects on the toxicity of insecticide to the genus Callosobruchus (Coleoptera: Bruchidae) , 2011 .

[49]  W. Huber,et al.  which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. MAnorm: a robust model for quantitative comparison of ChIP-Seq data sets , 2011 .

[50]  P. Somta,et al.  Inheritance of seed resistance to bruchids in cultivated mungbean (Vigna radiata, L. Wilczek) , 2007, Euphytica.

[51]  M. Ishimoto,et al.  Genetic localization of a bruchid resistance gene and its relationship to insecticidal cyclopeptide alkaloids, the vignatic acids, in mungbean (Vigna radiata L. Wilczek) , 1998, Molecular and General Genetics MGG.

[52]  T. Higgins,et al.  Transgenic Pea Seeds Expressing the α-Amylase Inhibitor of the Common Bean are Resistant to Bruchid Beetles , 1994, Bio/Technology.

[53]  N. Tomooka,et al.  Development of Bruchid-Resistant Mungbean Line Using Wild Mungbean Germplasm in Thailand , 1992 .

[54]  D. Raoult,et al.  Determination of the genome size of Ehrlichia spp., using pulsed field gel electrophoresis. , 1999, FEMS microbiology letters.

[55]  J. Greilhuber,et al.  Genome size determination in peronosporales (Oomycota) by Feulgen image analysis. , 1998, Fungal genetics and biology : FG & B.

[56]  H. Sung,et al.  A novel defensin encoded by a mungbean cDNA exhibits insecticidal activity against bruchid. , 2002, Journal of agricultural and food chemistry.

[57]  Dawn H. Nagel,et al.  The B73 Maize Genome: Complexity, Diversity, and Dynamics , 2009, Science.

[58]  J. Palmer,et al.  Deoxyribonucleic acid sequence organization in the mung bean genome. , 1979, Biochemistry.