Molecular analysis of the myosin gene family in Arabidopsis thaliana

Myosin is believed to act as the molecular motor for many actin-based motility processes in eukaryotes. It is becoming apparent that a single species may possess multiple myosin isoforms, and at least seven distinct classes of myosin have been identified from studies of animals, fungi, and protozoans. The complexity of the myosin heavy-chain gene family in higher plants was investigated by isolating and characterizing myosin genomic and cDNA clones from Arabidopsis thaliana. Six myosin-like genes were identified from three polymerase chain reaction (PCR) products (PCR1, PCR11, PCR43) and three cDNA clones (ATM2, MYA2, MYA3). Sequence comparisons of the deduced head domains suggest that these myosins are members of two major classes. Analysis of the overall structure of the ATM2 and MYA2 myosins shows that they are similar to the previously-identified ATM1 and MYA1 myosins, respectively. The MYA3 appears to possess a novel tail domain, with five IQ repeats, a six-member imperfect repeat, and a segment of unique sequence. Northern blot analyses indicate that some of the Arabidopsis myosin genes are preferentially expressed in different plant organs. Combined with previous studies, these results show that the Arabidopsis genome contains at least eight myosin-like genes representing two distinct classes.

Expression of HPV-11 L1 protein in transgenic Arabidopsis thaliana and Nicotiana tabacum

Background We have investigated the possibility and feasibility of producing the HPV-11 L1 major capsid protein in transgenic Arabidopsis thaliana ecotype Columbia and Nicotiana tabacum cv. Xanthi as potential sources for an inexpensive subunit vaccine. Results Transformation of plants was only achieved with the HPV-11 L1 gene with the C-terminal nuclear localization signal (NLS-) encoding region removed, and not with the full-length gene. The HPV-11 L1 NLS- gene was stably integrated and inherited through several generations of transgenic plants. Plant-derived HPV-11 L1 protein was capable of assembling into virus-like particles (VLPs), although resulting particles displayed a pleomorphic phenotype. Neutralising monoclonal antibodies binding both surface-linear and conformation-specific epitopes bound the A. thaliana-derived particles and - to a lesser degree - the N. tabacum-derived particles, suggesting that plant-derived and insect cell-derived VLPs displayed similar antigenic properties. Yields of up to 12 μg/g of HPV-11 L1 NLS- protein were harvested from transgenic A. thaliana plants, and 2 μg/g from N. tabacum plants - a significant increase over previous efforts. Immunization of New Zealand white rabbits with ∼50 μg of plant-derived HPV-11 L1 NLS- protein induced an antibody response that predominantly recognized insect cell-produced HPV-11 L1 NLS- and not NLS+ VLPs. Evaluation of the same sera concluded that none of them were able to neutralise pseudovirion in vitro. Conclusion We expressed the wild-type HPV-11 L1 NLS- gene in two different plant species and increased yields of HPV-11 L1 protein by between 500 and 1000-fold compared to previous reports. Inoculation of rabbits with extracts from both plant types resulted in a weak immune response, and antisera neither reacted with native HPV-11 L1 VLPs, nor did they neutralise HPV-11 pseudovirion infectivity. This has important and potentially negative implications for the production of HPV-11 vaccines in plants. © 2007 Kohl et al; licensee BioMed Central Ltd.

DRB2 is required for microRNA biogenesis in Arabidopsis thaliana

Background The Arabidopsis thaliana (Arabidopsis) DOUBLE-STRANDED RNA BINDING (DRB) protein family consists of five members, DRB1 to DRB5. The biogenesis of two developmentally important small RNA (sRNA) species, the microRNAs (miRNAs) and trans-acting small interfering RNAs (tasiRNAs) by DICER-LIKE (DCL) endonucleases requires the assistance of DRB1 and DRB4 respectively. The importance of miRNA-directed target gene expression in plant development is exemplified by the phenotypic consequence of loss of DRB1 activity (drb1 plants). Principal Findings Here we report that the developmental phenotype of the drb235 triple mutant plant is the result of deregulated miRNA biogenesis in the shoot apical meristem (SAM) region. The expression of DRB2, DRB3 and DRB5 in wild-type seedlings is restricted to the SAM region. Small RNA sequencing of the corresponding tissue of drb235 plants revealed altered miRNA accumulation. Approximately half of the miRNAs detected remained at levels equivalent to those of wild-type plants. However, the accumulation of the remaining miRNAs was either elevated or reduced in the triple mutant. Examination of different single and multiple drb mutants revealed a clear association between the loss of DRB2 activity and altered accumulation for both the elevated and reduced miRNA classes. Furthermore, we show that the constitutive over-expression of DRB2 outside of its wild-type expression domain can compensate for the loss of DRB1 activity in drb1 plants. Conclusions/Significance Our results suggest that in the SAM region, DRB2 is both antagonistic and synergistic to the role of DRB1 in miRNA biogenesis, adding an additional layer of gene regulatory complexity in this developmentally important tissue.

The Arabidopsis thaliana double-stranded RNA binding protein DRB1 directs guide strand selection from microRNA duplexes

In Arabidopsis thaliana (Arabidopsis), DICER-LIKE1 (DCL1) functions together with the double-stranded RNA binding protein (dsRBP), DRB1, to process microRNAs (miRNAs) from their precursor transcripts prior to their transfer to the RNA-induced silencing complex (RISC). miRNA-loaded RISC directs RNA silencing of cognate mRNAs via ARGONAUTE1 (AGO1)-catalyzed cleavage. Short interefering RNAs (siRNAs) are processed from viral-derived or transgene-encoded molecules of doublestranded RNA (dsRNA) by the DCL/dsRBP partnership, DCL4/DRB4, and are also loaded to AGO1-catalyzed RISC for cleavage of complementary mRNAs. Here, we use an artificial miRNA (amiRNA) technology, transiently expressed in Nicotiana benthamiana, to produce a series of amiRNA duplexes with differing intermolecular thermostabilities at the 5′ end of duplex strands. Analyses of amiRNA duplex strand accumulation and target transcript expression revealed that strand selection (amiRNA and amiRNA*) is directed by asymmetric thermostability of the duplex termini. The duplex strand possessing a lower 59 thermostability was preferentially retained by RISC to guide mRNA cleavage of the corresponding target transgene. In addition, analysis of endogenous miRNA duplex strand accumulation in Arabidopsis drb1 and drb2345 mutant plants revealed that DRB1 dictates strand selection, presumably by directional loading of the miRNA duplex onto RISC for passenger strand degradation. Bioinformatic and Northern blot analyses of DCL4/DRB4-dependent small RNAs (miRNAs and siRNAs) revealed that small RNAs produced by this DCL/dsRBP combination do not conform to the same terminal thermostability rules as those governing DCL1/DRB1-processed miRNAs. This suggests that small RNA processing in the DCL1/DRB1-directed miRNA and DCL4/DRB4-directed sRNA biogenesis pathways operates via different mechanisms.

DRB2, DRB3 and DRB5 function in a non-canonical microRNA pathway in Arabidopsis thaliana

DOUBLE-STRANDED RNA BIN DIN G (DRB) proteins have been functionally characterized in viruses, prokaryotes and eukaryotes and are involved in all aspects of RNA biology. Arabidopsis thaliana (Arabidopsis) encodes five closely related DRB proteins, DRB1 to DRB5. DRB1 and DRB4 are required by DICER-LIKE (DCL) proteins DCL1 and DCL4 to accurately and efficiently process structurally distinct double-stranded RNA (dsRNA) precursor substrates in the microRNA (miRNA) and trans-acting small-interfering RNA (tasiRNA) biogenesis pathways respectively. We recently reported that DRB2 is also involved in the biogenesis of specific miRNA subsets. Furthermore, the severity of the developmental phenotype displayed by the drb235 triple mutant plant, compared with those expressed by either drb2, drb3 and drb5 single mutants, or double mutant combinations thereof, indicates that DRB3 and DRB5 function in the same non-canonical miRNA pathway as DRB2. Through the use of our artificial miRNA (amiRNA) plant expression vector, pBlueGreen 2,3 we demonstrate here that unlike DRB2, DRB3 and DRB5 are not involved in the dsRNA processing stages of the miRNA biogenesis pathway, but are required to mediate RNA silencing of target genes of DRB2-associated miRNA s. © 2012 Landes Bioscience.

The use of artificial MicroRNA technology to control gene expression in Arabidopsis thaliana

In plants, double-stranded RNA (dsRNA) is an effective trigger of RNA silencing, and several classes of endogenous small RNA (sRNA), processed from dsRNA substrates by DICER-like (DCL) endonucleases, are essential in controlling gene expression. One such sRNA class, the microRNAs (miRNAs) control the expression of closely related genes to regulate all aspects of plant development, including the determination of leaf shape, leaf polarity, flowering time, and floral identity. A single miRNA sRNA silencing signal is processed from a long precursor transcript of nonprotein-coding RNA, termed the primary miRNA (pri-miRNA). A region of the pri-miRNA is partially self-complementary allowing the transcript to fold back onto itself to form a stem-loop structure of imperfectly dsRNA. Artificial miRNA (amiRNA) technology uses endogenous pri-miRNAs, in which the miRNA and miRNA*(passenger strand of the miRNA duplex) sequences have been replaced with corresponding amiRNA/ amiRNA*sequences that direct highly efficient RNA silencing of the targeted gene. Here, we describe the rules for amiRNA design, as well as outline the PCR and bacterial cloning procedures involved in the construction of an amiRNA plant expression vector to control target gene expression in Arabidopsis thaliana. © 2014 Springer Science+Business Media New York.

Efficient silencing of endogenous microRNAs using artificial microRNAs in Arabidopsis thaliana

We report here that the expression of endogenous microRNAs (miRNAs) can be efficiently silenced in Arabidopsis thaliana (Arabidopsis) using artificial miRNA (amiRNA) technology. We demonstrate that an amiRNA designed to target a mature miRNA directs silencing against all miRNA family members, whereas an amiRNA designed to target the stem-loop region of a miRNA precursor transcript directs silencing against only the individual family member targeted. Furthermore, our results indicate that amiRNAs targeting both the mature miRNA and stem-loop sequence direct RNA silencing through cleavage of the miRNA precursor transcript, which presumably occurs in the nucleus of a plant cell during the initial stages of miRNA biogenesis. This suggests that small RNA (sRNA)-guided RNA cleavage in plants occurs not only in the cytoplasm, but also in the nucleus. Many plant miRNA gene families have been identified via sequencing and bioinformatic analysis, but, to date, only a small tranche of these have been functionally characterized due to a lack of effective forward or reverse genetic tools. Our findings therefore provide a new and powerful reverse-genetic tool for the analysis of miRNA function in plants. © The Author 2010. Published by the Molecular Plant Shanghai Editorial Office in association with Oxford University Press on behalf of CSPP and IPPE, SIBS, CAS.

Effect of 5′UTR introns on gene expression in Arabidopsis thaliana

Background The majority of introns in gene transcripts are found within the coding sequences (CDSs). A small but significant fraction of introns are also found to reside within the untranslated regions (5′UTRs and 3′UTRs) of expressed sequences. Alignment of the whole genome and expressed sequence tags (ESTs) of the model plant Arabidopsis thaliana has identified introns residing in both coding and non-coding regions of the genome. Results A bioinformatic analysis revealed some interesting observations: (1) the density of introns in 5′UTRs is similar to that in CDSs but much higher than that in 3′UTRs; (2) the 5′UTR introns are preferentially located close to the initiating ATG codon; (3) introns in the 5′UTRs are, on average, longer than introns in the CDSs and 3′UTRs; and (4) 5′UTR introns have a different nucleotide composition to that of CDs and 3′UTR introns. Furthermore, we show that the 5′UTR intron of the A. thaliana EFIα-A3 gene affects the gene expression and the size of the 5′UTR intron influences the level of gene expression. Conclusion Introns within the 5′UTR show specific features that distinguish them from introns that reside within the coding sequence and the 3′UTR. In the EFIα-A3 gene, the presence of a long intron in the 5′UTR is sufficient to enhance gene expression in plants in a size dependent manner.

Environmental regulation of leaf colour in red 35S : PAP1 Arabidopsis thaliana

High-temperature, low-light (HTLL) treatment of 35S:PAP1 Arabidopsis thaliana over-expressing the PAP1 (Production of Anthocyanin Pigment 1) gene results in reversible reduction of red colouration, suggesting the action of additional anthocyanin regulators. High-performance liquid chromatography (HPLC), liquid chromatography mass spectrometry (LCMS) and Affimetrix®-based microarrays were used to measure changes in anthocyanin, flavonoids, and gene expression in response to HTLL. HTLL treatment of control and 35S:PAP1 A. thaliana resulted in a reversible reduction in the concentrations of major anthocyanins despite ongoing over-expression of the PAP1 MYB transcription factor. Twenty-one anthocyanins including eight cis-coumaryl esters were identified by LCMS. The concentrations of nine anthocyanins were reduced and those of three were increased, consistent with a sequential process of anthocyanin degradation. Analysis of gene expression showed down-regulation of flavonol and anthocyanin biosynthesis and of transport-related genes within 24 h of HTLL treatment. No catabolic genes up-regulated by HTLL were found. Reductions in the concentrations of anthocyanins and down-regulation of the genes of anthocyanin biosynthesis were achieved by environmental manipulation, despite ongoing over-expression of PAP1. Quantitative PCR showed reduced expression of three genes (TT8, TTG1 and EGL3) of the PAP1 transcriptional complex, and increased expression of the potential transcriptional repressors AtMYB3, AtMYB6 and AtMYBL2 coincided with HTLL-induced down-regulation of anthocyanin biosynthesis. HTLL treatment offers a model system with which to explore anthocyanin catabolism and to discover novel genes involved in the environmental control of anthocyanins.

Brca2 is involved in meiosis in Arabidopsis thaliana as suggested by its interaction with Dmc1

Two BRCA2-like sequences are present in the Arabidopsis genome. Both genes are expressed in flower buds and encode nearly identical proteins, which contain four BRC motifs. In a yeast two-hybrid assay, the Arabidopsis Brca2 proteins interact with Rad51 and Dmc1. RNAi constructs aimed at silencing the BRCA2 genes at meiosis triggered a reproducible sterility phenotype, which was associated with dramatic meiosis alterations. We obtained the same phenotype upon introduction of RNAi constructs aimed at silencing the RAD51 gene at meiosis in dmc1 mutant plants. The meiotic figures we observed strongly suggest that homologous recombination is highly disturbed in these meiotic cells, leaving aberrant recombination events to repair the meiotic double-strand breaks. The 'brca2' meiotic phenotype was eliminated in spo11 mutant plants. Our experiments point to an essential role of Brca2 at meiosis in Arabidopsis. We also propose a role for Rad51 in the dmc1 context.

Proteomic identification of putative microRNA394 target genes in Arabidopsis thaliana identifies MLP family members critical for normal development

Expression of the F-Box protein Leaf Curling Responsiveness (LCR) is regulated by microRNA, miR394, and alterations to this interplay in Arabidopsis thaliana produce defects in leaf polarity and shoot apical meristem (SAM) organisation. Although the miR394-LCR node has been documented in Arabidopsis, the identification of proteins targeted by LCR F-box itself has proven problematic. Here, a proteomic analysis of shoot apices from plants with altered LCR levels identified a member of the Major Latex Protein (MLP) family gene as a potential LCR F-box target. Bioinformatic and molecular analyses also suggested that other MLP family members are likely to be targets for this post-translational regulation. Direct interaction between LCR F-Box and MLP423 was validated. Additional MLP members had reduction in protein accumulation, in varying degrees, mediated by LCR F-Box. Transgenic Arabidopsis lines, in which MLP28 expression was reduced through an artificial miRNA technology, displayed severe developmental defects, including changes in leaf patterning and morphology, shoot apex defects, and eventual premature death. These phenotypic characteristics resemble those of Arabidopsis plants modified to over-express LCR. Taken together, the results demonstrate that MLPs are driven to degradation by LCR, and indicate that MLP gene family is target of miR394-LCR regulatory node, representing potential targets for directly post-translational regulation mediated by LCR F-Box. In addition, MLP28 family member is associated with the LCR regulation that is critical for normal Arabidopsis development.

An egg apparatus-specific enhancer of Arabidopsis, identified by enhancer detection

Despite a central role in angiosperm reproduction, few gametophyte-specific genes and promoters have been isolated, particularly for the inaccessible female gametophyte (embryo sac). Using the Ds-based enhancer-detector line ET253, we have cloned an egg apparatus-specific enhancer (EASE) from Arabidopsis (Arabidopsis thaliana). The genomic region flanking the Ds insertion site was further analyzed by examining its capability to control gusA and GFP reporter gene expression in the embryo sac in a transgenic context. Through analysis of a 5' and 3' deletion series in transgenic Arabidopsis, the sequence responsible for egg apparatus-specific expression was delineated to 77 bp. Our data showed that this enhancer is unique in the Arabidopsis genome, is conserved among different accessions, and shows an unusual pattern of sequence variation. This EASE works independently of position and orientation in Arabidopsis but is probably not associated with any nearby gene, suggesting either that it acts over a large distance or that a cryptic element was detected. Embryo-specific ablation in Arabidopsis was achieved by transactivation of a diphtheria toxin gene under the control of the EASE. The potential application of the EASE element and similar control elements as part of an open-source biotechnology toolkit for apomixis is discussed.

Interaction of NPR1 with basic leucine zipper protein transcription factors that bind sequences required for salicylic acid induction of the PR–1 gene

The Arabidopsis thaliana NPR1 has been shown to be a key regulator of gene expression during the onset of a plant disease-resistance response known as systemic acquired resistance. The npr1 mutant plants fail to respond to systemic acquired resistance-inducing signals such as salicylic acid (SA), or express SA-induced pathogenesis-related (PR) genes. Using NPR1 as bait in a yeast two-hybrid screen, we identified a subclass of transcription factors in the basic leucine zipper protein family (AHBP-1b and TGA6) and showed that they interact specifically in yeast and in vitro with NPR1. Point mutations that abolish the NPR1 function in A. thaliana also impair the interactions between NPR1 and the transcription factors in the yeast two-hybrid assay. Furthermore, a gel mobility shift assay showed that the purified transcription factor protein, AHBP-1b, binds specifically to an SA-responsive promoter element of the A. thaliana PR-1 gene. These data suggest that NPR1 may regulate PR-1 gene expression by interacting with a subclass of basic leucine zipper protein transcription factors.

A myosin from a higher plant has structural similarities to class V myosins

In plant cells, myosin is believed to be the molecular motor responsible for actin-based motility processes such as cytoplasmic streaming and directed vesicle transport. In an effort to characterize plant myosin, a cDNA encoding a myosin heavy chain was isolated from Arabidopsis thaliana. The predicted product of the MYA1 gene is 173 kDa and is structurally similar to the class V myosins. It is composed of the highly-conserved NH2-terminal "head" domain, a putative calmodulin-binding "neck" domain an alpha-helical coiled-coil domain, and a COOH-terminal domain. Northern blot analysis shows that the Arabidopsis MYA1 gene is expressed in all the major plant tissues (flower, leaf, root, and stem). We suggest that the MYA1 myosin may be involved in a general intracellular transport process in plant cells.

Multifractal characterisation and analysis of complex networks

Complex networks have been studied extensively due to their relevance to many real-world systems such as the world-wide web, the internet, biological and social systems. During the past two decades, studies of such networks in different fields have produced many significant results concerning their structures, topological properties, and dynamics. Three well-known properties of complex networks are scale-free degree distribution, small-world effect and self-similarity. The search for additional meaningful properties and the relationships among these properties is an active area of current research. This thesis investigates a newer aspect of complex networks, namely their multifractality, which is an extension of the concept of selfsimilarity. The first part of the thesis aims to confirm that the study of properties of complex networks can be expanded to a wider field including more complex weighted networks. Those real networks that have been shown to possess the self-similarity property in the existing literature are all unweighted networks. We use the proteinprotein interaction (PPI) networks as a key example to show that their weighted networks inherit the self-similarity from the original unweighted networks. Firstly, we confirm that the random sequential box-covering algorithm is an effective tool to compute the fractal dimension of complex networks. This is demonstrated on the Homo sapiens and E. coli PPI networks as well as their skeletons. Our results verify that the fractal dimension of the skeleton is smaller than that of the original network due to the shortest distance between nodes is larger in the skeleton, hence for a fixed box-size more boxes will be needed to cover the skeleton. Then we adopt the iterative scoring method to generate weighted PPI networks of five species, namely Homo sapiens, E. coli, yeast, C. elegans and Arabidopsis Thaliana. By using the random sequential box-covering algorithm, we calculate the fractal dimensions for both the original unweighted PPI networks and the generated weighted networks. The results show that self-similarity is still present in generated weighted PPI networks. This implication will be useful for our treatment of the networks in the third part of the thesis. The second part of the thesis aims to explore the multifractal behavior of different complex networks. Fractals such as the Cantor set, the Koch curve and the Sierspinski gasket are homogeneous since these fractals consist of a geometrical figure which repeats on an ever-reduced scale. Fractal analysis is a useful method for their study. However, real-world fractals are not homogeneous; there is rarely an identical motif repeated on all scales. Their singularity may vary on different subsets; implying that these objects are multifractal. Multifractal analysis is a useful way to systematically characterize the spatial heterogeneity of both theoretical and experimental fractal patterns. However, the tools for multifractal analysis of objects in Euclidean space are not suitable for complex networks. In this thesis, we propose a new box covering algorithm for multifractal analysis of complex networks. This algorithm is demonstrated in the computation of the generalized fractal dimensions of some theoretical networks, namely scale-free networks, small-world networks, random networks, and a kind of real networks, namely PPI networks of different species. Our main finding is the existence of multifractality in scale-free networks and PPI networks, while the multifractal behaviour is not confirmed for small-world networks and random networks. As another application, we generate gene interactions networks for patients and healthy people using the correlation coefficients between microarrays of different genes. Our results confirm the existence of multifractality in gene interactions networks. This multifractal analysis then provides a potentially useful tool for gene clustering and identification. The third part of the thesis aims to investigate the topological properties of networks constructed from time series. Characterizing complicated dynamics from time series is a fundamental problem of continuing interest in a wide variety of fields. Recent works indicate that complex network theory can be a powerful tool to analyse time series. Many existing methods for transforming time series into complex networks share a common feature: they define the connectivity of a complex network by the mutual proximity of different parts (e.g., individual states, state vectors, or cycles) of a single trajectory. In this thesis, we propose a new method to construct networks of time series: we define nodes by vectors of a certain length in the time series, and weight of edges between any two nodes by the Euclidean distance between the corresponding two vectors. We apply this method to build networks for fractional Brownian motions, whose long-range dependence is characterised by their Hurst exponent. We verify the validity of this method by showing that time series with stronger correlation, hence larger Hurst exponent, tend to have smaller fractal dimension, hence smoother sample paths. We then construct networks via the technique of horizontal visibility graph (HVG), which has been widely used recently. We confirm a known linear relationship between the Hurst exponent of fractional Brownian motion and the fractal dimension of the corresponding HVG network. In the first application, we apply our newly developed box-covering algorithm to calculate the generalized fractal dimensions of the HVG networks of fractional Brownian motions as well as those for binomial cascades and five bacterial genomes. The results confirm the monoscaling of fractional Brownian motion and the multifractality of the rest. As an additional application, we discuss the resilience of networks constructed from time series via two different approaches: visibility graph and horizontal visibility graph. Our finding is that the degree distribution of VG networks of fractional Brownian motions is scale-free (i.e., having a power law) meaning that one needs to destroy a large percentage of nodes before the network collapses into isolated parts; while for HVG networks of fractional Brownian motions, the degree distribution has exponential tails, implying that HVG networks would not survive the same kind of attack.