Genome-Wide Identification and Analysis of Genes Encoding Proteolytic Enzymes in Pineapple

Pineapple, Ananas comosus, is an economically important fruit crop. Recently its genome was completely sequenced and a total of 27,024 protein coding genes were predicted. Using a set of well evaluated bioinformatics tools we have predicted the protein subcellular locations and comparatively analyzed the protein conserved domains of the predicted proteomes in pineapple, Oryza sativa (rice), Sorghum bicolor (sorghum), and Brachypodium distachyson. Our analysis revealed that ~24–26 % of proteins were located in nucleus, 17–21 % in cytosol, 9–11 % in chloroplast, and 8–11 % proteins were secreted in these monocot plants. The secretomes in the four species were analyzed comparatively and a large number of secreted glycosyl hydrolases were identified. As pineapple proteolytic enzymes, knowns as bromelains, have been used for medical treatments, we focused on genome-wide identification and analysis of pineapple genes encoding proteases. A total of 512 pineapple genes encoding putative proteolytic enzymes were identified, with 152 secreted, 74 localized in cytosol, 67 in nucleus, 60 in chloroplast, 18 in mitochondria, and the remaining in other subcellular locations. The top large protease families in pineapple were papain family cysteine protease (62 genes), peptidase S8 family (56 genes), aspartyl protease family (38 genes), and serine carboxypeptidase (33 genes). Gene expression analysis revealed that among 512 protease genes 432 were expressed in various tissues and 72 genes were differentially expressed. The highly expressed protease genes were identified including 7 papain family cysteine proteases. The protease genes with the predicted protein subcellular locations will facilitate the efforts for examining their biological roles in pineapple growth and development and for expressing the recombinant proteases for medical use. The information of protein subcellular location of all plant species can be accessed at the PlantSecKB website (http://proteomics.ysu.edu/secretomes/plant.php).     

Characterization and the Expression Analysis of Nitrate Transporter (NRT) Gene Family in Pineapple

Nitrate is the preferred nitrogen source of higher plants and an essential nutrient for plant growth and development. Nitrate transporters (NRTs) play vital roles in the nitrate uptake and transportation. However, the NRT gene family in pineapple is still unexplored. In this study, we performed a genome-wide analysis of the pineapple genome and identified 48 NRT genes (AcNRTs) distributed unevenly across 9 chromosomes and 2 scaffolds. Phylogenetic analysis showed that these genes can be divided into three groups, namely, AcNRT1/PTR, AcNRT2 and AcNRT3/NAR1 with 44, 3 and 1 members, respectively. AcNRTs within the same phylogenetic group share similar gene structure and domain composition. In addition, syntenic and phylogenetic analyses identified 34 Arabidopsis NRT genes with 31 pineapple NRT genes as orthologs. By investigating the expression profiles of these genes in various tissues, we showed that the expression pattern of some AcNRTs genes is tissue-specific. Furthermore, we examined the expression of the AcNRT2s under nitrate starvation and found that AcNRT2.1 and AcNRT2.2 both have the strongest response in roots suggesting that AcNRTs may play a broad role in the pineapple in response to nitrate deficiency. Taken together, our data provide insights into the evolution and function of pineapple NRTs and pave a path for future functional investigation of pineapple NRTs genes.     

A preliminary integrated genetic map distinguishes every chromosome pair and locates essential genes related to abiotic adaptation of <Emphasis Type="Italic">Crassostrea angulata</Emphasis>/<Emphasis Type="Italic">gigas</Emphasis>

Background

The re-sequencing of C. angulata has revealed many polymorphisms in candidate genes related to adaptation to abiotic stress that are not present in C. gigas; these genes, therefore, are probably related to the ability of this oyster to retain high concentrations of toxic heavy metals. There is, in addition, an unresolved controversy as to whether or not C. angulata and C. gigas are the same species or subspecies. Both oysters have 20 metacentric chromosomes of similar size that are morphologically indistinguishable. From a genomic perspective, as a result of the great variation and selection for heterozygotes in C. gigas, the assembly of its draft genome was difficult: it is fragmented in more than seven thousand scaffolds.

Results

In this work sixty BAC sequences of C. gigas downloaded from NCBI were assembled in BAC-contigs and assigned to BACs that were used as probes for mFISH in C. angulata and C. gigas. In addition, probes of H3, H4 histone, 18S and 5S rDNA genes were also used. Hence we obtained markers identifying 8 out the 10 chromosomes constituting the karyotype. Chromosomes 1 and 9 can be distinguished morphologically. The bioinformatic analysis carried out with the BAC-contigs annotated 88 genes. As a result, genes associated with abiotic adaptation, such as metallothioneins, have been positioned in the genome. The gene ontology analysis has also shown many molecular functions related to metal ion binding, a phenomenon associated with detoxification processes that are characteristic in oysters. Hence the provisional integrated map obtained in this study is a useful complementary tool for the study of oyster genomes.

Conclusions

In this study 8 out of 10 chromosome pairs of Crassostrea angulata/gigas were identified using BAC clones as probes. As a result all chromosomes can now be distinguished. Moreover, FISH showed that H3 and H4 co-localized in two pairs of chromosomes different that those previously escribed. 88 genes were annotated in the BAC-contigs most of them related with Molecular Functions of protein binding, related to the resistance of the species to abiotic stress. An integrated genetic map anchored to the genome has been obtained in which the BAC-contigs structure were not concordant with the gene structure of the C. gigas scaffolds displayed in the Genomicus database.

     

Manipulating mRNA splicing by base editing in plants

Precursor-mRNAs(pre-mRNA) can be processed into one or more mature m RNA isoforms through constitutive or alternative splicing pathways. Constitutive splicing of pre-mRNA plays critical roles in gene expressional regulation, such as intronmediated enhancement(IME), whereas alternative splicing(AS) dramatically increases the protein diversity and gene functional regulation. However, the unavailability of mutants for individual spliced isoforms in plants has been a major limitation in studying the function of mRNA splicing. Here, we describe an efficient tool for manipulating the splicing of plant genes. Using a Cas9-directed base editor, we converted the 5′ splice sites in four Arabidopsis genes from the activated GT form to the inactive AT form. Silencing the AS of HAB 1.1(encoding a type 2 C phosphatase) validated its function in abscisic acid signaling, while perturbing the AS of RS31 A revealed its functional involvement in plant response to genotoxic treatment for the first time. Lastly,altering the constitutive splicing of Act2 via base editing facilitated the analysis of IME. This strategy provides an efficient tool for investigating the function and regulation of gene splicing in plants and other eukaryotes.