Comprehensive expression analysis of rice Armadillo gene family during abiotic stress and development

Genome-wide mapping of DNase I hypersensitive sites revealed differential chromatin accessibility and regulatory DNA elements under drought stress in rice cultivars

Drought stress (DS) is one of the major constraints limiting yield in crop plants including rice. Gene regulation under DS is largely governed by accessibility of the transcription factors (TFs) to their cognate cis-regulatory elements (CREs). In this study, we used DNase I hypersensitive assays followed by sequencing to identify the accessible chromatin regions under DS in a drought-sensitive (IR64) and a drought-tolerant (N22) rice cultivar. Our results indicated that DNase I hypersensitive sites (DHSs) were highly enriched at transcription start sites (TSSs) and numerous DHSs were detected in the promoter regions. DHSs were concurrent with epigenetic marks and the genes harboring DHSs in their TSS and promoter regions were highly expressed. In addition, DS induced changes in DHSs (∆DHSs) in TSS and promoter regions were positively correlated with upregulation of several genes involved in drought/abiotic stress response, those encoding TFs and located within drought-associated quantitative trait loci, much preferentially in the drought-tolerant cultivar. The CREs representing the binding sites of TFs involved in DS response were detected within the ∆DHSs, suggesting differential accessibility of TFs to their cognate sites under DS in different rice cultivars, which may be further deployed for enhancing drought tolerance in rice.

A systematic analysis of ARM genes revealed that GhARM144 regulates the resistance against Verticillium dahliae via interaction with GhOSM34

Proteins of the armadillo repeat gene family play important roles in plant pathogen response. Here, 169 armadillo (ARM) genes were identified in upland cotton (Gossypium hirsutum). Phylogenetic analysis grouped these into 11 subfamilies, with conserved protein structures within each subfamily. The results signify that the expansion of the gene family occurred via whole genome duplication and dispersed duplication. Expression profiling and network analysis suggest that GhARM144 may regulate cotton resistance to Verticillium dahliae. GhARM144 was upregulated in roots by V. dahliae infection or salicylic acid treatment. This upregulation indicates a negative regulatory role of GhARM144' in the cotton immune responses, potentially by manipulating salicylic acid biosynthesis. Protein interaction studies found that GhARM144 associates with an osmotin-like protein, GhOSM34, at the plasma membrane. Silencing GhOSM34 reduced the resistance to V. dahliae, suggesting it may play a positive regulatory role. The results demonstrate that GhARM144 modulates cotton immunity through interaction with GhOSM34 and salicylic acid signalling. Further study of these proteins may yield insights into disease resistance mechanisms in cotton and other plants.

Genome-wide identification and transcriptome profiling expression analysis of the U-box E3 ubiquitin ligase gene family related to abiotic stress in maize (Zea mays L.)

BACKGROUND

The U-box gene family encodes E3 ubiquitin ligases involved in plant hormone signaling pathways and abiotic stress responses. However, there has yet to be a comprehensive analysis of the U-box gene family in maize (Zea mays L.) and its responses to abiotic stress.

RESULTS

In this study, 85 U-box family proteins were identified in maize and were classified into four subfamilies based on phylogenetic analysis. In addition to the conserved U-box domain, we identified additional functional domains, including Pkinase, ARM, KAP and Tyr domains, by analyzing the conserved motifs and gene structures. Chromosomal localization and collinearity analysis revealed that gene duplications may have contributed to the expansion and evolution of the U-box gene family. GO annotation and KEGG pathway enrichment analysis identified a total of 105 GO terms and 21 KEGG pathways that were notably enriched, including ubiquitin-protein transferase activity, ubiquitin conjugating enzyme activity and ubiquitin-mediated proteolysis pathway. Tissue expression analysis showed that some ZmPUB genes were specifically expressed in certain tissues and that this could be due to their functions. In addition, RNA-seq data for maize seedlings under salt stress revealed 16 stress-inducible plant U-box genes, of which 10 genes were upregulated and 6 genes were downregulated. The qRT-PCR results for genes responding to abiotic stress were consistent with the transcriptome analysis. Among them, ZmPUB13, ZmPUB18, ZmPUB19 and ZmPUB68 were upregulated under all three abiotic stress conditions. Subcellular localization analysis showed that ZmPUB19 and ZmPUB59 were located in the nucleus.

CONCLUSIONS

Overall, our study provides a comprehensive analysis of the U-box gene family in maize and its responses to abiotic stress, suggesting that U-box genes play an important role in the stress response and providing insights into the regulatory mechanisms underlying the response to abiotic stress in maize.

Meta-QTL analysis in wheat: progress, challenges and opportunities

Wheat, an important cereal crop globally, faces major challenges due to increasing global population and changing climates. The production and productivity are challenged by several biotic and abiotic stresses. There is also a pressing demand to enhance grain yield and quality/nutrition to ensure global food and nutritional security. To address these multifaceted concerns, researchers have conducted numerous meta-QTL (MQTL) studies in wheat, resulting in the identification of candidate genes that govern these complex quantitative traits. MQTL analysis has successfully unraveled the complex genetic architecture of polygenic quantitative traits in wheat. Candidate genes associated with stress adaptation have been pinpointed for abiotic and biotic traits, facilitating targeted breeding efforts to enhance stress tolerance. Furthermore, high-confidence candidate genes (CGs) and flanking markers to MQTLs will help in marker-assisted breeding programs aimed at enhancing stress tolerance, yield, quality and nutrition. Functional analysis of these CGs can enhance our understanding of intricate trait-related genetics. The discovery of orthologous MQTLs shared between wheat and other crops sheds light on common evolutionary pathways governing these traits. Breeders can leverage the most promising MQTLs and CGs associated with multiple traits to develop superior next-generation wheat cultivars with improved trait performance. This review provides a comprehensive overview of MQTL analysis in wheat, highlighting progress, challenges, validation methods and future opportunities in wheat genetics and breeding, contributing to global food security and sustainable agriculture.

Armadillo repeat only protein GS10 negatively regulates brassinosteroid signaling to control rice grain size

Grain yield and grain quality are major determinants in modern breeding controlled by many quantitative traits loci (QTLs) in rice (Oryza sativa). However, the mechanisms underlying grain shape and quality are poorly understood. Here, we characterize a QTL for grain size and grain quality via map-based cloning from wild rice (W1943), GS10 (Grain Size on Chromosome 10), which encodes a protein with 6 tandem armadillo repeats. The null mutant gs10 shows slender and narrow grains with altered cell size, which has a pleiotropic effect on other agronomical traits. Functional analysis reveals that GS10 interacts with TUD1 (Taihu Dwarf1) and is epistatic to OsGSK2 (glycogen synthase kinase 2) through regulating grain shape and lamina joint inclination, indicating it is negatively involved in brassinosteroid (BR) signaling. Pyramiding gs10 and the grain size gene GW5 into cultivar GLA4 substantially improved grain shape and appearance quality. Natural variation analysis revealed that gs10 from the wild rice Oryza rufipogon W1943 is a rare allele across the rice population. Collectively, these findings advance our understanding of the underlying mechanism of grain shape and provide the beneficial allele of gs10 for future rice breeding and genetic improvement.

Genome-wide association mapping of genomic regions associated with drought stress tolerance at seedling and reproductive stages in bread wheat

Understanding the genetic architecture of drought stress tolerance in bread wheat at seedling and reproductive stages is crucial for developing drought-tolerant varieties. In the present study, 192 diverse wheat genotypes, a subset from the Wheat Associated Mapping Initiative (WAMI) panel, were evaluated at the seedling stage in a hydroponics system for chlorophyll content (CL), shoot length (SLT), shoot weight (SWT), root length (RLT), and root weight (RWT) under both drought and optimum conditions. Following that, a genome-wide association study (GWAS) was carried out using the phenotypic data recorded during the hydroponics experiment as well as data available from previously conducted multi-location field trials under optimal and drought stress conditions. The panel had previously been genotyped using the Infinium iSelect 90K SNP array with 26,814 polymorphic markers. Using single as well as multi-locus models, GWAS identified 94 significant marker-trait associations (MTAs) or SNPs associated with traits recorded at the seedling stage and 451 for traits recorded at the reproductive stage. The significant SNPs included several novel, significant, and promising MTAs for different traits. The average LD decay distance for the whole genome was approximately 0.48 Mbp, ranging from 0.07 Mbp (chromosome 6D) to 4.14 Mbp (chromosome 2A). Furthermore, several promising SNPs revealed significant differences among haplotypes for traits such as RLT, RWT, SLT, SWT, and GY under drought stress. Functional annotation and in silico expression analysis revealed important putative candidate genes underlying the identified stable genomic regions such as protein kinases, O-methyltransferases, GroES-like superfamily proteins, NAD-dependent dehydratases, etc. The findings of the present study may be useful for improving yield potential, and stability under drought stress conditions.

Enrichment of intrinsically disordered residues in ohnologs facilitates abiotic stress resilience in Brassica rapa

Arabidopsis thaliana and Brassica rapa are in the same evolutionary lineage, although the latter experienced an additional whole genome triplication event. Therefore, it would be intriguing to investigate the traits that gene duplication imposes to mediate plant stress tolerance. Here, we noticed that B. rapa abiotic stress resistance (ASR) genes which code at least one stress responsive domain have a significantly higher number of paralogs than A. thaliana. Analysing the disordered content of the ASR genes in both species, we found that intrinsically disordered residues (IDR) are specifically enriched in whole genome duplication (WGD) derived paralogs. Subsequently, domain similarity analysis between WGD pairs of both species has revealed that majority of WGD pairs in B. rapa did not share domains with each other. Furthermore, domain enrichment analysis has shown that B. rapa paralogs contain 36 distinct stress responsive enriched domains, significantly higher than A. thaliana paralogs. Next, we performed MSA to investigate the domain conservation between orthologs and ohnologs pairs, we found that 80.13% of B. rapa ohnologs acquire new domains, depicting the fact that ohnologs play a significant role in stress-related behaviours. The average IDR content of the ohnologs enriching new domains after gene duplication in B. rapa (0.19), is also significantly higher than A. thaliana (0.04). Interestingly, we also found that all of these attributes i.e., exhibiting higher number of WGD paralogs and enhancement of IDR in ASR genes of B. rapa compared to A. thaliana is exclusive for ASR genes only. No such significant differences were observed in randomly selected non-ASR genes between the two species. Together these results provide strong support for the hypothesis that augmentation of IDR content followed by a whole genome duplication event imposes the stress resistance potentiality in B. rapa. This research will shed light on the mechanism of how B. rapa is able to successfully adapt to stress over the evolutionary timescale.

The curious case of SARM1: Dr. Jekyll and Mr. Hyde in cell death and immunity?

Sterile alpha and toll/interleukin-1 receptor motif-containing protein 1 (SARM1) was first identified as a novel ortholog of Drosophila protein CG7915 and was subsequently placed as the fifth member of the human TIR-containing adaptor protein. SARM1 holds a unique position in this family where, unlike other members, it downregulates NFκB activity in response to immunogenic stimulation, interacts with another member of the family, TRIF, to negatively regulate its function, and it also mediates cell death responses. Over the past decade, SARM1 has emerged as one of the primary mediators of programmed axonal degeneration and this robust regulation of axonal degeneration-especially in models of peripheral neuropathy and traumatic injury-makes it an attractive target for therapeutic intervention. The TIR domain of SARM1 possesses an intrinsic NADase activity resulting in cellular energy deficits within the axons, a striking deviation from its other family members of human TLR adaptors. Interestingly, the TIR NADase activity, as seen in SARM1, is also observed in several prokaryotic TIR-containing proteins where they are involved in immune evasion once within the host. Although the immune function of SARM1 is yet to be conclusively discerned, this closeness in function with the prokaryotic TIR-domain containing proteins, places it at an interesting juncture of evolution raising questions about its origin and function in cell death and immunity. In this review, we discuss how a conserved immune adaptor protein like SARM1 switches to a pro-neurodegenerative function and the evolutionarily significance of the process.

Genome-wide association study in two-row spring barley landraces identifies QTL associated with plantlets root system architecture traits in well-watered and osmotic stress conditions

Water availability is undoubtedly one of the most important environmental factors affecting crop production. Drought causes a gradual deprivation of water in the soil from top to deep layers and can occur at diverse stages of plant development. Roots are the first organs that perceive water deficit in soil and their adaptive development contributes to drought adaptation. Domestication has contributed to a bottleneck in genetic diversity. Wild species or landraces represent a pool of genetic diversity that has not been exploited yet in breeding program. In this study, we used a collection of 230 two-row spring barley landraces to detect phenotypic variation in root system plasticity in response to drought and to identify new quantitative trait loci (QTL) involved in root system architecture under diverse growth conditions. For this purpose, young seedlings grown for 21 days in pouches under control and osmotic-stress conditions were phenotyped and genotyped using the barley 50k iSelect SNP array, and genome-wide association studies (GWAS) were conducted using three different GWAS methods (MLM GAPIT, FarmCPU, and BLINK) to detect genotype/phenotype associations. In total, 276 significant marker-trait associations (MTAs; p-value (FDR)< 0.05) were identified for root (14 and 12 traits under osmotic-stress and control conditions, respectively) and for three shoot traits under both conditions. In total, 52 QTL (multi-trait or identified by at least two different GWAS approaches) were investigated to identify genes representing promising candidates with a role in root development and adaptation to drought stress.

Genomics and transcriptomics to protect rice (Oryza sativa. L.) from abiotic stressors: -pathways to achieving zero hunger

More over half of the world's population depends on rice as a major food crop. Rice (Oryza sativa L.) is vulnerable to abiotic challenges including drought, cold, and salinity since it grown in semi-aquatic, tropical, or subtropical settings. Abiotic stress resistance has bred into rice plants since the earliest rice cultivation techniques. Prior to the discovery of the genome, abiotic stress-related genes were identified using forward genetic methods, and abiotic stress-tolerant lines have developed using traditional breeding methods. Dynamic transcriptome expression represents the degree of gene expression in a specific cell, tissue, or organ of an individual organism at a specific point in its growth and development. Transcriptomics can reveal the expression at the entire genome level during stressful conditions from the entire transcriptional level, which can be helpful in understanding the intricate regulatory network relating to the stress tolerance and adaptability of plants. Rice (Oryza sativa L.) gene families found comparatively using the reference genome sequences of other plant species, allowing for genome-wide identification. Transcriptomics via gene expression profiling which have recently dominated by RNA-seq complements genomic techniques. The identification of numerous important qtl,s genes, promoter elements, transcription factors and miRNAs involved in rice response to abiotic stress was made possible by all of these genomic and transcriptomic techniques. The use of several genomes and transcriptome methodologies to comprehend rice (Oryza sativa, L.) ability to withstand abiotic stress have been discussed in this review.

Classification and Expression Profile of the U-Box E3 Ubiquitin Ligase Enzyme Gene Family in Maize (Zea mays L.)

The U-box E3 (PUB) family genes encode the E3 ubiquitin ligase enzyme, which determines substrate specific recognition during protein ubiquitination. They are widespread in plants and are critical for plant growth, development, and response to external stresses. However, there are few studies on the functional characteristic of PUB gene family in the important staple crop, maize (Zea mays L.). In this study, the PUB gene in maize was aimed to identify and classify through whole-genome screening. Phylogenetic tree, gene structure, conserved motif, chromosome location, gene duplication (GD), synteny, and cis-acting regulatory element of PUB member were analyzed. The expression profiles of ZmPUB gene family in maize during development and under abiotic stress and hormones treatment were analyzed by the RNA-seq data. A total of 79 PUB genes were identified in maize genome, and they were stratified into seven categories. There were 25 pairs of segmental duplications (SD) and 1 pair of tandem duplication (TD) identified in the maize PUB gene family. A close relationship was observed between the monocot plant maize and rice in PUB gene family. There were 94 kinds of cis-acting elements identified in the maize PUB gene family, which included 46 biotic- and abiotic-responsive elements, 19 hormone-responsive elements, 13 metabolic and growth-related elements. The expression profiles of maize PUB gene family showed characteristics of tissue specificity and response to abiotic stress and hormones treatment. These results provided an extensive overview of the maize PUB gene family.

Identifying mutations in sd1, Pi54 and Pi-ta, and positively selected genes of TN1, the first semidwarf rice in Green Revolution

BACKGROUND

Taichung Native 1 (TN1) is the first semidwarf rice cultivar that initiated the Green Revolution. As TN1 is a direct descendant of the Dee-geo-woo-gen cultivar, the source of the sd1 semidwarf gene, the sd1 gene can be defined through TN1. Also, TN1 is susceptible to the blast disease and is described as being drought-tolerant. However, genes related to these characteristics of TN1 are unknown. Our aim was to identify and characterize TN1 genes related to these traits.

RESULTS

Aligning the sd1 of TN1 to Nipponbare sd1, we found a 382-bp deletion including a frameshift mutation. Sanger sequencing validated this deleted region in sd1, and we proposed a model of the sd1 gene that corrects errors in the literature. We also predicted the blast disease resistant (R) genes of TN1. Orthologues of the R genes in Tetep, a well-known resistant cultivar that is commonly used as a donor for breeding new blast resistant cultivars, were then sought in TN1, and if they were present, we looked for mutations. The absence of Pi54, a well-known R gene, in TN1 partially explains why TN1 is more susceptible to blast than Tetep. We also scanned the TN1 genome using the PosiGene software and identified 11 genes deemed to have undergone positive selection. Some of them are associated with drought-resistance and stress response.

CONCLUSIONS

We have redefined the deletion of the sd1 gene in TN1, a direct descendant of the Dee-geo-woo-gen cultivar, and have corrected some literature errors. Moreover, we have identified blast resistant genes and positively selected genes, including genes that characterize TN1's blast susceptibility and abiotic stress response. These new findings increase the potential of using TN1 to breed new rice cultivars.

Dissection of Maize Drought Tolerance at the Flowering Stage Using Genome-Wide Association Studies

Drought is one of the most critical environmental factors constraining maize production. When it occurs at the flowering stage, serious yield losses are caused, and often, the damage is irretrievable. In this study, anthesis to silk interval (ASI), plant height (PH), and ear biomass at the silking date (EBM) of 279 inbred lines were studied under both water-stress (WS) and well-water (WW) field conditions, for three consecutive years. Averagely, ASI was extended by 25.96%, EBM was decreased by 17.54%, and the PH was reduced by 12.47% under drought stress. Genome-wide association studies were carried out using phenotypic values under WS, WW, and drought-tolerance index (WS-WW or WS/WW) and applying a mixed linear model that controls both population structure and relative kinship. In total, 71, 159, and 21 SNPs, located in 32, 59, and 12 genes, were significantly (P < 10−5) associated with ASI, EBM, and PH, respectively. Only a few overlapped candidate genes were found to be associated with the same drought-related traits under different environments, for example, ARABIDILLO 1, glycoprotein, Tic22-like, and zinc-finger family protein for ASI; 26S proteasome non-ATPase and pyridoxal phosphate transferase for EBM; 11-ß-hydroxysteroid dehydrogenase, uncharacterised, Leu-rich repeat protein kinase, and SF16 protein for PH. Furthermore, most candidate genes were revealed to be drought-responsive in an association panel. Meanwhile, the favourable alleles/key variations were identified with a haplotype analysis. These candidate genes and their key variations provide insight into the genetic basis of drought tolerance, especially for the female inflorescence, and will facilitate drought-tolerant maize breeding.

Meta-QTLs, ortho-MQTLs and candidate genes for the traits contributing to salinity stress tolerance in common wheat (Triticum aestivum L.)

A meta-analysis of QTLs associated with the traits contributing to salinity tolerance was undertaken in wheat to detect consensus and robust meta-QTLs (MQTLs) using 844 known QTLs retrieved from 26 earlier studies. A consensus map with a total length of 4621.56 cM including 7710 markers was constructed using 21 individual linkage maps and three previously published integrated genetic maps. Out of 844 QTLs, 571 QTLs were projected on the consensus map which gave origin to 100 MQTLs. Interestingly, 49 MQTLs were co-located with marker-trait associations reported in wheat genome-wide association studies for the traits contributing to salinity stress tolerance. Five potential MQTLs associated with the major salinity-responsive traits were also identified to be utilized in the breeding programme. In the resulted MQTLs, the average confidence interval (CI, 3.58 cM) was reduced up to 4.16 folds compared to the mean CI of the initial QTLs. Furthermore, as many as 617 gene models including 81 most likely candidate genes (CGs) were identified in the high confidence MQTL regions. These most likely CGs encoded proteins mainly belonging to the following families: B-box-type zinc finger, cytochrome P450 protein, pentatricopeptide repeat, phospholipid/glycerol acyltransferase, F-box protein, small auxin-up RNA, UDP-glucosyltransferase, glutathione S-transferase protein, etc. In addition, ortho-MQTL analysis based on synteny among wheat, rice and barley was also performed which permitted the identification of six ortho-MQTLs among these three cereals. This meta-analysis defines a genome-wide landscape on the most stable and consistent loci associated with reliable molecular markers and candidate genes for salinity tolerance in wheat.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01112-0.

OsPUB41, a U-box E3 ubiquitin ligase, acts as a negative regulator of drought stress response in rice (Oryza Sativa L.)

OsPUB41 plays a negative role in drought stress response through the mediation of OsUBC25 and interacts with OsCLC6, suggesting a putative substrate. The notable expansion of Plant U-Box E3 ligases (PUB), compared with those in mammals, implies that PUB proteins have evolved to perform plant-specific functions. OsPUB41, a potential ortholog of CMPG1, was recently reported to regulate the cell wall degrading enzyme (CWDE)-induced innate immune response in rice. Here, we characterized the OsPUB41 gene, which encodes a dual-localized cytosolic and nuclear U-box E3 ligase in rice. OsPUB41 expression was specifically induced by dehydration among various abiotic stresses and abscisic acid (ABA) treatments. Furthermore, we revealed that the core U-box motif of OsPUB41 possesses the E3 ligase activity that can be activated by OsUBC25 in rice. The Ubi:RNAi-OsPUB41 knock-down and ospub41 suppression mutant plants exhibited enhanced tolerance to drought stress compared with the wild-type rice plants in terms of transpirational water loss, long-term dehydration response, and chlorophyll content. Moreover, the knock-down or suppression of the OsPUB41 gene did not cause adverse effect on rice yield-related traits. Yeast two-hybrid and an in vitro pull-down analyses revealed that OsCLC6, a chloride channel, is a putative substrate of OsPUB41. Overall, these results suggest that OsPUB41 acts as a negative regulator of dehydration conditions and interacts with OsCLC6, implying that it is a substrate of OsPUB41.

High-density mapping for gray leaf spot resistance using two related tropical maize recombinant inbred line populations

Gray leaf spot (GLS) caused by Cercospora zeae-maydis or Cercospora zeina is one of the devastating maize foliar diseases worldwide. Identification of GLS-resistant quantitative trait loci (QTL)/genes plays an urgent role in improving GLS resistance in maize breeding practice. Two groups of recombinant inbred line (RIL) populations derived from CML373 × Ye107 and Chang7-2 × Ye107 were generated and subjected to genotyping-by-sequencing (GBS). A total of 1,929,222,287 reads in CML373 × Ye107 (RIL-YCML) and 2,585,728,312 reads in Chang7-2 × Ye107 (RIL-YChang), with an average of 10,961,490 (RIL-YCML) and 13,609,096 (RIL-YChang) reads per individual, were got, which was roughly equal to 0.70-fold and 0.87-fold coverage of the maize B73 RefGen_V4 genome for each F₇ individual, respectively. 6418 and 5139 SNP markers were extracted to construct two high-density genetic maps. Comparative analysis using these physically mapped marker loci demonstrated a satisfactory colinear relationship with the reference genome. 11 GLS-resistant QTL have been detected. The individual QTL accounted for 1.53-24.00% of the phenotypic variance explained (PVE). The new consensus QTL (qYCM-DS3-3/qYCM-LT3-1/qYCM-LT3-2) with the largest effect was located in chromosome bin 3.05, with an interval of 2.7 Mb, representing 13.08 to 24.00% of the PVE. Further gene annotation indicated that there were four candidate genes (GRMZM2G032384, GRMZM2G041415, GRMZM2G041544, and GRMZM2G035992) for qYCM-LT3-1, which may be related to GLS resistance. Combining RIL populations and GBS-based high-density genetic maps, a new larger effect QTL was delimited to a narrow genomic interval, which will provide a new resistance source for maize breeding programs.

Molecular Characterization Revealed the Role of Thaumatin-Like Proteins of Bread Wheat in Stress Response

Thaumatin-like proteins (TLPs) are related to pathogenesis-related-5 (PR-5) family and involved in stress response. Herein, a total of 93 TLP genes were identified in the genome of Triticum aestivum. Further, we identified 26, 27, 39, and 37 TLP genes in the Brachypodium distachyon, Oryza sativa, Sorghum bicolor, and Zea mays genomes for comparative characterization, respectively. They could be grouped into small and long TLPs with conserved thaumatin signature motif. Tightly clustered genes exhibited conserved gene and protein structure. The physicochemical analyses suggested significant differences between small and long TLPs. Evolutionary analyses suggested the role of duplication events and purifying selection in the expansion of the TLP gene family. Expression analyses revealed the possible roles of TLPs in plant development and abiotic and fungal stress response. Recombinant expression of TaTLP2-B in Saccharomyces cerevisiae provided significant tolerance against cold, heat, osmotic, and salt stresses. The results depicted the importance of TLPs in cereal crops that would be highly useful in future crop improvement programs.

Fine mapping and candidate gene prediction of a major quantitative trait locus for tassel branch number in maize

Tassel branch number (TBN) is the principal component of tassel inflorescence architecture in the maize plant. TBN is believed to be controlled by a set of quantitative trait loci (QTLs). However, it is necessary to identify and genetically evaluate these QTLs before the TBN can be improved upon using a molecular breeding approach. Therefore, in this study, we developed the chromosome segment introgression line (CSIL) TBN1 with the Zong31 (Z31) background and a higher TBN, and then we utilized the CSIL-TBN1-derived populations and identified a major QTL, qTBN6a, by linkage analysis. Fine mapping of the qTBN6a QTL was validated using a set of sub-CSILs and located in a 240-kb genomic region (Bin6.07) in B73RefGen_v4. One allele included in the introgression fragment had a positive effect, noticeably increasing the TBN and demonstrating the potential to improve the TBN of Z31. Afterward, in the qTBN6a interval, gene expression, sequence alignment, functional analysis, and the analysis of motifs in the 5' UTR suggested that candidate genes of qTBN6a are important functional genes at the early stage of immature infected tassel development. Among these candidate genes, a long W22::Mu-insertion/deletion in exon one and an 11-bp insertion/deletion in the promoter region may affect the variation of the qTBN6a QTL observed between Z31 and TBN1. In addition, the candidate genes of qTBN6a were found to encode a pentatricopeptide repeat (PPR)-containing protein and a histone deacetylase (HDA), which are known to be closely associated with RNA editing and stability and chromatin state activity for the transcription of gene expression, respectively. Finally, a model of qTBN6a based on the synergistic regulation of PPR and HDA for the maintenance of inflorescence meristem (IM) identity and its differentiation to the branch meristem (BM) in TBN1 was suggested. Collectively, our results provide an available locus for the molecular improvement of TBN and the isolation of functional genes underlying this QTL.

Understanding salt tolerance mechanism using transcriptome profiling and de novo assembly of wild tomato Solanum chilense

Soil salinity affects the plant growth and productivity detrimentally, but Solanum chilense, a wild relative of cultivated tomato (Solanum lycopersicum L.), is known to have exceptional salt tolerance. It has precise adaptations against direct exposure to salt stress conditions. Hence, a better understanding of the mechanism to salinity stress tolerance by S. chilense can be accomplished by comprehensive gene expression studies. In this study 1-month-old seedlings of S. chilense and S. lycopersicum were subjected to salinity stress through application of sodium chloride (NaCl) solution. Through RNA-sequencing here we have studied the differences in the gene expression patterns. A total of 386 million clean reads were obtained through RNAseq analysis using the Illumina HiSeq 2000 platform. Clean reads were further assembled de novo into a transcriptome dataset comprising of 514,747 unigenes with N50 length of 578 bp and were further aligned to the public databases. Genebank non-redundant (Nr), Viridiplantae, Gene Ontology (GO), KOG, and KEGG databases classification suggested enrichment of these unigenes in 30 GO categories, 26 KOG, and 127 pathways, respectively. Out of 265,158 genes that were differentially expressed in response to salt treatment, 134,566 and 130,592 genes were significantly up and down-regulated, respectively. Upon placing all the differentially expressed genes (DEG) in known signaling pathways, it was evident that most of the DEGs involved in cytokinin, ethylene, auxin, abscisic acid, gibberellin, and Ca²⁺ mediated signaling pathways were up-regulated. Furthermore, GO enrichment analysis was performed using REVIGO and up-regulation of multiple genes involved in various biological processes in chilense under salinity were identified. Through pathway analysis of DEGs, “Wnt signaling pathway” was identified as a novel pathway for the response to the salinity stress. Moreover, key genes for salinity tolerance, such as genes encoding proline and arginine metabolism, ROS scavenging system, transporters, osmotic regulation, defense and stress response, homeostasis and transcription factors were not only salt-induced but also showed higher expression in S. chilense as compared to S. lycopersicum. Thus indicating that these genes may have an important role in salinity tolerance in S. chilense. Overall, the results of this study improve our understanding on possible molecular mechanisms underlying salt tolerance in plants in general and tomato in particular.

The Genome of Microthlaspi erraticum (Brassicaceae) Provides Insights Into the Adaptation to Highly Calcareous Soils

Microthlaspi erraticum is widely distributed in temperate Eurasia, but restricted to Ca²⁺-rich habitats, predominantly on white Jurassic limestone, which is made up by calcium carbonate, with little other minerals. Thus, naturally occurring Microthlaspi erraticum individuals are confronted with a high concentration of Ca²⁺ ions while Mg²⁺ ion concentration is relatively low. As there is a competitive uptake between these two ions, adaptation to the soil condition can be expected. In this study, it was the aim to explore the genomic consequences of this adaptation by sequencing and analysing the genome of Microthlaspi erraticum. Its genome size is comparable with other diploid Brassicaceae, while more genes were predicted. Two Mg²⁺ transporters known to be expressed in roots were duplicated and one showed a significant degree of positive selection. It is speculated that this evolved due to the pressure to take up Mg²⁺ ions efficiently in the presence of an overwhelming amount of Ca²⁺ ions. Future studies on plants specialized on similar soils and affinity tests of the transporters are needed to provide unequivocal evidence for this hypothesis. If verified, the transporters found in this study might be useful for breeding Brassicaceae crops for higher yield on Ca²⁺-rich and Mg²⁺ -poor soils.

The Effect of Drought on Transcriptome and Hormonal Profiles in Barley Genotypes With Contrasting Drought Tolerance

Like many cereal crops, barley is also negatively affected by drought stress. However, due to its simple genome as well as enhanced stress resilient nature compared to rice and wheat, barley has been considered as a model to decipher drought tolerance in cereals. In the present study, transcriptomic and hormonal profiles along with several biochemical features were compared between drought-tolerant (Otis) and drought-sensitive (Baronesse) barley genotypes subjected to drought to identify molecular and biochemical differences between the genotypes. The drought-induced decrease in the leaf relative water content, net photosynthesis, and biomass accumulation was relatively low in Otis compared to Baronesse. The hormonal profiles did not reveal significant differences for majority of the compounds other than the GA20 and the cis-zeatin-o-glucoside (c-ZOG), whose levels were greatly increased in Otis compared to Baronesse under drought. The major differences that emerged from the transcriptome analysis are; (1), the overall number of differentially expressed genes was relatively low in drought-tolerant Otis compared to drought-sensitive Baronesse; (2), a wax biosynthesis gene (CER1), and NAC transcription factors were specifically induced in Otis but not in Baronesse; (3), the degree of upregulation of betaine aldehyde dehydrogenase and a homeobox transcription factor (genes with proven roles in imparting drought tolerance), was greater in Otis compared to Baronesse; (4) the extent of downregulation of gene expression profiles for proteins of the reaction center photosystem II (PSII) (D1 and D2) was low in Otis compared to Baronesse; and, (5), alternative splicing (AS) was also found to differ between the genotypes under drought. Taken together, the overall transcriptional responses were low in drought-tolerant Otis but the genes that could confer drought tolerance were either specifically induced or greatly upregulated in the tolerant genotype and these differences could be important for drought tolerance in barley.

Thaumatin-like protein kinases: Molecular characterization and transcriptional profiling in five cereal crops

Thaumatin-like protein kinases (TLPKs) are defense related proteins having antimicrobial property. Herein, we identified two TLPKs in the genome of Brachypodium distachyon and Oryza sativa, four in Hordeum vulgare and Sorghum bicolor, and 16 in Triticum aestivum. All the TLPKs were located at only one chromosome in each plant except T. aestivum, where they were located on chromosome 2 and chromosome 3. Paralogous analysis suggested the occurrence of one duplication event (DE) in each B. distachyon and O. sativa, two in H. vulgare while four DEs in T. aestivum genome during the evolution of TLPKs. The majority of TLPKs were intron less, while a few contains one or two introns. The introns were found in each 0, 1 and 2 phase. Protein structure analysis suggested the occurrence of a thaumatin and a kinase domain with a transmembrane (TM) helix in each TLPK. Further, a thaumatin family signature motif "GX[GF]XCXT[GA]DCX(1,2)GX(2,3)C", a "REDDD" motif and 16 cysteine residues were found conserved in the majority of TLPKs. Expression analysis indicated variable expression of TLPKs in various tissues of different cereal crops. They were high expressing in reproductive tissues in B. distachyon, while in leaves in T. aestivum. Modulated expression of TaTLPKs in the presence of fungal pathogen, and heat, drought and salt stress in T. aestivum suggested their roles in stress response. Co-expression analysis showed interaction of TLPKs with various development and stress related genes. The results indicated diverse roles of TLPKs, which can be utilized for the development of eco-friendly pest resistant crops in future.

Thaumatin-like protein kinases: Molecular characterization and transcriptional profiling in five cereal crops

Thaumatin-like protein kinases (TLPKs) are defense related proteins having antimicrobial property. Herein, we identified two TLPKs in the genome of Brachypodium distachyon and Oryza sativa, four in Hordeum vulgare and Sorghum bicolor, and 16 in Triticum aestivum. All the TLPKs were located at only one chromosome in each plant except T. aestivum, where they were located on chromosome 2 and chromosome 3. Paralogous analysis suggested the occurrence of one duplication event (DE) in each B. distachyon and O. sativa, two in H. vulgare while four DEs in T. aestivum genome during the evolution of TLPKs. The majority of TLPKs were intron less, while a few contains one or two introns. The introns were found in each 0, 1 and 2 phase. Protein structure analysis suggested the occurrence of a thaumatin and a kinase domain with a transmembrane (TM) helix in each TLPK. Further, a thaumatin family signature motif "GX[GF]XCXT[GA]DCX(1,2)GX(2,3)C", a "REDDD" motif and 16 cysteine residues were found conserved in the majority of TLPKs. Expression analysis indicated variable expression of TLPKs in various tissues of different cereal crops. They were high expressing in reproductive tissues in B. distachyon, while in leaves in T. aestivum. Modulated expression of TaTLPKs in the presence of fungal pathogen, and heat, drought and salt stress in T. aestivum suggested their roles in stress response. Co-expression analysis showed interaction of TLPKs with various development and stress related genes. The results indicated diverse roles of TLPKs, which can be utilized for the development of eco-friendly pest resistant crops in future.

Identification of microRNAs and relative target genes in Moringa oleifera leaf and callus

MicroRNAs, a class of small, non-coding RNAs, play important roles in plant growth, development and stress response by negatively regulating gene expression. Moringa oleifera Lam. plant has many medical and nutritional uses; however, little attention has been dedicated to its potential for the bio production of active compounds. In this study, 431 conserved and 392 novel microRNA families were identified and 9 novel small RNA libraries constructed from leaf, and cold stress treated callus, using high-throughput sequencing technology. Based on the M. oleifera genome, the microRNA repertoire of the seed was re-evaluated. qRT-PCR analysis confirmed the expression pattern of 11 conserved microRNAs in all groups. MicroRNA159 was found to be the most abundant conserved microRNA in leaf and callus, while microRNA393 was most abundantly expressed in the seed. The majority of predicted microRNA target genes were transcriptional factors involved in plant reproduction, growth/development and abiotic/biotic stress response. In conclusion, this is the first comprehensive analysis of microRNAs in M. oleifera leaf and callus which represents an important addition to the existing M. oleifera seed microRNA database and allows for possible exploitation of plant microRNAs induced with abiotic stress, as a tool for bio-enrichment with pharmacologically important phytochemicals.

Communications Between the Endoplasmic Reticulum and Other Organelles During Abiotic Stress Response in Plants

To adapt to constantly changing environmental conditions, plants have evolved sophisticated tolerance mechanisms to integrate various stress signals and to coordinate plant growth and development. It is well known that inter-organellar communications play important roles in maintaining cellular homeostasis in response to environmental stresses. The endoplasmic reticulum (ER), extending throughout the cytoplasm of eukaryotic cells, is a central organelle involved in lipid metabolism, Ca²⁺ homeostasis, and synthesis and folding of secretory and transmembrane proteins crucial to perceive and transduce environmental signals. The ER communicates with the nucleus via the highly conserved unfolded protein response pathway to mitigate ER stress. Importantly, recent studies have revealed that the dynamic ER network physically interacts with other intracellular organelles and endomembrane compartments, such as the Golgi complex, mitochondria, chloroplast, peroxisome, vacuole, and the plasma membrane, through multiple membrane contact sites between closely apposed organelles. In this review, we will discuss the signaling and metabolite exchanges between the ER and other organelles during abiotic stress responses in plants as well as the ER-organelle membrane contact sites and their associated tethering complexes.

Molecular and functional characterization of Arabidopsis thaliana VPNB1 gene involved in plant vascular development

Armadillo (ARM) repeat containing proteins constitute a large family in plants and are involved in diverse cellular functions, like signal transduction, proliferation and differentiation. In animals, ARM repeat proteins have been implicated in cancer development. In this study, we aimed in characterizing the VPNB1 gene from Arabidopsis thaliana and its role in plant development, by implementing a number of genetic and molecular approaches. AtVPNB1 encodes for an ARM repeat protein of unknown function, exclusively expressed in the cambium as well as in the differentiating xylem and phloem cells of the vascular system. Subcellular localization experiments showed that VPNB is confined in nucleoplasmic speckle-like structures unrelated to cajal bodies. Transgenic VPNB-impaired plants exhibit a slower growing phenotype and a non-canonical pattern of xylem tissue. On the contrary, VPNB overexpression lines display an inverted phenotype of increased growth, accompanied by an increased deposition of phloem and xylem cell layers. In line with the above data, qPCR analysis revealed a deregulation of several key master genes of secondary wall biosynthesis, underlining the involvement of VPNB1 in the regulation and differentiation of the root and shoot vascular tissue.

Association mapping by aerial drone reveals 213 genetic associations for Sorghum bicolor biomass traits under drought

Background

Sorghum bicolor is the fifth most commonly grown cereal worldwide and is remarkable for its drought and abiotic stress tolerance. For these reasons and the large size of biomass varieties, it has been proposed as a bioenergy crop. However, little is known about the genes underlying sorghum’s abiotic stress tolerance and biomass yield.

Results

To uncover the genetic basis of drought tolerance in sorghum at a genome-wide level, we undertook a high-density phenomics genome wide association study (GWAS) in which 648 diverse sorghum lines were phenotyped at two locations in California once per week by drone over the course of a growing season. Biomass, height, and leaf area were measured by drone for individual field plots, subjected to two drought treatments and a well-watered control. The resulting dataset of ~ 171,000 phenotypic data-points was analyzed along with 183,989 genotype by sequence markers to reveal 213 high-quality, replicated, and conserved GWAS associations.

Conclusions

The genomic intervals defined by the associations include many strong candidate genes, including those encoding heat shock proteins, antifreeze proteins, and other domains recognized as important to plant stress responses. The markers identified by our study can be used for marker assisted selection for drought tolerance and biomass. In addition, our results are a significant step toward identifying specific sorghum genes controlling drought tolerance and biomass yield.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5055-5) contains supplementary material, which is available to authorized users.

Identification of induced mutations in hexaploid wheat genome using exome capture assay

Wheat is a staple food crop of many countries. Improving resilience to biotic and abiotic stresses remain key breeding targets. Among these, rust diseases are the most detrimental in terms of depressing wheat production. In the present study, chemical mutagenesis was used to induce mutations in the wheat variety NN-Gandum-1. This cultivar is moderately resistant to leaf and yellow rust. The aim of mutagenesis was to improve resistance to the disease as well as to study function of genes conferring resistance to the disease. In the present investigation, a 0.8% EMS dose was found optimum for supporting 45-55% germination of NN-Gandum-1. A total of 3,634 M2 fertile plants were produced from each of the M1 plant. Out of these, 33 (0.91%) and 20 plants (0.55%) showed absolute resistance to leaf and yellow rust, respectively. While 126 (3.46%) and 127 plants (3.49%) exhibited high susceptibility to the leaf and yellow rust, respectively. In the M4 generation, a total of 11 M4 lines (nine absolute resistant and two highly susceptible) and one wild type were selected for NGS-based exome capture assay. A total of 104,779 SNPs were identified that were randomly distributed throughout the wheat sub genomes (A, B and D). Induced mutations in intronic sequences predominated. The highest total number of SNPs detected in this assay were mapped to chr.2B (14,273 SNPs), which contains the highest number of targeted base pairs in the assay. The average mutation density across all regions interrogated was estimated to be one mutation per 20.91 Mb. The highest mutation frequency was found in chr.2D (1/11.7 kb) and the lowest in chr.7D (1/353.4 kb). Out of the detected mutations, 101 SNPs were filtered using analysis criteria aimed to enrich for mutations that may affect gene function. Out of these, one putative SNP detected in Lr21 were selected for further analysis. The SNP identified in chimeric allele (Lr21) of a resistant mutant (N1-252) was located in a NBS domain of chr.1BS at 3.4 Mb position. Through computational analysis, it was demonstrated that this identified SNP causes a substitution of glutamic acid with alanine, resulting in a predicted altered protein structure. This mutation, therefore, is a candidate for contributing to the resistance phenotype in the mutant line. Based on this work, we conclude that the wheat mutant resource developed is useful as a source of novel genetic variation for forward-genetic screens and also as a useful tool for gaining insights into the important biological circuits of different traits of complex genomes like wheat.

Genome wide association study for gray leaf spot resistance in tropical maize core

Gray leaf spot is a maize foliar disease with worldwide distribution and can drastically reduce the production in susceptible genotypes. Published works indicate that resistance to gray leaf spot is a complex trait controlled by multiple genes, with additive effect and influenced by environment. The aim of this study was to identify genomic regions, including putative genes, associated with resistance to gray leaf spot under natural conditions of disease occurrence. A genome wide association study was conducted with 355,972 single nucleotide polymorphism markers on a phenotypic data composed by 157 tropical maize inbred lines, evaluated at Maringá -Brazil. Seven single nucleotide polymorphisms were significantly associated with gray leaf spot, some of which were localized to previously reported quantitative trait loci regions. Three gene models linked to the associated single nucleotide polymorphism were expressed at flowering time and tissue related with gray leaf spot infection, explaining a considerable proportion of the phenotypic variance, ranging from 0.34 to 0.38. The gene model GRMZM2G073465 (bin 10.07) encodes a cysteine protease3 protein, gene model GRMZM2G007188 (bin 1.02) expresses a rybosylation factor-like protein and the gene model GRMZM2G476902 (bin 4.08) encodes an armadillo repeat protein. These three proteins are related with plant defense pathway. Once these genes are validated in next studies, they will be useful for marker-assisted selection and can help improve the understanding of maize resistance to gray leaf spot.

Transcriptome assembly and annotation of johnsongrass (Sorghum halepense) rhizomes identify candidate rhizome-specific genes

Rhizomes facilitate the wintering and vegetative propagation of many perennial grasses. Sorghum halepense (johnsongrass) is an aggressive perennial grass that relies on a robust rhizome system to persist through winters and reproduce asexually from its rootstock nodes. This study aimed to sequence and assemble expressed transcripts within the johnsongrass rhizome. A de novo transcriptome assembly was generated from a single johnsongrass rhizome meristem tissue sample. A total of 141,176 probable protein-coding sequences from the assembly were identified and assigned gene ontology terms using Blast2GO. Estimated expression analysis and BLAST results were used to reduce the assembly to 64,447 high-confidence sequences. The johnsongrass assembly was compared to Sorghum bicolor, a related nonrhizomatous species, along with an assembly of similar rhizome tissue from the perennial grain crop Thinopyrum intermedium. The presence/absence analysis yielded a set of 98 expressed johnsongrass contigs that are likely associated with rhizome development.

Subfunctionalization of duplicate MYB genes in Solanum commersonii generated the cold-induced ScAN2 and the anthocyanin regulator ScAN1

Wild potato species are useful sources of allelic diversity and loci lacking in the cultivated potato. In these species, the presence of anthocyanins in leaves has been associated with a greater tolerance to cold stress. However, the molecular mechanisms that allow potatoes to withstand cold exposure remain unclear. Here, we show that the expression of AN2, a MYB transcription factor, is induced by low temperatures in wild, cold-tolerant Solanum commersonii, and not in susceptible Solanum tuberosum varieties. We found that AN2 is a paralog of the potato anthocyanin regulator AN1, showing similar interaction ability with basic helix-loop-helix (bHLH) co-partners. Their sequence diversity resulted in a different capacity to promote accumulation of phenolics when tested in tobacco. Indeed, functional studies demonstrated that AN2 is less able to induce anthocyanins than AN1, but nevertheless it has a strong ability to induce accumulation of hydroxycinnamic acid derivatives. We propose that the duplication of R2R3 MYB genes resulted in subsequent subfunctionalization, where AN1 specialized in anthocyanin production and AN2 conserved the ability to respond to cold stress, inducing mainly the synthesis of hydroxycinnamic acid derivatives. These results contribute to understanding the evolutionary significance of gene duplication on phenolic compound regulation.

Genomewide identification of genes involved in the potato response to drought indicates functional evolutionary conservation with Arabidopsis plants

Potato is one of the four most important food crop plants worldwide and is strongly affected by drought. The following two pairs of potato cultivars, which are related in ancestry but show different drought tolerances, were chosen for comparative gene expression studies: Gwiazda/Oberon and Tajfun/Owacja. Comparative RNA-seq analyses of gene expression differences in the transcriptomes obtained from drought-tolerant versus drought-sensitive plants during water shortage conditions were performed. The 23 top-ranking genes were selected, 22 of which are described here as novel potato drought-responsive genes. Moreover, all but one of the potato genes selected have homologues in the Arabidopsis genome. Of the seven tested A. thaliana mutants with altered expression of the selected homologous genes, compared to the wild-type Arabidopsis plants, six showed an improved tolerance to drought. These genes encode carbohydrate transporter, mitogen-activated protein kinase kinase kinase 15 (MAPKKK15), serine carboxypeptidase-like 19 protein (SCPL19), armadillo/beta-catenin-like repeat-containing protein, high-affinity nitrate transporter 2.7 and nonspecific lipid transfer protein type 2 (nsLPT). The evolutionary conservation of the functions of the selected genes in the plant response to drought confirms the importance of these identified potato genes in the ability of plants to cope with water shortage conditions. Knowledge regarding these gene functions can be used to generate potato cultivars that are resistant to unfavourable conditions. The approach used in this work and the obtained results allowed for the identification of new players in the plant response to drought.

Identification, characterization, expression profiling, and virus-induced gene silencing of armadillo repeat-containing proteins in tomato suggest their involvement in tomato leaf curl New Delhi virus resistance

Armadillo repeat family is well-characterized in several plant species for their involvement in multiple regulatory processes including growth, development, and stress response. We have previously shown a three-fold higher expression of ARM protein-encoding in tomato cultivar tolerant to tomato leaf curl New Delhi virus (ToLCNDV) compared to susceptible cultivar upon virus infection. This suggests the putative involvement of ARM proteins in defense response against virus infection; however, no comprehensive investigation has been performed to address this inference. In the present study, we have identified a total of 46 ARM-repeat proteins (SlARMs), and 41 U-box-containing proteins (SlPUBs) in tomato. These proteins and their corresponding genes were studied for their physicochemical properties, gene structure, domain architecture, chromosomal localization, phylogeny, and cis-regulatory elements in the upstream promoter region. Expression profiling of candidate genes in response to ToLCNDV infection in contrasting tomato cultivars showed significant upregulation of SlARM18 in the tolerant cultivar. Virus-induced gene silencing of SlARM18 in the tolerant tomato cultivar conferred susceptibility, which suggests the involvement of this gene in resistance mechanism. Further studies are underway to functionally characterize SlARM18 to delineate its precise role in defense mechanism.

Genome-wide analysis of DUF221 domain-containing gene family in Oryza species and identification of its salinity stress-responsive members in rice

DUF221 domain-containing genes (DDP genes) play important roles in developmental biology, hormone signalling transduction, and responses to abiotic stress. Therefore to understand their structural and evolutionary relationship, we did a genome-wide analysis of this important gene family in rice. Further, through comparative genomics, DDP genes from Oryza sativa subsp. (indica), nine different wild species of rice and Arabidopsis were also identified. We also found an expansion of the DDP gene families in rice and Arabidopsis which is due to the segmental duplication events in some of the gene family members. In general, a highly purifying selection was found acting on all the deduced paralogous and orthologous DDP gene pairs. The data from microarray and subsequent qRT-PCR analysis revealed that although several OsDDPs were differentially regulated under salinity stress, yet OsDDP6 was upregulated at all the developmental stages in salt tolerant rice genotype, FL478. Interestingly, OsDDP6 was found to be involved in proline metabolism pathway as indicated by protein network analysis. The diverse gene structures, varied transmembrane topologies and the differential expression patterns implied the functional diversity in DDP genes. Therefore, the comprehensive evolutionary analysis of DDP genes from different Oryza species and Arabidopsis performed in this study will provide the basis for further functional validation studies vis-à-vis DDP genes of rice and other plant species.

Metazoan evolution of the armadillo repeat superfamily

The superfamily of armadillo repeat proteins is a fascinating archetype of modular-binding proteins involved in various fundamental cellular processes, including cell-cell adhesion, cytoskeletal organization, nuclear import, and molecular signaling. Despite their diverse functions, they all share tandem armadillo (ARM) repeats, which stack together to form a conserved three-dimensional structure. This superhelical armadillo structure enables them to interact with distinct partners by wrapping around them. Despite the important functional roles of this superfamily, a comprehensive analysis of the composition, classification, and phylogeny of this protein superfamily has not been reported. Furthermore, relatively little is known about a subset of ARM proteins, and some of the current annotations of armadillo repeats are incomplete or incorrect, often due to high similarity with HEAT repeats. We identified the entire armadillo repeat superfamily repertoire in the human genome, annotated each armadillo repeat, and performed an extensive evolutionary analysis of the armadillo repeat proteins in both metazoan and premetazoan species. Phylogenetic analyses of the superfamily classified them into several discrete branches with members showing significant sequence homology, and often also related functions. Interestingly, the phylogenetic structure of the superfamily revealed that about 30 % of the members predate metazoans and represent an ancient subset, which is gradually evolving to acquire complex and highly diverse functions.

Homologous U-box E3 Ubiquitin Ligases OsPUB2 and OsPUB3 Are Involved in the Positive Regulation of Low Temperature Stress Response in Rice (Oryza sativa L.)

Rice U-box E3 Ub ligases (OsPUBs) are implicated in biotic stress responses. However, their cellular roles in response to abiotic stress are poorly understood. In this study, we performed functional analyses of two homologous OsPUB2 and OsPUB3 in response to cold stress (4°C). OsPUB2 was up-regulated by high salinity, drought, and cold, whereas OsPUB3 was constitutively expressed. A subcellular localization assay revealed that OsPUB2 and OsPUB3 were localized to the exocyst positive organelle (EXPO)-like punctate structures. OsPUB2 was also localized to the nuclei. OsPUB2 and OsPUB3 formed a hetero-dimeric complex as well as homo-dimers in yeast cells and in vitro. OsPUB2/OsPUB3 exhibited self-ubiquitination activities in vitro and were rapidly degraded in the cell-free extracts with apparent half-lives of 150-160 min. This rapid degradation of OsPUB2/OsPUB3 was delayed in the presence of the crude extracts of cold-treated seedlings (apparent half-lives of 200-280 min). Moreover, a hetero-dimeric form of OsPUB2/OsPUB3 was more stable than the homo-dimers. These results suggested that OsPUB2 and OsPUB3 function coordinately in response to cold stress. OsPUB2- and OsPUB3-overexpressing transgenic rice plants showed markedly better tolerance to cold stress than did the wild-type plants in terms of survival rates, chlorophyll content, ion leakage, and expression levels of cold stress-inducible marker genes. Taken together, these results suggested that the two homologous rice U-box E3 Ub ligases OsPUB2 and OsPUB3 are positive regulators of the response to cold stress.

Rice Improvement Through Genome-Based Functional Analysis and Molecular Breeding in India

Rice is one of the main pillars of food security in India. Its improvement for higher yield in sustainable agriculture system is also vital to provide energy and nutritional needs of growing world population, expected to reach more than 9 billion by 2050. The high quality genome sequence of rice has provided a rich resource to mine information about diversity of genes and alleles which can contribute to improvement of useful agronomic traits. Defining the function of each gene and regulatory element of rice remains a challenge for the rice community in the coming years. Subsequent to participation in IRGSP, India has continued to contribute in the areas of diversity analysis, transcriptomics, functional genomics, marker development, QTL mapping and molecular breeding, through national and multi-national research programs. These efforts have helped generate resources for rice improvement, some of which have already been deployed to mitigate loss due to environmental stress and pathogens. With renewed efforts, Indian researchers are making new strides, along with the international scientific community, in both basic research and realization of its translational impact.

Molecular Basis of Transient Neonatal Zinc Deficiency: NOVEL ZnT2 MUTATIONS DISRUPTING ZINC BINDING AND PERMEATION

A gradually increasing number of transient neonatal zinc deficiency (TNZD) cases was recently reported, all of which were associated with inactivating ZnT2 mutations. Here we characterized the impact of three novel heterozygous ZnT2 mutations G280R, T312M, and E355Q, which cause TNZD in exclusively breastfed infants of Japanese mothers. We used the bimolecular fluorescence complementation (BiFC) assay to provide direct visual evidence for the in situ dimerization of these ZnT2 mutants, and to explore their subcellular localization. Moreover, using three complementary functional assays, zinc accumulation using BiFC-Zinquin and Zinpyr-1 fluorescence as well as zinc toxicity assay, we determined the impact of these ZnT2 mutations on vesicular zinc accumulation. Although all three mutants formed homodimers with the wild type (WT) ZnT2 and retained substantial vesicular localization, as well as vesicular zinc accumulation, they had no dominant-negative effect over the WT ZnT2. Furthermore, using advanced bioinformatics, structural modeling, and site-directed mutagenesis we found that these mutations localized at key residues, which play an important physiological role in zinc coordination (G280R and E355Q) and zinc permeation (T312M). Collectively, our findings establish that some heterozygous loss of function ZnT2 mutations disrupt zinc binding and zinc permeation, thereby suggesting a haploinsufficiency state for the unaffected WT ZnT2 allele in TNZD pathogenesis. These results highlight the burning need for the development of a suitable genetic screen for the early diagnosis of TNZD to prevent morbidity.

Presence of unique glyoxalase III proteins in plants indicates the existence of shorter route for methylglyoxal detoxification

Glyoxalase pathway, comprising glyoxalase I (GLY I) and glyoxalase II (GLY II) enzymes, is the major pathway for detoxification of methylglyoxal (MG) into D-lactate involving reduced glutathione (GSH). However, in bacteria, glyoxalase III (GLY III) with DJ-1/PfpI domain(s) can do the same conversion in a single step without GSH. Our investigations for the presence of DJ-1/PfpI domain containing proteins in plants have indicated the existence of GLY III-like proteins in monocots, dicots, lycopods, gymnosperm and bryophytes. A deeper in silico analysis of rice genome identified twelve DJ-1 proteins encoded by six genes. Detailed analysis has been carried out including their chromosomal distribution, genomic architecture and localization. Transcript profiling under multiple stress conditions indicated strong induction of OsDJ-1 in response to exogenous MG. A member of OsDJ-1 family, OsDJ-1C, showed high constitutive expression at all developmental stages and tissues of rice. MG depletion study complemented by simultaneous formation of D-lactate proved OsDJ-1C to be a GLY III enzyme that converts MG directly into D-lactate in a GSH-independent manner. Site directed mutagenesis of Cys-119 to Alanine significantly reduces its GLY III activity indicating towards the existence of functional GLY III enzyme in rice—a shorter route for MG detoxification.

Genome-wide identification and characterization of the apple (Malus domestica) HECT ubiquitin-protein ligase family and expression analysis of their responsiveness to abiotic stresses

The ubiquitin-protein ligases (E3s) directly participate in ubiquitin (Ub) transferring to the target proteins in the ubiquitination pathway. The HECT ubiquitin-protein ligase (UPL), one type of E3s, is characterized as containing a conserved HECT domain of approximately 350 amino acids in the C terminus. Some UPLs were found to be involved in trichome development and leaf senescence in Arabidopsis. However, studies on plant UPLs, such as characteristics of the protein structure, predicted functional motifs of the HECT domain, and the regulatory expression of UPLs have all been limited. Here, we present genome-wide identification of the genes encoding UPLs (HECT gene) in apple. The 13 genes (named as MdUPL1-MdUPL13) from ten different chromosomes were divided into four groups by phylogenetic analysis. Among these groups, the encoding genes in the intron-exon structure and the included additional functional domains were quite different. Notably, the F-box domain was first found in MdUPL7 in plant UPLs. The HECT domain in different MdUPL groups also presented different spatial features and three types of conservative motifs were identified. The promoters of each MdUPL member carried multiple stress-response related elements by cis-acting element analysis. Experimental results demonstrated that the expressions of several MdUPLs were quite sensitive to cold-, drought-, and salt-stresses by qRT-PCR assay. The results of this study helped to elucidate the functions of HECT proteins, especially in Rosaceae plants.

Role of Protein Tyrosine Phosphatases in Plants

Reversible protein phosphorylation is a crucial regulatory mechanism that controls many biological processes in eukaryotes. In plants, phosphorylation events primarily occur on serine (Ser) and threonine (Thr) residues, while in certain cases, it was also discovered on tyrosine (Tyr) residues. In contrary to plants, extensive reports on Tyr phosphorylation regulating a large numbers of biological processes exist in animals. Despite of such prodigious function in animals, Tyr phosphorylation is a least studied mechanism of protein regulation in plants. Recently, various chemical analytical procedures have strengthened the view that Tyr phosphorylation is equally prevalent in plants as in animals. However, regardless of Tyr phosphorylation events occuring in plants, no evidence could be found for the existence of gene encoding for Tyr phosphorylation i.e. the typical Tyr kinases. Various methodologies have suggested that plant responses to stress signals and developmental processes involved modifications in protein Tyr phosphorylation. Correspondingly, various reports have established the role of PTPs (Protein Tyrosine Phosphatases) in the dephosphorylation and inactivation of mitogen activated protein kinases (MAPKs) hence, in the regulation of MAPK signaling cascade. Besides this, many dual specificity protein phosphatases (DSPs) are also known to bind starch and regulate starch metabolism through reversible phosphorylation. Here, we are emphasizing the significant progress on protein Tyr phosphatases to understand the role of these enzymes in the regulation of post-translational modification in plant physiology and development.

Genome-Wide Association Mapping of Fertility Reduction upon Heat Stress Reveals Developmental Stage-Specific QTLs in Arabidopsis thaliana

For crops that are grown for their fruits or seeds, elevated temperatures that occur during flowering and seed or fruit set have a stronger effect on yield than high temperatures during the vegetative stage. Even short-term exposure to heat can have a large impact on yield. In this study, we used Arabidopsis thaliana to study the effect of short-term heat exposure on flower and seed development. The impact of a single hot day (35°C) was determined in more than 250 natural accessions by measuring the lengths of the siliques along the main inflorescence. Two sensitive developmental stages were identified, one before anthesis, during male and female meiosis, and one after anthesis, during fertilization and early embryo development. In addition, we observed a correlation between flowering time and heat tolerance. Genome-wide association mapping revealed four quantitative trait loci (QTLs) strongly associated with the heat response. These QTLs were developmental stage specific, as different QTLs were detected before and after anthesis. For a number of QTLs, T-DNA insertion knockout lines could validate assigned candidate genes. Our findings show that the regulation of complex traits can be highly dependent on the developmental timing.

Expansion and Function of Repeat Domain Proteins During Stress and Development in Plants

The recurrent repeats having conserved stretches of amino acids exists across all domains of life. Subsequent repetition of single sequence motif and the number and length of the minimal repeating motifs are essential characteristics innate to these proteins. The proteins with tandem peptide repeats are essential for providing surface to mediate protein-protein interactions for fundamental biological functions. Plants are enriched in tandem repeat containing proteins typically distributed into various families. This has been assumed that the occurrence of multigene repeats families in plants enable them to cope up with adverse environmental conditions and allow them to rapidly acclimatize to these conditions. The evolution, structure, and function of repeat proteins have been studied in all kingdoms of life. The presence of repeat proteins is particularly profuse in multicellular organisms in comparison to prokaryotes. The precipitous expansion of repeat proteins in plants is presumed to be through internal tandem duplications. Several repeat protein gene families have been identified in plants. Such as Armadillo (ARM), Ankyrin (ANK), HEAT, Kelch-like repeats, Tetratricopeptide (TPR), Leucine rich repeats (LRR), WD40, and Pentatricopeptide repeats (PPR). The structure and functions of these repeat proteins have been extensively studied in plants suggesting a critical role of these repeating peptides in plant cell physiology, stress and development. In this review, we illustrate the structural, functional, and evolutionary prospects of prolific repeat proteins in plants.