HIV Reverse Transcriptase Fidelity And Inhibition Are Modulated By Divalent Cations In A Concentration-Dependent Manner In Vitro

Human immunodeficiency virus (HIV) rapidly evolves through generation and selection of mutants that can escape drug therapy. This process is fueled, in part, by the presumably highly error prone polymerase reverse transcriptase (RT). Fidelity of polymerases can be influenced by cation co-factors. Physiologically, magnesium (Mg2+) is used as a co-factor by RT to perform catalysis, however, alternative cations including manganese (Mn2+), cobalt (Co2+), and zinc (Zn2+) can also be used. I demonstrate here that fidelity and inhibition of HIV RT can be influenced differently, in vitro, by divalent cations depending on their concentration. The reported mutation frequency for purified HIV RT in vitro is typically in the 10-4 range (per nucleotide addition), making the enzyme several-fold less accurate than most polymerases. Paradoxically, results examining HIV replication in cells indicate an error frequency that is ~10 times lower than the error rate obtained in the test tube. Here, I reconcile, at least in part, these discrepancies by showing that HIV RT fidelity in vitro is in the same range as cellular results, in physiological concentrations of free Mg2+ (~0.25 mM). At low Mg2+, mutation rates were 5-10 times lower compared to high Mg2+ conditions (5-10 mM). Alternative divalent cations also have a concentration-dependent effect on RT fidelity. Presumed promutagenic cations Mn2+ and Co2+ decreases the fidelity of RT only at elevated concentrations, and Zn2+, when present in low concentration, increases the fidelity of HIV-1 RT by ~2.5 fold compared to Mg2+. HIV-1 and HIV-2 RT inhibition by nucleoside (NRTIs) and non-nucleoside RT inhibitors (NNRTIs) in vitro is also affected by the Mg2+ concentration. NRTIs lacking 3'-OH group inhibited both enzymes less efficiently in low Mg2+ than in high Mg2+; whereas inhibition by the “translocation defective RT inhibitor”, which retains the 3ʹ-OH, was unaffected by Mg2+ concentration, suggesting that NRTIs with a 3ʹ-OH group may be more potent than other NRTIs. In contrast, NNRTIs were more effective in low vs. high Mg2+ conditions. Overall, the studies presented reveal strategies for designing novel RT inhibitors and strongly emphasize the need for studying HIV RT and RT inhibitors in physiologically relevant low Mg2+ conditions.

Dinoflagellate Genomic Organization and Phylogenetic Marker Discovery Utilizing Deep Sequencing Data

Dinoflagellates possess large genomes in which most genes are present in many copies. This has made studies of their genomic organization and phylogenetics challenging. Recent advances in sequencing technology have made deep sequencing of dinoflagellate transcriptomes feasible. This dissertation investigates the genomic organization of dinoflagellates to better understand the challenges of assembling dinoflagellate transcriptomic and genomic data from short read sequencing methods, and develops new techniques that utilize deep sequencing data to identify orthologous genes across a diverse set of taxa. To better understand the genomic organization of dinoflagellates, a genomic cosmid clone of the tandemly repeated gene Alchohol Dehydrogenase (AHD) was sequenced and analyzed. The organization of this clone was found to be counter to prevailing hypotheses of genomic organization in dinoflagellates. Further, a new non-canonical splicing motif was described that could greatly improve the automated modeling and annotation of genomic data. A custom phylogenetic marker discovery pipeline, incorporating methods that leverage the statistical power of large data sets was written. A case study on Stramenopiles was undertaken to test the utility in resolving relationships between known groups as well as the phylogenetic affinity of seven unknown taxa. The pipeline generated a set of 373 genes useful as phylogenetic markers that successfully resolved relationships among the major groups of Stramenopiles, and placed all unknown taxa on the tree with strong bootstrap support. This pipeline was then used to discover 668 genes useful as phylogenetic markers in dinoflagellates. Phylogenetic analysis of 58 dinoflagellates, using this set of markers, produced a phylogeny with good support of all branches. The Suessiales were found to be sister to the Peridinales. The Prorocentrales formed a monophyletic group with the Dinophysiales that was sister to the Gonyaulacales. The Gymnodinales was found to be paraphyletic, forming three monophyletic groups. While this pipeline was used to find phylogenetic markers, it will likely also be useful for finding orthologs of interest for other purposes, for the discovery of horizontally transferred genes, and for the separation of sequences in metagenomic data sets.

Shaping the Dicot Fruit: Molecular and Genomic Approaches to Fruit Development

The fruit is one of the most complex and important structures produced by flowering plants, and understanding the development and maturation process of fruits in different angiosperm species with diverse fruit structures is of immense interest. In the work presented here, molecular genetics and genomic analysis are used to explore the processes that form the fruit in two species: The model organism Arabidopsis and the diploid strawberry Fragaria vesca. One important basic question concerns the molecular genetic basis of fruit patterning. A long-standing model of Arabidopsis fruit (the gynoecium) patterning holds that auxin produced at the apex diffuses downward, forming a gradient that provides apical-basal positional information to specify different tissue types along the gynoecium’s length. The proposed gradient, however, has never been observed and the model appears inconsistent with a number of observations. I present a new, alternative model, wherein auxin acts to establish the adaxial-abaxial domains of the carpel primordia, which then ensures proper development of the final gynoecium. A second project utilizes genomics to identify genes that regulate fruit color by analyzing the genome sequences of Fragaria vesca, a species of wild strawberry. Shared and distinct SNPs among three F. vesca accessions were identified, providing a foundation for locating candidate mutations underlying phenotypic variations among different F. vesca accessions. Through systematic analysis of relevant SNP variants, a candidate SNP in FveMYB10 was identified that may underlie the fruit color in the yellow-fruited accessions, which was subsequently confirmed by functional assays. Our lab has previously generated extensive RNA-sequencing data that depict genome-scale gene expression profiles in F. vesca fruit and flower tissues at different developmental stages. To enhance the accessibility of this dataset, the web-based eFP software was adapted for this dataset, allowing visualization of gene expression in any tissues by user-initiated queries. Together, this thesis work proposes a well-supported new model of fruit patterning in Arabidopsis and provides further resources for F. vesca, including genome-wide variant lists and the ability to visualize gene expression. This work will facilitate future work linking traits of economic importance to specific genes and gaining novel insights into fruit patterning and development.

Exploring the Function and Regulation of Arabidopsis EIN2 in Ethylene Signaling

Ethylene is an essential plant hormone involved in nearly all stages of plant growth and development. EIN2 (ETHYLENE INSENSITIVE2) is a master positive regulator in the ethylene signaling pathway, consisting of an N-terminal domain and a C-terminal domain. The EIN2 N-terminal domain localizes to the endoplasmic reticulum (ER) membrane and shows sequence similarity to Nramp metal ion transporters. The cytosolic C-terminal domain is unique to plants and signals downstream. There have been several major gaps in our knowledge of EIN2 function. It was unknown how the ethylene signal gets relayed from the known upstream component CTR1 (CONSTITUTIVE RESPONSE1) a Ser/Thr kinase at the ER, to EIN2. How the ethylene signal was transduced from EIN2 to the next downstream component transcription factor EIN3 (ETHYLENE INSENSITIVE3) in the nucleus was also unknown. The N-terminal domain of EIN2 shows homology to Nramp metal ion transporters and whether EIN2 can also function as a metal transporter has been a question plaguing the ethylene field for almost two decades. Here, EIN2 was found to interact with the CTR1 protein kinase, leading to the discovery that CTR1 phosphorylates the C-terminal domain of EIN2 in Arabidopsis thaliana. Using tags at the termini of EIN2, it was deduced that in the presence of ethylene, the EIN2 C-terminal domain is cleaved and translocates into the nucleus, where it could somehow activate downstream ethylene responses. The EIN2 C-terminal domain interacts with nuclear proteins, RTE3 and EER5, which are components of the TREX-2 mRNA export complex, although the role of these interactions remains unclear. The EIN2 N-terminal domain was found to be capable of divalent metal transport when expressed in E. coli and S. cerevisiae leading to the hypothesis that metal transport plays a role in ethylene signaling. This hypothesis was tested using a novel missense allele, ein2 G36E, substituting a highly conserved residue that is required for metal transport in Nramp proteins. This G36E substitution did not disrupt metal ion transport of EIN2, but the ethylene insensitive phenotype of this mutant indicates that the EIN2 N-terminal domain is important for positively regulating the C-terminal domain. The defect of the ein2 G36E mutant does not prevent proper expression or subcellular localization, but might affect protein modifications. The ein2 G36E allele is partially dominant, mostly likely displaying haploinsufficiency. Overexpression of the EIN2 N-terminal domain in the ein2 G36E mutant did not rescue ethylene insensitivity, suggesting the N-terminal domain functions in cis to regulate the C-terminal domain. These findings advance our knowledge of EIN2, which is critical to understanding ethylene signaling.

Error, Bias, and Long-Branch Attraction in Data for Two Chloroplast Photosystem Genes in Seed Plants

Sequences of two chloroplast photosystem genes, psaA and psbB, together comprising about 3,500 bp, were obtained for all five major groups of extant seed plants and several outgroups among other vascular plants. Strongly supported, but significantly conflicting, phylogenetic signals were obtained in parsimony analyses from partitions of the data into first and second codon positions versus third positions. In the former, both genes agreed on a monophyletic gymnosperms, with Gnetales closely related to certain conifers. In the latter, Gnetales are inferred to be the sister group of all other seed plants, with gymnosperms paraphyletic. None of the data supported the modern ‘‘anthophyte hypothesis,’’ which places Gnetales as the sister group of flowering plants. A series of simulation studies were undertaken to examine the error rate for parsimony inference. Three kinds of errors were examined: random error, systematic bias (both properties of finite data sets), and statistical inconsistency owing to long-branch attraction (an asymptotic property). Parsimony reconstructions were extremely biased for third-position data for psbB. Regardless of the true underlying tree, a tree in which Gnetales are sister to all other seed plants was likely to be reconstructed for these data. None of the combinations of genes or partitions permits the anthophyte tree to be reconstructed with high probability. Simulations of progressively larger data sets indicate the existence of long-branch attraction (statistical inconsistency) for third-position psbB data if either the anthophyte tree or the gymnosperm tree is correct. This is also true for the anthophyte tree using either psaA third positions or psbB first and second positions. A factor contributing to bias and inconsistency is extremely short branches at the base of the seed plant radiation, coupled with extremely high rates in Gnetales and nonseed plant outgroups. M. J. Sanderson,* M. F. Wojciechowski,*† J.-M. Hu,* T. Sher Khan,* and S. G. Brady

Biotechnology education : topics of interest to students and teachers

This paper presents the findings of a survey that investigates the biotechnology topics of interest according to students and teachers for inclusion in biology lessons and reports on the similarities and differences in teachers’ and students’ biotechnology topics of interest. The study is of significance as biotechnology has been identified as a key area of technological and economic importance worldwide yet there is scant literature relating to teachers’ and students’ interests concerning biotechnology education topics. 500 students and their 15 teachers completed the survey. Interviews were conducted with 3 teachers and 60 students. Responses indicate there is a mismatch in the interests of students and teachers, and what they perceive as being possible topics for inclusion in biology and biotechnology lessons. Where teachers are provided with the freedom to design and assess their own units of work, this mismatch of interests causes problems. The study found students withdrawing from biology courses in post compulsory settings due to lack of interest, and perceived lack of relevance of the course. It is possible that this lack of agreement on topics of interest is a factor in the world wide decline of enrolments in the sciences.

Age-related maculopathy : linking aetiology and pathophysiological changes to the ischaemia hypothesis

Age-related maculopathy (ARM) has remained a challenging topic with respect to its aetiology, pathomechanisms, early detection and treatment since the late 19th century when it was first described as its own entity. ARM was previously considered an inflammatory disease, a degenerative disease, a tumor and as the result of choroidal hemodynamic disturbances and ischaemia. The latter processes have been repeatedly suggested to have a key role in its development and progression. In vivo experiments under hypoxic conditions could be models for the ischaemic deficits in ARM. Recent research has also linked ARM with gene polymorphisms. It is however unclear what triggers a person's gene susceptibility. In this manuscript, a linking hypothesis between aetiological factors including ischaemia and genetics and the development of early clinicopathological changes in ARM is proposed. New clinical psychophysical and electrophysiological tests are introduced that can detect ARM at an early stage. Models of early ARM based upon hemodynamic, photoreceptor and post-receptoral deficits are described and the mechanisms by which ischaemia may be involved as a final common pathway are considered. In neovascular age-related macular degeneration (neovascular AMD), ischaemia is thought to promote release of vascular endothelial growth factor (VEGF) which induces chorioretinal neovascularisation. VEGF is critical in the maintenance of the healthy choriocapillaris. In the final section of the manuscript the documentation of the effect of new anti-VEGF treatments on retinal function in neovascular AMD is critically viewed.

In vitro cultivation of adult human chondrocytes : importance of culture system, oxygen and zonal differences

Articular cartilage exhibits limited intrinsic regenerative capacity and focal tissue defects can lead to the development of osteoarthritis (OA), a painful and debilitating loss of cartilage tissue. In Australia, 1.4 million people are affected by OA and its prevalence is increasing in line with current demographics. As treatment options are limited, new therapeutic approaches are being investigated including biological resurfacing of joints with tissue-engineered cartilage. Despite some progress in the field, major challenges remain to be addressed for large scale clinical success. For example, large numbers of chondrogenic cells are required for cartilage formation, but chondrocytes lose their chondrogenic phenotype (dedifferentiate) during in vitro propagation. Additionally, the zonal organization of articular cartilage is critical for normal cartilage function, but development of zonal structure has been largely neglected in cartilage repair strategies. Therefore, we hypothesised that culture conditions for freshly isolated human articular chondrocytes from non-OA and OA sources can be improved by employing microcarrier cultures and a reduced oxygen environment and that oxygen is a critical factor in the maintenance of the zonal chondrocyte phenotype. Microcarriers have successfully been used to cultivate bovine chondrocytes, and offer a potential alternative for clinical expansion of human chondrocytes. We hypothesised that improved yields can be achieved by propagating human chondrocytes on microcarriers. We found that cells on microcarriers acquired a flattened, polygonal morphology and initially proliferated faster than monolayercultivated cells. However, microcarrier cultivation over four weeks did not improve growth rates or the chondrogenic potential of non-OA and OA human articular chondrocytes over conventional monolayer cultivation. Based on these observations, we aimed to optimise culture conditions by modifying oxygen tension, to more closely reflect the in vivo environment. We found that propagation at 5% oxygen tension (moderate hypoxia) did not improve proliferation or redifferentiation capacity of human osteoarthritic chondrocytes. Moderate hypoxia increased the expression of chondrogenic markers during redifferentiation. However, osteoarthritic chondrocytes cultivated on microcarriers exhibited lower expression levels of chondrogenic surface marker proteins and had at best equivalent redifferentiation capacities compared to monolayer-cultured cells. This suggests that monolayer culture with multiple passaging potentially selects for a subpopulation of cells with higher differentiation capacity, which are otherwise rare in osteoarthritic, aged cartilage. However, fibroblastic proteins were found to be highly expressed in all cultures of human osteoarthritic chondrocytes indicating the presence of a high proportion of dedifferentiated, senescent cells with a chondrocytic phenotype that was not rescued by moderate hypoxia. The different zones of cartilage support chondrocyte subpopulations, which exhibit characteristic protein expression and experience varying oxygen tensions. We, therefore, hypothesised that oxygen tension affects the zonal marker expression of human articular chondrocytes isolated from the different cartilage layers. We found that zonal chondrocytes maintained these phenotypic differences during in vitro cultivation. Low oxygen environments favoured the expression of the zonal marker proteoglycan 4 in superficial cells, most likely through the promotion of chondrogenesis. The putative zonal markers clusterin and cartilage intermediate layer protein were found to be expressed by all subpopulations of human osteoarthritic chondrocytes ex vivo and, thus, may not be reliable predictors of in vitro stratification using these clinically relevant cells. The findings in this thesis underline the importance of considering low oxygen conditions and zonal stratification when creating native-like cartilaginous constructs. We have not yet found the right cues to successfully cultivate clinically-relevant human osteoarthritic chondrocytes in vitro. A more thorough understanding of chondrocyte biology and the processes of chondrogenesis are required to ensure the clinical success of cartilage tissue engineering.

Novel analytical and numerical methods for solving fractional dynamical systems

During the past three decades, the subject of fractional calculus (that is, calculus of integrals and derivatives of arbitrary order) has gained considerable popularity and importance, mainly due to its demonstrated applications in numerous diverse and widespread fields in science and engineering. For example, fractional calculus has been successfully applied to problems in system biology, physics, chemistry and biochemistry, hydrology, medicine, and finance. In many cases these new fractional-order models are more adequate than the previously used integer-order models, because fractional derivatives and integrals enable the description of the memory and hereditary properties inherent in various materials and processes that are governed by anomalous diffusion. Hence, there is a growing need to find the solution behaviour of these fractional differential equations. However, the analytic solutions of most fractional differential equations generally cannot be obtained. As a consequence, approximate and numerical techniques are playing an important role in identifying the solution behaviour of such fractional equations and exploring their applications. The main objective of this thesis is to develop new effective numerical methods and supporting analysis, based on the finite difference and finite element methods, for solving time, space and time-space fractional dynamical systems involving fractional derivatives in one and two spatial dimensions. A series of five published papers and one manuscript in preparation will be presented on the solution of the space fractional diffusion equation, space fractional advectiondispersion equation, time and space fractional diffusion equation, time and space fractional Fokker-Planck equation with a linear or non-linear source term, and fractional cable equation involving two time fractional derivatives, respectively. One important contribution of this thesis is the demonstration of how to choose different approximation techniques for different fractional derivatives. Special attention has been paid to the Riesz space fractional derivative, due to its important application in the field of groundwater flow, system biology and finance. We present three numerical methods to approximate the Riesz space fractional derivative, namely the L1/ L2-approximation method, the standard/shifted Gr¨unwald method, and the matrix transform method (MTM). The first two methods are based on the finite difference method, while the MTM allows discretisation in space using either the finite difference or finite element methods. Furthermore, we prove the equivalence of the Riesz fractional derivative and the fractional Laplacian operator under homogeneous Dirichlet boundary conditions – a result that had not previously been established. This result justifies the aforementioned use of the MTM to approximate the Riesz fractional derivative. After spatial discretisation, the time-space fractional partial differential equation is transformed into a system of fractional-in-time differential equations. We then investigate numerical methods to handle time fractional derivatives, be they Caputo type or Riemann-Liouville type. This leads to new methods utilising either finite difference strategies or the Laplace transform method for advancing the solution in time. The stability and convergence of our proposed numerical methods are also investigated. Numerical experiments are carried out in support of our theoretical analysis. We also emphasise that the numerical methods we develop are applicable for many other types of fractional partial differential equations.

A functional screen identifies lateral transfer of beta-glucuronidase (gus) from bacteria to fungi

Lateral gene transfer (LGT) from prokaryotes to microbial eukaryotes is usually detected by chance through genome-sequencing projects. Here, we explore a different, hypothesis-driven approach. We show that the fitness advantage associated with the transferred gene, typically invoked only in retrospect, can be used to design a functional screen capable of identifying postulated LGT cases. We hypothesized that beta-glucuronidase (gus) genes may be prone to LGT from bacteria to fungi (thought to lack gus) because this would enable fungi to utilize glucuronides in vertebrate urine as a carbon source. Using an enrichment procedure based on a glucose-releasing glucuronide analog (cellobiouronic acid), we isolated two gus(+) ascomycete fungi from soils (Penicillium canescens and Scopulariopsis sp.). A phylogenetic analysis suggested that their gus genes, as well as the gus genes identified in genomic sequences of the ascomycetes Aspergillus nidulans and Gibberella zeae, had been introgressed laterally from high-GC gram(+) bacteria. Two such bacteria (Arthrobacter spp.), isolated together with the gus(+) fungi, appeared to be the descendants of a bacterial donor organism from which gus had been transferred to fungi. This scenario was independently supported by similar substrate affinities of the encoded beta-glucuronidases, the absence of introns from fungal gus genes, and the similarity between the signal peptide-encoding 5' extensions of some fungal gus genes and the Arthrobacter sequences upstream of gus. Differences in the sequences of the fungal 5' extensions suggested at least two separate introgression events after the divergence of the two main Euascomycete classes. We suggest that deposition of glucuronides on soils as a result of the colonization of land by vertebrates may have favored LGT of gus from bacteria to fungi in soils.

Candidate metastasis suppressor genes uncovered by array comparative genomic hybridization in a mouse allograft model of prostate cancer.

Background The purpose of this study was to identify candidate metastasis suppressor genes from a mouse allograft model of prostate cancer (NE-10). This allograft model originally developed metastases by twelve weeks after implantation in male athymic nude mice, but lost the ability to metastasize after a number of in vivo passages. We performed high resolution array comparative genomic hybridization on the metastasizing and non-metastasizing allografts to identify chromosome imbalances that differed between the two groups of tumors. Results This analysis uncovered a deletion on chromosome 2 that differed between the metastasizing and non-metastasizing tumors. Bioinformatics filters were employed to mine this region of the genome for candidate metastasis suppressor genes. Of the 146 known genes that reside within the region of interest on mouse chromosome 2, four candidate metastasis suppressor genes (Slc27a2, Mall, Snrpb, and Rassf2) were identified. Quantitative expression analysis confirmed decreased expression of these genes in the metastasizing compared to non-metastasizing tumors. Conclusion This study presents combined genomics and bioinformatics approaches for identifying potential metastasis suppressor genes. The genes identified here are candidates for further studies to determine their functional role in inhibiting metastases in the NE-10 allograft model and human prostate cancer.

Multiple evolutionary rate classes in animal genome evolution

The proportion of functional sequence in the human genome is currently a subject of debate. The most widely accepted figure is that approximately 5% is under purifying selection. In Drosophila, estimates are an order of magnitude higher, though this corresponds to a similar quantity of sequence. These estimates depend on the difference between the distribution of genomewide evolutionary rates and that observed in a subset of sequences presumed to be neutrally evolving. Motivated by the widening gap between these estimates and experimental evidence of genome function, especially in mammals, we developed a sensitive technique for evaluating such distributions and found that they are much more complex than previously apparent. We found strong evidence for at least nine well-resolved evolutionary rate classes in an alignment of four Drosophila species and at least seven classes in an alignment of four mammals, including human. We also identified at least three rate classes in human ancestral repeats. By positing that the largest of these ancestral repeat classes is neutrally evolving, we estimate that the proportion of nonneutrally evolving sequence is 30% of human ancestral repeats and 45% of the aligned portion of the genome. However, we also question whether any of the classes represent neutrally evolving sequences and argue that a plausible alternative is that they reflect variable structure-function constraints operating throughout the genomes of complex organisms.

Comparison of koala LPCoLN and human strains of Chlamydia pneumoniae highlights extended genetic diversity in the species

Background Chlamydia pneumoniae is a widespread pathogen causing upper and lower respiratory tract infections in addition to a range of other diseases in humans and animals. Previous whole genome analyses have focused on four essentially clonal (> 99% identity) C. pneumoniae human genomes (AR39, CWL029, J138 and TW183), providing relatively little insight into strain diversity and evolution of this species. Results We performed individual gene-by-gene comparisons of the recently sequenced C. pneumoniae koala genome and four C. pneumoniae human genomes to identify species-specific genes, and more importantly, to gain an insight into the genetic diversity and evolution of the species. We selected genes dispersed throughout the chromosome, representing genes that were specific to C. pneumoniae, genes with a demonstrated role in chlamydial biology and/or pathogenicity (n = 49), genes encoding nucleotide salvage or amino acid biosynthesis proteins (n = 6), and extrachromosomal elements (9 plasmid and 2 bacteriophage genes). Conclusions We have identified strain-specific differences and targets for detection of C. pneumoniae isolates from both human and animal origin. Such characterisation is necessary for an improved understanding of disease transmission and intervention.