trans gene regulation in adaptive evolution: A genetic algorithm model

This is a continuation of earlier studies on the evolution of infinite populations of haploid genotypes within a genetic algorithm framework. We had previously explored the evolutionary consequences of the existence of indeterminate—“plastic”—loci, where a plastic locus had a finite probability in each generation of functioning (being switched “on”) or not functioning (being switched “off”). The relative probabilities of the two outcomes were assigned on a stochastic basis. The present paper examines what happens when the transition probabilities are biased by the presence of regulatory genes. We find that under certain conditions regulatory genes can improve the adaptation of the population and speed up the rate of evolution (on occasion at the cost of lowering the degree of adaptation). Also, the existence of regulatory loci potentiates selection in favour of plasticity. There is a synergistic effect of regulatory genes on plastic alleles: the frequency of such alleles increases when regulatory loci are present. Thus, phenotypic selection alone can be a potentiating factor in a favour of better adaptation.

Positional Preferences of Polypurine/Polypyrimidine Tracts in Saccharomyces cerevisiae Genome: Implications for cis Regulation of Gene Expression

The complete genome of the baker's yeast S. cerevisiae was analyzed for the presence of polypurine/polypyrimidine (poly[pu/py]) repeats and their occurrences were classified on the basis of their location within and outside open reading frames (ORFs). The analysis reveals that such sequence motifs are present abundantly both in coding as well as noncoding regions. Clear positional preferences are seen when these tracts occur in noncoding regions. These motifs appear to occur predominantly at a unit nucleosomal length both upstream and downstream of ORFs. Moreover, there is a biased distribution of polypurines in the coding strands when these motifs occur within open reading frames. The significance of the biased distribution is discussed with reference to the occurrence of these motifs in other known mRNA sequences and expressed sequence tags. A model for cis regulation of gene expression is proposed based on the ability of these motifs to form an intermolecular triple helix structure when present within the coding region and/or to modulate nucleosome positioning via enhanced histone affinity when present outside coding regions.

Identification and functional characterization of a cis-acting positive DNA element regulating CYP 2B1/B2 gene transcription in rat liver

A positive cis-acting DNA element in the near 5'-upstream region of the CYP2B1/B2 genes in rat liver was found to play an important role in the transcription of these genes. An oligonucleotide covering -69 to -98 nt mimicked the gel mobility shift pattern given by the fragment -179 to +29 nt, which was earlier found adequate to confer the regulatory features of this gene. Two major complexes were seen, of which the slower and faster moving complexes became intense under uninduced and Phenobarbitone-induced conditions respectively. Minigene cloned DNA plasmid covering -179 to +181 nt in pUC 19 and Bal 31 mutants derived from this parent were transcribed in whole nuclei and cell free transcription extracts and mutants containing only upto -75 nt of the upstream were poorly transcribed. Transcription extracts from phenobarbitone-injected rat liver nuclei were significantly more active than extracts from uninduced rats in transcribing the minigene constructs. Addition of the oligonucleotide (-69 to -98nt) specifically inhibited the transcription of the minigene construct (-179 to +181 nt) in the cell free transcription system. It is therefore, concluded that the region -69 to -98 nt acts as a positive cis-acting element in the transcription of the CYP2B1/B2 genes and in mediating the inductive effects of phenobarbitone.

Use of protein A gene fusions for the analysis of structure-function relationship of the transactivator protein C of bacteriophage Mu

A sensitive dimerization assay for DNA binding proteins has been developed using gene fusion technology. For this purpose, we have engineered a gene fusion using protein A gene of Staphylococcus aureus and C gene, the late gene transactivator of bacteriophage Mu. The C gene was fused to the 3' end of the gene for protein A to generate an A- C fusion. The overexpressed fusion protein was purified in a single step using immunoglobulin affinity chromatography. Purified fusion protein exhibits DNA binding activity as demonstrated by electrophoretic mobility shift assays. When the fusion protein A-C was mixed with C and analyzed for DNA binding, in addition to C and A-C specific complexes, a single intermediate complex comprising of a heterodimer of C and A-C fusion proteins was observed. Further, the protein A moiety in the fusion protein A-C does not contribute to DNA binding as demonstrated by proteolytic cleavage and circular dichroism (CD) analysis. The assay has also been applied to analyze the DNA binding domain of C protein by generating fusions between protein A and N- and C-terminal deletion mutants of C. The results indicate a role for the region towards the carboxy terminal of the protein in DNA binding. The general applicability of this method is discussed.

Involvement of Synthesis and Phosphorylation of Nuclear Protein Factors That Bind to the Positivecis-Acting Element in the Transcriptional Activation of the CYP2B1/B2 Gene by Phenobarbitonein Vivo

The synthesis and phosphorylation of protein factor(s) that bind to the positivecis-acting element (−69 to −98 nt) of the CYP2B1/B2 gene have been examinedin vivoin the rat. Treatment of rats with cycloheximide, a protein synthetic inhibitor, suppresses basal as well as phenobarbitone-induced levels of CYP2B1/B2 mRNA and its run-on transcription. Under these conditions, complex formation of the nuclear extract with the positive element is also inhibited, as judged by gel shift assays. Treatment of rats with 2-aminopurine, a general protein kinase inhibitor, blocks the phenobarbitone-mediated increase in CYP2B1/B2 mRNA, cell-free transcription of a minigene construct containing the positive element, pP450e179DNA, and binding of nuclear proteins to the positive element. Treatment of rats with okadaic acid, a protein phosphatase inhibitor, mimics the effects of phenobarbitone, but only partially. Thus, both phenobarbitone and okadaic acid individually enhance binding of the nuclear protein(s) to the positive element, cell-free transcription of the minigene construct, and phosphorylation of the not, vert, similar26- and 94-kDa proteins binding to the positive element. But unlike phenobarbitone, okadaic acid is not an inducer of CYP2B1/B2 mRNA or its run-on transcription. Thus, phenobarbitone-responsive positive element interactions constitute only a minimal requirement, and okadaic acid is perhaps not able to bring about the total requirement for activation of CYP2B1/B2 gene transcription that should include interaction between the minimal promoter and further upstream elements. An intriguing feature is the antagonistic effect of okadaic acid on phenobarbitone-mediated effects on CYP2B1/B2 mRNA levels, cell-free and run-on transcription, and nuclear protein binding to the positive element. The reason for this antagonism is not clear. It is concluded that phenobarbitone treatment enhancesin vivothe synthesis and phosphorylation of protein factors binding to the positive element and these constitute a minimal requirement for the transcriptional activation of the CYP2B1/B2 gene.

Genetic studies of the PRP17 gene of Saccharomyces cerevisiae: a domain essential for function maps to a nonconserved region of the protein

The PRP17 gene product is required for the second step of pre-mRNA splicing reactions. The C-terminal half of this protein bears four repeat units with homology to the beta transducin repeat. Missense mutations in three temperature-sensitive prp17 mutants map to a region in the N-terminal half of the protein. We have generated, in vitro, 11 missense alleles at the beta transducin repeat units and find that only one affects function in vivo. A phenotypically silent missense allele at the fourth repeat unit enhances the slow-growing phenotype conferred by an allele at the third repeat, suggesting an interaction between these domains. Although many missense mutations in highly conserved amino acids lack phenotypic effects, deletion analysis suggests an essential role for these units. Only mutations in the N-terminal nonconserved domain of PRP17 are synthetically lethal in combination with mutations in PRP16 and PRP18, two other gene products required for the second splicing reaction. A mutually allele-specific interaction between Prp17 and snr7, with mutations in U5 snRNA, was observed. We therefore suggest that the functional region of Prp17p that interacts with Prp18p, Prp16p, and U5 snRNA is the N terminal region of the protein.

Gene Expression : From Microarrays to Functional Genomics

The time of the large sequencing projects has enabled unprecedented possibilities of investigating more complex aspects of living organisms. Among the high-throughput technologies based on the genomic sequences, the DNA microarrays are widely used for many purposes, including the measurement of the relative quantity of the messenger RNAs. However, the reliability of microarrays has been strongly doubted as robust analysis of the complex microarray output data has been developed only after the technology had already been spread in the community. An objective of this study consisted of increasing the performance of microarrays, and was measured by the successful validation of the results by independent techniques. To this end, emphasis has been given to the possibility of selecting candidate genes with remarkable biological significance within specific experimental design. Along with literature evidence, the re-annotation of the probes and model-based normalization algorithms were found to be beneficial when analyzing Affymetrix GeneChip data. Typically, the analysis of microarrays aims at selecting genes whose expression is significantly different in different conditions followed by grouping them in functional categories, enabling a biological interpretation of the results. Another approach investigates the global differences in the expression of functionally related groups of genes. Here, this technique has been effective in discovering patterns related to temporal changes during infection of human cells. Another aspect explored in this thesis is related to the possibility of combining independent gene expression data for creating a catalog of genes that are selectively expressed in healthy human tissues. Not all the genes present in human cells are active; some involved in basic activities (named housekeeping genes) are expressed ubiquitously. Other genes (named tissue-selective genes) provide more specific functions and they are expressed preferably in certain cell types or tissues. Defining the tissue-selective genes is also important as these genes can cause disease with phenotype in the tissues where they are expressed. The hypothesis that gene expression could be used as a measure of the relatedness of the tissues has been also proved. Microarray experiments provide long lists of candidate genes that are often difficult to interpret and prioritize. Extending the power of microarray results is possible by inferring the relationships of genes under certain conditions. Gene transcription is constantly regulated by the coordinated binding of proteins, named transcription factors, to specific portions of the its promoter sequence. In this study, the analysis of promoters from groups of candidate genes has been utilized for predicting gene networks and highlighting modules of transcription factors playing a central role in the regulation of their transcription. Specific modules have been found regulating the expression of genes selectively expressed in the hippocampus, an area of the brain having a central role in the Major Depression Disorder. Similarly, gene networks derived from microarray results have elucidated aspects of the development of the mesencephalon, another region of the brain involved in Parkinson Disease.

Type III Secretion System of Phytopathogenic Bacterium Pseudomonas syringae:From Gene to Function

The type III secretion system (T3SS) is an essential requirement for the virulence of many Gram-negative bacteria which infect plants, animals and men. Pathogens use the T3SS to deliver effector proteins from the bacterial cytoplasm to the eukaryotic host cells, where the effectors subvert host defenses. The best candidates for directing effector protein traffic are the bacterial type III-associated appendages, called needles or pili. In plant pathogenic bacteria, the best characterized example of a T3SS-associated appendage is the HrpA pilus of the plant pathogen Pseudomonas syringae pv. tomato DC3000. The components of the T3SS in plant pathogens are encoded by a cluster of hrp (hypersensitive reaction and pathogenicity) genes. Two major classes of T3SS-secreted proteins are: harpin proteins such as HrpZ which are exported into extracellular space, and avirulence (Avr) proteins such as AvrPto which are translocated directly to the plant cytoplasm. This study deals with the structural and functional characterization of the T3SS-associated HrpA pilus and the T3SS-secreted harpins. By insertional mutagenesis analysis of HrpA, we located the optimal epitope insertion site in the amino-terminus of HrpA, and revealed the potential application of the HrpA pilus as a carrier of antigenic determinants for vaccination. By pulse-expression of proteins combined with immuno-electron microscopy, we discovered the Hrp pilus assembly strategy as addition of HrpA subunits to the distal end of the growing pilus, and we showed for the first time that secretion of HrpZ occurs at the tip of the pilus. The pilus thus functions as a conduit delivering proteins to the extracellular milieu. By using phage-display and scanning-insertion mutagenesis methods we identified a conserved HrpZ-binding peptide and localized the peptide-binding site to the central domain of HrpZ. We also found that the HrpZ specifically interacts with a host bean protein. Taken together, the current results provide deeper insight into the molecular mechanism of T3SS-associated pilus assembly and effector protein translocation, which will be helpful for further studies on the pathogenic mechanisms of Gram-negative bacteria and for developing new strategies to prevent bacterial infection.

Phylogenetics of Myrmica ants and their social parasites

Understanding the overwhelming diversity of life calls for complex organisational schemes. The field of systematics may thus be seen as the cornerstone of evolutionary biology. In the last few decades, systematics has been rejuvenated through the introduction of molecular methods such as DNA barcoding and multi-gene phylogenetic approaches. These methods may shed new light on established taxonomic ideas and problems. For example, the classification of ants has aroused much debate due to reinterpretation of morphological characters or contradictions between molecular data and morphology. Only in the last few years a consensus was reached regarding the phylogeny of ant subfamilies. However, the situation remains deplorable for lower taxonomic ranks such as subfamilies, tribes and genera. This thesis describes the systematics and evolution of the Holarctic ant genus Myrmica and the tribe to which it belongs, Myrmicini. Using barcoding, molecular-phylogenetic data and divergence time estimations, it addresses questions regarding the taxonomy, morphology and biogeography of this group. Furthermore, the interrelationships between socially parasitic Myrmica species and their hosts (other species in the genus) were inferred. The phylogeny suggests that social parasitism evolved several times in Myrmica. Finally, this thesis investigated whether coevolution shaped the phylogeny of socially parasitic Maculinea butterflies that live inside Myrmica colonies. No evidence was found for coevolution.

FGFR1 regulated gene-expression, cell proliferation and differentiation in the developing midbrain and hindbrain

The neuroectodermal tissue close to the midbrain hindbrain boundary (MHB) is an important secondary organizer in the developing neural tube. This so-called isthmic organizer (IsO) regulates cellular survival, patterning and proliferation in the midbrain (Mb) and rhombomere 1 (R1) of the hindbrain. Signaling molecules of the IsO, such as fibroblast growth factor 8 (FGF8) and WNT1 are expressed in distinct bands of cells around the MHB. It has been previously shown that FGF-receptor 1 (FGFR1) is required for the normal development of this brain region in the mouse embryo. In the present study, we have compared the gene expression profiles of wild-type and Fgfr1 mutant embryos. We show that the loss of Fgfr1 results in the downregulation of several genes expressed close to the MHB and in the disappearance of gene expression gradients in the midbrain and R1. Our microarray screen identified several previously uncharacterized genes which may participate in the development of midbrain R1 region. Our results also show altered neurogenesis in the midbrain and R1 of the Fgfr1 mutants. Interestingly, the neuronal progenitors in midbrain and R1 show different responses to the loss of signaling through FGFR1. As Wnt1 expression at the MHB region requires the FGF signaling pathway, WNT target genes, including Drapc1, were also identified in our screen. The microarray data analysis also suggested that the cells next to the midbrain hindbrain boundary express distinct cell cycle regulators. We showed that the cells close to the border appeared to have unique features. These cells proliferate less rapidly than the surrounding cells. Unlike the cells further away from the boundary, these cells express Fgfr1 but not the other FGF receptors. The slowly proliferating boundary cells are necessary for development of the characteristic isthmic constriction. They may also contribute to compartmentalization of this brain region.

Dung beetle radiations in Madagascar

Speciation on islands is affected by island size and the range of habitats and resources available and often also by limited interactions with other taxa. An ancestral population may evolve into a large number of species via an adaptive radiation. In Madagascar, most groups of animals and plants have radiated on the island, having arrived via oceanic dispersal during the long isolation of Madagascar. Characteristic features of Malagasy biota are exceptionally high level of endemism, high species richness as well as lack of many higher taxa that are dominant on the African mainland. Malagasy dung beetles are dominated by two tribes, Canthonini and Helictopleurina, with more than 250 endemic species. In this thesis I have reconstructed molecular phylogenies for the two tribes using several gene regions and different phylogenetic methods. Evolution of closely related species and among populations of the same species was examined with haplotype networks. The Malagasy Canthonini consists of three large lineages, while Helictopleurina forms a monophyletic group. The ancestors of each of the four clades colonised Madagascar at different times during Cenozoic. The subsequent radiations differ in terms of the number of extant species (from 37 to more than 100) and the level of ecological differentiation. In addition, Onthophagini (6 species) and Scarabaeini (3) have colonised Madagascar several times, but they have not radiated and the few species have not entered forests where Canthonini and Helictopleurina mostly occur. Among the three Canthonini radiations, speciation appears to have been mostly allopatric in the oldest and the youngest clades, while in the Epactoides clade sister species have diverged in their ecologies but have similar geographical distributions, indicating that speciation may have occurred in regional sympatry. The most likely isolating mechanisms have been rivers and forest refugia during dry and cool geological periods. Most species are generalists feeding on both carrion and dung, and competition among ecologically similar species may prevent their coexistence in the same communities. Some species have evolved to forage in the canopy and a few species have shifted to use cattle dung, a new resource in the open habitats following the introduction of cattle 1500 years ago. The latter shift has allowed species to expand their geographical ranges.

Ecological consequences of genetic modifications - an invasion analysis approach

Biological invasions are considered as one of the greatest threats to biodiversity, as they may lead to disruption and homogenization of natural communities, and in the worst case, to native species extinctions. The introduction of gene modified organisms (GMOs) to agricultural, fisheries and forestry practices brings them into contact with natural populations. GMOs may appear as new invasive species if they are able to (1) invade into natural habitats or (2) hybridize with their wild relatives. The benefits of GMOs, such as increased yield or decreased use of insecticides or herbicides in cultivation, may thus be reduced due the potential risks they may cause. A careful ecological risk analysis therefore has to precede any responsible GMO introduction. In this thesis I study ecological invasion in relation to GMOs, and what kind of consequences invasion may have in natural populations. A set of theoretical models that combine life-history evolution, population dynamics, and population genetics were developed for the hazard identification part of ecological risks assessment of GMOs. In addition, the potential benefits of GMOs in management of an invasive pest were analyzed. In the first study I showed that a population that is fluctuating due to scramble-type density dependence (due to, e.g., nutrient competition in plants) may be invaded by a population that is relatively more limited by a resource (e.g., light in plants) that is a cause of contest-type density dependence. This result emphasises the higher risk of invasion in unstable environments. The next two studies focused on escape of a growth hormone (GH) transgenic fish into a natural population. The results showed that previous models may have given too pessimistic a view of the so called Trojan gene -effect, where the invading genotype is harmful for the population as a whole. The previously suggested population extinctions did not occur in my studies, since the changes in mating preferences caused by the GH-fish were be ameliorated by decreased level of competition. The GH-invaders may also have to exceed a threshold density before invasion can be successful. I also showed that the prevalence of mature parr (aka. sneaker) strategy among GH-fish may have clear effect on invasion outcome. The fourth study assessed the risks and developed methods against the invasion of the Colorado Potato Beetle (CPB, Leptinotarsa decemlineata). I showed that the eradication of CPB is most important for the prevention of their establishment, but the cultivation of transgenic Bt-potato could also be effective. In general, my results emphasise that invasion of transgenic species or genotypes to be possible under certain realistic conditions and resulting in competitive exclusion, population decline through outbreeding depression and genotypic displacement of native species. Ecological risk assessment should regard the decline and displacement of the wild genotype by an introduced one as a consequence that is as serious as the population extinction. It will also be crucial to take into account different kinds of behavioural differences among species when assessing the possible hazards that GMOs may cause if escaped. The benefits found of GMO crops effectiveness in pest management may also be too optimistic since CPB may evolve resistance to Bt-toxin. The models in this thesis could be further applied in case specific risk assessment of GMOs by supplementing them with detailed data of the species biology, the effect of the transgene introduced to the species, and also the characteristics of the populations or the environments in the risk of being invaded.

ERAP2 functional knockout in humans does not alter surface heavy chains or HLA-B27, inflammatory cytokines or endoplasmic reticulum stress markers

Introduction Single nucleotide polymorphisms in ERAP2 are strongly associated with ankylosing spondylitis (AS). One AS-associated single nucleotide polymorphism, rs2248374, causes a truncated ERAP2 protein that is degraded by nonsense-mediated decay. Approximately 25% of the populations of European ancestry are therefore natural ERAP2 knockouts. We investigated the effect of this associated variant on HLA class I allele presentation, surface heavy chains, endoplasmic reticulum (ER) stress markers and cytokine gene transcription in AS. Methods Patients with AS and healthy controls with either AA or GG homozygous status for rs2248374 were studied. Antibodies to CD14, CD19-ECD, HLA-A-B-C, Valpha7.2, CD161, anti-HC10 and anti-HLA-B27 were used to analyse peripheral blood mononuclear cells. Expression levels of ER stress markers (GRP78 and CHOP) and proinflammatory genes (tumour necrosis factor (TNF), IL6, IL17 and IL22) were assessed by qPCR. Results There was no significant difference in HLAclass I allele presentation or major histocompatibility class I heavy chains or ER stress markers GRP78 and CHOP or proinflammatory gene expression between genotypes for rs2248374 either between cases, between cases and controls, and between controls. Discussion Large differences were not seen in HLAB27 expression or cytokine levels between subjects with and without ERAP2 in AS cases and controls. This suggests that ERAP2 is more likely to influence AS risk through other mechanisms.

Genetic aspects of outcome in kidney transplantation: cytokine and thrombosis associated candidate genes and gene expression biomarkers

Kidney transplantation (Tx) is the treatment of choice for end stage renal disease. Immunosuppressive medications are given to prevent an immunological rejection of the transplant. However, immunosuppressive drugs increase e.g. the risk of infection, cancer or nephrotoxicity. A major genetic contributors to immunological acceptance of the graft are human leukocyte antigen (HLA) genes. Also other non-HLA gene polymorphisms may predict the future risk of complications before Tx, possibly enabling individualised immunotherapy. Graft function after Tx is monitored using non-specific clinical symptoms and laboratory markers. The definitive diagnosis of graft rejection however relies on a biopsy of the graft. In the acute rejection (AR) diagnostics there is a need for an alternative to biopsy that would be an easily repeatable and simple method for regular use. Frequent surveillance of acute or subclinical rejection (SCR) may improve long-term function. In this thesis, associations between cytokine and thrombosis associated candidate genes and the outcome of kidney Tx were studied. Cytotoxic and co-stimulatory T lymphocyte molecule gene expression biomarkers for the diagnosis of the AR and the SCR were also investigated. We found that polymorphisms in the cytokine genes tumor necrosis factor and interleukin 10 (IL10) of the recipients were associated with AR. In addition, certain IL10 gene polymorphisms of the donors were associated with the incidence of cytomegalovirus infection and occurrence of later infection in a subpopulation of recipients. Further, polymorphisms in genes related to the risk of thrombosis and those of certain cytokines were not associated with the occurrence of thrombosis, infarction, AR or graft survival. In the study of biomarkers for AR, whole blood samples were prospectively collected from adult kidney Tx patients. With real-time quantitative PCR (RT-QPCR) gene expression quantities of CD154 and ICOS differentiated the patients with AR from those without, but not from the patients with other causes of graft dysfunction. Biomarkers for SCR were studied in paediatric kidney Tx patients. We used RT-QPCR to quantify the gene expression of immunological candidate genes in a low-density array format. In addition, we used RT-QPCR to validate the results of the microarray analysis. No gene marker differentiated patients with SCR from those without SCR. This research demonstrates the lack of robust markers among polymorphisms or biomarkers in investigated genes that could be included in routine analysis in a clinical laboratory. In genetic studies, kidney Tx can be regarded as a complex trait, i.e. several environmental and genetic factors may determine its outcome. A number of currently unknown genetic factors probably influence the results of Tx.