DNA methylation plays important roles in genome protection, regulation of gene expression and is associated with plants development. Plant DNA methylation pattern was mediated by cytosine-5 DNA methyltransferase and demethylase. Although the genomes of AA and BB wild peanuts have been fully sequenced, these two gene families have not been studied. In this study we report the identification and analysis of putative cytosine-5 DNA methyltransferases (C5-MTases) and demethylases in AA and BB wild peanuts. Cytosine-5 DNA methyltransferases in AA and BB wild peanuts could be classified in MET, CMT, and DRM2 groups based on their domain organization. This result was supported by the gene and protein structural characteristics and phylogenetic analysis. We found that some wild peanut DRM2 members didn't contain UBA domain which was different from other plants such as Arabidopsis, maize and soybean. Five DNA demethylase encoding genes were found in AA genome and five in BB genome. The selective pressure analysis showed that wild peanut C5-MTase genes mainly underwent purifying selection but many positive selection sites can be detected. Conversely, DNA demethylase genes mainly underwent positive selection during evolution. Additionally, the expression dynamic of cytosine-5 DNA methyltransferase and demethylase genes in different cultivated peanut tissues were analyzed. Expression result showed that cold, heat or PEG stress could influence the expression level of C5-MTase and DNA demethylase genes in cultivated peanut. These results are useful for better understanding the complexity of these two gene families, and will facilitate epigenetic studies in peanut in the future.

Legumes play an important role as food and forage crops in international agriculture especially in developing countries. Legumes have a unique biological process called nitrogen fixation (NF) by which they convert atmospheric nitrogen to ammonia. Although legume genomes have undergone polyploidization, duplication and divergence, NF-related genes, because of their essential functional role for legumes, might have remained conserved. To understand the relationship of divergence and evolutionary processes in legumes, this study analyzes orthologs and paralogs for selected 20 NF-related genes by using comparative genomic approaches in six legumes i.e., Medicago truncatula (Mt), Cicer arietinum, Lotus japonicus, Cajanus cajan (Cc), Phaseolus vulgaris (Pv), and Glycine max (Gm). Subsequently, sequence distances, numbers of synonymous substitutions per synonymous site (Ks) and non-synonymous substitutions per non-synonymous site (Ka) between orthologs and paralogs were calculated and compared across legumes. These analyses suggest the closest relationship between Gm and Cc and the highest distance between Mt and Pv in six legumes. Ks proportional plots clearly showed ancient genome duplication in all legumes, whole genome duplication event in Gm and also speciation pattern in different legumes. This study also reports some interesting observations e.g., no peak at Ks 0.4 in Gm-Gm, location of two independent genes next to each other in Mt and low Ks values for outparalogs for three genes as compared to other 12 genes. In summary, this study underlines the importance of NF-related genes and provides important insights in genome organization and evolutionary aspects of six legume species analyzed.

Plants are exposed to different external conditions that affect growth, development, and productivity. Water deficit is one of these adverse conditions caused by drought, salinity, and extreme temperatures. Plants have developed different responses to prevent, ameliorate or repair the damage inflicted by these stressful environments. One of these responses is the activation of a set of genes encoding a group of hydrophilic proteins that typically accumulate to high levels during seed dehydration, at the last stage of embryogenesis, hence named Late Embryogenesis Abundant (LEA) proteins. LEA proteins also accumulate in response to water limitation in vegetative tissues, and have been classified in seven groups based on their amino acid sequence similarity and on the presence of distinctive conserved motifs. These proteins are widely distributed in the plant kingdom, from ferns to angiosperms, suggesting a relevant role in the plant response to this unfavorable environmental condition. In this review, we analyzed the LEA proteins from those legumes whose complete genomes have been sequenced such as Phaseolus vulgaris, Glycine max, Medicago truncatula, Lotus japonicus, Cajanus cajan, and Cicer arietinum. Considering their distinctive motifs, LEA proteins from the different groups were identified, and their sequence analysis allowed the recognition of novel legume specific motifs. Moreover, we compile their transcript accumulation patterns based on publicly available data. In spite of the limited information on these proteins in legumes, the analysis and data compiled here confirm the high correlation between their accumulation and water deficit, reinforcing their functional relevance under this detrimental conditions.

The study identified 9045 high-quality SNPs employing both genome-wide GBS- and candidate gene-based SNP genotyping assays in 172, including 93 cultivated (desi and kabuli) and 79 wild chickpea accessions. The GWAS in a structured population of 93 sequenced accessions detected 15 major genomic loci exhibiting significant association with seed coat color. Five seed color-associated major genomic loci underlying robust QTLs mapped on a high-density intra-specific genetic linkage map were validated by QTL mapping. The integration of association and QTL mapping with gene haplotype-specific LD mapping and transcript profiling identified novel allelic variants (non-synonymous SNPs) and haplotypes in a MATE secondary transporter gene regulating light/yellow brown and beige seed coat color differentiation in chickpea. The down-regulation and decreased transcript expression of beige seed coat color-associated MATE gene haplotype was correlated with reduced proanthocyanidins accumulation in the mature seed coats of beige than light/yellow brown seed colored desi and kabuli accessions for their coloration/pigmentation. This seed color-regulating MATE gene revealed strong purifying selection pressure primarily in LB/YB seed colored desi and wild Cicer reticulatum accessions compared with the BE seed colored kabuli accessions. The functionally relevant molecular tags identified have potential to decipher the complex transcriptional regulatory gene function of seed coat coloration and for understanding the selective sweep-based seed color trait evolutionary pattern in cultivated and wild accessions during chickpea domestication. The genome-wide integrated approach employed will expedite marker-assisted genetic enhancement for developing cultivars with desirable seed coat color types in chickpea.

The impact of deleterious variation on both plant fitness and crop productivity is not completely understood and is a hot topic of debates. The deleterious mutations in plants have been solely predicted using sequence conservation methods rather than function-based classifiers due to lack of well-annotated mutational datasets in these organisms. Here, we developed a machine learning classifier based on a dataset of deleterious and neutral mutations in Arabidopsis thaliana by extracting 18 informative features that discriminate deleterious mutations from neutral, including 9 novel features not used in previous studies. We examined linear SVM, Gaussian SVM, and Random Forest classifiers, with the latter performing best. Random Forest classifiers exhibited a markedly higher accuracy than the popular PolyPhen-2 tool in the Arabidopsis dataset. Additionally, we tested whether the Random Forest, trained on the Arabidopsis dataset, accurately predicts deleterious mutations in Orýza sativa and Pisum sativum and observed satisfactory levels of performance accuracy (87% and 93%, respectively) higher than obtained by the PolyPhen-2. Application of Transfer learning in classifiers did not improve their performance. To additionally test the performance of the Random Forest classifier across different angiosperm species, we applied it to annotate deleterious mutations in Cicer arietinum and validated them using population frequency data. Overall, we devised a classifier with the potential to improve the annotation of putative functional mutations in QTL and GWAS hit regions, as well as for the evolutionary analysis of proliferation of deleterious mutations during plant domestication; thus optimizing breeding improvement and development of new cultivars.

This work was designed to evaluate whether external application of nitric oxide (NO) in the form of its donor S-nitroso-N-acetylpenicillamine (SNAP) could mitigate the deleterious effects of NaCl stress on chickpea (Cicer arietinum L.) plants. SNAP (50 μM) was applied to chickpea plants grown under non-saline and saline conditions (50 and 100 mM NaCl). Salt stress inhibited growth and biomass yield, leaf relative water content (LRWC) and chlorophyll content of chickpea plants. High salinity increased electrolyte leakage, carotenoid content and the levels of osmolytes (proline, glycine betaine, soluble proteins and soluble sugars), hydrogen peroxide (H2O2) and malondialdehyde (MDA), as well as the activities of antioxidant enzymes, such as superoxide dismutase (SOD), catalase (CAT), ascorbate peroxidase (APX), and glutathione reductase in chickpea plants. Expression of the representative SOD, CAT and APX genes examined was also up-regulated in chickpea plants by salt stress. On the other hand, exogenous application of NO to salinized plants enhanced the growth parameters, LRWC, photosynthetic pigment production and levels of osmolytes, as well as the activities of examined antioxidant enzymes which is correlated with up-regulation of the examined SOD, CAT and APX genes, in comparison with plants treated with NaCl only. Furthermore, electrolyte leakage, H2O2 and MDA contents showed decline in salt-stressed plants supplemented with NO as compared with those in NaCl-treated plants alone. Thus, the exogenous application of NO protected chickpea plants against salt stress-induced oxidative damage by enhancing the biosyntheses of antioxidant enzymes, thereby improving plant growth under saline stress. Taken together, our results demonstrate that NO has capability to mitigate the adverse effects of high salinity on chickpea plants by improving LRWC, photosynthetic pigment biosyntheses, osmolyte accumulation and antioxidative defense system.

Isoflavone reductase (IFR) is an enzyme involved in the biosynthetic pathway of isoflavonoid phytoalexin in plants. IFRs are unique to the plant kingdom and are considered to have crucial roles in plant response to various biotic and abiotic environmental stresses. Here, we report the characterization of a novel member of the soybean isoflavone reductase gene family GmIFR. Overexpression of GmIFR transgenic soybean exhibited enhanced resistance to Phytophthora sojae. Following stress treatments, GmIFR was significantly induced by P. sojae, ethephon (ET), abscisic acid (placeCityABA), salicylic acid (SA). It is located in the cytoplasm when transiently expressed in soybean protoplasts. The daidzein levels reduced greatly for the seeds of transgenic plants, while the relative content of glyceollins in transgenic plants was significantly higher than that of non-transgenic plants. Furthermore, we found that the relative expression levels of reactive oxygen species (ROS) of transgenic soybean plants were significantly lower than those of non-transgenic plants after incubation with P. sojae, suggesting an important role of GmIFR might function as an antioxidant to reduce ROS in soybean. The enzyme activity assay suggested that GmIFR has isoflavone reductase activity.

Transcription factors (TFs) are the key players in gene expression and their study is highly significant for shedding light on the molecular mechanisms and evolutionary history of organisms. During host-pathogen interaction, extensive reprogramming of gene expression facilitated by TFs is likely to occur in both host and pathogen. To date, the knowledge about TF repertoire in filamentous fungi is in infancy. The necrotrophic fungus Ascochyta rabiei, that causes destructive Ascochyta blight (AB) disease of chickpea (Cicer arietinum), demands more comprehensive study for better understanding of Ascochyta-legume pathosystem. In the present study, we performed the genome-wide identification and analysis of TFs in A. rabiei. Taking advantage of A. rabiei genome sequence, we used a bioinformatic approach to predict the TF repertoire of A. rabiei. For identification and classification of A. rabiei TFs, we designed a comprehensive pipeline using a combination of BLAST and InterProScan software. A total of 381 A. rabiei TFs were predicted and divided into 32 fungal specific families of TFs. The gene structure, domain organization and phylogenetic analysis of abundant families of A. rabiei TFs were also carried out. Comparative study of A. rabiei TFs with that of other necrotrophic, biotrophic, hemibiotrophic, symbiotic, and saprotrophic fungi was performed. It suggested presence of both conserved as well as unique features among them. Moreover, cis-acting elements on promoter sequences of earlier predicted A. rabiei secretome were also identified. With the help of published A. rabiei transcriptome data, the differential expression of TF and secretory protein coding genes was analyzed. Furthermore, comprehensive expression analysis of few selected A. rabiei TFs using quantitative real-time polymerase chain reaction revealed variety of expression patterns during host colonization. These genes were expressed in at least one of the time points tested post infection. Overall, this study illustrates the first genome-wide identification and analysis of TF repertoire of A. rabiei. This work would provide the basis for further studies to dissect role of TFs in the molecular mechanisms during A. rabiei-chickpea interactions.

Ascochyta blight (AB) is a fungal disease that can significantly reduce chickpea production in Australia and other regions of the world. In this study, 69 chickpea genotypes were sequenced using whole genome re-sequencing (WGRS) methods. They included 48 Australian varieties differing in their resistance ranking to AB, 16 advanced breeding lines from the Australian chickpea breeding program, four landraces, and one accession representing the wild chickpea species Cicer reticulatum. More than 800,000 single nucleotide polymorphisms (SNPs) were identified. Population structure analysis revealed relatively narrow genetic diversity amongst recently released Australian varieties and two groups of varieties separated by the level of AB resistance. Several regions of the chickpea genome were under positive selection based on Tajima's D test. Both Fst genome- scan and genome-wide association studies (GWAS) identified a 100 kb region (AB4.1) on chromosome 4 that was significantly associated with AB resistance. The AB4.1 region co-located to a large QTL interval of 7 Mb∼30 Mb identified previously in three different mapping populations which were genotyped at relatively low density with SSR or SNP markers. The AB4.1 region was validated by GWAS in an additional collection of 132 advanced breeding lines from the Australian chickpea breeding program, genotyped with approximately 144,000 SNPs. The reduced level of nucleotide diversity and long extent of linkage disequilibrium also suggested the AB4.1 region may have gone through selective sweeps probably caused by selection of the AB resistance trait in breeding. In total, 12 predicted genes were located in the AB4.1 QTL region, including those annotated as: NBS-LRR receptor-like kinase, wall-associated kinase, zinc finger protein, and serine/threonine protein kinases. One significant SNP located in the conserved catalytic domain of a NBS-LRR receptor-like kinase led to amino acid substitution. Transcriptional analysis using qPCR showed that some predicted genes were significantly induced in resistant lines after inoculation compared to non-inoculated plants. This study demonstrates the power of combining WGRS data with relatively simple traits to rapidly develop "functional makers" for marker-assisted selection and genomic selection.

Background: Chickpea (Cicer arietinum L.) contributes 75% of total pulse production. Being cheaper than animal protein, makes it important in dietary requirement of developing countries. Weed not only competes with chickpea resulting into drastic yield reduction but also creates problem of harboring fungi, bacterial diseases and insect pests. Chemical approach having new herbicide discovery has constraint of limited lead molecule options, statutory regulations and environmental clearance. Through genetic approach, transgenic herbicide tolerant crop has given successful result but led to serious concern over ecological safety thus non-transgenic approach like marker assisted selection is desirable. Since large variability in tolerance limit of herbicide already exists in chickpea varieties, thus the genes offering herbicide tolerance can be introgressed in variety improvement programme. Transcriptome studies can discover such associated key genes with herbicide tolerance in chickpea. Results: This is first transcriptomic studies of chickpea or even any legume crop using two herbicide susceptible and tolerant genotypes exposed to imidazoline (Imazethapyr). Approximately 90 million paired-end reads generated from four samples were processed and assembled into 30,803 contigs using reference based assembly. We report 6,310 differentially expressed genes (DEGs), of which 3,037 were regulated by 980 miRNAs, 1,528 transcription factors associated with 897 DEGs, 47 Hub proteins, 3,540 putative Simple Sequence Repeat-Functional Domain Marker (SSR-FDM), 13,778 genic Single Nucleotide Polymorphism (SNP) putative markers and 1,174 Indels. Randomly selected 20 DEGs were validated using qPCR. Pathway analysis suggested that xenobiotic degradation related gene, glutathione S-transferase (GST) were only up-regulated in presence of herbicide. Down-regulation of DNA replication genes and up-regulation of abscisic acid pathway genes were observed. Study further reveals the role of cytochrome P450, xyloglucan endotransglucosylase/hydrolase, glutamate dehydrogenase, methyl crotonoyl carboxylase and of thaumatin-like genes in herbicide resistance. Conclusion: Reported DEGs can be used as genomic resource for future discovery of candidate genes associated with herbicide tolerance. Reported markers can be used for future association studies in order to develop marker assisted selection (MAS) for refinement. In endeavor of chickpea variety development programme, these findings can be of immense use in improving productivity of chickpea germplasm.

Chickpea (Cicer arietinum); the second largest legume grown worldwide is prone to drought and various pathogen infections. These drought and pathogen stresses often occur concurrently in the field conditions. However, the molecular events in response to that are largely unknown. The present study examines the transcriptome dynamics in chickpea plants exposed to a combination of water-deficit stress and Ralstonia solanacearum infection. R. solanacearum is a potential wilt disease causing pathogen in chickpea. Drought stressed chickpea plants were infected with this pathogen and the plants were allowed to experience progressive drought with 2 and 4 days of R. solanacearum infection called short duration stress (SD stresses) and long duration stress (LD stresses), respectively. Our study showed that R. solanacearum multiplication decreased under SD-combined stress compared to SD-pathogen but there was no significant change in LD-combined stress compared to LD-pathogen. The microarray analysis during these conditions showed that 821 and 1039 differentially expressed genes (DEGs) were unique to SD- and LD-combined stresses, respectively, when compared with individual stress conditions. Three and fifteen genes were common among all the SD-stress treatments and LD-stress treatments, respectively. Genes involved in secondary cell wall biosynthesis, alkaloid biosynthesis, defense related proteins, and osmo-protectants were up-regulated during combined stress. The expression of genes involved in lignin and cellulose biosynthesis were specifically up-regulated in SD-combined, LD-combined, and LD-pathogen stress. A close transcriptomic association of LD-pathogen stress with SD-combined stress was observed in this study which indicates that R. solanacearum infection also exerts drought stress along with pathogen stress thus mimics combined stress effect. Furthermore the expression profiling of candidate genes using real-time quantitative PCR validated the microarray data. The study showed that down-regulation of defense-related genes during LD-combined stress resulted in an increased bacterial multiplication as compared to SD-combined stress. Overall, our study highlights a sub-set of DEGs uniquely expressed in response to combined stress, which serve as potential candidates for further functional characterization to delineate the molecular response of the plant to concurrent drought-pathogen stress.

Chickpea is an economically and nutritionally important grain legume globally, however, cold stress has adverse effects on its growth. In cold countries, like Canada where the growing season is short, having cold stress-tolerant varieties is crucial. Crop wild relatives of chickpea, especially Cicer reticulatum, can survive in suboptimal environments and are an important resource for crop improvement. In this study, we explored the performance of eleven C. reticulatum wild accessions and two chickpea cultivars, CDC Leader and CDC Consul, together with a cold sensitive check ILC533 under freezing stress. Freezing tolerance was scored based on a 1-9 scale. The wild relatives, particularly Kesen_075 and CudiA_152, had higher frost tolerance compared to the cultivars, which all died after frost treatment. We completed transcriptome analysis via mRNA sequencing to assess changes in gene expression in response to freezing stress and identified 6,184 differentially expressed genes (DEGs) in CDC Consul, and 7,842 DEGs in Kesen_075. GO (gene ontology) analysis of the DEGs revealed that those related to stress responses, endogenous and external stimuli responses, secondary metabolite processes, and photosynthesis were significantly over-represented in CDC Consul, while genes related to endogenous stimulus responses and photosynthesis were significantly over-represented in Kesen_075. These results are consistent with Kesen_075 being more tolerant to freezing stress than CDC Consul. Moreover, our data revealed that the expression of CBF pathway-related genes was impacted during freezing conditions in Kesen_075, and expression of these genes is believed to alleviate the damage caused by freezing stress. We identified genomic regions associated with tolerance to freezing stress in an F2 population derived from a cross between CDC Consul and Kesen_075 using QTL-seq analysis. Eight QTLs (P<0.05) on chromosomes Ca3, Ca4, Ca6, Ca7, Ca8, and two QTLs (P<0.01) on chromosomes Ca4 and Ca8, were associated with tolerance to freezing stress. Interestingly, 58 DEGs co-located within these QTLs. To our knowledge, this is the first study to explore the transcriptome and QTLs associated with freezing tolerance in wild relatives of chickpea under controlled conditions. Altogether, these findings provide comprehensive information that aids in understanding the molecular mechanism of chickpea adaptation to freezing stress and further provides functional candidate genes that can assist in breeding of freezing-stress tolerant varieties.

Iron deficiency currently affects over two billion people worldwide despite significant advances in technology and society aimed at mitigating this global health problem. Biofortification of food staples with iron (Fe) represents a sustainable approach for alleviating human Fe deficiency in developing countries, however, biofortification efforts have focused extensively on cereal staples while pulses have been largely overlooked. In this study we describe a genetic engineering (GE) approach to biofortify the pulse crop, chickpea (Cicer arietinum L.), with Fe using a combination of the chickpea nicotianamine synthase 2 (CaNAS2) and soybean (Glycine max) ferritin (GmFER) genes which function in Fe transport and storage, respectively. This study consists of three main components: (1) the establishment for baseline Fe concentration of existing germplam, (2) the isolation and study of expression pattern of the novel CaNAS2 gene, and (3) the generation of GE chickpea overexpressing the CaNAS2 and GmFER genes. Seed of six commercial chickpea cultivars was collected from four different field locations in Australia and assessed for seed Fe concentration. The results revealed little difference between the cultivars assessed, and that chickpea seed Fe was negatively affected where soil Fe bioavailability is low. The desi cultivar HatTrick was then selected for further study. From it, the CaNAS2 gene was cloned and its expression in different tissues examined. The gene was found to be expressed in multiple vegetative tissues under Fe-sufficient conditions, suggesting that it may play a housekeeping role in systemic translocation of Fe. Two GE chickpea events were then generated and the overexpression of the CaNAS2 and GmFER transgenes confirmed. Analysis of nicotianamine (NA) and Fe levels in the GE seeds revealed that NA was nearly doubled compared to the null control while Fe concentration was not changed. Increased NA content in chickpea seed is likely to translate into increased Fe bioavailability and may thus overcome the effect of the bioavailability inhibitors found in pulses; however, further study is required to confirm this. This is the first known example of GE Fe biofortified chickpea; information gleaned from this study can feed into future pulse biofortification work to help alleviate global Fe deficiency.

Chickpea (Cicer arietinum L.) is an important cool season food legume, however, its production is severely constrained by the foliar disease Ascochyta blight caused by the fungus Ascochyta rabiei (syn. Phoma rabiei). Several disease management options have been developed to control the pathogen, including breeding for host plant resistance. However, the pathogen population is evolving to produce more aggressive isolates. For host resistance to be effective, the plant must quickly recognize the pathogen and instigate initial defense mechanisms, optimally at the point of contact. Given that the most resistant host genotypes display rapid pathogen recognition and response, the approach taken was to assess the type, speed and pattern of recognition via Resistance Gene Analog (RGA) transcription among resistant and susceptible cultivated chickpea varieties. RGAs are key factors in the recognition of plant pathogens and the signaling of inducible defenses. In this study, a suite of RGA loci were chosen for further investigation from both published literature and from newly mined homologous sequences within the National Center for Biotechnology Information (NCBI) database. Following their validation in the chickpea genome, 10 target RGAs were selected for differential expression analysis in response to A. rabiei infection. This was performed in a set of four chickpea varieties including two resistant cultivars (ICC3996 and PBA Seamer), one moderately resistant cultivar (PBA HatTrick) and one susceptible cultivar (Kyabra). Gene expression at each RGA locus was assessed via qPCR at 2, 6, and 24 h after A. rabiei inoculation with a previously characterized highly aggressive isolate. As a result, all loci were differentially transcribed in response to pathogen infection in at least one genotype and at least one time point after inoculation. Among these, the differential expression of four RGAs was significant and consistently increased in the most resistant genotype ICC3996 immediately following inoculation, when spore germination began and ahead of penetration into the plant's epidermal tissues. Further in silico analyses indicated that the differentially transcribed RGAs function through ADP-binding within the pathogen recognition pathway. These represent clear targets for future functional validation and potential for selective resistance breeding for introgression into elite cultivars.

Red clover (Trifolium pratense) is an important forage plant worldwide. This study was directed to broadening current knowledge of red clover's coding regions and enhancing its utilization in practice by specific reanalysis of previously published assembly. A total of 42,996 genes were characterized using Illumina paired-end sequencing after manual revision of Blast2GO annotation. Genes were classified into metabolic and biosynthetic pathways in response to biological processes, with 7,517 genes being assigned to specific pathways. Moreover, 17,727 enzymatic nodes in all pathways were described. We identified 6,749 potential microsatellite loci in red clover coding sequences, and we characterized 4,005 potential simple sequence repeat (SSR) markers as generating polymerase chain reaction products preferentially within 100-350 bp. Marker density of 1 SSR marker per 12.39 kbp was achieved. Aligning reads against predicted coding sequences resulted in the identification of 343,027 single nucleotide polymorphism (SNP) markers, providing marker density of one SNP marker per 144.6 bp. Altogether, 95 SSRs in coding sequences were analyzed for 50 red clover varieties and a collection of 22 highly polymorphic SSRs with pooled polymorphism information content >0.9 was generated, thus obtaining primer pairs for application to diversity studies in T. pratense. A set of 8,623 genome-wide distributed SNPs was developed and used for polymorphism evaluation in individual plants. The polymorphic information content ranged from 0 to 0.375. Temperature switch PCR was successfully used in single-marker SNP genotyping for targeted coding sequences and for heterozygosity or homozygosity confirmation in validated five loci. Predicted large sets of SSRs and SNPs throughout the genome are key to rapidly implementing genome-based breeding approaches, for identifying genes underlying key traits, and for genome-wide association studies. Detailed knowledge of genetic relationships among breeding material can also be useful for breeders in planning crosses or for plant variety protection. Single-marker assays are useful for diagnostic applications.

Sucrose transporters (SUTs) are essential for the export and efficient movement of sucrose from source leaves to sink organs in plants. The angiosperm SUT family was previously classified into three or four distinct groups, Types I, II (subgroup IIB), and III, with dicot-specific Type I and monocot-specific Type IIB functioning in phloem loading. To shed light on the underlying drivers of SUT evolution, Bayesian phylogenetic inference was undertaken using 41 sequenced plant genomes, including seven basal lineages at key evolutionary junctures. Our analysis supports four phylogenetically and structurally distinct SUT subfamilies, originating from two ancient groups (AG1 and AG2) that diverged early during terrestrial colonization. In both AG1 and AG2, multiple intron acquisition events in the progenitor vascular plant established the gene structures of modern SUTs. Tonoplastic Type III and plasmalemmal Type II represent evolutionarily conserved descendants of AG1 and AG2, respectively. Type I and Type IIB were previously thought to evolve after the dicot-monocot split. We show, however, that divergence of Type I from Type III SUT predated basal angiosperms, likely associated with evolution of vascular cambium and phloem transport. Type I SUT was subsequently lost in monocots along with vascular cambium, and independent evolution of Type IIB coincided with modified monocot vasculature. Both Type I and Type IIB underwent lineage-specific expansion. In multiple unrelated taxa, the newly-derived SUTs exhibit biased expression in reproductive tissues, suggesting a functional link between phloem loading and reproductive fitness. Convergent evolution of Type I and Type IIB for SUT function in phloem loading and reproductive organs supports the idea that differential vascular development in dicots and monocots is a strong driver for SUT family evolution in angiosperms.

Heat shock transcription factors (Hsfs) are important transcription factors (TFs) in protecting plants from damages caused by various stresses. The released whole genome sequences of wild peanuts make it possible for genome-wide analysis of Hsfs in peanut. In this study, a total of 16 and 17 Hsf genes were identified from Arachis duranensis and A. ipaensis, respectively. We identified 16 orthologous Hsf gene pairs in both peanut species; however HsfXs was only identified from A. ipaensis. Orthologous pairs between two wild peanut species were highly syntenic. Based on phylogenetic relationship, peanut Hsfs were divided into groups A, B, and C. Selection pressure analysis showed that group B Hsf genes mainly underwent positive selection and group A Hsfs were affected by purifying selection. Small scale segmental and tandem duplication may play important roles in the evolution of these genes. Cis-elements, such as ABRE, DRE, and HSE, were found in the promoters of most Arachis Hsf genes. Five AdHsfs and two AiHsfs contained fungal elicitor responsive elements suggesting their involvement in response to fungi infection. These genes were differentially expressed in cultivated peanut under abiotic stress and Aspergillus flavus infection. AhHsf2 and AhHsf14 were significantly up-regulated after inoculation with A. flavus suggesting their possible role in fungal resistance.

Jasmonates (JA) are well-known phytohormones which play important roles in plant development and defense against pathogens. Jasmonate ZIM domain (JAZ) proteins are plant-specific proteins and act as transcriptional repressors of JA-responsive genes. JA regulates both biotic and abiotic stress responses in plants; however, its role in nutrient deficiency responses is very elusive. Although, JA is well-known for root growth inhibition, little is known about behavior of JAZ genes in response to nutrient deficiencies, under which root architectural alteration is an important adaptation. Using protein sequence homology and a conserved-domains approach, here we identify 10 novel JAZ genes from the recently sequenced Chickpea genome, which is one of the most nutrient efficient crops. Both rice and chickpea JAZ genes express in tissue- and stimuli-specific manners. Many of which are preferentially expressed in root. Our analysis further showed differential expression of JAZ genes under macro (NPK) and micronutrients (Zn, Fe) deficiency in rice and chickpea roots. While both rice and chickpea JAZ genes showed a certain level of specificity toward type of nutrient deficiency, generally majority of them showed induction under K deficiency. Generally, JAZ genes showed an induction at early stages of stress and expression declined at later stages of macro-nutrient deficiency. Our results suggest that JAZ genes might play a role in early nutrient deficiency response both in monocot and dicot roots, and information generated here can be further used for understanding the possible roles of JA in root architectural alterations for nutrient deficiency adaptations.

Lupins, like other legumes, have a unique biosynthesis scheme of 5-deoxy-type flavonoids and isoflavonoids. A key enzyme in this pathway is chalcone isomerase (CHI), a member of CHI-fold protein family, encompassing subfamilies of CHI1, CHI2, CHI-like (CHIL), and fatty acid-binding (FAP) proteins. Here, two Lupinus angustifolius (narrow-leafed lupin) CHILs, LangCHIL1 and LangCHIL2, were identified and characterized using DNA fingerprinting, cytogenetic and linkage mapping, sequencing and expression profiling. Clones carrying CHIL sequences were assembled into two contigs. Full gene sequences were obtained from these contigs, and mapped in two L. angustifolius linkage groups by gene-specific markers. Bacterial artificial chromosome fluorescence in situ hybridization approach confirmed the localization of two LangCHIL genes in distinct chromosomes. The expression profiles of both LangCHIL isoforms were very similar. The highest level of transcription was in the roots of the third week of plant growth; thereafter, expression declined. The expression of both LangCHIL genes in leaves and stems was similar and low. Comparative mapping to reference legume genome sequences revealed strong syntenic links; however, LangCHIL2 contig had a much more conserved structure than LangCHIL1. LangCHIL2 is assumed to be an ancestor gene, whereas LangCHIL1 probably appeared as a result of duplication. As both copies are transcriptionally active, questions arise concerning their hypothetical functional divergence. Screening of the narrow-leafed lupin genome and transcriptome with CHI-fold protein sequences, followed by Bayesian inference of phylogeny and cross-genera synteny survey, identified representatives of all but one (CHI1) main subfamilies. They are as follows: two copies of CHI2, FAPa2 and CHIL, and single copies of FAPb and FAPa1. Duplicated genes are remnants of whole genome duplication which is assumed to have occurred after the divergence of Lupinus, Arachis, and Glycine.