Comprehensive sequence and expression profile analysis of PEX11 gene family in rice

Inferring Historical Introgression with Deep Learning

Resolving phylogenetic relationships among taxa remains a challenge in the era of big data due to the presence of genetic admixture in a wide range of organisms. Rapidly developing sequencing technologies and statistical tests enable evolutionary relationships to be disentangled at a genome-wide level, yet many of these tests are computationally intensive and rely on phased genotypes, large sample sizes, restricted phylogenetic topologies, or hypothesis testing. To overcome these difficulties, we developed a deep learning-based approach, named ERICA, for inferring genome-wide evolutionary relationships and local introgressed regions from sequence data. ERICA accepts sequence alignments of both population genomic data and multiple genome assemblies, and efficiently identifies discordant genealogy patterns and exchanged regions across genomes when compared with other methods. We further tested ERICA using real population genomic data from Heliconius butterflies that have undergone adaptive radiation and frequent hybridization. Finally, we applied ERICA to characterize hybridization and introgression in wild and cultivated rice, revealing the important role of introgression in rice domestication and adaptation. Taken together, our findings demonstrate that ERICA provides an effective method for teasing apart evolutionary relationships using whole genome data, which can ultimately facilitate evolutionary studies on hybridization and introgression.

E3 ligase, the Oryza sativa salt-induced RING finger protein 4 (OsSIRP4), negatively regulates salt stress responses via degradation of the OsPEX11-1 protein

OsSIRP4 is an E3 ligase that acts as a negative regulator in the plant response to salt stress via the 26S proteasomal system regulation of substrate proteins, OsPEX11-1, which it provides important information for adaptation and regulation in rice. Plants are sessile organisms that can be exposed to environmental stress. Plants alter their cellular processes to survive under potentially unfavorable conditions. Protein ubiquitination is an important post-translational modification that has a crucial role in various cellular signaling processes in abiotic stress response. In this study, we characterized Oryza sativa salt-induced RING finger protein 4, OsSIRP4, a membrane and cytosol-localized RING E3 ligase in rice. OsSIRP4 transcripts were highly induced under salt stress in rice. We found that OsSIRP4 possesses E3 ligase activity; however, no E3 ligase activity was observed with a single amino acid substitution (OsSIRP4^C269A). The results of the yeast two hybrid system, in vitro pull-down assay, BiFC analysis, in vitro ubiquitination assay, and in vitro degradation assay indicate that OsSIRP4 regulates degradation of a substrate protein, OsPEX11-1 (Oryza sativa peroxisomal biogenesis factor 11-1) via the 26S proteasomal system. Phenotypic analysis of OsSIRP4-overexpressing plants demonstrated hypersensitivity to salt response compared to that of the wild type and mutated OsSIRP4^C269A plants. In addition, OsSIRP4-overexpressing plants exhibited significant low enzyme activities of superoxide dismutase, catalase, and peroxidase, and accumulation of proline and soluble sugar, but a high level of H₂O₂. Furthermore, qRT data on transgenic plants suggest that OsSIRP4 acted as a negative regulator of salt response by diminishing the expression of genes related to Na⁺/K⁺ homeostasis (AtSOS1, AtAKT1, AtNHX1, and AtHKT1;1) in transgenic plants under salt stress. These results suggest that OsSIRP4 plays a negative regulatory role in response to salt stress by modulating the target protein levels.

Differential Dynamic Changes of Reduced Trait Model for Analyzing the Plastic Response to Drought Phases: A Case Study in Spring Wheat

Current limited water availability due to climate changes results in severe drought stress and desiccation in plants. Phenotyping drought tolerance remains challenging. In particular, our knowledge about the discriminating power of traits for capturing a plastic phenotype in high-throughput settings is scant. The study is designed to investigate the differential performance and broad-sense heritability of a battery set of morphological, physiological, and cellular traits to understand the adaptive phenotypic response to drought in spring wheat during the tillering stage. The potential of peroxisome abundance to predict the adaptive response under severe drought was assessed using a high-throughput technique for peroxisome quantification in plants. The research dissected the dynamic changes of some phenological traits during three successive phases of drought using two contrasting genotypes of adaptability to drought. The research demonstrates 5 main findings: (1) a reduction of the overall dimension of the phenological traits for robust phenotyping of the adaptive performance under drought; (2) the abundance of peroxisomes in response to drought correlate negatively with grain yield; (3) the efficiency of ROS homeostasis through peroxisome proliferation which seems to be genetically programmed; and (4) the dynamics of ROS homeostasis seems to be timing dependent mechanism, the tolerant genotype response is earlier than the susceptible genotype. This work will contribute to the identification of robust plastic phenotypic tools and the understanding of the mechanisms for adaptive behavior under drought conditions.

SUMMARY STATEMENT

This study presents the estimated broad-sense heritability of 24 phenological traits under drought compared with non-stressed conditions. The results demonstrated a reduced model of the overall dimension of the phenological traits for phenotyping drought tolerant response including a novel trait (peroxisome abundance). Also, it displays that the adaptive mechanism through peroxisomes proliferation that is a genetic-dependent manner and related to the stress phase, since tolerant plants can sense the stress and maintain the cellular balance earlier than the sensitive plants.

Using in silico techniques: Isolation and characterization of an insect cuticle-degrading-protease gene from Beauveria bassiana

Cuticle-degrading-proteases (CDPs) secreted by Beauveria spp. are pivotal biocontrol substances, possessing commercial potential for developing bio-pesticides. Therefore, a thoughtful and contemplative understanding and assessment of the structural and functional features of these proteases would markedly assist the development of biogenic pesticides. Computational molecular biology is a new facile alternative approach to the tedious experimental molecular biology; therefore, by using bioinformatics tools, we isolated and characterized an insect CDP gene from Beauveria bassiana 70 s.l. genomic DNA. The CDP gene (1240 bp with GeneBank accession no. KT804651.1) consisted of three introns and four CDS exons, and shared 74-100% sequence identity to the reference CDP genes. Its phylogenetic tree results showed a unique evolution pattern, and the predicted amino acid peptide (PAAP) consisted of 344 amino acid residues with pI, molecular weight, instability index, grand average hydropathicity value and aliphatic index of 7.2, 35.4 kDa, 24.45, -0.149, and 76.63, respectively. The gene possessed 74-89% amino acid sequence similarity to the 12 reference strains. Three motifs (Peptidase_S8 subtilase family) were detected in the PAAP, and the computed 3D structure possessed 79.09% structural identity to alkaline serine proteases. The PAAP had four (three serine proteases and one Pyridoxal-dependent decarboxylase) conserved domains, a disulfide bridge, two calcium binding sites, MY domain, and three predicted active sites in the serine family domains. These results will set the groundwork for further exploitation of proteases and understanding the mechanism of disease caused by cuticle-degrading-serine-proteases from entomopathogenic fungi.

Giant peroxisomes in a moss (Physcomitrella patens) peroxisomal biogenesis factor 11 mutant

Peroxisomal biogenesis factor 11 (PEX11) proteins are found in yeasts, mammals and plants, and play a role in peroxisome morphology and regulation of peroxisome division. The moss Physcomitrella patens has six PEX11 isoforms which fall into two subfamilies, similar to those found in monocots and dicots. We carried out targeted gene disruption of the Phypa_PEX11-1 gene and compared the morphological and cellular phenotypes of the wild-type and mutant strains. The mutant grew more slowly and the development of gametophores was retarded. Mutant chloronemal filaments contained large cellular structures which excluded all other cellular organelles. Expression of fluorescent reporter proteins revealed that the mutant strain had greatly enlarged peroxisomes up to 10 μm in diameter. Expression of a vacuolar membrane marker confirmed that the enlarged structures were not vacuoles, or peroxisomes sequestered within vacuoles as a result of pexophagy. Phypa_PEX11 targeted to peroxisome membranes could rescue the knock out phenotype and interacted with Fission1 on the peroxisome membrane. Moss PEX11 functions in peroxisome division similar to PEX11 in other organisms but the mutant phenotype is more extreme and environmentally determined, making P. patens a powerful system in which to address mechanisms of peroxisome proliferation and division.

An ancestral role in peroxisome assembly is retained by the divisional peroxin Pex11 in the yeast Yarrowia lipolytica

The peroxin Pex11 has a recognized role in peroxisome division. Pex11p remodels and elongates peroxisomal membranes prior to the recruitment of dynamin-related GTPases that act in membrane scission to divide peroxisomes. We performed a comprehensive comparative genomics survey to understand the significance of the evolution of the Pex11 protein family in yeast and other eukaryotes. Pex11p is highly conserved and ancestral, and has undergone numerous lineage-specific duplications, whereas other Pex11 protein family members are fungal-specific innovations. Functional characterization of the in-silico-predicted Pex11 protein family members of the yeast Yarrowia lipolytica, i.e. Pex11p, Pex11Cp and Pex11/25p, demonstrated that Pex11Cp and Pex11/25p have a role in the regulation of peroxisome size and number characteristic of Pex11 protein family members. Unexpectedly, deletion of PEX11 in Y. lipolytica produces cells that lack morphologically identifiable peroxisomes, mislocalize peroxisomal matrix proteins and preferentially degrade peroxisomal membrane proteins, i.e. they exhibit the classical pex mutant phenotype, which has not been observed previously in cells deleted for the PEX11 gene. Our results are consistent with an unprecedented role for Pex11p in de novo peroxisome assembly.

Expression profile of maize microRNAs corresponding to their target genes under drought stress

Microarray assay of four inbred lines was used to identify 303 microRNAs differentially expressed under drought stress. The microRNAs were used for bioinformatics prediction of their target genes. The majority of the differentially expressed microRNA families showed different expression profiles at different time points of the stress process among the four inbred lines. Digital gene expression profiling revealed 54 genes targeted by 128 of the microRNAs differentially expressed under the same stress conditions. The differential expression of miR159 and miR168 was further validated by locked nucleic acid northern hybridization. These results indicated that miR159 and miR168, as well as numerous other microRNAs, play critical roles in signaling pathways of maize response to drought stress. However, the level of the post-transcriptional regulation mediated by microRNAs had different responses among genotypes, and the gene expression related to signaling pathways under drought stress is also regulated, possibly by multiple mechanisms.

Peroxisomes and photomorphogenesis

In higher plants, light-grown seedlings exhibit photomorphogenesis, a developmental program controlled by a complex web of interactions between photoreceptors, central repressors, and downstream effectors that leads to changes in gene expression and physiological changes. Light induces peroxisomal proliferation through a phytochrome A-mediated pathway, in which the transcription factor HYH activates the peroxisomal proliferation factor gene PEX11b. Microarray analysis revealed that light activates the expression of a number of peroxisomal genes, especially those involved in photorespiration, a process intimately associated with photosynthesis. In contrast, light represses the expression of genes involved in β-oxidation and the glyoxylate cycle, peroxisomal pathways essential for seedling establishment before photosynthesis begins. Furthermore, the peroxisome is a source of signaling molecules, notably nitric oxide, which promotes photomorphogenesis. Lastly, a gain-of-function mutant of the peroxisomal membrane-tethered RING-type E3 ubiquitin ligase PEX2 partially suppresses the phenotype of the photomorphogenic mutant det1. Possible mechanisms underlying this phenomenon are discussed.

Expression profile analysis of the polygalacturonase-inhibiting protein genes in rice and their responses to phytohormones and fungal infection

Polygalacturonase-inhibiting proteins (PGIPs) are typically leucine-rich repeat (LRR) proteins that can inhibit the activity of fungal polygalacturonases (PGs). In this study, two new Ospgip genes, named Ospgip6 and Ospgip7 with consensus sequence of ten imperfect LRR motif located on rice chromosomes 8 and 9, were identified using BLAST analysis. Both of them appear to be extracellular glycoproteins. To have a global view of the dynamic gene expression pattern, seven Ospgip genes were first analyzed using the Affymetrix rice genome array data from online resource. All of these seven Ospgip genes showed variable expression patterns among tissues/organs. In order to further investigate the potential function of these Ospgip genes, the responses of Ospgip genes to the treatment of various phytohormones (abscisic acid, brassinosteroid, gibberellic acid, 3-indole acetic acid, jasmonic acid, kinetin, naphthalene acetic acid and salicylic acid) as well as fungal infection were analyzed by real-time PCR using time course array. Generally, all the Ospgip genes were slightly up-regulated in the indica rice cultivar Minghui 63 under GA(3), KT and NAA treatments (except Ospgip2, which was down-regulated under KT treatment). In the japonica rice cultivar Zhonghua 11, Ospgip genes were regulated by most treatments with the response time variability. We also analyzed putative cis-elements in the promoter regions of Ospgip genes. This dataset provided a versatile resource to understand the regulatory network of Ospgip genes during the process of phytohormones treatment and fungal infection in the model monocotyledonous plant, rice, and could aid in the transgenic breeding against rice fungal diseases.

KEY MESSAGE

All the seven Ospgip genes showed variable expression patterns in Minghui 63 and their expressions were regulated by different phytohormone treatments or fungal infection in Minghui 63 and Zhonghua 11.

[Bioinformatic research of the family of PEX11, peroxisome proliferous factor in fungus]

The family members of PEX11 are key factors involved in regulation of peroxisome proliferation. Sixty-six PEX11p candidates of PEX11 gene family from 26 representative fungal species were obtained and analyzed by bioinformatic strategies. In most filamentous fungi, 2 or 3 potential PEX11ps were found, in contrast with 1 or 2 in yeast species. Compared with other fungal species, the Ascomycetes tend to have more PEX11ps, and even 5 in several individuals. The data of phylogenetic analysis and protein structure indicated that all of the PEX11ps were divided into 3 groups: I, II, and III. The members of group I and group III existed in most species, while those in group II were found only in Pezizomycotina. By MEME analysis, 5-6 conserved motifs were found in each PEX11ps. Among them,motif 8 in C-terminal had the most conservation, indicating that this motif probably plays a key role in maintaining the proper function of PEX11p.

Membrane elongation factors in organelle maintenance: the case of peroxisome proliferation

Separation of metabolic pathways in organelles is critical for eukaryotic life. Accordingly, the number, morphology and function of organelles have to be maintained through processes linked with membrane remodeling events. Despite their acknowledged significance and intense study many questions remain about the molecular mechanisms by which organellar membranes proliferate. Here, using the example of peroxisome proliferation, we give an overview of how proteins elongate membranes. Subsequent membrane fission is achieved by dynamin-related proteins shared with mitochondria. We discuss basic criteria that membranes have to fulfill for these fission factors to complete the scission. Because peroxisome elongation is always associated with unequal distribution of matrix and membrane proteins, we propose peroxisomal division to be non-stochastic and asymmetric. We further show that these organelles need not be functional to carry on membrane elongation and present the most recent findings concerning members of the Pex11 protein family as membrane elongation factors. These factors, beside known proteins such as BAR-domain proteins, represent another family of proteins containing an amphipathic α-helix with membrane bending activity.

A dynamic gene expression atlas covering the entire life cycle of rice

Growth and development of a plant are controlled by programmed expression of suits of genes at the appropriate time, tissue and abundance. Although genomic resources have been developed rapidly in recent years in rice, a model plant for cereal genome research, data of gene expression profiling are still insufficient to relate the developmental processes to transcriptomes, leaving a large gap between the genome sequence and phenotype. In this study, we generated genome-wide expression data by hybridizing 190 Affymetrix GeneChip Rice Genome Arrays with RNA from 39 tissues collected throughout the life cycle of the rice plant from two varieties, Zhenshan 97 and Minghui 63. Analyses of the global transcriptomes revealed many interesting features of dynamic patterns of gene expression across the tissues and stages. In total, 38 793 probe sets were detected as expressed and 69% of the expressed transcripts showed significantly variable expression levels among tissues/organs. We found that similarity of transcriptomes among organs corresponded well to their developmental relatedness. About 5.2% of the expressed transcripts showed tissue-specific expression in one or both varieties and 22.7% of the transcripts exhibited constitutive expression including 19 genes with high and stable expression in all the tissues. This dataset provided a versatile resource for plant genomic research, which can be used for associating the transcriptomes to the developmental processes, understanding the regulatory network of these processes, tracing the expression profile of individual genes and identifying reference genes for quantitative expression analyses.

Salt stress causes peroxisome proliferation, but inducing peroxisome proliferation does not improve NaCl tolerance in Arabidopsis thaliana

The PEX11 family of peroxisome membrane proteins have been shown to be involved in regulation of peroxisome size and number in plant, animals, and yeast cells. We and others have previously suggested that peroxisome proliferation as a result of abiotic stress may be important in plant stress responses, and recently it was reported that several rice PEX11 genes were up regulated in response to abiotic stress. We sought to test the hypothesis that promoting peroxisome proliferation in Arabidopsis thaliana by over expression of one PEX11 family member, PEX11e, would give increased resistance to salt stress. We could demonstrate up regulation of PEX11e by salt stress and increased peroxisome number by both PEX11e over expression and salt stress, however our experiments failed to find a correlation between PEX11e over expression and increased peroxisome metabolic activity or resistance to salt stress. This suggests that although peroxisome proliferation may be a consequence of salt stress, it does not affect the ability of Arabidopsis plants to tolerate saline conditions.

Molecular characterization, expression pattern, and function analysis of the OsBC1L family in rice

COBRA-like proteins play important roles in cell expansion and cell wall biosynthesis in Arabidopsis. In rice, a COBRA-like gene, BRITTLE CULM1 (BC1), has been identified as a regulator controlling the culm mechanical strength. Analysis of the rice genome indicated that BC1 belongs to an 11-member multigene family, termed the OsBC1L family in this study. Based on sequence comparisons and phylogenetic analysis, the OsBC1L family comprises two main subgroups. Expression patterns examined by microarray and reverse transcription polymerase chain reaction revealed that OsBC1L genes exhibit universal or specific expression patterns. Through T-DNA or Tos17 insertion mutant lines, the functions of six OsBC1L family members have been examined by investigating the phenotype variations of knockout mutants under normal growth conditions. Results suggest that the OsBC1L genes perform a range of functions and participate in various developmental processes in rice.

The ankyrin repeat gene family in rice: genome-wide identification, classification and expression profiling

Ankyrin repeat (ANK) containing proteins comprise a large protein family. Although many members of this family have been implicated in plant growth, development and signal transduction, only a few ANK genes have been reported in rice. In this study, we analyzed the structures, phylogenetic relationship, genome localizations and expression profiles of 175 ankyrin repeat genes identified in rice (OsANK). Domain composition analysis suggested OsANK proteins can be classified into ten subfamilies. Chromosomal localizations of OsANK genes indicated nine segmental duplication events involving 17 genes and 65 OsANK genes were involved in tandem duplications. The expression profiles of 158 OsANK genes were analyzed in 24 tissues covering the whole life cycle of two rice genotypes, Minghui 63 and Zhenshan 97. Sixteen genes showed preferential expression in given tissues compared to all the other tissues in Minghui 63 and Zhenshan 97. Nine genes were preferentially expressed in stamen of 1 day before flowering, suggesting that these genes may play important roles in pollination and fertilization. Expression data of OsANK genes were also obtained with tissues of seedlings subjected to three phytohormone (NAA, GA3 and KT) and light/dark treatments. Eighteen genes showed differential expression with at least one phytohormone treatment while under light/dark treatments, 13 OsANK genes showed differential expression. Our data provided a very useful reference for cloning and functional analysis of members of this gene family in rice.

Comprehensive sequence and expression profile analysis of Hsp20 gene family in rice

The Hsp20 genes represent the most abundant small heat shock proteins (sHSPs) in plants. Hsp20 gene family has been shown to be involved in preventing heat shock and promoting resistance to environmental stress factors, but very little is known about this gene family in rice. Here, we report the identification and characterization of 39 OsHsp20 genes in rice, describing the gene structure, gene expression, genome localization, and phylogenetic relationship of each member. We have used RT-PCR to perform a characterization of the normal and heat shock-induced expression of selective OsHsp20 genes. A genome-wide microarray based gene expression analysis involving 25 stages of vegetative and reproductive development in three rice cultivars has revealed that 36 OsHsp20 genes were expressed in at least one of the experimental stages studied. Among these, transcripts of OsHsp20 were accumulated differentially during vegetative and reproductive developmental stages and preferentially down-regulated in Shanyou 63. In addition, OsHsp20 genes were identified as showing prominent heterosis in family-level expression. Our results suggest that the expression patterns of the OsHsp20 genes are diversified not only in developmental stages but also in variety level.