Genome-wide survey and comprehensive expression profiling of Aux/IAA gene family in chickpea and soybean

Characterization of the Soluble NSF Attachment Protein gene family identifies two members involved in additive resistance to a plant pathogen

Proteins with Tetratricopeptide-repeat (TPR) domains are encoded by large gene families and distributed in all plant lineages. In this study, the Soluble NSF-Attachment Protein (SNAP) subfamily of TPR containing proteins is characterized. In soybean, five members constitute the SNAP gene family: GmSNAP18, GmSNAP11, GmSNAP14, GmSNAP02, and GmSNAP09. Recently, GmSNAP18 has been reported to mediate resistance to soybean cyst nematode (SCN). Using a population of recombinant inbred lines from resistant and susceptible parents, the divergence of the SNAP gene family is analysed over time. Phylogenetic analysis of SNAP genes from 22 diverse plant species showed that SNAPs were distributed in six monophyletic clades corresponding to the major plant lineages. Conservation of the four TPR motifs in all species, including ancestral lineages, supports the hypothesis that SNAPs were duplicated and derived from a common ancestor and unique gene still present in chlorophytic algae. Syntenic analysis of regions harbouring GmSNAP genes in soybean reveals that this family expanded from segmental and tandem duplications following a tetraploidization event. qRT-PCR analysis of GmSNAPs indicates a co-regulation following SCN infection. Finally, genetic analysis demonstrates that GmSNAP11 contributes to an additive resistance to SCN. Thus, GmSNAP11 is identified as a novel minor gene conferring resistance to SCN.

Evolution and structural diversification of Nictaba-like lectin genes in food crops with a focus on soybean (Glycine max)

Background and Aims

The Nictaba family groups all proteins that show homology to Nictaba, the tobacco lectin. So far, Nictaba and an Arabidopsis thaliana homologue have been shown to be implicated in the plant stress response. The availability of more than 50 sequenced plant genomes provided the opportunity for a genome-wide identification of Nictaba -like genes in 15 species, representing members of the Fabaceae, Poaceae, Solanaceae, Musaceae, Arecaceae, Malvaceae and Rubiaceae. Additionally, phylogenetic relationships between the different species were explored. Furthermore, this study included domain organization analysis, searching for orthologous genes in the legume family and transcript profiling of the Nictaba -like lectin genes in soybean.

Methods

Using a combination of BLASTp, InterPro analysis and hidden Markov models, the genomes of Medicago truncatula , Cicer arietinum , Lotus japonicus , Glycine max , Cajanus cajan , Phaseolus vulgaris , Theobroma cacao , Solanum lycopersicum , Solanum tuberosum , Coffea canephora , Oryza sativa , Zea mays, Sorghum bicolor , Musa acuminata and Elaeis guineensis were searched for Nictaba -like genes. Phylogenetic analysis was performed using RAxML and additional protein domains in the Nictaba-like sequences were identified using InterPro. Expression analysis of the soybean Nictaba -like genes was investigated using microarray data.

Key Results

Nictaba -like genes were identified in all studied species and analysis of the duplication events demonstrated that both tandem and segmental duplication contributed to the expansion of the Nictaba gene family in angiosperms. The single-domain Nictaba protein and the multi-domain F-box Nictaba architectures are ubiquitous among all analysed species and microarray analysis revealed differential expression patterns for all soybean Nictaba-like genes.

Conclusions

Taken together, the comparative genomics data contributes to our understanding of the Nictaba -like gene family in species for which the occurrence of Nictaba domains had not yet been investigated. Given the ubiquitous nature of these genes, they have probably acquired new functions over time and are expected to take on various roles in plant development and defence.

Genetic Analysis of NBS-LRR Gene Family in Chickpea and Their Expression Profiles in Response to Ascochyta Blight Infection

Ascochyta blight is one of the major diseases of chickpea worldwide. The genetic resistance to ascochyta blight in chickpea is complex and governed by multiple QTLs. However, the molecular mechanism of quantitative disease resistance to ascochyta blight and the genes underlying these QTLs are still unknown. Most often disease resistance is determined by resistance (R) genes. The most predominant R-genes contain nucleotide binding site and leucine rich repeat (NBS-LRR) domains. A total of 121 NBS-LRR genes were identified in the chickpea genome. Ninety-eight of these genes contained all essential conserved domains while 23 genes were truncated. The NBS-LRR genes were grouped into eight distinct classes based on their domain architecture. Phylogenetic analysis grouped these genes into two major clusters based on their structural variation, the first cluster with toll or interleukin-1 like receptor (TIR) domain and the second cluster either with or without a coiled-coil domain. The NBS-LRR genes are distributed unevenly across the eight chickpea chromosomes and nearly 50% of the genes are present in clusters. Thirty of the NBS-LRR genes were co-localized with nine of the previously reported ascochyta blight QTLs and were tested as potential candidate genes for ascochyta blight resistance. Expression pattern of these genes was studied in two resistant (CDC Corinne and CDC Luna) and one susceptible (ICCV 96029) genotypes at different time points after ascochyta blight infection using real-time quantitative PCR. Twenty-seven NBS-LRR genes showed differential expression in response to ascochyta blight infection in at least one genotype at one time point. Among these 27 genes, the majority of the NBS-LRR genes showed differential expression after inoculation in both resistant and susceptible genotypes which indicates the involvement of these genes in response to ascochyta blight infection. Five NBS-LRR genes showed genotype specific expression. Our study provides a new insight of NBS-LRR gene family in chickpea and the potential involvement of NBS-LRR genes in response to ascochyta blight infection.

Characterization of the FAD2 Gene Family in Soybean Reveals the Limitations of Gel-Based TILLING in Genes with High Copy Number

Soybean seed oil typically contains 18-20% oleic acid. Increasing the content of oleic acid is beneficial for health and biodiesel production. Mutations in FAD2-1 genes have been reported to increase seed oleic acid content. A subset of 1,037 mutant families from a mutagenized soybean cultivar (cv.) Forrest population was screened using reverse genetics (TILLING) to identify mutations within FAD2 genes. Although no fad2 mutants were identified using gel-based TILLING, four fad2-1A and one fad2-1B mutants were identified to have high seed oleic acid content using forward genetic screening and subsequent target sequencing. TILLING has been successfully used as a non-transgenic reverse genetic approach to identify mutations in genes controlling important agronomic traits. However, this technique presents limitations in traits such as oil composition due to gene copy number and similarities within the soybean genome. In soybean, FAD2 are present as two copies, FAD2-1 and FAD2-2. Two FAD2-1 members: FAD2-1A and FAD2-1B; and three FAD2-2 members: FAD2-2A, FAD2-2B, and FAD2-2C have been reported. Syntenic, phylogenetic, and in silico analysis revealed two additional members constituting the FAD2 gene family: GmFAD2-2D and GmFAD2-2E, located on chromosomes 09 and 15, respectively. They are presumed to have diverged from other FAD2-2 members localized on chromosomes 19 (GmFAD2-2A and GmFAD2-2B) and 03 (GmFAD2-2C). This work discusses alternative solutions to the limitations of gel-based TILLING in functional genomics due to high copy number and multiple paralogs of the FAD2 gene family in soybean.

Genome-Wide Analysis of Soybean JmjC Domain-Containing Proteins Suggests Evolutionary Conservation Following Whole-Genome Duplication

Histone modifications, such as methylation and demethylation, play an important role in regulating chromatin structure and gene expression. The JmjC domain-containing proteins, an important family of histone lysine demethylases (KDMs), play a key role in maintaining homeostasis of histone methylation in vivo. In this study, we performed a comprehensive analysis of the jumonji C (JmjC) gene family in the soybean genome and identified 48 JmjC genes (GmJMJs) distributed unevenly across 18 chromosomes. Phylogenetic analysis showed that these JmjC domain-containing genes can be divided into eight groups. GmJMJs within the same phylogenetic group share similar exon/intron organization and domain composition. In addition, 16 duplicated gene pairs were formed by a Glycine-specific whole-genome duplication (WGD) event approximately 13 million years ago (Mya). By investigating the expression profiles of these gene pairs in various tissues, we showed that the expression pattern is conserved in the polyploidy-derived JmjC duplicates, demonstrating that the majority of GmJMJs were preferentially retained after the most recent WGD event and suggesting important roles for demethylase duplications in soybean evolution. These results shed light on the evolutionary history of this family in soybean and provide insights into the JmjCs which will be helpful to reveal their functions in controlling soybean development.

The SOD Gene Family in Tomato: Identification, Phylogenetic Relationships, and Expression Patterns

Superoxide dismutases (SODs) are critical antioxidant enzymes that protect organisms from reactive oxygen species (ROS) caused by adverse conditions, and have been widely found in the cytoplasm, chloroplasts, and mitochondria of eukaryotic and prokaryotic cells. Tomato (Solanum lycopersicum L.) is an important economic crop and is cultivated worldwide. However, abiotic and biotic stresses severely hinder growth and development of the plant, which affects the production and quality of the crop. To reveal the potential roles of SOD genes under various stresses, we performed a systematic analysis of the tomato SOD gene family and analyzed the expression patterns of SlSOD genes in response to abiotic stresses at the whole-genome level. The characteristics of the SlSOD gene family were determined by analyzing gene structure, conserved motifs, chromosomal distribution, phylogenetic relationships, and expression patterns. We determined that there are at least nine SOD genes in tomato, including four Cu/ZnSODs, three FeSODs, and one MnSOD, and they are unevenly distributed on 12 chromosomes. Phylogenetic analyses of SOD genes from tomato and other plant species were separated into two groups with a high bootstrap value, indicating that these SOD genes were present before the monocot-dicot split. Additionally, many cis-elements that respond to different stresses were found in the promoters of nine SlSOD genes. Gene expression analysis based on RNA-seq data showed that most genes were expressed in all tested tissues, with the exception of SlSOD6 and SlSOD8, which were only expressed in young fruits. Microarray data analysis showed that most members of the SlSOD gene family were altered under salt- and drought-stress conditions. This genome-wide analysis of SlSOD genes helps to clarify the function of SlSOD genes under different stress conditions and provides information to aid in further understanding the evolutionary relationships of SOD genes in plants.

EcoTILLING-Based Association Mapping Efficiently Delineates Functionally Relevant Natural Allelic Variants of Candidate Genes Governing Agronomic Traits in Chickpea

The large-scale mining and high-throughput genotyping of novel gene-based allelic variants in natural mapping population are essential for association mapping to identify functionally relevant molecular tags governing useful agronomic traits in chickpea. The present study employs an alternative time-saving, non-laborious and economical pool-based EcoTILLING approach coupled with agarose gel detection assay to discover 1133 novel SNP allelic variants from diverse coding and regulatory sequence components of 1133 transcription factor (TF) genes by genotyping in 192 diverse desi and kabuli chickpea accessions constituting a seed weight association panel. Integrating these SNP genotyping data with seed weight field phenotypic information of 192 structured association panel identified eight SNP alleles in the eight TF genes regulating seed weight of chickpea. The associated individual and combination of all SNPs explained 10-15 and 31% phenotypic variation for seed weight, respectively. The EcoTILLING-based large-scale allele mining and genotyping strategy implemented for association mapping is found much effective for a diploid genome crop species like chickpea with narrow genetic base and low genetic polymorphism. This optimized approach thus can be deployed for various genomics-assisted breeding applications with optimal expense of resources in domesticated chickpea. The seed weight-associated natural allelic variants and candidate TF genes delineated have potential to accelerate marker-assisted genetic improvement of chickpea.