Construct design for efficient, effective and high-throughput gene silencing in plants

Post-transcriptional silencing of plant genes using anti-sense or co-suppression constructs usually results in only a modest proportion of silenced individuals. Recent work has demonstrated the potential for constructs encoding self-complementary 'hairpin' RNA (hpRNA) to efficiently silence genes. In this study we examine design rules for efficient gene silencing, in terms of both the proportion of independent transgenic plants showing silencing, and the degree of silencing. Using hpRNA constructs containing sense/anti-sense arms ranging from 98 to 853 nt gave efficient silencing in a wide range of plant species, and inclusion of an intron in these constructs had a consistently enhancing effect. Intron-containing constructs (ihpRNA) generally gave 90-100% of independent transgenic plants showing silencing. The degree of silencing with these constructs was much greater than that obtained using either co-suppression or anti-sense constructs. We have made a generic vector, pHANNIBAL, that allows a simple, single PCR product from a gene of interest to be easily converted into a highly effective ihpRNA silencing construct. We have also created a high-throughput vector, pHELLSGATE, that should facilitate the cloning of gene libraries or large numbers of defined genes, such as those in EST collections, using an in vitro recombinase system. This system may facilitate the large-scale determination and discovery of plant gene functions in the same way as RNAi is being used to examine gene function in Caenorhabditis elegans.

The Anadara trapezia transcriptome : a resource for molluscan physiological genomics

In this study we undertook deep sequencing of the blood cockle, Anadara trapezia, transcriptome to generate genomic resources for future functional genomics analyses. Over 27 million high quality paired end reads were assembled into 75 024 contigs. Of these contigs, 29 013 (38.7%) received significant BLASTx hits and gene ontology (GO) terms were assigned to 13 718 of these sequences. This resourcewill facilitate physiological genomic studies to test the gene expression response of A. trapezia to various environmental stresses.

targetTB: A target identification pipeline for Mycobacterium tuberculosis through an interactome, reactome and genome-scale structural analysis

Background: Tuberculosis still remains one of the largest killer infectious diseases, warranting the identification of newer targets and drugs. Identification and validation of appropriate targets for designing drugs are critical steps in drug discovery, which are at present major bottle-necks. A majority of drugs in current clinical use for many diseases have been designed without the knowledge of the targets, perhaps because standard methodologies to identify such targets in a high-throughput fashion do not really exist. With different kinds of 'omics' data that are now available, computational approaches can be powerful means of obtaining short-lists of possible targets for further experimental validation. Results: We report a comprehensive in silico target identification pipeline, targetTB, for Mycobacterium tuberculosis. The pipeline incorporates a network analysis of the protein-protein interactome, a flux balance analysis of the reactome, experimentally derived phenotype essentiality data, sequence analyses and a structural assessment of targetability, using novel algorithms recently developed by us. Using flux balance analysis and network analysis, proteins critical for survival of M. tuberculosis are first identified, followed by comparative genomics with the host, finally incorporating a novel structural analysis of the binding sites to assess the feasibility of a protein as a target. Further analyses include correlation with expression data and non-similarity to gut flora proteins as well as 'anti-targets' in the host, leading to the identification of 451 high-confidence targets. Through phylogenetic profiling against 228 pathogen genomes, shortlisted targets have been further explored to identify broad-spectrum antibiotic targets, while also identifying those specific to tuberculosis. Targets that address mycobacterial persistence and drug resistance mechanisms are also analysed. Conclusion: The pipeline developed provides rational schema for drug target identification that are likely to have high rates of success, which is expected to save enormous amounts of money, resources and time in the drug discovery process. A thorough comparison with previously suggested targets in the literature demonstrates the usefulness of the integrated approach used in our study, highlighting the importance of systems-level analyses in particular. The method has the potential to be used as a general strategy for target identification and validation and hence significantly impact most drug discovery programmes.

Host-specificity of Salmonella enterica serovar Gallinarum: Insights from comparative genomics

In this study, we have identified the possible genetic factors responsible for fowl-adaptation of Salmonella enterica serovar Gallinarum (S. Gallinarum). By comparing the genes related to Salmonella pathogenicity islands (SPI) of S. Gallinarum with those of Salmonella enterica serovar Enteritidis (S. Enteritidis) we have identified twenty-four positively selected genes. Our results suggest that the genes encoding the structural components of SPI-2 encoded type three secretion apparatus (TTSS) and the effector proteins that are secreted via SPI-1 encoded TTSS have evolved under positive selection pressure in these serovars. We propose that these positively selected genes play important roles in conferring different host-specificities to S. Gallinarum and S. Enteritidis.

Relationship Between the Levels of Non-structural Carbohydrates, Digging Date, Nursery-growing Environment, and Chilling in Strawberry Transplants in a Subtropical Environment

Experiments were conducted to study the effect of time of digging and nursery-growing environment on the levels of non-structural carbohydrates in 'Festival' strawberry transplants (Fragaria xananassa) over 2 years in southeastern Queensland, Australia. We were interested in determining whether there was a strong relationship between the potential productivity of this material and reserves in the plants. First, bare-rooted plants were obtained from Stanthorpe in southern Queensland from early March to mid-April/late April. Second, bare-rooted plants were sourced from Stanthorpe (a warm-growing area) or from Toolangi in Victoria (a cool-growing area). In Year 1 of the experiments, the nursery material from the different treatments was grown at Nambour in southeastern Queensland and fruit yield determined. The total weight of nonstructural carbohydrates/plant increased as digging was delayed and was higher in the plants from Stanthorpe than the plants from Toolangi. Plants dug on 17 Mar. in Year 1 had higher weights of non-structural carbohydrates [292 mg/plant dry weight (DW)] than plants dug on 3 Mar. (224 mg/plant) and higher early yield to the end of June or to the end of July and higher total yield to mid-October adjusted by the length of the growing season for the different treatments. Plants dug on 1 Apr. (408 mg/plant) or on 13 Apr. (445 mg/plant) had higher reserves than the plants dug on 17 Mar. but lower yields. Only the differences in yields between the plants dug on 3 Mar. and 17 Mar. reflected the differences in carbohydrates. The stock from Stanthorpe had greater reserves (408 mg/plant) than the stock from Toolangi (306 mg/plant) but similar yields in Year 1 possibly because of poorer flowering in the nursery plants. It was concluded that carbohydrate reserves in transplants only partially reflect their productivity in this environment.

Tripogon loliiformis elicits a rapid physiological and structural response to dehydration for desiccation tolerance

Resurrection plants can withstand extreme dehydration to an air-dry state and then recover upon receiving water. Tripogon loliiformis (F.Muell.) C.E.Hubb. is a largely uncharacterised native Australian desiccation-tolerant grass that resurrects from the desiccated state within 72 h. Using a combination of structural and physiological techniques the structural and physiological features that enable T. loliiformis to tolerate desiccation were investigated. These features include: - (i) a myriad of structural changes such as leaf folding, cell wall folding and vacuole fragmentation that mitigate desiccation stress; - (ii) potential role of sclerenchymatous tissue within leaf folding and radiation protection; - (iii) retention of ~70% chlorophyll in the desiccated state; - (iv) early response of photosynthesis to dehydration by 50% reduction and ceasing completely at 80 and 70% relative water content, respectively; - (v) a sharp increase in electrolyte leakage during dehydration, and; - (vi) confirmation of membrane integrity throughout desiccation and rehydration. Taken together, these results demonstrate that T. loliiformis implements a range of structural and physiological mechanisms that minimise mechanical, oxidative and irradiation stress. These results provide powerful insights into tolerance mechanisms for potential utilisation in the enhancement of stress-tolerance in crop plants.

Arabidopsis thaliana J-class heat shock proteins: cellular stress sensors

Plants are sessile organisms that have evolved a variety of mechanisms to maintain their cellular homeostasis under stressful environmental conditions. Survival of plants under abiotic stress conditions requires specialized group of heat shock protein machinery, belonging to Hsp70:J-protein family. These heat shock proteins are most ubiquitous types of chaperone machineries involved in diverse cellular processes including protein folding, translocation across cell membranes, and protein degradation. They play a crucial role in maintaining the protein homeostasis by reestablishing functional native conformations under environmental stress conditions, thus providing protection to the cell. J-proteins are co-chaperones of Hsp70 machine, which play a critical role by stimulating Hsp70s ATPase activity, thereby stabilizing its interaction with client proteins. Using genome-wide analysis of Arabidopsis thaliana, here we have outlined identification and systematic classification of J-protein co-chaperones which are key regulators of Hsp70s function. In comparison with Saccharomyces cerevisiae model system, a comprehensive domain structural organization, cellular localization, and functional diversity of A. thaliana J-proteins have also been summarized. Electronic supplementary material The online version of this article (doi:10.1007/s10142-009-0132-0) contains supplementary material, which is available to authorized users.

Production of F4 Fimbrial Adhesin in Plants: a Model for Oral Porcine Vaccine against Enterotoxigenic Escherichia coli

F4 fimbriae of enterotoxigenic Escherichia coli (ETEC) are highly stable multimeric structures with a capacity to evoke mucosal immune responses. With these characters F4 offer a unique model system to study oral vaccination against ETEC-induced porcine postweaning diarrhea. Postweaning diarrhea is a major problem in piggeries worldwide and results in significant economic losses. No vaccine is currently available to protect weaned piglets against ETEC infections. Transgenic plants provide an economically feasible platform for large-scale production of vaccine antigens for animal health. In this study, the capacity of transgenic plants to produce FaeG protein, the major structural subunit and adhesin of F4 fimbria, was evaluated. Using the model plant tobacco, the optimal subcellular location for FaeG accumulation was examined. Targeting of FaeG into chloroplasts offered a superior accumulation level of 1% of total soluble proteins (TSP) over the other investigated subcellular locations, namely, the endoplasmic reticulum and the apoplast. Moreover, we determined whether the FaeG protein, when isolated from its fimbrial background and produced in a plant cell, would retain the key properties of an oral vaccine, i.e. stability in gastrointestinal conditions, binding to porcine intestinal F4 receptors (F4R), and inhibition of the F4-possessing (F4+) ETEC attachment to F4R. The chloroplast-derived FaeG protein did show resistance against low pH and proteolysis in the simulated gastrointestinal conditions and was able to bind to the F4R, subsequently inhibiting the F4+ ETEC binding in a dose-dependent manner. To investigate the oral immunogenicity of FaeG protein, the edible crop plant alfalfa was transformed with the chloroplast-targeting construct and equally to tobacco plants, a high-yield FaeG accumulation of 1% of TSP was obtained. A similar yield was also obtained in the seeds of barley, a valuable crop plant, when the FaeG-encoding gene was expressed under an endosperm-specific promoter and subcellularly targeted into the endoplasmic reticulum. Furthermore, desiccated alfalfa plants and barley grains were shown to have a capacity to store FaeG protein in a stable form for years. When the transgenic alfalfa plants were administred orally to weaned piglets, slight F4-specific systemic and mucosal immune responses were induced. Co-administration of the transgenic alfalfa and the mucosal adjuvant cholera toxin enhanced the F4-specific immune response; the duration and number of F4+ E. coli excretion following F4+ ETEC challenge were significantly reduced as compared with pigs that had received nontransgenic plant material. In conclusion, the results suggest that transgenic plants producing the FaeG subunit protein could be used for production and delivery of oral vaccines against porcine F4+ ETEC infections. The findings here thus present new approaches to develop the vaccination strategy against porcine postweaning diarrhea.

Structural studies on tobacco streak virus coat protein: Insights into the pleomorphic nature of ilarviruses

Tobacco streak virus (TSV), the type member of Ilarvirus genus, is a major plant pathogen. TSV purified from infected plants consists of a ss-RNA genome encapsidated in spheroidal particles with diameters of 27, 30 and 33 nm constructed from multiple copies of a single species of coat protein (CP) subunits. Apart from protecting the viral genome, CPs of ilarviruses play several key roles in the life cycle of these viruses. Unlike the related bromo and cucumoviruses, ilarvirus particles are labile and pleomorphic, which has posed difficulties in their crystallization and structure determination. In the current study, a truncated TSV-CP was crystallized in two distinct forms and their structures were determined at resolutions of 2.4 angstrom and 2.1 angstrom, respectively. The core of TSV CP was found to possess the canonical beta-barrel jelly roll tertiary structure observed in several other viruses. Dimers of CP with swapped C-terminal arms (C-arm) were observed in both the crystal forms. The C-arm was found to be flexible and is likely to be responsible for the polymorphic and pleomorphic nature of TSV capsids. Consistent with this observation, mutations in the hinge region of the C-arm that reduce the flexibility resulted in the formation of more uniform particles. TSV CP was found to be structurally similar to that of Alfalfa mosaic virus (AMV) accounting for similar mechanism of genome activation in alfamo and ilar viruses. This communication represents the first report on the structure of the CP from an ilarvirus. (C) 2015 Elsevier Inc. All rights reserved.

Functional genomic studies of the structure and regulation of eukaryotic transcriptomes

The main focus of this thesis is the use of high-throughput sequencing technologies in functional genomics (in particular in the form of ChIP-seq, chromatin immunoprecipitation coupled with sequencing, and RNA-seq) and the study of the structure and regulation of transcriptomes. Some parts of it are of a more methodological nature while others describe the application of these functional genomic tools to address various biological problems. A significant part of the research presented here was conducted as part of the ENCODE (ENCyclopedia Of DNA Elements) Project.

The first part of the thesis focuses on the structure and diversity of the human transcriptome. Chapter 1 contains an analysis of the diversity of the human polyadenylated transcriptome based on RNA-seq data generated for the ENCODE Project. Chapter 2 presents a simulation-based examination of the performance of some of the most popular computational tools used to assemble and quantify transcriptomes. Chapter 3 includes a study of variation in gene expression, alternative splicing and allelic expression bias on the single-cell level and on a genome-wide scale in human lymphoblastoid cells; it also brings forward a number of critical to the practice of single-cell RNA-seq measurements methodological considerations.

The second part presents several studies applying functional genomic tools to the study of the regulatory biology of organellar genomes, primarily in mammals but also in plants. Chapter 5 contains an analysis of the occupancy of the human mitochondrial genome by TFAM, an important structural and regulatory protein in mitochondria, using ChIP-seq. In Chapter 6, the mitochondrial DNA occupancy of the TFB2M transcriptional regulator, the MTERF termination factor, and the mitochondrial RNA and DNA polymerases is characterized. Chapter 7 consists of an investigation into the curious phenomenon of the physical association of nuclear transcription factors with mitochondrial DNA, based on the diverse collections of transcription factor ChIP-seq datasets generated by the ENCODE, mouseENCODE and modENCODE consortia. In Chapter 8 this line of research is further extended to existing publicly available ChIP-seq datasets in plants and their mitochondrial and plastid genomes.

The third part is dedicated to the analytical and experimental practice of ChIP-seq. As part of the ENCODE Project, a set of metrics for assessing the quality of ChIP-seq experiments was developed, and the results of this activity are presented in Chapter 9. These metrics were later used to carry out a global analysis of ChIP-seq quality in the published literature (Chapter 10). In Chapter 11, the development and initial application of an automated robotic ChIP-seq (in which these metrics also played a major role) is presented.

The fourth part presents the results of some additional projects the author has been involved in, including the study of the role of the Piwi protein in the transcriptional regulation of transposon expression in Drosophila (Chapter 12), and the use of single-cell RNA-seq to characterize the heterogeneity of gene expression during cellular reprogramming (Chapter 13).

The last part of the thesis provides a review of the results of the ENCODE Project and the interpretation of the complexity of the biochemical activity exhibited by mammalian genomes that they have revealed (Chapters 15 and 16), an overview of the expected in the near future technical developments and their impact on the field of functional genomics (Chapter 14), and a discussion of some so far insufficiently explored research areas, the future study of which will, in the opinion of the author, provide deep insights into many fundamental but not yet completely answered questions about the transcriptional biology of eukaryotes and its regulation.

Structural and Biophysical Characterization of Variants of the Mechanosensitive Channel of Large Conductance (MSCL)

The ability to sense mechanical force is vital to all organisms to interact with and respond to stimuli in their environment. Mechanosensation is critical to many physiological functions such as the senses of hearing and touch in animals, gravitropism in plants and osmoregulation in bacteria. Of these processes, the best understood at the molecular level involve bacterial mechanosensitive channels. Under hypo-osmotic stress, bacteria are able to alleviate turgor pressure through mechanosensitive channels that gate directly in response to tension in the membrane lipid bilayer. A key participant in this response is the mechanosensitive channel of large conductance (MscL), a non-selective channel with a high conductance of ~3 nS that gates at tensions close to the membrane lytic tension.

It has been appreciated since the original discovery by C. Kung that the small subunit size (~130 to 160 residues) and the high conductance necessitate that MscL forms a homo-oligomeric channel. Over the past 20 years of study, the proposed oligomeric state of MscL has ranged from monomer to hexamer. Oligomeric state has been shown to vary between MscL homologues and is influenced by lipid/detergent environment. In this thesis, we report the creation of a chimera library to systematically survey the correlation between MscL sequence and oligomeric state to identify the sequence determinants of oligomeric state. Our results demonstrate that although there is no combination of sequences uniquely associated with a given oligomeric state (or mixture of oligomeric states), there are significant correlations. In the quest to characterize the oligomeric state of MscL, an exciting discovery was made about the dynamic nature of the MscL complex. We found that in detergent solution, under mild heating conditions (37 °C – 60 °C), subunits of MscL can exchange between complexes, and the dynamics of this process are sensitive to the protein sequence.

Extensive efforts were made to produce high diffraction quality crystals of MscL for the determination of a high resolution X-ray crystal structure of a full length channel. The surface entropy reduction strategy was applied to the design of S. aureus MscL variants and while the strategy appears to have improved the crystallizability of S. aureus MscL, unfortunately the diffraction qualities of these crystals were not significantly improved. MscL chimeras were also screened for crystallization in various solubilization detergents, but also failed to yield high quality crystals.

MscL is a fascinating protein and continues to serve as a model system for the study of the structural and functional properties of mechanosensitive channels. Further characterization of the MscL chimera library will offer more insight into the characteristics of the channel. Of particular interest are the functional characterization of the chimeras and the exploration of the physiological relevance of intercomplex subunit exchange.

Structural analysis of peptides of PSII light-harvesting complexes in siphonous green algae, Codium fragile

peptide composition and arrangement of 4 major light-harvesting complexes LHCP1-3 and LHCP3, isolated from siphonous green algae (Codium fragile (Sur.) Hariot.) were investigated. LHCP1 showed five main peptides, 34.4, 31.5, 29.5, 28.2 and 26.5 kD in SDS-PAGE, the 34.4 and 31.5 kD peptides were never found in higher plants. LHCP3 contained the other four kinds of LHCP1 peptides except 34.4 kD, while LHCP3, consisted of only 28.2 and 26.5 kD peptides. We found that 34.4, 28.2 and 26.5 kD peptides were easy to decompose from LHCP1 when subjected to SDS-PACE without pretreatment. They might be located at the exterior of LHCP1, while the 31.5 and 29.5 kD peptides were at the central part. The 28.2 and 26.5 kD peptides often occurred in CPa, the center complex of PS II. They are possibly the LHC II peptides tightly associated with CC II. According to the results described above, a peptide map of LHCP1 was sketched.

Brachypodium distachyon. A new model system for functional genomics in grasses

John Draper, Luis A.J. Mur, Glyn Jenkins, Gadab C. Ghosh-Biswas, Pauline Bablak, Robert Hasterok,and Andrew P.M. Routledge (2001). Brachypodium distachyon. A new model system for functional genomics in grasses. Plant Physiology, 127 (4), 1539-1555. Sponsorship: BBSRC / Gatsby Foundation RAE2008

Bayesian variable selection in structured high-dimensional covariate spaces with applications in genomics

We consider the problem of variable selection in regression modeling in high-dimensional spaces where there is known structure among the covariates. This is an unconventional variable selection problem for two reasons: (1) The dimension of the covariate space is comparable, and often much larger, than the number of subjects in the study, and (2) the covariate space is highly structured, and in some cases it is desirable to incorporate this structural information in to the model building process. We approach this problem through the Bayesian variable selection framework, where we assume that the covariates lie on an undirected graph and formulate an Ising prior on the model space for incorporating structural information. Certain computational and statistical problems arise that are unique to such high-dimensional, structured settings, the most interesting being the phenomenon of phase transitions. We propose theoretical and computational schemes to mitigate these problems. We illustrate our methods on two different graph structures: the linear chain and the regular graph of degree k. Finally, we use our methods to study a specific application in genomics: the modeling of transcription factor binding sites in DNA sequences. © 2010 American Statistical Association.

Investigation into mercury bound to biothiols:structural identification using ESI-ion-trap MS and introduction of a method for their HPLC separation with simultaneous detection by ICP-MS and ESI-MS

Mercury in plants or animal tissue is supposed to occur in the form of complexes formed with biologically relevant thiols (biothiols), rather than as free cation. We describe a technique for the separation and molecular identification of mercury and methylmercury complexes derived from their reactions with cysteine (Cys) and glutathione (GS): Hg(Cys)(2), Hg(GS)(2), MeHgCys, MeHgGS. Complexes were characterised by electrospray mass spectrometry (MS) equipped with an ion trap and the fragmentation pattern of MeHgCys was explained by using MP2 and B3LYP calculations, showing the importance of mercury-amine interactions in the gas phase. Chromatographic baseline separation was performed within 10 min with formic acid as the mobile phase on a reversed-phase column. Detection was done by online simultaneous coupling of ES-MS and inductively coupled plasma MS. When the mercury complexes were spiked in real samples (plant extracts), no perturbation of the separation and detection conditions was observed, suggesting that this method is capable of detecting mercury biothiol complexes in plants.