Comprehensive analysis of cystatin family genes suggests their putative functions in sexual reproduction, embryogenesis, and seed formation

Mutation of Arabidopsis SME1 and Sm core assembly improves oxidative stress resilience

Alternative splicing is a key posttranscriptional gene regulatory process, acting in diverse adaptive and basal plant processes. Splicing of precursor-messenger RNA (pre-mRNA) is catalyzed by a dynamic ribonucleoprotein complex, designated the spliceosome. In a suppressor screen, we identified a nonsense mutation in the Smith (Sm) antigen protein SME1 to alleviate photorespiratory H₂O₂-dependent cell death in catalase deficient plants. Similar attenuation of cell death was observed upon chemical inhibition of the spliceosome, suggesting pre-mRNA splicing inhibition to be responsible for the observed cell death alleviation. Furthermore, the sme1-2 mutants showed increased tolerance to the reactive oxygen species inducing herbicide methyl viologen. Both an mRNA-seq and shotgun proteomic analysis in sme1-2 mutants displayed a constitutive molecular stress response, together with extensive alterations in pre-mRNA splicing of transcripts encoding metabolic enzymes and RNA binding proteins, even under unstressed conditions. Using SME1 as a bait to identify protein interactors, we provide experimental evidence for almost 50 homologs of the mammalian spliceosome-associated protein to reside in the Arabidopsis thaliana spliceosome complexes and propose roles in pre-mRNA splicing for four uncharacterized plant proteins. Furthermore, as for sme1-2, a mutant in the Sm core assembly protein ICLN resulted in a decreased sensitivity to methyl viologen. Taken together, these data show that both a perturbed Sm core composition and assembly results in the activation of a defense response and in enhanced resilience to oxidative stress.

Discovering novel reproductive genes in a non-model fly using de novo GridION transcriptomics

Gene discovery has important implications for investigating phenotypic trait evolution, adaptation, and speciation. Male reproductive tissues, such as accessory glands (AGs), are hotspots for recruitment of novel genes that diverge rapidly even among closely related species/populations. These genes synthesize seminal fluid proteins that often affect post-copulatory sexual selection-they can mediate male-male sperm competition, ejaculate-female interactions that modify female remating and even influence reproductive incompatibilities among diverging species/populations. Although de novo transcriptomics has facilitated gene discovery in non-model organisms, reproductive gene discovery is still challenging without a reference database as they are often novel and bear no homology to known proteins. Here, we use reference-free GridION long-read transcriptomics, from Oxford Nanopore Technologies (ONT), to discover novel AG genes and characterize their expression in the widespread dung fly, Sepsis punctum. Despite stark population differences in male reproductive traits (e.g.: Body size, testes size, and sperm length) as well as female re-mating, the male AG genes and their secretions of S. punctum are still unknown. We implement a de novo ONT transcriptome pipeline incorporating quality-filtering and rigorous error-correction procedures, and we evaluate gene sequence and gene expression results against high-quality Illumina short-read data. We discover highly-expressed reproductive genes in AG transcriptomes of S. punctum consisting of 40 high-quality and high-confidence ONT genes that cross-verify against Illumina genes, among which 26 are novel and specific to S. punctum. Novel genes account for an average of 81% of total gene expression and may be functionally relevant in seminal fluid protein production. For instance, 80% of genes encoding secretory proteins account for 74% total gene expression. In addition, median sequence similarities of ONT nucleotide and protein sequences match within-Illumina sequence similarities. Read-count based expression quantification in ONT is congruent with Illumina's Transcript per Million (TPM), both in overall pattern and within functional categories. Rapid genomic innovation followed by recruitment of de novo genes for high expression in S. punctum AG tissue, a pattern observed in other insects, could be a likely mechanism of evolution of these genes. The study also demonstrates the feasibility of adapting ONT transcriptomics for gene discovery in non-model systems.

Review: Unraveling the origin of the structural and functional diversity of plant cystatins

The regulation of protease activity is a critical factor for the physiological balance during plant growth and development. Among the proteins involved in controlling protease activity are the cystatins, well-described inhibitors of cysteine proteases present in viruses, bacteria and most Eukaryotes. Plant cystatins, commonly called phytocystatins, display unique structural and functional diversity and are classified according to their molecular weight as type-I, -II, and -III. Their gene structure is highly conserved across Viridiplantae and provides insights into their evolutionary relationships. Many type-I phytocystatins with introns share sequence similarities with type-II phytocystatins. New data shows that they could have originated from recent losses of the carboxy-terminal extension present in type-II phytocystatins. Intronless type-I phytocystatins originated from a single event shared by flowering plants. Pieces of evidence show multiple events of gene duplications, intron losses, and gains throughout the expansion and diversity of the phytocystatin family. Gene duplication events in Gymnosperms and Eudicots resulted in inhibitors with amino acid substitutions that may modify their interaction with target proteases and other proteins. This review brings a phylogenomic analysis of plant cystatin evolution and contributes to a broader understanding of their origins. A complete functional genomic analysis among phytocystatins and their roles in plant development and responses to abiotic and biotic stresses remains a question to be fully solved.

Armet from whitefly saliva acts as an effector to suppress plant defences by targeting tobacco cystatin

Arginine rich, mutated in early stage of tumours (Armet), is a well-characterized bifunctional protein as an unfolded protein response component intracellularly and a neurotrophic factor extracellularly in mammals. Recently, a new role of Armet as an effector protein mediating insect-plant interactions has been reported; however, its molecular mechanisms underlying the regulation of plant defences remain unclear. We investigated the molecular mechanisms underlying whitefly-secreted Armet-mediated regulation of insect-plant interaction by agrobacterium-mediated transient expression, RNA interference, electrical penetration graph, protein-protein interaction studies, virus-induced gene silencing assay, phytohormone analysis and whitefly bioassays. Armet, secreted by Bemisia tabaci whitefly, is highly expressed in the primary salivary gland and is delivered into tobacco plants during feeding. Overexpression of the BtArmet gene in tobacco enhanced whitefly performance, while silencing the BtArmet gene in whitefly interrupted whitefly feeding and suppressed whitefly performance on tobacco plants. BtArmet was shown to interact with NtCYS6, a cystatin protein essential for tobacco anti-whitefly resistance, and counteract the negative effects of NtCYS6 on whitefly. These results indicate that BtArmet is a salivary effector and acts to promote whitefly performance on tobacco plants through binding to the tobacco cystatin NtCYS6. Our findings provide novel insight into whitefly-plant interactions.

Genome-Wide Identification and Characterization of the Cystatin Gene Family in Bread Wheat (Triticum aestivum L.)

Cystatins, as reversible inhibitors of papain-like and legumain proteases, have been identified in several plant species. Although the cystatin family plays crucial roles in plant development and defense responses to various stresses, this family in wheat (Triticum aestivum L.) is still poorly understood. In this study, 55 wheat cystatins (TaCystatins) were identified. All TaCystatins were divided into three groups and both the conserved gene structures and peptide motifs were relatively conserved within each group. Homoeolog analysis suggested that both homoeolog retention percentage and gene duplications contributed to the abundance of the TaCystatin family. Analysis of duplication events confirmed that segmental duplications played an important role in the duplication patterns. The results of codon usage pattern analysis showed that TaCystatins had evident codon usage bias, which was mainly affected by mutation pressure. TaCystatins may be regulated by cis-acting elements, especially abscisic acid and methyl jasmonate responsive elements. In addition, the expression of all selected TaCystatins was significantly changed following viral infection and cold stress, suggesting potential roles in response to biotic and abiotic challenges. Overall, our work provides new insights into TaCystatins during wheat evolution and will help further research to decipher the roles of TaCystatins under diverse stress conditions.

In silico approach on sequential and structural variability in oryzacystatin and its interaction with cysteine protease enzymes of insect

Phytoprotease inhibitors (PI) are important defence compounds produced by plants against microbes and insect herbivory. Oryzacystatins (OCs) are a group of protease inhibitors from Oryza sativa L. that are specific against cysteine protease enzymes. This study revealed the evolutionary relationship of eleven different oryzacystatins from rice and their interaction with cysteine protease enzymes from brown planthopper (BPH) and striped stem borer (SSB). Three-dimensional structure of eleven different oryzacystatins and six cysteine protease enzymes were homologically modelled, and their interaction was analysed to explore the sequence heterogeneity, structural variability and functional significance. OC XI and OC V showed higher docking score and hydrogen bond interaction with all the six tested cysteine protease enzymes. N terminal glycine residue, central conserved QVVXG, C terminal AVVXXXPW regions are involved in interacting with the active site residue of protease enzymes. Substitution of N terminal glycine by any other residue in OC VI significantly reduced the interaction efficiency with cysteine proteases. In OC XI, glutamine in the fourth position of QVVXG showed higher interaction efficiency with all cysteine proteases than serine at the fourth amino acid position. N terminal glycine plays a vital role in OC XI for interacting with active sites of cysteine protease enzymes whereas, in the OC V central conserved region QVVSG and C terminal PW region plays a major part in the interaction. However, either N terminal or C terminal region along with the central conserved region of oryzacystatin is involved in mediating an efficient interaction with the active site residues of cysteine proteases. Molecular dynamic (MD) simulation study revealed the stability of the OC X1 - cathepsin O2 like and OC V - cathepsin F like complexes during a simulation for 20 ns. Insilico results of the present study predict the potential of oryzacystatins interaction with cysteine protease enzymes of insects.

CRISPR/Cas9 Guided Mutagenesis of Grain Size 3 Confers Increased Rice (Oryza sativa L.) Grain Length by Regulating Cysteine Proteinase Inhibitor and Ubiquitin-Related Proteins

Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)/CRISPR-associated protein (Cas9)-mediated genome editing has become an important way for molecular breeding in crop plants. To promote rice breeding, we edited the Grain Size 3 (GS3) gene for obtaining valuable and stable long-grain rice mutants. Furthermore, isobaric tags for the relative and absolute quantitation (iTRAQ)-based proteomic method were applied to determine the proteome-wide changes in the GS3 mutants compared with wild type (WT). Two target sites were designed to construct the vector, and the Agrobacterium-mediated method was used for rice transformation. Specific mutations were successfully introduced, and the grain length (GL) and 1000-grain weight (GWT) of the mutants were increased by 31.39% and 27.15%, respectively, compared with WT. The iTRAQ-based proteomic analysis revealed that a total of 31 proteins were differentially expressed in the GS3 mutants, including 20 up-regulated and 11 down-regulated proteins. Results showed that differentially expressed proteins (DEPs) were mainly related to cysteine synthase, cysteine proteinase inhibitor, vacuolar protein sorting-associated, ubiquitin, and DNA ligase. Furthermore, functional analysis revealed that DEPs were mostly enriched in cellular process, metabolic process, binding, transmembrane, structural, and catalytic activities. Pathway enrichment analysis revealed that DEPs were mainly involved in lipid metabolism and oxylipin biosynthesis. The protein-to-protein interaction (PPI) network found that proteins related to DNA damage-binding, ubiquitin-40S ribosomal, and cysteine proteinase inhibitor showed a higher degree of interaction. The homozygous mutant lines featured by stable inheritance and long-grain phenotype were obtained using the CRISPR/Cas9 system. This study provides a convenient and effective way of improving grain yield, which could significantly accelerate the breeding process of long-grain japonica parents and promote the development of high-yielding rice.

Genome-wide identification and analysis of cystatin family genes in Sorghum (Sorghum bicolor (L.) Moench)

To set a systematic study of the Sorghum cystatins (SbCys) gene family, a genome-wide analysis of the SbCys family genes was performed by bioinformatics-based methods. In total, 18 SbCys genes were identified in Sorghum, which were distributed unevenly on chromosomes, and two genes were involved in a tandem duplication event. All SbCys genes had similar exon/intron structure and motifs, indicating their high evolutionary conservation. Transcriptome analysis showed that 16 SbCys genes were expressed in different tissues, and most genes displayed higher expression levels in reproductive tissues than in vegetative tissues, indicating that the SbCys genes participated in the regulation of seed formation. Furthermore, the expression profiles of the SbCys genes revealed that seven cystatin family genes were induced during Bipolaris sorghicola infection and only two genes were responsive to aphid infestation. In addition, quantitative real-time polymerase chain reaction (qRT-PCR) confirmed that 17 SbCys genes were induced by one or two abiotic stresses (dehydration, salt, and ABA stresses). The interaction network indicated that SbCys proteins were associated with several biological processes, including seed development and stress responses. Notably, the expression of SbCys4 was up-regulated under biotic and abiotic stresses, suggesting its potential roles in mediating the responses of Sorghum to adverse environmental impact. Our results provide new insights into the structural and functional characteristics of the SbCys gene family, which lay the foundation for better understanding the roles and regulatory mechanism of Sorghum cystatins in seed development and responses to different stress conditions.

A new Piper nigrum cysteine proteinase inhibitor, PnCPI, with antifungal activity: molecular cloning, recombinant expression, functional analyses and molecular modeling

A new Piper nigrum cysteine proteinase inhibitor, PnCPI, belonging to group I of phytocystatins, with inhibitory activity against papain and growth of Fusarium solani f. sp. piperis, was isolated and characterized. Previous studies (de Souza et al. 2011) have identified a partial cDNA sequence of putative cysteine proteinase inhibitor differentially expressed in roots of black pepper (P. nigrum L.) infected by F. solani f. sp. piperis. Here, we aimed to isolate the full-length cDNA and genomic sequences of the P. nigrum cysteine proteinase inhibitor gene, named PnCPI. Sequence analyses showed that the PnCPI gene encodes a deduced protein of 108 amino acid residues with a predicted molecular mass of 12.3 kDa and isoelectric point of 6.51. Besides the LARFAV-like sequence, common to all phytocystatins, PnCPI contains three conserved motifs of the superfamily cystatin: a glycine residue at the N-terminal region, the QxVxG reactive site more centrally positioned, and one tryptophan in the C-terminal region. PnCPI, belonging to group I of phytocystatins, showed high identity with cystatins isolated from several plant species. Sequence analyses also revealed no putative signal peptide at the N-terminal of PnCPI, as well as no introns within the genomic sequence corresponding to the PnCPI coding region. Molecular modeling showed the ability of PnCPI to interact with papain, while its inhibitory activity against this protease was confirmed after heterologous expression in Escherichia coli. The effects of heat treatments on the inhibitory activity of recombinant PnCPI, rPnCPI, were evaluated. In addition, rPnCPI exhibited in vitro activity against F. solani f. sp. piperis, revealing a new cystatin with the potential antifungal application. The identification of PnCPI as a functional cystatin able to inhibit the in vitro growth of F. solani f. sp. piperis indicates other factors contributing to in vivo susceptibility of black pepper to root rot disease.

Parasite Cystatin: Immunomodulatory Molecule with Therapeutic Activity against Immune Mediated Disorders

The use of parasites or their products for treating chronic inflammation associated diseases (CIADs) has generated significant attention recently. Findings from basic and clinical research have provided valuable information on strengthening the notion that parasites' molecules can be developed as biotherapeutic agents. Completion of the genome, secreotome, and proteome of the parasites has provided an excellent platform for screening and identifying several host immunomodulatory molecules from the parasites and evaluate their therapeutic potential for CIADs. One of the widely studied host immunomodulatory molecules of the parasites is the cysteine protease inhibitor (cystatin), which is primarily secreted by the parasites to evade host immune responses. In this review, we have attempted to summarize the findings to date on the use of helminth parasite-derived cystatin as a therapeutic agent against CIADs. Although several studies suggest a role for alternatively activated macrophages, other regulatory cells, and immunosuppressive molecules, in this immunoregulatory activity of the parasite-derived cystatin, there is still no clear demonstration as to how cystatin induces its anti-inflammatory effect in suppressing CIADs.

Recombinant cystatins in plants

Dozens of studies have assessed the practical value of plant cystatins as ectopic inhibitors of Cys proteases in biological systems. The potential of these proteins in crop protection to control herbivorous pests and pathogens has been documented extensively over the past 25 years. Their usefulness to regulate endogenous Cys proteases in planta has also been considered recently, notably to implement novel traits of agronomic relevance in crops or to generate protease activity-depleted environments in plants or plant cells used as bioreactors for recombinant proteins. After a brief update on the basic structural characteristics of plant cystatins, we summarize recent advances on the use of these proteins in plant biotechnology. Attention is also paid to the molecular improvement of their structural properties for the improvement of their protease inhibitory effects or the fine-tuning of their biological target range.

Nodule-specific PLAT domain proteins are expanded in the Medicago lineage and required for nodulation

Symbiotic nitrogen fixation in legumes is mediated by an interplay of signaling processes between plant hosts and rhizobial symbionts. In legumes, several secreted protein families have undergone expansions and play key roles in nodulation. Thus, identifying lineage-specific expansions (LSEs) of nodulation-associated genes can be a strategy to discover candidate gene families. Using bioinformatic tools, we identified 13 LSEs of nodulation-related secreted protein families, each unique to either Glycine, Arachis or Medicago lineages. In the Medicago lineage, nodule-specific Polycystin-1, Lipoxygenase, Alpha Toxin (PLAT) domain proteins (NPDs) expanded to five members. We examined NPD function using CRISPR/Cas9 multiplex genome editing to create Medicago truncatula NPD knockout lines, targeting one to five NPD genes. Mutant lines with differing combinations of NPD gene inactivations had progressively smaller nodules, earlier onset of nodule senescence, or ineffective nodules compared to the wild-type control. Double- and triple-knockout lines showed dissimilar nodulation phenotypes but coincided in upregulation of a DHHC-type zinc finger and an aspartyl protease gene, possible candidates for the observed disturbance of proper nodule function. By postulating that gene family expansions can be used to detect candidate genes, we identified a family of nodule-specific PLAT domain proteins and confirmed that they play a role in successful nodule formation.

Plant proteases during developmental programmed cell death

Proteases are among the key regulators of most forms of programmed cell death (PCD) in animals. Many PCD processes have also been associated with protease expression or activation in plants, However, functional evidence for the roles and actual modes of action of plant proteases in PCD remains surprisingly limited. In this review, we provide an update on protease involvement in the context of developmentally regulated plant PCD. To illustrate the diversity of protease functions, we focus on several prominent developmental PCD processes, including xylem and tapetum maturation, suspensor elimination, endosperm degradation, and seed coat formation, as well as plant senescence processes. Despite the substantial advances in the field, protease functions are often only correlatively linked to developmental PCD, and the specific molecular roles of proteases in many developmental PCD processes remain to be elucidated.

Comparative Analysis of WUSCHEL-Related Homeobox Genes Revealed Their Parent-of-Origin and Cell Type-Specific Expression Pattern During Early Embryogenesis in Tobacco

WUSCHEL-related homeobox (WOX) gene is a plant-specific clade of homeobox transcription factors. Increasing evidences reveal that WOXs play critical roles in early embryogenesis, which involves zygote development, initiation of zygote division, and apical or basal cell lineage establishment. However, how WOXs regulate these developmental events remains largely unknown, and even detailed expression pattern in gametes and early proembryos is not yet available. Here, 13 WOX family genes were identified in Nicotiana tabacum genome. Comparative analysis of 13 WOX family genes with their homologs in Arabidopsis thaliana reveals relatively conserved expression pattern of WUS and WOX5 in shoot/root apical meristem. Whereas variations were also found, e.g., lacking homolog of WOX8 (a marker for suspensor cell) in tobacco genome and the expression of WOX2/WOX9 in both apical cell and basal cell. Transient transcriptional activity analysis revealed that WOXs in WUS clade have repressive activities for their target's transcription, whereas WOXs in ancient and intermediate clade have activation activities, giving a molecular basis for the phylogenetic classification of tobacco WOXs into three major clades. Expression pattern analysis revealed that some WOXs (e.g., WOX 13a) expressed in both male and female gametes and some WOXs (e.g., WOX 11 and WOX 13b) displayed the characteristics of parent-of-origin genes. Interestingly, some WOXs (e.g., WOX2 and WOX9), which are essential for early embryo patterning, were de novo transcribed in zygote, indicating relevant mechanism for embryo pattern formation is only established in zygote right after fertilization and not carried in by gametes. We also found that most WOXs displayed a stage-specific and cell type-specific expression pattern. Taken together, this work provides a detailed landscape of WOXs in tobacco during fertilization and early embryogenesis, which will facilitate the understanding of their specific roles in these critical developmental processes of embryogenesis.

Ribosomal protein NtRPL17 interacts with kinesin-12 family protein NtKRP and functions in the regulation of embryo/seed size and radicle growth

We previously reported that a novel motor protein belonging to the kinesin-12 family, NtKRP, displays critical roles in regulating embryo and seed size establishment. However, it remains unknown exactly how NtKRP contributes to this developmental process. Here, we report that a 60S ribosomal protein NtRPL17 directly interacts with NtKRP. The phenotypes of NtRPL17 RNAi lines show notable embryo and seed size reduction. Structural observations of the NtRPL17-silenced embryos/seeds reveal that the embryo size reduction is due to a decrease in cell number. In these embryos, cell division cycle progression is delayed at the G2/M transition. These phenotypes are similar to that in NtKRP-silenced embryos/seeds, indicating that NtKRP and NtRPL17 function as partners in the same regulatory pathway during seed development and specifically regulate cell cycle progression to control embryo/seed size. This work reveals that NtRPL17, as a widely distributed ribosomal protein, plays a critical role in seed development and provides a new clue in the regulation of seed size. Confirmation of the interaction between NtKRP and NtRPL17 and their co-function in the control of the cell cycle also suggests that the mechanism might be conserved in both plants and animals.

Defense-related proteins involved in sugarcane responses to biotic stress

Sugarcane is one of the most important agricultural crops in the world. However, pathogen infection and herbivore attack cause constant losses in yield. Plants respond to pathogen infection by inducing the expression of several protein types, such as glucanases, chitinases, thaumatins, peptidase inhibitors, defensins, catalases and glycoproteins. Proteins induced by pathogenesis are directly or indirectly involved in plant defense, leading to pathogen death or inducing other plant defense responses. Several of these proteins are induced in sugarcane by different pathogens or insects and have antifungal or insecticidal activity. In this review, defense-related proteins in sugarcane are described, with their putative mechanisms of action, pathogen targets and biotechnological perspectives.

Search for Nodulation and Nodule Development-Related Cystatin Genes in the Genome of Soybean (Glycine max)

Nodulation, nodule development and senescence directly affects nitrogen fixation efficiency, and previous studies have shown that inhibition of some cysteine proteases delay nodule senescence, so their nature inhibitors, cystatin genes, are very important in nodulation, nodule development, and senescence. Although several cystatins are actively transcribed in soybean nodules, their exact roles and functional diversities in legume have not been well explored in genome-wide survey studies. In this report, we performed a genome-wide survey of cystatin family genes to explore their relationship to nodulation and nodule development in soybean and identified 20 cystatin genes that encode peptides with 97-245 amino acid residues, different isoelectric points (pI) and structure characteristics, and various putative plant regulatory elements in 3000 bp putative promoter fragments upstream of the 20 soybean cystatins in response to different abiotic/biotic stresses, hormone signals, and symbiosis signals. The expression profiles of these cystatin genes in soybean symbiosis with rhizobium strain Bradyrhizobium japonicum strain 113-2 revealed that 7 cystatin family genes play different roles in nodulation as well as nodule development and senescence. However, these genes were not root nodule symbiosis (RNS)-specific and did not encode special clade cystatin protein with structures related to nodulation and nodule development. Besides, only two of these soybean cystatins were not upregulated in symbiosis after ABA treatment. The functional analysis showed that a candidate gene Glyma.15G227500 (GmCYS16) was likely to play a positive role in soybean nodulation. Besides, evolutionary relationships analysis divided the cystatin genes from Arabidopsis thaliana, Nicotiana tabacum, rice, barley and four legume plants into three groups. Interestingly, Group A cystatins are special in legume plants, but only include one of the above-mentioned 7 cystatin genes related to nodulation and nodule development. Overall, our results provide useful information or clues for our understanding of the functional diversity of legume cystatin family proteins in soybean nodulation and nodule development and for finding nodule-specific cysteine proteases in soybean.

Programmed cell death in seeds of angiosperms

During the diversification of angiosperms, seeds have evolved structural, chemical, molecular and physiologically developing changes that specially affect the nucellus and endosperm. All through seed evolution, programmed cell death (PCD) has played a fundamental role. However, examples of PCD during seed development are limited. The present review examines PCD in integuments, nucellus, suspensor and endosperm in those representative examples of seeds studied to date.

Genome-wide identification and characterization of cystatin family genes in rice (Oryza sativa L.)

11 Cystatin genes in rice were identified, and their expression patterns were comprehensively analyzed, which reveals multiple roles in both seed development and plant response to environmental variations. Cystatin is a group of small proteins and known to inhibit the activities of cysteine proteases in the papain C1A and legumain C13 peptidase families in plants. Cystatin family genes have only been well characterized recently in a few plant species such as Hordeum vulgare and Nicotiana tabacum, which show their critical roles in programmed cell death and responses to biotic stresses. Up to now, little is known about cystatin family genes and their roles in Oryza sativa, a model plant for cereal biology study. Here, we identified 11 cystatin genes in rice genome. Comprehensive expression profile analysis reveals that cystatin family genes in rice display diverse expression pattern. They are temporally regulated at different developmental stages during the process of seed production and germination. Our experiments also reveal that the majority of cystatin genes are responsive to plant hormones and different environmental cues including cold, drought and other abiotic stresses, while some others are very stable under different stresses, indicating their fundamental roles in normal plant development. In addition, their distribution in rice chromosomes and their evolutionary relation to the members of Cystatin family in A. thaliana and N. tabacum have also been analyzed. These works suggest multiple roles of cystatin family genes in both seed development and plant response to environmental variations.

A comprehensive, genome-wide analysis of autophagy-related genes identified in tobacco suggests a central role of autophagy in plant response to various environmental cues

Autophagy is an evolutionarily conserved mechanism in both animals and plants, which has been shown to be involved in various essential developmental processes in plants. Nicotiana tabacum is considered to be an ideal model plant and has been widely used for the study of the roles of autophagy in the processes of plant development and in the response to various stresses. However, only a few autophagy-related genes (ATGs) have been identified in tobacco up to now. Here, we identified 30 ATGs belonging to 16 different groups in tobacco through a genome-wide survey. Comprehensive expression profile analysis reveals an abroad expression pattern of these ATGs, which could be detected in all tissues tested under normal growth conditions. Our series tests further reveal that majority of ATGs are sensitive and responsive to different stresses including nutrient starvation, plant hormones, heavy metal and other abiotic stresses, suggesting a central role of autophagy, likely as an effector, in plant response to various environmental cues. This work offers a detailed survey of all ATGs in tobacco and also suggests manifold functions of autophagy in both normal plant growth and plant response to environmental stresses.