The GP problem: quantifying gene-to-phenotype relationships

In this paper we refer to the gene-to-phenotype modeling challenge as the GP problem. Integrating information across levels of organization within a genotype-environment system is a major challenge in computational biology. However, resolving the GP problem is a fundamental requirement if we are to understand and predict phenotypes given knowledge of the genome and model dynamic properties of biological systems. Organisms are consequences of this integration, and it is a major property of biological systems that underlies the responses we observe. We discuss the E(NK) model as a framework for investigation of the GP problem and the prediction of system properties at different levels of organization. We apply this quantitative framework to an investigation of the processes involved in genetic improvement of plants for agriculture. In our analysis, N genes determine the genetic variation for a set of traits that are responsible for plant adaptation to E environment-types within a target population of environments. The N genes can interact in epistatic NK gene-networks through the way that they influence plant growth and development processes within a dynamic crop growth model. We use a sorghum crop growth model, available within the APSIM agricultural production systems simulation model, to integrate the gene-environment interactions that occur during growth and development and to predict genotype-to-phenotype relationships for a given E(NK) model. Directional selection is then applied to the population of genotypes, based on their predicted phenotypes, to simulate the dynamic aspects of genetic improvement by a plant-breeding program. The outcomes of the simulated breeding are evaluated across cycles of selection in terms of the changes in allele frequencies for the N genes and the genotypic and phenotypic values of the populations of genotypes.

Trait physiology and crop modelling as a framework to link phenotypic complexity to underlying genetic systems

New tools derived from advances in molecular biology have not been widely adopted in plant breeding for complex traits because of the inability to connect information at gene level to the phenotype in a manner that is useful for selection. In this study, we explored whether physiological dissection and integrative modelling of complex traits could link phenotype complexity to underlying genetic systems in a way that enhanced the power of molecular breeding strategies. A crop and breeding system simulation study on sorghum, which involved variation in 4 key adaptive traits-phenology, osmotic adjustment, transpiration efficiency, stay-green-and a broad range of production environments in north-eastern Australia, was used. The full matrix of simulated phenotypes, which consisted of 547 location-season combinations and 4235 genotypic expression states, was analysed for genetic and environmental effects. The analysis was conducted in stages assuming gradually increased understanding of gene-to-phenotype relationships, which would arise from physiological dissection and modelling. It was found that environmental characterisation and physiological knowledge helped to explain and unravel gene and environment context dependencies in the data. Based on the analyses of gene effects, a range of marker-assisted selection breeding strategies was simulated. It was shown that the inclusion of knowledge resulting from trait physiology and modelling generated an enhanced rate of yield advance over cycles of selection. This occurred because the knowledge associated with component trait physiology and extrapolation to the target population of environments by modelling removed confounding effects associated with environment and gene context dependencies for the markers used. Developing and implementing this gene-to-phenotype capability in crop improvement requires enhanced attention to phenotyping, ecophysiological modelling, and validation studies to test the stability of candidate genetic regions.

Trait physiology and crop modelling to link phenotypic complexity to underlying genetic systems

New tools derived from advances in molecular biology have not been widely adopted in plant breeding because of the inability to connect information at gene level to the phenotype in a manner that is useful for selection. We explore whether a crop growth and development modelling framework can link phenotype complexity to underlying genetic systems in a way that strengthens molecular breeding strategies. We use gene-to-phenotype simulation studies on sorghum to consider the value to marker-assisted selection of intrinsically stable QTLs that might be generated by physiological dissection of complex traits. The consequences on grain yield of genetic variation in four key adaptive traits – phenology, osmotic adjustment, transpiration efficiency, and staygreen – were simulated for a diverse set of environments by placing the known extent of genetic variation in the context of the physiological determinants framework of a crop growth and development model. It was assumed that the three to five genes associated with each trait, had two alleles per locus acting in an additive manner. The effects on average simulated yield, generated by differing combinations of positive alleles for the traits incorporated, varied with environment type. The full matrix of simulated phenotypes, which consisted of 547 location-season combinations and 4235 genotypic expression states, was analysed for genetic and environmental effects. The analysis was conducted in stages with gradually increased understanding of gene-to-phenotype relationships, which would arise from physiological dissection and modelling. It was found that environmental characterisation and physiological knowledge helped to explain and unravel gene and environment context dependencies. We simulated a marker-assisted selection (MAS) breeding strategy based on the analyses of gene effects. When marker scores were allocated based on the contribution of gene effects to yield in a single environment, there was a wide divergence in rate of yield gain over all environments with breeding cycle depending on the environment chosen for the QTL analysis. It was suggested that knowledge resulting from trait physiology and modelling would overcome this dependency by identifying stable QTLs. The improved predictive power would increase the utility of the QTLs in MAS. Developing and implementing this gene-to-phenotype capability in crop improvement requires enhanced attention to phenotyping, ecophysiological modelling, and validation studies to test the stability of candidate QTLs.

An antisense approach to phenotype-based gene cloning in Tetrahymena

We report a pioneering approach using Tetrahymena thermophila that permits rapid identification of genes based on their null or hypomorphic phenotypes. This technique involves cell transformation with a library of plasmids that encode 26S ribosomal subunits containing short insertions. The insertions correspond to antisense sequences for a large number of genes. The majority of cells each acquires a single antisense sequence, which silences a single genomic locus. Because the insertion site within the ribosomal sequence is known, the silenced gene is easily amplified. We demonstrate that this approach can be used to identify genes required for dense core granule exocytosis.

Recurrent mutations in the CDKL5 gene: Genotype-phenotype relationships.

Mutations in the cyclin-dependent kinase-like 5 gene (CDKL5) have been described in epileptic encephalopathies in females with infantile spasms with features that overlap with Rett syndrome. With more than 80 reported patients, the phenotype of CDKL5-related encephalopathy is well-defined. The main features consist of seizures starting before 6 months of age, severe intellectual disability with absent speech and hand stereotypies and deceleration of head growth, which resembles Rett syndrome. However, some clinical discrepancies suggested the influence of genetics and/or environmental factors. No genotype-phenotype correlation has been defined and thus there is a need to examine individual mutations. In this study, we analyzed eight recurrent CDKL5 mutations to test whether the clinical phenotype of patients with the same mutation is similar and whether patients with specific CDKL5 mutations have a milder phenotype than those with other CDKL5 mutations. Patients bearing missense mutations in the ATP binding site such as the p.Ala40Val mutation typically walked unaided, had normocephaly, better hand use ability, and less frequent refractory epilepsy when compared to girls with other CDKL5 mutations. In contrast, patients with mutations in the kinase domain (such as p.Arg59X, p.Arg134X, p.Arg178Trp/Pro/Gln, or c.145 + 2T > C) and frameshift mutations in the C-terminal region (such as c.2635_2636delCT) had a more severe phenotype with infantile spasms, refractory epileptic encephalopathy, absolute microcephaly, and inability to walk. It is important for clinicians to have this information when such patients are diagnosed. © 2012 Wiley Periodicals, Inc.

Inducible caspase 9 suicide gene to improve the safety of allodepleted T cells after haploidentical stem cell transplantation.

Addback of donor T cells following T cell-depleted stem cell transplantation (SCT) can accelerate immune reconstitution and be effective against relapsed malignancy. After haploidentical SCT, a high risk of graft-versus-host disease (GVHD) essentially precludes this option, unless the T cells are first depleted of alloreactive precursor cells. Even then, the risks of severe GVHD remain significant. To increase the safety of the approach and thereby permit administration of larger T cell doses, we used a suicide gene, inducible caspase 9 (iCasp9), to transduce allodepleted T cells, permitting their destruction should administration have adverse effects. We made a retroviral vector encoding iCasp9 and a selectable marker (truncated CD19). Even after allodepletion (using anti-CD25 immunotoxin), donor T cells could be efficiently transduced, expanded, and subsequently enriched by CD19 immunomagnetic selection to >90% purity. These engineered cells retained antiviral specificity and functionality, and contained a subset with regulatory phenotype and function. Activating iCasp9 with a small-molecule dimerizer rapidly produced >90% apoptosis. Although transgene expression was downregulated in quiescent T cells, iCasp9 remained an efficient suicide gene, as expression was rapidly upregulated in activated (alloreactive) T cells. We have demonstrated the clinical feasibility of this approach after haploidentical transplantation by scaling up production using clinical grade materials.

From gene to protein: Investigating the disease etiology of retinitis pigmentosa

Retinitis pigmentosa (RP) is a name given to a group of inherited retinal dystrophies that lead to progressive photoreceptor degeneration, and thus, visual impairment. It is evident at both the clinical and the molecular level that these are heterogeneous disorders, with wide variation in severity, mode of inheritance, and phenotype. The genetics of RP are not simple; the disease can be inherited in dominant, recessive, X-linked, and digenic modes. Autosomal dominant RP (adRP) results from mutations in at least ten mapped loci, but there may be dozens of genetic loci where mutations can cause RP. To date, there are over a hundred genes known to cause retinal degenerative diseases, and less than half of these have been cloned (RetNet). Among the dozens of retinitis pigmentosa loci known to exist, only a few have been identified and the remainders are inferred from linkage studies. Today, the genes for seven of the twelve-adRP loci have been identified, and these are rhodopsin, peripherin/RDS, NRL, ROM1, CRX, RP13 and RP1. My research projects involved a combination of the continued search for genes involved in retinal dystrophies, as well the investigation into the role of peripherin/RDS and RP1 in the disease etiology of autosomal dominant RP. ^ Most of the mutations leading to inherited retinal disorders have been identified in predominately retina expressed genes like rhodopsin, peripherin/RDS, and RP1. Expressed sequence tags (ESTs) that were retina-specific were culled from sequence databases and, together with laboratory analysis, were analyzed as potential candidate genes for retinal dystrophies. Thirteen of the fifty-five identified retina-specific ESTs mapped to within candidate regions for inherited retinopathies. One of these is RP1L1, a homologue of RP1 and a potential cause of adRP. ^ Once a disease-associated gene has been identified, elucidating the role of that gene in the visual process is essential for understanding what happens when the process is defective as it is in adRP. My next projects involved investigating the role of a novel 5′ donor +3 splice site mutation on the mRNA of peripherin/RDS in adRP affected individuals, and comparative sequencing in RP1 to define conserved regions of the protein. Comparative sequencing is a powerful way to delineate critical regions of a sequence because different regions of a gene have different functions, and each region is subject to different levels of functional or structural constraints. Establishing a framework of conserved domains is beneficial not only for structural or functional studies, but can also aid in determining the potential effects of mutations. With the completion of sequencing of human genome, and other organisms such as Saccharomyces cerevisiae, Caenorhabditis elegans , and Drosophila, the facility of comparative sequencing will only increase in the future. Comparative sequencing has already become an established procedure for pinpointing conserved regions of a protein, and is an efficient way to target regions of a protein for experimental and/or evolutionary analysis. ^

A simple approach to numerical modelling of propane combustion in fluidized beds

A mathematical model is proposed for the evolution of temperature, chemical composition, and energy release in bubbles, clouds, and emulsion phase during combustion of gaseous premixtures of air and propane in a bubbling fluidized bed. The analysis begins as the bubbles are formed at the orifices of the distributor, until they explode inside the bed or emerge at the free surface of the bed. The model also considers the freeboard region of the fluidized bed until the propane is thoroughly burned. It is essentially built upon the quasi-global mechanism of Hautman et al. (1981) and the mass and heat transfer equations from the two-phase model of Davidson and Harrison (1963). The focus is not on a new modeling approach, but on combining the classical models of the kinetics and other diffusional aspects to obtain a better insight into the events occurring inside a fluidized bed reactor. Experimental data are obtained to validate the model by testing the combustion of commercial propane, in a laboratory-scale fluidized bed, using four sand particle sizes: 400–500, 315–400, 250–315, and 200–250 µm. The mole fractions of CO2, CO, and O2 in the flue gases and the temperature of the fluidized bed are measured and compared with the numerical results.

Transformation of Medicago truncatula with the arginine decarboxylase gene to modify polyamine metabolism toward water deficit resistance

Dissertation submitted to obtain a Ph.D. (Doutoramento) degree in Biology at the Instituto de Tecnologia Química e Biológica da Universidade Nova de Lisboa

Leishmania cysteine proteinases: from gene to subunit vaccine.

Whole genome sequences of microbial pathogens present new opportunities for clinical application. Presently, genome sequencing of the human protozoan parasite Leishmania major is in progress. The driving forces behind the genome project are to identify genes with key cellular functions and new drug targets, to increase knowledge on mechanisms of drug resistance and to favor technology transfer to scientists from endemic countries. Sequencing of the genome is also aimed at the identification of genes that are expressed in the infectious stages of the parasite and in particular in the intracellular form of the parasite. Several protective antigens of Leishmania have been identified. In addition to these antigens, lysosomal cysteine proteinases (CPs) have been characterized in different strains of Leishmania and Trypanosoma, as new target molecules. Recently, we have isolated and characterized Type I (CPB) and Type II (CPA) cysteine proteinase encoding genes from L. major. The exact function of cysteine proteinases of Leishmania is not completely understood, although there are a few reports describing their role as virulence factors. One specific feature of CPB in Leishmania and other trypanosomatids, is the presence of a Cterminal extension (CTE) which is possibly indicative of conserved structure and function. Recently, we demonstrated that DNA immunization of genetically susceptible BALB / c mice, using a cocktail of CPB and CPA genes, induced long lasting protection against L. major infection. This review intends to give an overview of the current knowledge on genetic vaccination used against leishmaniasis and the importance of CP genes for such an approach.

Evolution of C(4) phosphoenolpyruvate carboxykinase in grasses, from genotype to phenotype.

C(4) photosynthesis is an adaptation over the classical C(3) pathway that has evolved multiple times independently. These convergences are accompanied by strong variations among the independent C(4) lineages. The decarboxylating enzyme used to release CO(2) around Rubisco particularly differs between C(4) species, a criterion used to distinguish three distinct biochemical C(4) subtypes. The phosphoenolpyruvate carboxykinase (PCK) serves as a primary decarboxylase in a minority of C(4) species. This enzyme is also present in C(3) plants, where it is responsible for nonphotosynthetic functions. The genetic changes responsible for the evolution of C(4)-specific PCK are still unidentified. Using phylogenetic analyses on PCK sequences isolated from C(3) and C(4) grasses, this study aimed at resolving the evolutionary history of C(4)-specific PCK enzymes. Four independent evolutions of C(4)-PCK were shown to be driven by positive selection, and nine C(4)-adaptive sites underwent parallel genetic changes in different C(4) lineages. C(4)-adaptive residues were also observed in C(4) species from the nicotinamide adenine dinucleotide phosphate-malic enzyme (NADP-ME) subtype and particularly in all taxa where a PCK shuttle was previously suggested to complement the NADP-ME pathway. Acquisitions of C(4)-specific PCKs were mapped on a species tree, which revealed that the PCK subtype probably appeared at the base of the Chloridoideae subfamily and was then recurrently lost and secondarily reacquired at least three times. Linking the genotype to subtype phenotype shed new lights on the evolutionary transitions between the different C(4) subtypes.

Computational drug design strategies applied to the modelling of human immunodeficiency virus-1 reverse transcriptase inhibitors

Reverse transcriptase (RT) is a multifunctional enzyme in the human immunodeficiency virus (HIV)-1 life cycle and represents a primary target for drug discovery efforts against HIV-1 infection. Two classes of RT inhibitors, the nucleoside RT inhibitors (NRTIs) and the nonnucleoside transcriptase inhibitors are prominently used in the highly active antiretroviral therapy in combination with other anti-HIV drugs. However, the rapid emergence of drug-resistant viral strains has limited the successful rate of the anti-HIV agents. Computational methods are a significant part of the drug design process and indispensable to study drug resistance. In this review, recent advances in computer-aided drug design for the rational design of new compounds against HIV-1 RT using methods such as molecular docking, molecular dynamics, free energy calculations, quantitative structure-activity relationships, pharmacophore modelling and absorption, distribution, metabolism, excretion and toxicity prediction are discussed. Successful applications of these methodologies are also highlighted.

Reaction of sweet orange cultivars expressing the attacin A gene to 'Candidatus Liberibacter asiaticus' infection

The objective of this work was to evaluate the reaction of four sweet orange cultivars expressing the attacin A gene to 'Candidatus Liberibacter asiaticus' (Las) infection, a bacterium associated to huanglongbing (HLB) disease. Transgenic sweet orange plants of Hamlin, Natal, Pêra, and Valência cultivars, as well as nontransgenic controls received inocula by grafting budwood sections of HLB-infected branches. Disease progression was evaluated through observations of leaf symptoms and by polymerase chain reaction (PCR) analysis, eight months after inoculation. A completely randomized design was used, with four experiments (one for each cultivar) performed simultaneously. Bacteria title was estimated by quantitative PCR (qPCR). HLB symptoms and Las titers were present in nontransgenic and transgenic plants expressing the attacin A gene of the four sweet orange cultivars, eight months after bacteria inoculation. Five transgenic lines (transformation events) of 'Pêra' sweet orange expressing the attacin A gene have significantly lower Las titers in comparison with nontransgenic plants of this cultivar.