Introns and protein revolution--an analysis of the exon/intron organisation of actin genes

Genome-wide identification of GhRLCK-VII subfamily genes in Gossypium hirsutum and investigation of their functions in resistance to Verticillium wilt

BACKGROUND

The receptor-like cytoplasmic kinases subfamily VII (RLCK-VII) is critical in regulating plant growth, development, and pattern-triggered immunity. However, a comprehensive exploration of these genes in the allotetraploid Gossypium hirsutum is still lacking. This study aimed to identify RLCK-VII genes in G. hirsutum and investigate their evolutionary history, structural features, expression patterns, and role in plant defense.

RESULTS

Seventy-two RLCK-VII genes in the G. hirsutum genome were unveiled and classified into nine groups following their phylogenetic analysis with Arabidopsis thaliana. Group VII-1 was the largest, accounting for 28%, while Groups VII-2 and VII-3 had only one member each. The analysis using MCScanX revealed that these 72 genes formed 166 collinear gene pairs and were resided on 26 chromosomes of G. hirsutum, suggesting that they were derived from whole genome segmental duplication events. Their calculated Ka/Ks values were below one, implying the occurrence of purification selection during the evolution and inhibition of gene function differentiation/loss. All members of the RLCK-VII subfamily possessed two conserved domains, PKinase-Tyr and PKinase, and several conserved PBS1 kinase subdomains, individually included in one of the ten motifs identified using MEME. The RNA-Seq results showed that RLCK-VII genes exhibited different spatiotemporal expression, indicating their involvement in cotton growth, development, and defense responses to Verticillium dahliae. The transcription patterns of RLCK-VII genes found by RNA-Seq were further validated using qRT-PCR assays after inoculating "20B12" (cotton cultivar) with "V991" (V. dahliae). The virus-induced gene silencing (VIGS) assays uncovered that two RLCK-VII genes (Gohir.A13G227248 and Gohir.A10G219900) were essential to G. hirsutum resistance to Verticillium wilt.

CONCLUSIONS

These observations offer valuable insight into the attributes and roles of RLCK-VII genes in G. hirsutum, potentially enable the breeding of new cotton cultivars with enhanced resistance to Verticillium wilt.

The role of Munc18-1 and its orthologs in modulation of cortical F-actin in chromaffin cells

Munc18-1 was originally described as an essential docking factor in chromaffin cells. Recent findings showed that Munc18-1 has an additional role in the regulation of the cortical F-actin network, which is thought to function as a physical barrier preventing secretory vesicles from access to their release sites under resting conditions. In our review, we discuss whether this function is evolutionarily conserved in all Sec1/Munc18-like (SM) proteins. In addition, we introduce a new quantification method that improves the analysis of cortical filamentous actin (F-actin) in comparison with existing methods. Since the docking process is highly evolutionarily conserved in the SM protein superfamily, we use our novel quantification method to investigate whether the F-actin-regulating function is similarly conserved among SM proteins. Our preliminary data suggest that the regulation of cortical F-actin is a shared function of SM proteins, and we propose a way to gain more insight in the molecular mechanism underlying the Munc18-1-mediated cortical F-actin regulation.

Genome-wide expansion and expression divergence of the basic leucine zipper transcription factors in higher plants with an emphasis on sorghum

Plant bZIP transcription factors play crucial roles in multiple biological processes. However, little is known about the sorghum bZIP gene family although the sorghum genome has been completely sequenced. In this study, we have carried out a genome-wide identification and characterization of this gene family in sorghum. Our data show that the genome encodes at least 92 bZIP transcription factors. These bZIP genes have been expanded mainly by segmental duplication. Such an expansion mechanism has also been observed in rice, arabidopsis and many other plant organisms, suggesting a common expansion mode of this gene family in plants. Further investigation shows that most of the bZIP members have been present in the most recent common ancestor of sorghum and rice and the major expansion would occur before the sorghum-rice split era. Although these bZIP genes have been duplicated with a long history, they exhibited limited functional divergence as shown by nonsynonymous substitutions (Ka)/synonymous substitutions (Ks) analyses. Their retention was mainly due to the high percentages of expression divergence. Our data also showed that this gene family might play a role in multiple developmental stages and tissues and might be regarded as important regulators of various abiotic stresses and sugar signaling.

The actin gene in Paracoccidioides brasiliensis: organization, expression and phylogenetic analyses

PbrACT1, the gene responsible for the synthesis of actin in Paracoccidioides brasiliensis, was found as a single copy, organized into six exons and five introns. Its open reading frame (ORF) codes for a putative protein of 375 amino acids, with a molecular mass of 41.5 kDa and an isoelectric point of 5.6. Analysis of the nucleotide sequence revealed a high homology to other fungal actins, the presence of characteristic fungal actin sequences, and heat shock elements at the 5' untranslated region (UTR). Phylogenetic analyses with deduced amino acid sequences of fungal actins grouped P. brasiliensis within the phylum Ascomycota, order Onygenales, in concordance with a few previous reports. Patterns of expression through the temperature-induced morphological transitions from mycelial to yeast-like shapes and reverse, suggests that PbrACT1 is regulated in this process. The PbrACT1 gene sequence is available at the GenBank database under accession number AY383732.

Tempo and Mode of Spliceosomal Intron Evolution in Actin of Foraminifera

Spliceosomal introns are present in almost all eukaryotic genes, yet little is known about their origin and turnover in the majority of eukaryotic phyla. There is no agreement whether most introns are ancestral and have been lost in some lineage or have been gained recently. We addressed this question by analyzing the spatial and temporal distribution of introns in actins of foraminifera, a group of testate protists whose exceptionally rich fossil record permits the calibration of molecular phylogenies to date intron origins. We identified 24 introns dispersed along the sequence of two foraminiferan actin paralogues and actin deviating proteins, an unconventional type of fast-evolving actin found in some foraminifera. Comparison of intron positions indicates that 20 of 24 introns are specific to foraminifera. Four introns shared between foraminifera and other eukaryotes were interpreted as parallel gains because they have been found only in single species belonging to phylogenetically distinctive lineages. Moreover, additional recent intron gain due to the transfer between the actin paralogues was observed in two cultured species. Based on a relaxed molecular clock timescale, we conclude that intron gains in actin took place throughout the evolution of foraminifera, with the oldest introns inserted between 550 and 500 million years ago and the youngest ones acquired less than 100 million years ago.

Evolution of intron/exon structure of DEAD helicase family genes in Arabidopsis, Caenorhabditis, and Drosophila

The DEAD box RNA helicase (RH) proteins are homologs involved in diverse cellular functions in all of the organisms from prokaryotes to eukaryotes. Nevertheless, there is a lack of conservation in the splicing pattern in the 53 Arabidopsis thaliana (AtRHs), the 32 Caenorhabditis elegans (CeRHs) and the 29 Drosophila melanogaster (DmRHs) genes. Of the 153 different observed intron positions, 4 are conserved between AtRHs, CeRHs, and DmRHs, and one position is also found in RHs from yeast and human. Of the 27 different AtRH structures with introns, 20 have at least one predicted ancient intron in the regions coding for the catalytic domain. In all of the organisms examined, we found at least one gene with most of its intron predicted to be ancient. In A. thaliana, the large diversity in RH structures suggests that duplications of the ancestral RH were followed by a high number of intron deletions and additions. The very high bias toward phase 0 introns is in favor of intron addition, preferentially in phase 0. Results from this comparative study of the same gene family in a plant and in two animals are discussed in terms of the general mechanisms of gene family evolution.

The evolution of controlled multitasked gene networks: the role of introns and other noncoding RNAs in the development of complex organisms

Eukaryotic phenotypic diversity arises from multitasking of a core proteome of limited size. Multitasking is routine in computers, as well as in other sophisticated information systems, and requires multiple inputs and outputs to control and integrate network activity. Higher eukaryotes have a mosaic gene structure with a dual output, mRNA (protein-coding) sequences and introns, which are released from the pre-mRNA by posttranscriptional processing. Introns have been enormously successful as a class of sequences and comprise up to 95% of the primary transcripts of protein-coding genes in mammals. In addition, many other transcripts (perhaps more than half) do not encode proteins at all, but appear both to be developmentally regulated and to have genetic function. We suggest that these RNAs (eRNAs) have evolved to function as endogenous network control molecules which enable direct gene-gene communication and multitasking of eukaryotic genomes. Analysis of a range of complex genetic phenomena in which RNA is involved or implicated, including co-suppression, transgene silencing, RNA interference, imprinting, methylation, and transvection, suggests that a higher-order regulatory system based on RNA signals operates in the higher eukaryotes and involves chromatin remodeling as well as other RNA-DNA, RNA-RNA, and RNA-protein interactions. The evolution of densely connected gene networks would be expected to result in a relatively stable core proteome due to the multiple reuse of components, implying that cellular differentiation and phenotypic variation in the higher eukaryotes results primarily from variation in the control architecture. Thus, network integration and multitasking using trans-acting RNA molecules produced in parallel with protein-coding sequences may underpin both the evolution of developmentally sophisticated multicellular organisms and the rapid expansion of phenotypic complexity into uncontested environments such as those initiated in the Cambrian radiation and those seen after major extinction events.

EID: the Exon-Intron Database-an exhaustive database of protein-coding intron-containing genes

To aid studies of molecular evolution and to assist in gene prediction research, we have constructed an Exon-Intron Database (EID) in FASTA format. Currently, the database is derived from GenBank release 112, and it contains 51 289 protein-coding genes (287 209 exons) that harbor introns, along with extensive descriptions of each gene and its DNA and protein sequences, as well as splice motif information. There is 17% redundancy inherited from GenBank-a purge at the 99% identity level reduced the database to 42 460 genes (243 589 exons). We have created subdatabases of genes whose intron positions have been experimentally determined. One such database, constructed by comparing genomic and mRNA sequences, contains 11 242 genes (62 474 exons). A larger database of 22 196 genes (105 595 exons) was constructed by selecting on keywords to eliminate computer-predicted genes. By examining the two nucleotides adjacent to the intron boundary, we infer that there is a 2% rate of errors or other deviations from the standard GTellipsisAG motif in nuclear genes. This criterion can be used to eliminate 4921 genes from the overall database. Various tools are provided to enable generation of user-specific subsets of the EID. The EID distribution can be obtained from http://mcb.harvard.edu/gilbert/EID

Type 2 rice metallothionein-like gene has two introns

A type 2 rice metallothionein-like gene was isolated from root by PCR and sequenced. The PCR fragment was designated as pcr1460, which overlaps with OsMT-2, a cDNA sequence previously characterized, with the presence of two additional segments, 583 and 613 bp in length. These segments are recognized as introns which divide the coding sequence into three exons, 65, 78 and 106 base pairs in length. The sequences flanking the introns conform with the GT/AG rule for splice junctions, and one exonic open reading frame can be identified in each of the introns. The observation that MT-like gene has two introns is the first of such a finding obtained from monocotyledonous plants.