Gene expression atlas of pigeonpea and its application to gain insights into genes associated with pollen fertility implicated in seed formation

Pigeonpea (Cajanus cajan) is an important grain legume of the semi-arid tropics, mainly used for its protein rich seeds. To link the genome sequence information with agronomic traits resulting from specific developmental processes, a Cajanus cajan gene expression atlas (CcGEA) was developed using the Asha genotype. Thirty tissues/organs representing developmental stages from germination to senescence were used to generate 590.84 million paired-end RNA-Seq data. The CcGEA revealed a compendium of 28 793 genes with differential, specific, spatio-temporal and constitutive expression during various stages of development in different tissues. As an example to demonstrate the application of the CcGEA, a network of 28 flower-related genes analysed for cis-regulatory elements and splicing variants has been identified. In addition, expression analysis of these candidate genes in male sterile and male fertile genotypes suggested their critical role in normal pollen development leading to seed formation. Gene network analysis also identified two regulatory genes, a pollen-specific SF3 and a sucrose-proton symporter, that could have implications for improvement of agronomic traits such as seed production and yield. In conclusion, the CcGEA provides a valuable resource for pigeonpea to identify candidate genes involved in specific developmental processes and to understand the well-orchestrated growth and developmental process in this resilient crop.

Genomic Regions Associated with Root Traits under Drought Stress in Tropical Maize (Zea mays L.)

An association mapping panel, named as CIMMYT Asia association mapping (CAAM) panel, involving 396 diverse tropical maize lines were phenotyped for various structural and functional traits of roots under drought and well-watered conditions. The experiment was conducted during Kharif (summer-rainy) season of 2012 and 2013 in root phenotyping facility at CIMMYT-Hyderabad, India. The CAAM panel was genotyped to generate 955, 690 SNPs through GBS v2.7 using Illumina Hi-seq 2000/2500 at Institute for Genomic Diversity, Cornell University, Ithaca, NY, USA. GWAS analysis was carried out using 331,390 SNPs filtered from the entire set of SNPs revealed a total of 50 and 67 SNPs significantly associated for root functional (transpiration efficiency, flowering period water use) and structural traits (rooting depth, root dry weight, root length, root volume, root surface area and root length density), respectively. In addition to this, 37 SNPs were identified for grain yield and shoot biomass under well-watered and drought stress. Though many SNPs were found to have significant association with the traits under study, SNPs that were common for more than one trait were discussed in detail. A total 18 SNPs were found to have common association with more than one trait, out of which 12 SNPs were found within or near the various gene functional regions. In this study we attempted to identify the trait specific maize lines based on the presence of favorable alleles for the SNPs associated with multiple traits. Two SNPs S3_128533512 and S7_151238865 were associated with transpiration efficiency, shoot biomass and grain yield under well-watered condition. Based on favorable allele for these SNPs seven inbred lines were identified. Similarly, four lines were identified for transpiration efficiency and shoot biomass under drought stress based on the presence of favorable allele for the common SNPs S1_211520521, S2_20017716, S3_57210184 and S7_130878458 and three lines were identified for flowering period water-use, transpiration efficiency, root dry weight and root volume based on the presence of favorable allele for the common SNPs S3_162065732 and S3_225760139.

Two key genomic regions harbour QTLs for salinity tolerance in ICCV 2 × JG 11 derived chickpea (Cicer arietinum L.) recombinant inbred lines

BACKGROUND

Although chickpea (Cicer arietinum L.), an important food legume crop, is sensitive to salinity, considerable variation for salinity tolerance exists in the germplasm. To improve any existing cultivar, it is important to understand the genetic and physiological mechanisms underlying this tolerance.

RESULTS

In the present study, 188 recombinant inbred lines (RILs) derived from the cross ICCV 2 × JG 11 were used to assess yield and related traits in a soil with 0 mM NaCl (control) and 80 mM NaCl (salinity) over two consecutive years. Salinity significantly (P < 0.05) affected almost all traits across years and yield reduction was in large part related to a reduction in seed number but also a reduction in above ground biomass. A genetic map was constructed using 56 polymorphic markers (28 simple sequence repeats; SSRs and 28 single nucleotide polymorphisms; SNPs). The QTL analysis revealed two key genomic regions on CaLG05 (28.6 cM) and on CaLG07 (19.4 cM), that harboured QTLs for six and five different salinity tolerance associated traits, respectively, and imparting either higher plant vigour (on CaLG05) or higher reproductive success (on CaLG07). Two major QTLs for yield in the salinity treatment (explaining 12 and 17% of the phenotypic variation) were identified within the two key genomic regions. Comparison with already published chickpea genetic maps showed that these regions conferred salinity tolerance across two other populations and the markers can be deployed for enhancing salinity tolerance in chickpea. Based on the gene ontology annotation, forty eight putative candidate genes responsive to salinity stress were found on CaLG05 (31 genes) and CaLG07 (17 genes) in a distance of 11.1 Mb and 8.2 Mb on chickpea reference genome. Most of the genes were known to be involved in achieving osmoregulation under stress conditions.

CONCLUSION

Identification of putative candidate genes further strengthens the idea of using CaLG05 and CaLG07 genomic regions for marker assisted breeding (MAB). Further fine mapping of these key genomic regions may lead to novel gene identification for salinity stress tolerance in chickpea.

Selection and Validation of Housekeeping Genes as Reference for Gene Expression Studies in Pigeonpea (Cajanus cajan) under Heat and Salt Stress Conditions

To identify stable housekeeping genes as a reference for expression analysis under heat and salt stress conditions in pigeonpea, the relative expression variation for 10 commonly used housekeeping genes (EF1α, UBQ10, GAPDH, 18Sr RNA, 25Sr RNA, TUB6, ACT1, IF4α, UBC, and HSP90) was studied in root, stem, and leaves tissues of Asha (ICPL 87119), a leading pigeonpea variety. Three statistical algorithms geNorm, NormFinder, and BestKeeper were used to define the stability of candidate genes. Under heat stress, UBC, HSP90, and GAPDH were found to be the most stable reference genes. In the case of salinity stress, GAPDH followed by UBC and HSP90 were identified to be the most stable reference genes. Subsequently, the above identified genes were validated using qRT-PCR based gene expression analysis of two universal stress-resposive genes namely uspA and uspB. The relative quantification of these two genes varied according to the internal controls (most stable, least stable, and combination of most stable and least stable housekeeping genes) and thus confirmed the choice as well as validation of internal controls in such experiments. The identified and validated housekeeping genes will facilitate gene expression studies under heat and salt stress conditions in pigeonpea.

Identification and Validation of Selected Universal Stress Protein Domain Containing Drought-Responsive Genes in Pigeonpea (Cajanus cajan L.)

Pigeonpea is a resilient crop, which is relatively more drought tolerant than many other legume crops. To understand the molecular mechanisms of this unique feature of pigeonpea, 51 genes were selected using the Hidden Markov Models (HMM) those codes for proteins having close similarity to universal stress protein domain. Validation of these genes was conducted on three pigeonpea genotypes (ICPL 151, ICPL 8755, and ICPL 227) having different levels of drought tolerance. Gene expression analysis using qRT-PCR revealed 6, 8, and 18 genes to be ≥2-fold differentially expressed in ICPL 151, ICPL 8755, and ICPL 227, respectively. A total of 10 differentially expressed genes showed ≥2-fold up-regulation in the more drought tolerant genotype, which encoded four different classes of proteins. These include plant U-box protein (four genes), universal stress protein A-like protein (four genes), cation/H(+) antiporter protein (one gene) and an uncharacterized protein (one gene). Genes C.cajan_29830 and C.cajan_33874 belonging to uspA, were found significantly expressed in all the three genotypes with ≥2-fold expression variations. Expression profiling of these two genes on the four other legume crops revealed their specific role in pigeonpea. Therefore, these genes seem to be promising candidates for conferring drought tolerance specifically to pigeonpea.

Identification and Validation of Selected Universal Stress Protein Domain Containing Drought-Responsive Genes in Pigeonpea (Cajanus cajan L.)

Pigeonpea is a resilient crop, which is relatively more drought tolerant than many other legume crops. To understand the molecular mechanisms of this unique feature of pigeonpea, 51 genes were selected using the Hidden Markov Models (HMM) those codes for proteins having close similarity to universal stress protein domain. Validation of these genes was conducted on three pigeonpea genotypes (ICPL 151, ICPL 8755, and ICPL 227) having different levels of drought tolerance. Gene expression analysis using qRT-PCR revealed 6, 8, and 18 genes to be ≥2-fold differentially expressed in ICPL 151, ICPL 8755, and ICPL 227, respectively. A total of 10 differentially expressed genes showed ≥2-fold up-regulation in the more drought tolerant genotype, which encoded four different classes of proteins. These include plant U-box protein (four genes), universal stress protein A-like protein (four genes), cation/H(+) antiporter protein (one gene) and an uncharacterized protein (one gene). Genes C.cajan_29830 and C.cajan_33874 belonging to uspA, were found significantly expressed in all the three genotypes with ≥2-fold expression variations. Expression profiling of these two genes on the four other legume crops revealed their specific role in pigeonpea. Therefore, these genes seem to be promising candidates for conferring drought tolerance specifically to pigeonpea.

Advances in genetics and molecular breeding of three legume crops of semi-arid tropics using next-generation sequencing and high-throughput genotyping technologies

Molecular markers are the most powerful genomic tools to increase the efficiency and precision of breeding practices for crop improvement. Progress in the development of genomic resources in the leading legume crops of the semi-arid tropics (SAT), namely, chickpea (Cicer arietinum), pigeonpea (Cajanus cajan) and groundnut (Arachis hypogaea), as compared to other crop species like cereals, has been very slow. With the advances in next-generation sequencing (NGS) and high-throughput (HTP) genotyping methods, there is a shift in development of genomic resources including molecular markers in these crops. For instance, 2,000 to 3,000 novel simple sequence repeats (SSR) markers have been developed each for chickpea, pigeonpea and groundnut. Based on Sanger, 454/FLX and Illumina transcript reads, transcriptome assemblies have been developed for chickpea (44,845 transcript assembly contigs, or TACs) and pigeonpea (21,434 TACs). Illumina sequencing of some parental genotypes of mapping populations has resulted in the development of 120 million reads for chickpea and 128.9 million reads for pigeonpea. Alignment of these Illumina reads with respective transcriptome assemblies have provided more than 10,000 SNPs each in chickpea and pigeonpea. A variety of SNP genotyping platforms including GoldenGate, VeraCode and Competitive Allele Specific PCR (KASPar) assays have been developed in chickpea and pigeonpea. By using above resources, the first-generation or comprehensive genetic maps have been developed in the three legume speciesmentioned above. Analysis of phenotyping data together with genotyping data has provided candidate markers for drought-tolerance-related root traits in chickpea, resistance to foliar diseases in groundnut and sterility mosaic disease (SMD) and fertility restoration in pigeonpea. Together with these traitassociated markers along with those already available, molecular breeding programmes have been initiated for enhancing drought tolerance, resistance to fusarium wilt and ascochyta blight in chickpea and resistance to foliar diseases in groundnut. These trait-associated robust markers along with other genomic resources including genetic maps and genomic resources will certainly accelerate crop improvement programmes in the SAT legumes.

Quantitative trait locus analysis and construction of consensus genetic map for drought tolerance traits based on three recombinant inbred line populations in cultivated groundnut (Arachis hypogaea L.)

Groundnut (Arachis hypogaea L.) is an important food and cash crop grown mainly in semi-arid tropics (SAT) regions of the world where drought is the major constraint on productivity. With the aim of understanding the genetic basis and identification of quantitative trait loci (QTL) for drought tolerance, two new recombinant inbred line (RIL) mapping populations, namely ICGS 76 × CSMG 84-1 (RIL-2) and ICGS 44 × ICGS 76 (RIL-3), were used. After screening of 3,215 simple sequence repeat (SSR) markers on the parental genotypes of these populations, two new genetic maps were developed with 119 (RIL-2) and 82 (RIL-3) SSR loci. Together with these maps and the reference map with 191 SSR loci based on TAG 24 × ICGV 86031 (RIL-1), a consensus map was constructed with 293 SSR loci distributed over 20 linkage groups, spanning 2,840.8 cM. As all these three populations segregate for drought-tolerance-related traits, a comprehensive QTL analysis identified 153 main effect QTL (M-QTL) and 25 epistatic QTL (E-QTL) for drought-tolerance-related traits. Localization of these QTL on the consensus map provided 16 genomic regions that contained 125 QTL. A few key genomic regions were selected on the basis of the QTL identified in each region, and their expected role in drought adaptation is also discussed. Given that no major QTL for drought adaptation were identified, novel breeding approaches such as marker-assisted recurrent selection (MARS) and genomic selection (GS) approaches are likely to be the preferred approaches for introgression of a larger number of QTL in order to breed drought-tolerant groundnut genotypes. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s11032-011-9660-0) contains supplementary material, which is available to authorized users.

Phenotyping chickpeas and pigeonpeas for adaptation to drought

The chickpea and pigeonpea are protein-rich grain legumes used for human consumption in many countries. Grain yield of these crops is low to moderate in the semi-arid tropics with large variation due to high GxE interaction. In the Indian subcontinent chickpea is grown in the post-rainy winter season on receding soil moisture, and in other countries during the cool and dry post winter or spring seasons. The pigeonpea is sown during rainy season which flowers and matures in post-rainy season. The rainy months are hot and humid with diurnal temperature varying between 25 and 35°C (maximum) and 20 and 25°C (minimum) with an erratic rainfall. The available soil water during post-rainy season is about 200-250 mm which is bare minimum to meet the normal evapotranspiration. Thus occurrence of drought is frequent and at varying degrees. To enhance productivity of these crops cultivars tolerant to drought need to be developed. ICRISAT conserves a large number of accessions of chickpea (>20,000) and pigeonpea (>15,000). However only a small proportion (<1%) has been used in crop improvement programs mainly due to non-availability of reliable information on traits of economic importance. To overcome this, core and mini core collections (10% of core, 1% of entire collection) have been developed. Using the mini core approach, trait-specific donor lines were identified for agronomic, quality, and stress related traits in both crops. Composite collections were developed both in chickpea (3000 accessions) and pigeonpea (1000 accessions), genotyped using SSR markers and genotype based reference sets of 300 accessions selected for each crop. Screening methods for different drought-tolerant traits such as early maturity (drought escape), large and deep root system, high water-use efficiency, smaller leaflets, reduced canopy temperature, carbon isotope discrimination, high leaf chlorophyll content (drought avoidance), and breeding strategies for improving drought tolerance have been discussed.

Large-scale transcriptome analysis in chickpea (Cicer arietinum L.), an orphan legume crop of the semi-arid tropics of Asia and Africa

Chickpea (Cicer arietinum L.) is an important legume crop in the semi-arid regions of Asia and Africa. Gains in crop productivity have been low however, particularly because of biotic and abiotic stresses. To help enhance crop productivity using molecular breeding techniques, next generation sequencing technologies such as Roche/454 and Illumina/Solexa were used to determine the sequence of most gene transcripts and to identify drought-responsive genes and gene-based molecular markers. A total of 103,215 tentative unique sequences (TUSs) have been produced from 435,018 Roche/454 reads and 21,491 Sanger expressed sequence tags (ESTs). Putative functions were determined for 49,437 (47.8%) of the TUSs, and gene ontology assignments were determined for 20,634 (41.7%) of the TUSs. Comparison of the chickpea TUSs with the Medicago truncatula genome assembly (Mt 3.5.1 build) resulted in 42,141 aligned TUSs with putative gene structures (including 39,281 predicted intron/splice junctions). Alignment of ∼37 million Illumina/Solexa tags generated from drought-challenged root tissues of two chickpea genotypes against the TUSs identified 44,639 differentially expressed TUSs. The TUSs were also used to identify a diverse set of markers, including 728 simple sequence repeats (SSRs), 495 single nucleotide polymorphisms (SNPs), 387 conserved orthologous sequence (COS) markers, and 2088 intron-spanning region (ISR) markers. This resource will be useful for basic and applied research for genome analysis and crop improvement in chickpea.

Assessment of transpiration efficiency in peanut (Arachis hypogaea L.) under drought using a lysimetric system

Transpiration efficiency (TE) is an important trait for drought tolerance in peanut (Arachis hypogaea L.). The variation in TE was assessed gravimetrically using a long time interval in nine peanut genotypes (Chico, ICGS 44, ICGV 00350, ICGV 86015, ICGV 86031, ICGV 91114, JL 24, TAG 24 and TMV 2) grown in lysimeters under well-watered or drought conditions. Transpiration was measured by regularly weighing the lysimeters, in which the soil surface was mulched with a 2-cm layer of polythene beads. TE in the nine genotypes used varied from 1.4 to 2.9 g kg(-1) under well-watered and 1.7 to 2.9 g kg(-1) under drought conditions, showing consistent variation in TE among genotypes. A higher TE was found in ICGV 86031 in both well-watered and drought conditions and lower TE was found in TAG-24 under both water regimes. Although total water extraction differed little across genotypes, the pattern of water extraction from the soil profile varied among genotypes. High water extraction within 24 days following stress imposition was negatively related to pod yield (r(2) = 0.36), and negatively related to water extraction during a subsequent period of 32 days (r(2) = 0.73). By contrast, the latter, i.e. water extraction during a period corresponding to grain filling (24 to 56 days after flowering) was positively related to pod yield (r(2) = 0.36). TE was positively correlated with pod weight (r(2) = 0.30) under drought condition. Our data show that under an intermittent drought regime, TE and water extraction from the soil profile during a period corresponding to pod filling were the most important components.

The first SSR-based genetic linkage map for cultivated groundnut (Arachis hypogaea L.)

Molecular markers and genetic linkage maps are pre-requisites for molecular breeding in any crop species. In case of peanut or groundnut (Arachis hypogaea L.), an amphidiploid (4X) species, not a single genetic map is, however, available based on a mapping population derived from cultivated genotypes. In order to develop a genetic linkage map for tetraploid cultivated groundnut, a total of 1,145 microsatellite or simple sequence repeat (SSR) markers available in public domain as well as unpublished markers from several sources were screened on two genotypes, TAG 24 and ICGV 86031 that are parents of a recombinant inbred line mapping population. As a result, 144 (12.6%) polymorphic markers were identified and these amplified a total of 150 loci. A total of 135 SSR loci could be mapped into 22 linkage groups (LGs). While six LGs had only two SSR loci, the other LGs contained 3 (LG_AhXV) to 15 (LG_AhVIII) loci. As the mapping population used for developing the genetic map segregates for drought tolerance traits, phenotyping data obtained for transpiration, transpiration efficiency, specific leaf area and SPAD chlorophyll meter reading (SCMR) for 2 years were analyzed together with genotyping data. Although, 2-5 QTLs for each trait mentioned above were identified, the phenotypic variation explained by these QTLs was in the range of 3.5-14.1%. In addition, alignment of two linkage groups (LGs) (LG_AhIII and LG_AhVI) of the developed genetic map was shown with available genetic maps of AA diploid genome of groundnut and Lotus and Medicago. The present study reports the construction of the first genetic map for cultivated groundnut and demonstrates its utility for molecular mapping of QTLs controlling drought tolerance related traits as well as establishing relationships with diploid AA genome of groundnut and model legume genome species. Therefore, the map should be useful for the community for a variety of applications.

Pathways database system: an integrated system for biological pathways

MOTIVATION

During the next phase of the Human Genome Project, research will focus on functional studies of attributing functions to genes, their regulatory elements, and other DNA sequences. To facilitate the use of genomic information in such studies, a new modeling perspective is needed to examine and study genome sequences in the context of many kinds of biological information. Pathways are the logical format for modeling and presenting such information in a manner that is familiar to biological researchers.

RESULTS

In this paper we present an integrated system, called Pathways Database System, with a set of software tools for modeling, storing, analyzing, visualizing, and querying biological pathways data at different levels of genetic, molecular, biochemical and organismal detail. The novel features of the system include: (a) genomic information integrated with other biological data and presented from a pathway, rather than from the DNA sequence, perspective; (b) design for biologists who are possibly unfamiliar with genomics, but whose research is essential for annotating gene and genome sequences with biological functions; (c) database design, implementation and graphical tools which enable users to visualize pathways data in multiple abstraction levels, and to pose predetermined queries; and (d) an implementation that allows for web(XML)-based dissemination of query outputs (i.e. pathways data) to researchers in the community, giving them control on the use of pathways data.

AVAILABILITY

Available on request from the authors.

Efficacy of intrathecal methotrexate with and without cranial radiotherapy in preventing central nervous system relapses in acute lymphocytic leukemia

Introduction of CNS chemoprophylaxis was a major milestone in the development of current therapy for acute lymphocytic leukemia. However, controversies are still existing for ideal form of CNS chemoprophylaxis. The present study was conducted to determine the efficacy of intrathecal methotrexate (IT-MTX) with and without cranial radiotherapy in preventing CNS relapses in Indian children. CNS chemoprophylaxis comprising of six injections of intrathecal methotrexate (12 mg/M2) was administered alone or along with cranial radiotherapy (2000 GY) in 76 children each after successful induction remission. Cranial radiotherapy (RT) with intrathecal methotrexate (IT-MTX) was observed to be more effective as CNS relapses were seen in 11.8% of children as compared to 16.8% of children receiving IT-MTX alone. IT-MTX along with cranial RT delayed the occurrence of CNS relapses and prolonged the event free survival periods.

Causes of mortality in children with acute lymphocytic leukemia

Fifty five deaths between January, 1982 to September, 1989 in children with acute lymphoblastic leukemia (ALL) were evaluated to determine the cause of mortality. Fifty cases died during remission. Infection alone was responsible for death in 26 of 55 (47.3%) cases while hemorrhage was seen in 7 (12.7%) children. Infection and hemorrhage together were responsible in another 13 cases. Gastrointestinal tract and pulmonary system were the major sites of bleeding. Infections either alone or in combination with other factors were responsible for death in 42 of 55 (76.5%) of children. Septicemia (n = 11), gastrointestinal (n = 15) and pulmonary infections (n = 10) and meningitis in 2 cases were the major sites of infections. Pseudomonas and Klebsiella in 6 cases each accounted for 54.5% of isolates.

Subacute hepatic failure; is it a distinct entity?

We prospectively studied patients with subacute hepatic failure due to subacute hepatitis to find out 1) its relative prevalence compared to acute liver failure due to fulminant hepatitis and chronic liver failure due to chronic active hepatitis; 2) its clinical, biochemical, and morphological hepatitis; 3) the role of virus B in its etiology; and 4) its prognosis and whether there were any predictors of bad prognosis. Thirty-three patients with subacute hepatic failure were registered during a 3-year period. Persistent or progressively deepening jaundice of 8 weeks duration and development of moderate to sever ascites in patients starting otherwise typical features of acute viral hepatitis, defined the subacute hepatic failure group. The characteristic features included moderate to deep icterus, ascites, and peripheral edema; encephalopathy and gastrointestinal bleeding were infrequent. Liver function tests were abnormal but not diagnostic. Submassive and bridging necrosis of the liver were the main histological findings. Virus B etiology was recorded in 42% of the patients. Mortality was 66%. This condition is highly fatal and not infrequent in India.

Study of an epidemic of jaundice, presumably due to toxic hepatitis, in Northwest India

An epidemic of jaundice probably due to toxic hepatitis occurred in three adjoining districts of Northwest India during the period November and December, 1974. The dogs of the villages were affected first, then the human beings. Detailed clinical features, appropriate laboratory tests, and liver biopsies were studied. A retrospective epidemiological survey was carried out. The disease had a subacute onset starting with high fever, followed by rapidly progressive jaundice. Ascites appeared simultaneously and soon became quite massive. Hepatomegaly was recorded when ascites decreased. Liver function tests suggested cholestatic jaundice. The mortality rate in the hospital was 10%. Clinical features in dogs were similar, but mortality was almost 100%. Liver histology was characterized by (1) edema and collagenization of the central veins, never with thrombosis, (2) cholangiolar proliferation, (3) moderate to severe ballooning of the hepatocytes, (4) perisinusoidal fibrosis, (5) cholestasis, and finally, (6) cirrhosis with reverse lobulation. Etiology of this epidemic of hepatitis could not be unequivocally established. Critical analysis of the data suggests that some food toxin may have been a factor in the outbreak of this unusual epidemic of toxic hepatitis.