Bunyamin Tar'an

Draft genome sequence of chickpea (Cicer arietinum) provides a resource for trait improvement

Chickpea (Cicer arietinum) is the second most widely grown legume crop after soybean, accounting for a substantial proportion of human dietary nitrogen intake and playing a crucial role in food security in developing countries. We report the ∼738-Mb draft whole genome shotgun sequence of CDC Frontier, a kabuli chickpea variety, which contains an estimated 28,269 genes. Resequencing and analysis of 90 cultivated and wild genotypes from ten countries identifies targets of both breeding-associated genetic sweeps and breeding-associated balancing selection. Candidate genes for disease resistance and agronomic traits are highlighted, including traits that distinguish the two main market classes of cultivated chickpea—desi and kabuli. These data comprise a resource for chickpea improvement through molecular breeding and provide insights into both genome diversity and domestication.

Quantitative trait loci for lodging resistance, plant height and partial resistance to mycosphaerella blight in field pea ( Pisum sativum L.)

With the development of genetic maps and the identification of the most-likely positions of quantitative trait loci (QTLs) on these maps, molecular markers for lodging resistance can be identified. Consequently, marker-assisted selection (MAS) has the potential to improve the efficiency of selection for lodging resistance in a breeding program. This study was conducted to identify genetic loci associated with lodging resistance, plant height and reaction to mycosphaerella blight in pea. A population consisting of 88 recombinant inbred lines (RILs) was developed from a cross between Carneval and MP1401. The RILs were evaluated in 11 environments across the provinces of Manitoba, Saskatchewan and Alberta, Canada in 1998, 1999 and 2000. One hundred and ninety two amplified fragment length polymorphism (AFLP) markers, 13 random amplified polymorphic DNA (RAPD) markers and one sequence tagged site (STS) marker were assigned to ten linkage groups (LGs) that covered 1,274 centi Morgans (cM) of the pea genome. Six of these LGs were aligned with the previous pea map. Two QTLs were identified for lodging resistance that collectively explained 58% of the total phenotypic variation in the mean environment. Three QTLs were identified each for plant height and resistance to mycosphaerella blight, which accounted for 65% and 36% of the total phenotypic variation, respectively, in the mean environment. These QTLs were relatively consistent across environments. The AFLP marker that was associated with the major locus for lodging resistance was converted into the sequence-characterized amplified-region (SCAR) marker. The presence or absence of the SCAR marker corresponded well with the lodging reaction of 50 commercial pea varieties.

Comprehensive transcriptome assembly of Chickpea (Cicer arietinum L.) using sanger and next generation sequencing platforms: development and applications.

A comprehensive transcriptome assembly of chickpea has been developed using 134.95 million Illumina single-end reads, 7.12 million single-end FLX/454 reads and 139,214 Sanger expressed sequence tags (ESTs) from >17 genotypes. This hybrid transcriptome assembly, referred to as Cicer arietinum Transcriptome Assembly version 2 (CaTA v2, available at http://data.comparative-legumes.org/transcriptomes/cicar/lista_cicar-201201), comprising 46,369 transcript assembly contigs (TACs) has an N50 length of 1,726 bp and a maximum contig size of 15,644 bp. Putative functions were determined for 32,869 (70.8%) of the TACs and gene ontology assignments were determined for 21,471 (46.3%). The new transcriptome assembly was compared with the previously available chickpea transcriptome assemblies as well as to the chickpea genome. Comparative analysis of CaTA v2 against transcriptomes of three legumes - Medicago, soybean and common bean, resulted in 27,771 TACs common to all three legumes indicating strong conservation of genes across legumes. CaTA v2 was also used for identification of simple sequence repeats (SSRs) and intron spanning regions (ISRs) for developing molecular markers. ISRs were identified by aligning TACs to the Medicago genome, and their putative mapping positions at chromosomal level were identified using transcript map of chickpea. Primer pairs were designed for 4,990 ISRs, each representing a single contig for which predicted positions are inferred and distributed across eight linkage groups. A subset of randomly selected ISRs representing all eight chickpea linkage groups were validated on five chickpea genotypes and showed 20% polymorphism with average polymorphic information content (PIC) of 0.27. In summary, the hybrid transcriptome assembly developed and novel markers identified can be used for a variety of applications such as gene discovery, marker-trait association, diversity analysis etc., to advance genetics research and breeding applications in chickpea and other related legumes.

Using molecular markers to pyramid genes for resistance to ascochyta blight and anthracnose in lentil (Lens culinaris Medik)

Ascochyta blight caused by the fungus Ascochyta lentis Vassilievsky and anthracnose caused by Colletotrichum truncatum [(Schwein.) Andrus & W.D. Moore] are the most destructive diseases of lentil in Canada. The diseases reduce both seed yield and seed quality. Previous studies demonstrated that two genes, ral1 and AbR1, confer resistance toA. lentis and a major gene controls the resistance to 95B36 isolate of C. truncatum. Molecular markers linked to each gene have been identified. The current study was conducted to pyramid the two genes for resistance to ascochyta blight and the gene for resistance to anthracnose into lentil breeding lines. A population (F6:7) consisting of 156 recombinant inbred lines (RILs) was developed from across between ‘CDC Robin’ and a breeding line ‘964a-46’. The RILs were screened for reaction to two isolates (A1 and 3D2) ofA. lentis and one isolate (95B36) ofC. truncatum. χ2 analysis of disease reactions demonstrated that the observed segregation ratios of resistant versus susceptible fit the two gene model for resistance to ascochyta blight and a single gene model for resistance to anthracnose. Using markers linked to ral1 (UBC 2271290), to AbR1(RB18680) and to the major gene for resistance to anthracnose (OPO61250),respectively, we confirmed that 11 RILs retained all the three resistance genes. More than 82% of the lines that had either or both RB18680 and UBC2271290markers were resistant to 3D2 isolate and had a mean disease score lower than 2.5. By contrast, 80% of the lines that had none of the RAPD markers were susceptible and had a mean disease score of 5.8. For the case of A1 isolate of A. lentis, more than 74% of the lines that carriedUBC2271290 were resistant, whereas more than 79% of the lines that do not have the marker were susceptible. The analysis of the RILs usingOPO61250 marker demonstrated that 11out of 72 resistant lines carried the marker, whereas 66 out of 84 susceptible lines had the marker present. Therefore, selecting materials with both markers for resistance to ascochyta blight and a marker for resistance to anthracnose can clearly make progress toward resistance in the population.

A chromosomal genomics approach to assess and validate the desi and kabuli draft chickpea genome assemblies

With the expansion of next-generation sequencing technology and advanced bioinformatics, there has been a rapid growth of genome sequencing projects. However, while this technology enables the rapid and cost-effective assembly of draft genomes, the quality of these assemblies usually falls short of gold standard genome assemblies produced using the more traditional BAC by BAC and Sanger sequencing approaches. Assembly validation is often performed by the physical anchoring of genetically mapped markers, but this is prone to errors and the resolution is usually low, especially towards centromeric regions where recombination is limited. New approaches are required to validate reference genome assemblies. The ability to isolate individual chromosomes combined with next-generation sequencing permits the validation of genome assemblies at the chromosome level. We demonstrate this approach by the assessment of the recently published chickpea kabuli and desi genomes. While previous genetic analysis suggests that these genomes should be very similar, a comparison of their chromosome sizes and published assemblies highlights significant differences. Our chromosomal genomics analysis highlights short defined regions that appear to have been misassembled in the kabuli genome and identifies large-scale misassembly in the draft desi genome. The integration of chromosomal genomics tools within genome sequencing projects has the potential to significantly improve the construction and validation of genome assemblies. The approach could be applied both for new genome assemblies as well as published assemblies, and complements currently applied genome assembly strategies.

Genotype-Specific Variation in the Structure of Root Fungal Communities Is Related to Chickpea Plant Productivity

ABSTRACT Increasing evidence supports the existence of variations in the association of plant roots with symbiotic fungi that can improve plant growth and inhibit pathogens. However, it is unclear whether intraspecific variations in the symbiosis exist among plant cultivars and if they can be used to improve crop productivity. In this study, we determined genotype-specific variations in the association of chickpea roots with soil fungal communities and evaluated the effect of root mycota on crop productivity. A 2-year field experiment was conducted in southwestern Saskatchewan, the central zone of the chickpea growing region of the Canadian prairie. The effects of 13 cultivars of chickpea, comprising a wide range of phenotypes and genotypes, were tested on the structure of root-associated fungal communities based on internal transcribed spacer (ITS) and 18S rRNA gene markers using 454 amplicon pyrosequencing. Chickpea cultivar significantly influenced the structure of the root fungal community. The magnitude of the effect varied with the genotypes evaluated, and effects were consistent across years. For example, the roots of CDC Corrine, CDC Cory, and CDC Anna hosted the highest fungal diversity and CDC Alma and CDC Xena the lowest. Fusarium sp. was dominant in chickpea roots but was less abundant in CDC Corrine than the other cultivars. A bioassay showed that certain of these fungal taxa, including Fusarium species, can reduce the productivity of chickpea, whereas Trichoderma harzianum can increase chickpea productivity. The large variation in the profile of chickpea root mycota, which included growth-promoting and -inhibiting species, supports the possibility of improving the productivity of chickpea by improving its root mycota in chickpea genetic improvement programs using traditional breeding techniques.

Fast track genetic improvement of ascochyta blight resistance and double podding in chickpea by marker-assisted backcrossing

Ascochyta blight (AB) caused by the fungus Ascochyta rabiei Pass. Lab. is one of the major diseases of chickpea worldwide and a constraint to production in western Canada. The use of varieties with high levels of resistance is considered the most economical solution for long-term ascochyta blight management in chickpea. QTL for resistance to ascochyta blight have been identified in chickpea. The availability of molecular markers associated with QTL for ascochyta blight resistant and double podding provides an opportunity to apply marker-assisted backcrossing to introgress the traits into adapted chickpea cultivars. In the present study, molecular markers that were linked to the QTL for ascochyta blight resistance and the double podding trait, and those unlinked to the resistance were used in foreground and background selection, respectively, in backcrosses between moderately resistant donors (CDC Frontier and CDC 425-14) and the adapted varieties (CDC Xena, CDC Leader and FLIP98-135C). The strategy included two backcrosses and selection for two QTL for ascochyta blight resistance and a locus associated with double podding. The fixation of the elite genetic background was monitored with 16–22 SSR markers to accelerate restoration of the genetic background at each backcross. By the BC2F1 generation, plants with improved ascochyta blight resistance and double podding were identified. The selected plants possessed the majority of elite parental type SSR alleles on all fragments analyzed except the segment of LG 4, LG 6 and LG 8 that possessed the target QTL. The results showed that the adapted variety could be efficiently converted into a variety with improved resistance in two backcross generations.

Variation in chickpea germplasm for tolerance to imazethapyr and imazamox herbicides.

Tolerance to the imidazolinone class of herbicides would be a desirable agronomic trait for chickpea (Cicer arietinum L.) grown in western Canada. Identification of germplasm tolerant to imidazolinones and incorporation of this tolerance into future varieties will allow an integrated weed management strategy in chickpea. The current study evaluated the variation of diverse chickpea germplasm and cultivars available in Canada for tolerance to the imidazolinone class of herbicides under greenhouse conditions. Large differences among the genotypes in response to a mixture of imazethapyr and imazamox were observed. Several accessions were identified with tolerance to a mixture of imazethapyr and imazamox. Conventional breeding for imazethapyr/imazamox tolerance in chickpea is feasible. The simple screening used in the current study allows for rapid progress towards the development of herbicide-tolerant cultivars.Key words: Chickpea, germplasm, imazethapyr, imazamox, tolerance

Improved sources of resistance to ascochyta blight in chickpea

Successful chickpea production in western Canada typically requires multiple applications of fungicides to minimize the severity of ascochyta blight (AB) caused by Ascochyta rabiei (Pass.) Lab. Although planting resistant cultivars could be economical and environmentally safer than fungicide use, varieties with a high level of resistance are not available. The objective of this research was to identify potentially useful parents for breeding programs aimed at the northern Great Plains by assessing the AB reaction of 12 desi and 12 kabuli chickpea varieties for their AB reaction on leaves, stems and pods under two fungicide regimes. The experiment was conducted at Swift Current and Shaunavon, Saskatchewan, in 2004 and 2005. Differences in AB severity on leaves, stems and pods, seed yield and 1000-seed weight occurred among varieties at all site-years tested. The variation was greater among kabuli varieties than desi varieties. Ascochyta blight severity was generally lower under the high fungicide regime. A...

Genome-Wide Analysis of the Aquaporin Gene Family in Chickpea (Cicer arietinum L.)

Aquaporins (AQPs) are essential membrane proteins that play critical role in the transport of water and many other solutes across cell membranes. In this study, a comprehensive genome-wide analysis identified 40 AQP genes in chickpea (Cicer arietinum L.). A complete overview of the chickpea AQP (CaAQP) gene family is presented, including their chromosomal locations, gene structure, phylogeny, gene duplication, conserved functional motifs, gene expression, and conserved promoter motifs. To understand AQPs evolution, a comparative analysis of chickpea AQPs with AQP orthologs from soybean, Medicago, common bean, and Arabidopsis was performed. The chickpea AQP genes were found on all of the chickpea chromosomes, except chromosome 7, with a maximum of six genes on chromosome 6, and a minimum of one gene on chromosome 5. Gene duplication analysis indicated that the expansion of chickpea AQP gene family might have been due to segmental and tandem duplications. CaAQPs were grouped into four subfamilies including 15 NOD26-like intrinsic proteins (NIPs), 13 tonoplast intrinsic proteins (TIPs), eight plasma membrane intrinsic proteins (PIPs), and four small basic intrinsic proteins (SIPs) based on sequence similarities and phylogenetic position. Gene structure analysis revealed a highly conserved exon-intron pattern within CaAQP subfamilies supporting the CaAQP family classification. Functional prediction based on conserved Ar/R selectivity filters, Frogers residues, and specificity-determining positions suggested wide differences in substrate specificity among the subfamilies of CaAQPs. Expression analysis of the AQP genes indicated that some of the genes are tissue-specific, whereas few other AQP genes showed differential expression in response to biotic and abiotic stresses. Promoter profiling of CaAQP genes for conserved cis-acting regulatory elements revealed enrichment of cis-elements involved in circadian control, light response, defense and stress responsiveness reflecting their varying pattern of gene expression and potential involvement in biotic and abiotic stress responses. The current study presents the first detailed genome-wide analysis of the AQP gene family in chickpea and provides valuable information for further functional analysis to infer the role of AQP in the adaptation of chickpea in diverse environmental conditions.