The genetic basis of morphological change in convergent evolution of natural populations: identifying candidate genes behind convergent evolution in blind cavefish astyanax mexicanus

Dissertation presented to obtain the Ph.D. degree in Biology at the Instituto de Tecnologia Química e Biológica, Universidade Nova de Lisboa.

Mini-review: Specificity and expression of CIITA, the master regulator of MHC class II genes.

The class II transactivator (CIITA) has been referred to as the "master control factor" for the expression of MHC class II (MHCII) genes. As our knowledge on the specificity and function of CIITA grows, it is becoming increasingly evident that this sobriquet is entirely justified. First, despite extensive investigations, the major target genes of CIITA remain those implicated in the presentation of antigenic peptides by MHCII molecules. Although other putative target genes have been reported, the contribution of CIITA to their expression remains indirect, controversial or comparatively minor relative to its decisive role as a regulator of MHCII and related genes. Second, the most important parameter dictating MHCII expression is by far the expression pattern of the gene encoding CIITA (MHC2TA). The vast majority of signals that activate or repress MHCII expression under physiological and pathological situations converge on one or more of the three alternative promoters that drive transcription of the MHC2TA gene. In short, with respect to its specificity and its exquisitely controlled pattern of expression, CIITA is by a long stretch the single most important transcription factor for the regulation of genes required for MHCII-restricted antigen-presentation.

Attenuated poxviruses expressing a synthetic HIV protein stimulate HLA-A2-restricted cytotoxic T-cell responses.

Efficient HIV vaccines have to trigger cell-mediated immunity directed against various viral antigens. However little is known about the breadth of the response induced by vaccines carrying multiple proteins. Here, we report on the immunogenicity of a construct harbouring a fusion of the HIV-1 IIIB gag, pol and nef genes (gpn) designed for optimal safety and equimolar expression of the HIV proteins. The attenuated poxviruses, MVA and NYVAC, harbouring the gpn construct, induced potent immune responses in conventional mice characterised by stimulation of Gpn-specific IFN-gamma-producing cells and cytotoxic T cells. In HLA-A2 transgenic mice, recombinant MVA elicited cytotoxic responses against epitopes recognised in most HLA-A2+ HIV-1-infected individuals. We also found that the MVA vaccine triggered the in vitro expansion of peripheral blood cells isolated from a HIV-1-seropositive patient and with similar specificity as found in immunised HLA-A2 transgenic mice. In conclusion, the synthetic HIV polyantigen Gpn delivered by MVA is immunogenic, efficiently processed and presented by human MHC class I molecules.

Chromosomal divergence and evolutionary inferences in Rhodniini based on the chromosomal location of ribosomal genes

In this study, we used fluorescence in situ hybridisation to determine the chromosomal location of 45S rDNA clusters in 10 species of the tribe Rhodniini (Hemiptera: Reduviidae: Triatominae). The results showed striking inter and intraspecific variability, with the location of the rDNA clusters restricted to sex chromosomes with two patterns: either on one (X chromosome) or both sex chromosomes (X and Y chromosomes). This variation occurs within a genus that has an unchanging diploid chromosome number (2n = 22, including 20 autosomes and 2 sex chromosomes) and a similar chromosome size and genomic DNA content, reflecting a genome dynamic not revealed by these chromosome traits. The rDNA variation in closely related species and the intraspecific polymorphism in Rhodnius ecuadoriensis suggested that the chromosomal position of rDNA clusters might be a useful marker to identify recently diverged species or populations. We discuss the ancestral position of ribosomal genes in the tribe Rhodniini and the possible mechanisms involved in the variation of the rDNA clusters, including the loss of rDNA loci on the Y chromosome, transposition and ectopic pairing. The last two processes involve chromosomal exchanges between both sex chromosomes, in contrast to the widely accepted idea that the achiasmatic sex chromosomes of Heteroptera do not interchange sequences.

Identification of new MHC class I-restricted human tumor antigens for immunotherapy

Summary Cancer is a leading cause of morbidity and mortality in Western countries (as an example, colorectal cancer accounts for about 300'000 new cases and 200'000 deaths each year in Europe and in the USA). Despite that many patients with cancer have complete macroscopic clearance of their disease after resection, radiotherapy and/or chemotherapy, many of these patients develop fatal recurrence. Vaccination with immunogenic peptide tumor antigens has shown encouraging progresses in the last decade; immunotherapy might therefore constitute a fourth therapeutic option in the future. We dissect here and critically evaluate the numerous steps of reverse immunology, a forecast procedure to identify antigenic peptides from the sequence of a gene of interest. Bioinformatic algorithms were applied to mine sequence databases for tumor-specific transcripts. A quality assessment of publicly available sequence databanks allowed defining strengths and weaknesses of bioinformatics-based prediction of colon cancer-specific alternative splicing: new splice variants could be identified, however cancer-restricted expression could not be significantly predicted. Other sources of target transcripts were quantitatively investigated by polymerase chain reactions, as cancer-testis genes or reported overexpressed transcripts. Based on the relative expression of a defined set of housekeeping genes in colon cancer tissues, we characterized a precise procedure for accurate normalization and determined a threshold for the definition of significant overexpression of genes in cancers versus normal tissues. Further steps of reverse immunology were applied on a splice variant of the Melan¬A gene. Since it is known that the C-termini of antigenic peptides are directly produced by the proteasome, longer precursor and overlapping peptides encoded by the target sequence were synthesized chemically and digested in vitro with purified proteasome. The resulting fragments were identified by mass spectroscopy to detect cleavage sites. Using this information and based on the available anchor motifs for defined HLA class I molecules, putative antigenic peptides could be predicted. Their relative affinity for HLA molecules was confirmed experimentally with functional competitive binding assays and they were used to search patients' peripheral blood lymphocytes for the presence of specific cytolytic T lymphocytes (CTL). CTL clones specific for a splice variant of Melan-A could be isolated; although they recognized peptide-pulsed cells, they failed to lyse melanoma cells in functional assays of antigen recognition. In the conclusion, we discuss advantages and bottlenecks of reverse immunology and compare the technical aspects of this approach with the more classical procedure of direct immunology, a technique introduced by Boon and colleagues more than 10 years ago to successfully clone tumor antigens.

Evolutionary simulations to detect functional lineage-specific genes.

MOTIVATION: Supporting the functionality of recent duplicate gene copies is usually difficult, owing to high sequence similarity between duplicate counterparts and shallow phylogenies, which hamper both the statistical and experimental inference. RESULTS: We developed an integrated evolutionary approach to identify functional duplicate gene copies and other lineage-specific genes. By repeatedly simulating neutral evolution, our method estimates the probability that an ORF was selectively conserved and is therefore likely to represent a bona fide coding region. In parallel, our method tests whether the accumulation of non-synonymous substitutions reveals signatures of selective constraint. We show that our approach has high power to identify functional lineage-specific genes using simulated and real data. For example, a coding region of average length (approximately 1400 bp), restricted to hominoids, can be predicted to be functional in approximately 94-100% of cases. Notably, the method may support functionality for instances where classical selection tests based on the ratio of non-synonymous to synonymous substitutions fail to reveal signatures of selection. Our method is available as an automated tool, ReEVOLVER, which will also be useful to systematically detect functional lineage-specific genes of closely related species on a large scale. AVAILABILITY: ReEVOLVER is available at http://www.unil.ch/cig/page7858.html.

Development of a T Cell Receptor Targeting an HLA-A*0201 Restricted Epitope from the Cancer-Testis Antigen SSX2 for Adoptive Immunotherapy of Cancer.

The clinical success of adoptive immunotherapy of cancer relies on the selection of target antigens that are highly expressed in tumor cells but absent in essential normal tissues. A group of genes that encode the cancer/testis or cancer germline antigens have been proposed as ideal targets for immunotherapy due to their high expression in multiple cancer types and their restricted expression in immunoprivileged normal tissues. In the present work we report the isolation and characterization of human T cell receptors (TCRs) with specificity for synovial sarcoma X breakpoint 2 (SSX2), a cancer/testis antigen expressed in melanoma, prostate cancer, lymphoma, multiple myeloma and pancreatic cancer, among other tumors. We isolated seven HLA-A2 restricted T cell receptors from natural T cell clones derived from tumor-infiltrated lymph nodes of two SSX2-seropositive melanoma patients, and selected four TCRs for cloning into retroviral vectors. Peripheral blood lymphocytes (PBL) transduced with three of four SSX2 TCRs showed SSX241-49 (KASEKIFYV) peptide specific reactivity, tumor cell recognition and tetramer binding. One of these, TCR-5, exhibited tetramer binding in both CD4 and CD8 cells and was selected for further studies. Antigen-specific and HLA-A*0201-restricted interferon-γ release, cell lysis and lymphocyte proliferation was observed following culture of TCR engineered human PBL with relevant tumor cell lines. Codon optimization was found to increase TCR-5 expression in transduced T cells, and this construct has been selected for development of clinical grade viral vector producing cells. The tumor-specific pattern of expression of SSX2, along with the potent and selective activity of TCR-5, makes this TCR an attractive candidate for potential TCR gene therapy to treat multiple cancer histologies.

Evolution of mRNA expression levels of autosomal and sex-linked genes in amniotes

In addition to differences in protein-coding gene sequences, changes in expression resulting from mutations in regulatory sequences have long been hypothesized to be responsible for phenotypic differences between species. However, unlike comparison of genome sequences, few studies, generally restricted to pairwise comparisons of closely related mammalian species, have assessed between-species differences at the transcriptome level. They reported that gene expression evolves at different rates in various organs and in a pattern that is overall consistent with neutral models of evolution. In the first part of my thesis, I investigated the evolution of gene expression in therian mammals (i.e.7 placental and marsupials), based on microarray data from human, mouse and the gray short-tailed opossum (Monodelphis domestica). In addition to autosomal genes, a special focus was given to the evolution of X-linked genes. The therian X chromosome was recently shown to be younger than previously thought and to harbor a specific gene content (e.g., genes involved in brain or reproductive functions) that is thought to have been shaped by specific sex-related evolutionary forces. Sex chromosomes derive from ordinary autosomes and their differentiation led to the degeneration of the Y chromosome (in mammals) or W chromosome (in birds). Consequently, X- or Z-linked genes differ in gene dose between males and females such that the heterogametic sex has half the X/Z gene dose compared to the ancestral state. To cope with this dosage imbalance, mammals have been reported to have evolved mechanisms of dosage compensation.¦In the first project, I could first show that transcriptomes evolve at different rates in different organs. Out of the five tissues I investigated, the testis is the most rapidly evolving organ at the gene expression level while the brain has the most conserved transcriptome. Second, my analyses revealed that mammalian gene expression evolution is compatible with a neutral model, where the rates of change in gene expression levels is linked to the efficiency of purifying selection in a given lineage, which, in turn, is determined by the long-term effective population size in that lineage. Thus, the rate of DNA sequence evolution, which could be expected to determine the rate of regulatory sequence change, does not seem to be a major determinant of the rate of gene expression evolution. Thus, most gene expression changes seem to be (slightly) deleterious. Finally, X-linked genes seem to have experienced elevated rates of gene expression change during the early stage of X evolution. To further investigate the evolution of mammalian gene expression, we generated an extensive RNA-Seq gene expression dataset for nine mammalian species and a bird. The analyses of this dataset confirmed the patterns previously observed with microarrays and helped to significantly deepen our view on gene expression evolution.¦In a specific project based on these data, I sought to assess in detail patterns of evolution of dosage compensation in amniotes. My analyses revealed the absence of male to female dosage compensation in monotremes and its presence in marsupials and, in addition, confirmed patterns previously described for placental mammals and birds. I then assessed the global level of expression of X/Z chromosomes and contrasted this with its ancestral gene expression levels estimated from orthologous autosomal genes in species with non-homologous sex chromosomes. This analysis revealed a lack of up-regulation for placental mammals, the level of expression of X-linked genes being proportional to gene dose. Interestingly, the ancestral gene expression level was at least partially restored in marsupials as well as in the heterogametic sex of monotremes and birds. Finally, I investigated alternative mechanisms of dosage compensation and found that gene duplication did not seem to be a widespread mechanism to restore the ancestral gene dose. However, I could show that placental mammals have preferentially down-regulated autosomal genes interacting with X-linked genes which underwent gene expression decrease, and thus identified a novel alternative mechanism of dosage compensation.

Restricted transgene expression in the brain with cell-type specific neuronal promoters.

Tissue-targeted expression is of major interest for studying the contribution of cellular subpopulations to neurodegenerative diseases. However, in vivo methods to investigate this issue are limited. Here, we report an analysis of the cell specificity of expression of fluorescent reporter genes driven by six neuronal promoters, with the ubiquitous phosphoglycerate kinase 1 (PGK) promoter used as a reference. Quantitative analysis of AcGFPnuc expression in the striatum and hippocampus of rodents showed that all lentiviral vectors (LV) exhibited a neuronal tropism; however, there was substantial diversity of transcriptional activity and cell-type specificity of expression. The promoters with the highest activity were those of the 67 kDa glutamic acid decarboxylase (GAD67), homeobox Dlx5/6, glutamate receptor 1 (GluR1), and preprotachykinin 1 (Tac1) genes. Neuron-specific enolase (NSE) and dopaminergic receptor 1 (Drd1a) promoters showed weak activity, but the integration of an amplification system into the LV overcame this limitation. In the striatum, the expression profiles of Tac1 and Drd1a were not limited to the striatonigral pathway, whereas in the hippocampus, Drd1a and Dlx5/6 showed the expected restricted pattern of expression. Regulation of the Dlx5/6 promoter was observed in a disease condition, whereas Tac1 activity was unaffected. These vectors provide safe tools that are more selective than others available, for the administration of therapeutic molecules in the central nervous system (CNS). Nevertheless, additional characterization of regulatory elements in neuronal promoters is still required.

Molecular biology of baculovirus and its use in biological control in Brazil

Baculoviruses are insect viruses found mainly in Lepidoptera. The family Baculoviridae is taxonomically divided in two genera, Nucleopolyhedrovirus and Granulovirus, which differ by occlusion body morphology. NPVs (Nucleopolyhedroviruses) have polyhedrical inclusion bodies (PIBs) containing multiple viral particles, while GVs (Granuloviruses) appear to be generally single particles occluded in oval shaped occlusion bodies. During the life cycle, two different viral progenies are produced: BV (Budded Virus) and PDV (Polyhedra Derived Virus), which are essential for the infectious process and virus propagation in host cells. Baculoviruses are being used for pest control and they are especially safe due to their specificity and invertebrate-restricted host range. Baculoviruses have been used as vectors for high level protein expression ofheterologous genes from prokaryotic and eukaryotic organisms. Also, recombinant DNA techniques have allowed the production of genetically modified viral insecticides. This study is a review on the taxonomy, structure, replication and molecular biology of baculoviruses, as well as their use as bioinsecticides in Brazil.

A comparative phenotypic study of kallmann syndrome patients carrying monoallelic and biallelic mutations in the prokineticin 2 or prokineticin receptor 2 genes.

Context: Both biallelic and monoallelic mutations in PROK2 or PROKR2 have been found in Kallmann syndrome (KS). Objective: The objective of the study was to compare the phenotypes of KS patients harboring monoallelic and biallelic mutations in these genes. Design and Patients: We studied clinical and endocrine features that reflect the functioning of the pituitary-gonadal axis, and the nonreproductive phenotype, in 55 adult KS patients (42 men and 13 women), of whom 41 had monoallelic mutations and 14 biallelic mutations in PROK2 or PROKR2. Results: Biallelic mutations were associated with more frequent cryptorchidism (70% vs. 34%, P < 0.05) and microphallus (90% vs. 28%, P < 0.001) and lower mean testicular volume (1.2 +/- 0.4 vs. 4.5 +/- 6.0 ml; P < 0.01) in male patients. Likewise, the testosterone level as well as the basal FSH level and peak LH level under GnRH-stimulation were lower in males with biallelic mutations (0.2 +/- 0.1 vs. 0.7 +/- 0.8 ng/ml; P = 0.05, 0.3 +/- 0.1 vs. 1.8 +/- 3.0 IU/liter; P < 0.05, and 0.8 +/- 0.8 vs. 5.2 +/- 5.5 IU/liter; P < 0.05, respectively). Nonreproductive, nonolfactory anomalies were rare in both sexes and were never found in patients with biallelic mutations. The mean body mass index of the patients (23.9 +/- 4.2 kg/m(2) in males and 26.3 +/- 6.6 kg/m(2) in females) did not differ significantly from that of gender-, age-, and treatment-matched KS individuals who did not carry a mutation in PROK2 or PROKR2. Finally, circadian cortisol levels evaluated in five patients, including one with biallelic PROKR2 mutations, were normal in all cases. Conclusion: Male patients carrying biallelic mutations in PROK2 or PROKR2 have a less variable and on average a more severe reproductive phenotype than patients carrying monoallelic mutations in these genes. Nonreproductive, nonolfactory clinical anomalies associated with KS seem to be restricted to patients with monoallelic mutations.

Alps, genes, and chromosomes: their role in the formation of species in the Sorex araneus group (Mammalia, Insectivora), as inferred from two hybrid zones.

During the Pleistocene glaciations, the Alps were an efficient barrier to gene flow between isolated populations, often leading to allopatric speciation. Afterwards, the Alps strongly influenced the post-glacial recolonization of Europe and represent a major suture zone between differentiated populations. Two hybrid zones in the Swiss and French Alps between genetically and chromosomally well-differentiated species-the Valais shrew, Sorex antinorii, and the common shrew, S. araneus-were studied karyotypically and by analyzing the distribution of seven microsatellite loci. In the center of the Haslital hybrid zone the two species coexist over a distance of 900 m. Hybrid karyotypes, among them the most complex known in Sorex, are rare. F-statistics based on microsatellite data revealed a strong heterozygote deficit only in the center of the zone, due to the sympatric distribution of the two species with little hybridization between them. Structuring within the species (both F(IS) and F(ST)) was low. An hierarchical analysis showed a high level of interspecific differentiation. Results were compared with those previously reported in another hybrid zone located at Les Houches in the French Alps. Genetic structuring within and between species was comparable in both hybrid zones, although chromosomal incompatibilities are more important in Haslital, where a linkage block of the race-specific chromosomes should additionally impede gene flow. Evidence for a more restricted gene flow in Haslital comes from the genetically intermediate hybrid karyotypes, whereas in Les Houches, hybrid karyotypes are genetically identical to individuals of the pure karyotypic races. Genic and chromosomal introgression was observed in Les Houches, but not in Haslital. The possible influence of a river, separating the two species at Les Houches, on gene flow is discussed.

Differential regulation of horizontally acquired and core genome genes by the bacterial modulator H-NS

Horizontal acquisition of DNA by bacteria dramatically increases genetic diversity and hence successful bacterial colonization of several niches, including the human host. A relevant issue is how this newly acquired DNA interacts and integrates in the regulatory networks of the bacterial cell. The global modulator H-NS targets both core genome and HGT genes and silences gene expression in response to external stimuli such as osmolarity and temperature. Here we provide evidence that H-NS discriminates and differentially modulates core and HGT DNA. As an example of this, plasmid R27-encoded H-NS protein has evolved to selectively silence HGT genes and does not interfere with core genome regulation. In turn, differential regulation of both gene lineages by resident chromosomal H-NS requires a helper protein: the Hha protein. Tight silencing of HGT DNA is accomplished by H-NS-Hha complexes. In contrast, core genes are modulated by H-NS homoligomers. Remarkably, the presence of Hha-like proteins is restricted to the Enterobacteriaceae. In addition, conjugative plasmids encoding H-NS variants have hitherto been isolated only from members of the family. Thus, the H-NS system in enteric bacteria presents unique evolutionary features. The capacity to selectively discriminate between core and HGT DNA may help to maintain horizontally transmitted DNA in silent form and may give these bacteria a competitive advantage in adapting to new environments, including host colonization.

Differential regulation of horizontally acquired and core genome genes by the bacterial modulator H-NS

Horizontal acquisition of DNA by bacteria dramatically increases genetic diversity and hence successful bacterial colonization of several niches, including the human host. A relevant issue is how this newly acquired DNA interacts and integrates in the regulatory networks of the bacterial cell. The global modulator H-NS targets both core genome and HGT genes and silences gene expression in response to external stimuli such as osmolarity and temperature. Here we provide evidence that H-NS discriminates and differentially modulates core and HGT DNA. As an example of this, plasmid R27-encoded H-NS protein has evolved to selectively silence HGT genes and does not interfere with core genome regulation. In turn, differential regulation of both gene lineages by resident chromosomal H-NS requires a helper protein: the Hha protein. Tight silencing of HGT DNA is accomplished by H-NS-Hha complexes. In contrast, core genes are modulated by H-NS homoligomers. Remarkably, the presence of Hha-like proteins is restricted to the Enterobacteriaceae. In addition, conjugative plasmids encoding H-NS variants have hitherto been isolated only from members of the family. Thus, the H-NS system in enteric bacteria presents unique evolutionary features. The capacity to selectively discriminate between core and HGT DNA may help to maintain horizontally transmitted DNA in silent form and may give these bacteria a competitive advantage in adapting to new environments, including host colonization.