Nagendra Kumar Singh

From QTL to variety-harnessing the benefits of QTLs for drought, flood and salt tolerance in mega rice varieties of India through a multi-institutional network.

Rice is a staple cereal of India cultivated in about 43.5Mha area but with relatively low average productivity. Abiotic factors like drought, flood and salinity affect rice production adversely in more than 50% of this area. Breeding rice varieties with inbuilt tolerance to these stresses offers an economically viable and sustainable option to improve rice productivity. Availability of high quality reference genome sequence of rice, knowledge of exact position of genes/QTLs governing tolerance to abiotic stresses and availability of DNA markers linked to these traits has opened up opportunities for breeders to transfer the favorable alleles into widely grown rice varieties through marker-assisted backcross breeding (MABB). A large multi-institutional project, From QTL to variety: marker-assisted breeding of abiotic stress tolerant rice varieties with major QTLs for drought, submergence and salt tolerance was initiated in 2010 with funding support from Department of Biotechnology, Government of India, in collaboration with International Rice Research Institute, Philippines. The main focus of this project is to improve rice productivity in the fragile ecosystems of eastern, northeastern and southern part of the country, which bear the brunt of one or the other abiotic stresses frequently. Seven consistent QTLs for grain yield under drought, namely, qDTY1.1, qDTY2.1, qDTY2.2, qDTY3.1, qDTY3.2, qDTY9.1 and qDTY12.1 are being transferred into submergence tolerant versions of three high yielding mega rice varieties, Swarna-Sub1, Samba Mahsuri-Sub1 and IR 64-Sub1. To address the problem of complete submergence due to flash floods in the major river basins, the Sub1 gene is being transferred into ten highly popular locally adapted rice varieties namely, ADT 39, ADT 46, Bahadur, HUR 105, MTU 1075, Pooja, Pratikshya, Rajendra Mahsuri, Ranjit, and Sarjoo 52. Further, to address the problem of soil salinity, Saltol, a major QTL for salt tolerance is being transferred into seven popular locally adapted rice varieties, namely, ADT 45, CR 1009, Gayatri, MTU 1010, PR 114, Pusa 44 and Sarjoo 52. Genotypic background selection is being done after BC2F2 stage using an in-house designed 50K SNP chip on a set of twenty lines for each combination, identified with phenotypic similarity in the field to the recipient parent. Near-isogenic lines with more than 90% similarity to the recipient parent are now in advanced generation field trials. These climate smart varieties are expected to improve rice productivity in the adverse ecologies and contribute to the farmers livelihood.

Improvement of Basmati rice varieties for resistance to blast and bacterial blight diseases using marker assisted backcross breeding.

Marker assisted backcross breeding was employed to incorporate the blast resistance genes, Pi2 and Pi54 and bacterial blight (BB) resistance genes xa13 and Xa21 into the genetic background of Pusa Basmati 1121 (PB1121) and Pusa Basmati 6. Foreground selection for target gene(s) was followed by arduous phenotypic and background selection which fast-tracked the recovery of recurrent parent genome (RPG) to an extent of 95.8% in one of the near-isogenic lines (NILs) namely, Pusa 1728-23-33-31-56, which also showed high degree of resemblance to recurrent parent, PB6 in phenotype. The phenotypic selection prior to background selection provided an additional opportunity for identifying the novel recombinants viz., Pusa 1884-9-12-14 and Pusa 1884-3-9-175, superior to parental lines in terms of early maturity, higher yield and improved quality parameters. There was no significant difference between the RPG recovery estimated based on SSR or SNP markers, however, the panel of SNPs markers was considered as the better choice for background selection as it provided better genome coverage and included SNPs in the genic regions. Multi-location evaluation of NILs depicted their stable and high mean performance in comparison to the respective recurrent parents. The Pi2+Pi54 carrying NILs were effective in combating a pan-India panel of Magnaporthe oryzae isolates with high level of field resistance in northern, eastern and southern parts of India. Alongside, the PB1121-NILs and PB6-NILs carrying BB resistance genes xa13+Xa21 were resistant against Xanthomonas oryzae pv. oryzae races of north-western, southern and eastern parts of the country. Three of NILs developed in this study, have been promoted to final stage of testing during the u200bKharif 2015 in the Indian National Basmati Trial.

Pathways associated with lignin biosynthesis in lignomaniac jute fibres.

We generated the bast transcriptomes of a deficient lignified phloem fibre mutant and its wild-type jute (Corchorus capsularis) using Illumina paired-end sequencing. A total of 34,163 wild-type and 29,463 mutant unigenes, with average lengths of 1442 and 1136 xa0bp, respectively, were assembled de novo,xa0~77–79xa0% of which were functionally annotated. These annotated unigenes were assigned to COG (~37–40xa0%) and GO (~22–28xa0%) classifications and mapped to 189 KEGG pathways (~19–21xa0%). We discovered 38 and 43 isoforms of 16 and 10 genes of the upstream shikimate-aromatic amino acid and downstream monolignol biosynthetic pathways, respectively, rendered their sequence similarities, confirmed the identities of 22 of these candidate gene families by phylogenetic analyses and reconstructed the pathway leading to lignin biosynthesis in jute fibres. We also identified major genes and bast-related transcription factors involved in secondary cell wall (SCW) formation. The quantitative RT-PCRs revealed that phenylalanine ammonia-lyase 1 (CcPAL1) was co-down-regulated with several genes of the upstream shikimate pathway in mutant bast tissues at an early growth stage, although its expression relapsed to the normal level at the later growth stage. However, cinnamyl alcohol dehydrogenase 7 (CcCAD7) was strongly down-regulated in mutant bast tissues irrespective of growth stages. CcCAD7 disruption at an early growth stage was accompanied by co-up-regulation of SCW-specific genes cellulose synthase A7 (CcCesA7) and fasciclin-like arabinogalactan 6 (CcFLA6), which was predicted to be involved in coordinating the S-layers’ deposition in the xylan-type jute fibres. Our results identified CAD as a promising target for developing low-lignin jute fibres using genomics-assisted molecular approaches.

A restriction-site-associated DNA (RAD) linkage map, comparative genomics and identification of QTL for histological fibre content coincident with those for retted bast fibre yield and its major components in jute (Corchorus olitorius L., Malvaceae s. l.)

We used RAD (restriction-site-associated DNA) sequencing to detect genome-wide SNPs and construct a dense linkage map using an intercross F2 population in jute (Corchorus olitorius). The linkage map comprising a total of 503 RAD markers in seven linkage groups spanned 358.5xa0cM with an average marker interval of 0.72xa0cM and covered 87.0xa0% of the genome. Genome-wide segregation distortion of the mapped loci (34.4xa0%) was non-random across the linkage map, with a directional bias mostly towards the female genotypes. Jute had maximum syntenic relationships with cocoa (47.5xa0% homology) and diploid cotton (29.2xa0% homology). However, synteny and collinearity were not conserved. Histological fibre content (FC; total number of fibre cell bundles in a stem cross section) was positively correlated with fibre yield (FY), plant height (PH), root weight (RW) and stem-base diameter (SBD). Broad-sense heritability estimates were high for all traits, with FC and FY showing maximum heritability (~93xa0%). QTL mapping based on the F2:3 phenotypes detected nine QTL across the two environments. The QTL for FC was coincident with one QTL each for FY, PH, RW and SBD on top of a single-SNP (C/T) marker at 40.2xa0cM on LG1, each accounting for ~7–11xa0% of the phenotypic variance. Two QTL linked in repulsion one each for PH and SBD, with varying degrees of overdominance, were associated with two single-SNP (C/T) markers on LG2, each accounting for ~17–18xa0% of the phenotypic variance. Few candidate genes were identified within the QTL regions. Our results would enable development of tools for marker-assisted selection in jute.

Mapping quantitative trait loci responsible for resistance to Bakanae disease in rice.

BackgroundBakanae or foot rot disease caused by Fusarium fujikuroi [teleomorph: Gibberella fujikuroi (Sawada) Ito] is emerging as a serious disease in rice. The disease causes both quantitative and qualitative losses to the grains under the field conditions. Breeding for resistance to Bakanae disease is a promising strategy to manage this emerging disease. In this study, we used a population of 168xa0F14 recombinant inbred lines (RILs) derived from two indica rice parents Pusa 1342, a highly resistant variety and Pusa Basmati 1121, a highly susceptible variety to map quantitative trait loci (QTLs) governing resistance against Bakanae disease.ResultsThe disease reaction of 168xa0F14 RILs were measured on the seedlings inoculated using Fusarium fujikuroi culture using high-throughput screening protocol under glasshouse conditions. Utilizing inclusive composite interval mapping, three QTLs governing resistance to Bakanae were identified, namely qBK1.1, qBK1.2 and qBK1.3 which accounted 4.76, 24.74 and 6.49xa0% of phenotypic variation, respectively. The major effect QTL designated qBK1.2 was mapped in 0.26xa0Mb region between RM5336 and RM10153. A total of 55 annotated genes were identified within the identified QTL region qBK1.2.ConclusionsThe novel QTLs identified in this study are useful resource for efficiently breeding rice cultivars resistantxa0to Bakanae disease. This is the first report on identification of QTLs governing resistance against Bakanae in rice using inclusive composite interval mapping strategy in a RIL population.

Draft Genome of the Wheat Rust Pathogen (Puccinia triticina) Unravels Genome-Wide Structural Variations during Evolution

Abstract Leaf rust is one of the most important diseases of wheat and is caused by Puccinia triticina, a highly variable rust pathogen prevalent worldwide. Decoding the genome of this pathogen will help in unraveling the molecular basis of its evolution and in the identification of genes responsible for its various biological functions. We generated high quality draft genome sequences (approximately 100- 106 Mb) of two races of P. triticina; the variable and virulent Race77 and the old, avirulent Race106. The genomes of races 77 and 106 had 33X and 27X coverage, respectively. We predicted 27678 and 26384 genes, with average lengths of 1,129 and 1,086 bases in races 77 and 106, respectively and found that the genomes consisted of 37.49% and 39.99% repetitive sequences. Genome wide comparative analysis revealed that Race77 differs substantially from Race106 with regard to segmental duplication (SD), repeat element, and SNP/InDel characteristics. Comparative analyses showed that Race 77 is a recent, highly variable and adapted Race compared with Race106. Further sequence analyses of 13 additional pathotypes of Race77 clearly differentiated the recent, active and virulent, from the older pathotypes. Average densities of 2.4 SNPs and 0.32 InDels per kb were obtained for all P. triticina pathotypes. Secretome analysis demonstrated that Race77 has more virulence factors than Race 106, which may be responsible for the greater degree of adaptation of this pathogen. We also found that genes under greater selection pressure were conserved in the genomes of both races, and may affect functions crucial for the higher levels of virulence factors in Race77. This study provides insights into the genome structure, genome organization, molecular basis of variation, and pathogenicity of P. triticina. The genome sequence data generated in this study have been submitted to public domain databases and will be an important resource for comparative genomics studies of the more than 4000 existing Puccinia species.

High Resolution Mapping of QTLs for Heat Tolerance in Rice Using a 5K SNP Array

BackgroundHeat stress is one of the major abiotic threats to rice production, next to drought and salinity stress. Incidence of heat stress at reproductive phase of the crop results in abnormal pollination leading to floret sterility, low seed set and poor grain quality. Identification of QTLs and causal genes for heat stress tolerance at flowering will facilitate breeding for improved heat tolerance in rice. In the present study, we used 272xa0F8 recombinant inbred lines derived from a cross between Nagina22, a well-known heat tolerant Aus cultivar and IR64, a heat sensitive popular Indica rice variety to map the QTLs for heat tolerance.ResultsTo enable precise phenotyping for heat stress tolerance, we used a controlled phenotyping facility available at ICAR-Indian Institute of Wheat and Barley Research, Karnal, India. Based on ‘days to 50% flowering’ data of the RILs, we followed staggered sowing to synchronize flowering to impose heat stress at uniform stage. Using the Illumina infinium 5K SNP array for genotyping the parents and the RILs, and stress susceptibility and stress tolerance indices (SSI and STI) of percent spikelet sterility and yield per plant (g), we identified five QTLs on chromosomes 3, 5, 9 and 12. The identified QTLs explained phenotypic variation in the range of 6.27 to 21. 29%. Of these five QTLs, two high effect QTLs, one novel (qSTIPSS9.1) and one known (qSTIY5.1/qSSIY5.2), were mapped in less than 400 Kbp genomic regions, comprising of 65 and 54 genes, respectively.ConclusionsThe present study identified two major QTLs for heat tolerance in rice in narrow physical intervals, which can be employed for crop improvement by marker assisted selection (MAS) after development of suitable scorable markers for breeding of high yielding heat tolerant rice varieties. This is the first report of a major QTL for heat tolerance on chromosome 9 of rice. Further, a known QTL for heat tolerance on chromosome 5 was narrowed down from 23xa0Mb to 331 Kbp in this study.

Marker Aided Incorporation of Saltol, a Major QTL Associated with Seedling Stage Salt Tolerance, into Oryza sativa ‘Pusa Basmati 1121’

Pusa Basmati 1121 (PB1121), an elite Basmati rice cultivar is vulnerable to salinity at seedling stage. A study was undertaken to impart seedling-stage salt tolerance into PB1121 by transferring a quantitative trait locus (QTL), Saltol, using FL478 as donor, through marker assisted backcrossing. Sequence tagged microsatellite site (STMS) marker RM 3412, tightly linked to Saltol was used for foreground selection. Background recovery was estimated using 90 genome-wide STMS markers. Systematic phenotypic selection helped in accelerated recovery of recurrent parent phenome (RPP). A set of 51 BC3F2 lines homozygous for Saltol were advanced to develop four improved near isogenic lines (NILs) of PB1121 with seedling stage salt tolerance. The background genome recovery in the NILs ranged from 93.3 to 99.4%. The improved NILs were either similar or better than the recurrent parent PB1121 for yield, grain and cooking quality and duration. Biochemical analyses revealed significant variation in shoot and root Na+ and K+ concentrations. Correlation between shoot and root Na+ concentration was stronger than that between root and shoot K+ concentration. The effect of QTL integration into the NILs was studied through expression profiling of OsHKT1;5, one of the genes present in the Saltol region. The NILs had significantly higher OsHKT1;5 expression than the recurrent parent PB1121, but lower than FL478 on salt exposure validating the successful introgression of Saltol in the NILs. This was also confirmed under agronomic evaluation, wherein the NILs showed greater salt tolerance at seedling stage. One of the NILs, Pusa1734-8-3-3 (NIL3) showed comparable yield and cooking quality to the recurrent parent PB1121, with high field level seedling stage salinity tolerance and shorter duration. This is the first report of successful introgression of Saltol into a Basmati rice cultivar.

Chloroplast Genome Sequence of Pigeonpea (Cajanus cajan (L.) Millspaugh) and Cajanus scarabaeoides (L.) Thouars: Genome Organization and Comparison with Other Legumes

Pigeonpea (Cajanus cajan (L.) Millspaugh), a diploid (2n = 22) legume crop with a genome size of 852 Mbp, serves as an important source of human dietary protein especially in South East Asian and African regions. In this study, the draft chloroplast genomes of Cajanus cajan and Cajanus scarabaeoides (L.) Thouars were generated. Cajanus scarabaeoides is an important species of the Cajanus gene pool and has also been used for developing promising CMS system by different groups. A male sterile genotype harboring the C. scarabaeoides cytoplasm was used for sequencing the plastid genome. The cp genome of C. cajan is 152,242bp long, having a quadripartite structure with LSC of 83,455 bp and SSC of 17,871 bp separated by IRs of 25,398 bp. Similarly, the cp genome of C. scarabaeoides is 152,201bp long, having a quadripartite structure in which IRs of 25,402 bp length separates 83,423 bp of LSC and 17,854 bp of SSC. The pigeonpea cp genome contains 116 unique genes, including 30 tRNA, 4 rRNA, 78 predicted protein coding genes and 5 pseudogenes. A 50 kb inversion was observed in the LSC region of pigeonpea cp genome, consistent with other legumes. Comparison of cp genome with other legumes revealed the contraction of IR boundaries due to the absence of rps19 gene in the IR region. Chloroplast SSRs were mined and a total of 280 and 292 cpSSRs were identified in C. scarabaeoides and C. cajan respectively. RNA editing was observed at 37 sites in both C. scarabaeoides and C. cajan, with maximum occurrence in the ndh genes. The pigeonpea cp genome sequence would be beneficial in providing informative molecular markers which can be utilized for genetic diversity analysis and aid in understanding the plant systematics studies among major grain legumes.

De novo transcriptome sequencing facilitates genomic resource generation in Tinospora cordifolia

Tinospora cordifolia is known for its medicinal properties owing to the presence of useful constituents such as terpenes, glycosides, steroids, alkaloids, and flavonoids belonging to secondary metabolism origin. However, there is little information available pertaining to critical genomic elements (ESTs, molecular markers) necessary for judicious exploitation of its germplasm. We employed 454 GS-FLX pyrosequencing of entire transcripts and altogether ∼25xa0K assembled transcripts or Expressed sequence tags (ESTs) were identified. As the interest in T. cordifolia is primarily due to its secondary metabolite constituents, the ESTs pertaining to terpenoids biosynthetic pathway were identified in the present study. Additionally, several ESTs were assigned to different transcription factor families. To validate our transcripts dataset, the novel EST-SSR markers were generated to assess the genetic diversity among germplasm of T. cordifolia. These EST-SSR markers were found to be polymorphic and the dendrogram based on dice similarity index revealed three distinct clustering of accessions. The present study demonstrates effectiveness in using both NEWBLER and MIRA sequence read assembler software for enriching transcript-dataset and thus enables better exploitation of EST resources for mining candidate genes and designing molecular markers.