Genome reconstruction and comparative genomics of pneumocystis jirovecii

Pneumocystis jirovecii is a fungus belonging to a basal lineage of the Ascomycotina, the Taphrinomycotina subphylum. It is a parasite specific to humans that dwells primarily in the lung and can cause severe pneumonia in individuals with debilitated immune system. Despite its clinical importance, many aspects of its biology remain poorly understood, at least in part because of the lack of a continuous in vitro cultivation system. The present thesis consists in the genome reconstruction and comparative genomics of P. jirovecii. It is made of three parts: (i) the de novo sequencing of P. jirovecii genome starting from a single broncho- alveolar lavage fluid of a single patient (ii) the de novo sequencing of the genome of the plant pathogen Taphrina deformans, a fungus closely related to P. jirovecii, and (iii) the genome scale comparison of P. jirovecii to other Taphrinomycotina members. Enrichment in P. jirovecii cells by immuno-precipitation, whole DNA random amplification, two complementary high throughput DNA sequencing methods, and in silico sorting and assembly of sequences were used for the de novo reconstruction of P. jirovecii genome from the microbiota of a single clinical specimen. An iterative ad hoc pipeline as well as numerical simulations was used to recover P. jirovecii sequences while purging out contaminants and assembly or amplification chimeras. This strategy produced a 8.1 Mb assembly, which encodes 3,898 genes. Homology searches, mapping on biochemical pathways atlases, and manual validations revealed that this genome lacks (i) most of the enzymes dedicated to the amino acids biosyntheses, and (ii) most virulence factors observed in other fungi, e.g. the glyoxylate shunt pathway and specific peptidases involved in the degradation of the host cell membrane. The same analyses applied to the available genomic sequences from Pneumocystis carinii the species infecting rats and Pneumocystis murina the species infecting mice revealed the same deficiencies. The genome sequencing of T. deformans yielded a 13 Mb assembly, which encodes 5,735 genes. T. deformans possesses enzymes involved plant cell wall degradation, secondary metabolism, the glyoxylate cycle, detoxification, sterol biosynthesis, as well as the biosyntheses of plant hormones such as abscisic acid or indole-3-acetic acid. T. deformans also harbors gene subsets that have counterparts in plant saprophytes or pathogens, which is consistent with its alternate saprophytic and pathogenic lifestyles. Mating genes were also identified. The homothallism of this fungus suggests a mating-type switching mechanism. Comparative analyses indicated that 81% of P. jirovecii genes are shared with eight other Taphrinomycotina members, including T. deformans, P. carinii and P. murina. These genes are mostly involved in housekeeping activities. The genes specific to the Pneumocystis genus represent 8%, and are involved in RNA metabolism and signaling. The signaling is known to be crucial for interaction of Pneumocystis spp with their environment. Eleven percent are unique to P. jirovecii and encode mostly proteins of unknown function. These genes in conjunction with other ones (e.g. the major surface glycoproteins) might govern the interaction of P. jirovecii with its human host cells, and potentially be responsible of the host specificity. P. jirovecii exhibits a reduced genome in size with a low GC content, and most probably scavenges vital compounds such as amino acids and cholesterol from human lungs. Consistently, its genome encodes a large set of transporters (ca. 22% of its genes), which may play a pivotal role in the acquisition of these compounds. All these features are generally observed in obligate parasite of various kingdoms (bacteria, protozoa, fungi). Moreover, epidemiological studies failed to evidence a free-living form of the fungus and Pneumocystis spp were shown to co-evolved with their hosts. Given also the lack of virulence factors, our observations strongly suggest that P. jirovecii is an obligate parasite specialized in the colonization of human lungs, and which causes disease only in individuals with compromised immune system. The same conclusion is most likely true for all other Pneumocystis spp in their respective mammalian host. - Pneumocystis jirovecii est un champignon appartenant à ine branche basale des Ascomycotina, le sous-embranchement des Taphrinomycotina. C'est un parasite spécifique aux humains qui réside principalement dans les poumons, et qui peut causer des pneumonies sévères chez des individus ayant un système immunitaire déficient. En dépit de son importance clinique, de nombreux aspects de sa biologie demeurent,largement méconnus, au moins en partie à cause de l'absence d'un système de culture in vitro continu. Cette thèse traite de la reconstruction du génome et de la génomique comparative de P. jirovecii. Elle comporte trois parties: (i) le séquençage de novo du génome de P. jirovecii à partir d'un lavage broncho-alvéolaire provenant d'un seul patient, (ii) le séquençage de novo du génome d'un champignon pathogène de plante Taphrina deformans qui est phylogénétiquement proche de P. jirovecii, et (iii) la comparaison du génome de P. jirovecii à celui d'autres membres du sous-embranchement des Taphrinomycotina. Un enrichissement en cellules de P. jirovecii par immuno-précipitation, une amplification aléatoire des molécules d'ADN, deux méthodes complémentaires de séquençage à haut débit, un tri in silico et un assemblage des séquences ont été utilisés pour reconstruire de novo le génome de P. jirovecii à partir du microbiote d'un seul échantillon clinique. Un pipeline spécifique ainsi que des simulations numériques ont été utilisés pour récupérer les séquences de P. jirovecii tout en éliminant les séquences contaminants et les chimères d'amplification ou d'assemblage. Cette stratégie a produit un assemblage de 8.1 Mb, qui contient 3898 gènes. Les recherches d'homologies, de cartographie des voies métaboliques et des validations manuelles ont révélé que ce génome est dépourvu (i) de la plupart des enzymes dédiées à la biosynthèse des acides aminés, et (ii) de la plupart des facteurs de virulence observés chez d'autres champignons, par exemple, le cycle du glyoxylate ainsi que des peptidases spécifiques impliquées dans la dégradation de la membrane de la cellule hôte. Les analyses appliquées aux données génomiques disponibles de Pneumocystis carinii, l'espèce infectant les rats, et de Pneumocystis murina, l'espèce infectant les souris, ont révélé les mêmes déficiences. Le séquençage du génome de T. deformans a généré un assemblage de 13.3 Mb qui contient 5735 gènes. T. deformans possède les gènes codant pour les enzymes impliquées dans la dégradation des parois cellulaires des plantes, le métabolisme secondaire, le cycle du glyoxylate, la détoxification, la biosynthèse des stérols ainsi que la biosynthèse d'hormones de plantes telles que l'acide abscissique ou l'acide indole 3-acétique. T. deformans possède également des sous-ensembles de gènes présents exclusivement chez des saprophytes ou des pathogènes de plantes, ce qui est consistent avec son mode de vie alternatif saprophyte et pathogène. Des gènes impliqués dans la conjugaison ont été identifiés. L'homothallisme de ce champignon suggère mécanisme de permutation du type conjuguant. Les analyses comparatives ont démontré que 81% des gènes de P. jirovecii sont présent chez les autres membres du sous-embranchement des Taphrinomycotina. Ces gènes sont essentiellement impliqués dans le métabolisme basai. Les gènes spécifiques au genre Pneumocystis représentent 8%, et sont impliqués dans le métabolisme de l'ARN et la signalisation. La signalisation est connue pour être cruciale pour l'interaction des espèces de Pneumocystis avec leur environnement. Les gènes propres à P. jirovecii représentent 11% et codent en majorité pour des protéines dont la fonction est inconnue. Ces gènes en conjonction avec d'autres (par exemple, les glycoprotéines de surface), pourraient être déterminants dans l'interaction de P. jirovecii avec les cellules de l'hôte humain, et être potentiellement responsable de la spécificité d'hôte. P. jirovecii possède un génome de taille réduite à faible pourcentage en GC et récupère très probablement des composés vitaux comme les acides aminés et le cholestérol à partir des poumons humains. De manière consistante, son génome code pour de nombreux transporteurs (22% de ses gènes), qui pourraient jouer un rôle essentiel dans l'acquisition de ces composés. Ces caractéristiques sont généralement observées chez les parasites obligatoires de plusieurs règnes (bactéries, protozoaires, champignons). De plus, les études épidémiologiques n'ont pas réussi à prouver l'existence d'ime forme vivant librement du champignon. Etant donné également l'absence de facteurs de virulence, nos observations suggèrent que P. jirovecii est un parasite obligatoire spécialisé dans la colonisation des poumons humains, ne causant une maladie que chez des individus ayant un système immunitaire compromis. La même conclusion est très probablement applicable à toutes les autres espèces de Pneumocystis dans leur hôte mammifère respectif.

Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

Fractional hepatic de novo lipogenesis in healthy subjects during near-continuous oral nutrition and bed rest: a comparison with published data in artificially fed, critically ill patients.

BACKGROUND AND AIMS: In critically ill patients, fractional hepatic de novo lipogenesis increases in proportion to carbohydrate administration during isoenergetic nutrition. In this study, we sought to determine whether this increase may be the consequence of continuous enteral nutrition and bed rest. We, therefore, measured fractional hepatic de novo lipogenesis in a group of 12 healthy subjects during near-continuous oral feeding (hourly isoenergetic meals with a liquid formula containing 55% carbohydrate). In eight subjects, near-continuous enteral nutrition and bed rest were applied over a 10 h period. In the other four subjects, it was extended to 34 h. Fractional hepatic de novo lipogenesis was measured by infusing(13) C-labeled acetate and monitoring VLDL-(13)C palmitate enrichment with mass isotopomer distribution analysis. Fractional hepatic de novo lipogenesis was 3.2% (range 1.5-7.5%) in the eight subjects after 10 h of near continuous nutrition and 1.6% (range 1.3-2.0%) in the four subjects after 34 h of near-continuous nutrition and bed rest. This indicates that continuous nutrition and physical inactivity do not increase hepatic de novo lipogenesis. Fractional hepatic de novo lipogenesis previously reported in critically ill patients under similar nutritional conditions (9.3%) (range 5.3-15.8%) was markedly higher than in healthy subjects (P<0.001). These data from healthy subjects indicate that fractional hepatic de novo lipogenesis is increased in critically ill patients.

Carbohydrate metabolism and de novo lipogenesis in human obesity.

Respiratory exchange was measured during 14 consecutive hours in six lean and six obese individuals after ingestion of 500 g of dextrin maltose to investigate and compare their capacity for net de novo lipogenesis. After ingestion of the carbohydrate load, metabolic rates rose similarly in both groups but fell earlier and more rapidly in the obese. RQs also rose rapidly and remained in the range of 0.95 to 1.00 for approximately 8 h in both groups. During this time, RQ exceeded 1.00 for only short periods of time with the result that 4 +/- 1 g and 5 +/- 3 g (NS) of fat were synthesized via de novo lipogenesis in excess of concomitant fat oxidation in the lean and obese subjects, respectively. Results demonstrate that net de novo lipid synthesis from an unusually large carbohydrate load is not greater in obese than in lean individuals.

De novo mutation in the KCNQ1 gene causal to Jervell and Lange-Nielsen syndrome.

Jervell and Lange-Nielsen syndrome (JLNS) is an autosomal recessive disorder, clinically characterized by severe cardiac arrhythmias [due to prolonged QTc interval in electrocardiogram (ECG)] and bilateral sensory neural deafness. Molecular defects causal to JLNS are either homozygous or compound heterozygous mutations, predominantly in the KCNQ1 gene and occasionally in the KCNE1 gene. As the molecular defect is bi-allelic, JLNS patients inherit one pathogenic mutation causal to the disorder from each parent. In this report, we show for the first time that such a disorder could also occur due to a spontaneous de novo mutation in the affected individual, not inherited from the parent, which makes this case unique unlike the previously reported JLNS cases.

Development and testing of novel chloroplast microsatellite markers for Lolium perenne and other grasses (Poaceae) from de novo sequencing and in silico sequences

Twelve primers to amplify microsatellite markers from the chloroplast genome of Lolium perenne were designed and optimized using de novo sequencing and in silico sequences. With one exception, each locus was polymorphic with a range from two to nine alleles in L. perenne. The newly developed primer pairs cross-amplified in different species of Lolium and in 50 other grass species representing nine grass subfamilies.

Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

Early occurrence of lung adenocarcinoma and breast cancer after radiotherapy of a chest wall sarcoma in a patient with a de novo germline mutation in TP53.

We report a 26-year-old female patient who was diagnosed within 4 years with chest sarcoma, lung adenocarcinoma, and breast cancer. While her family history was unremarkable, DNA sequencing of TP53 revealed a germline de novo non-sense mutation in exon 6 p.Arg213X. One year later, she further developed a contralateral ductal carcinoma in situ, and 18 months later a jaw osteosarcoma. This case illustrates the therapeutic pitfalls in the care of a young cancer patient with TP53 de novo germline mutations and the complications related to her first-line therapy. Suggestion is made to use the less stringent Chompret criteria for germline TP53 mutation screening. Our observation underlines the possibly negative effect of radiotherapy in generating second tumors in patients with a TP53 mutation. We also present a review of six previously reported cases, comparing their cancer phenotypes with those generally produced by TP53 mutations.

Recurrent de novo mitochondrial DNA mutations in respiratory chain deficiency.

Starting from a cohort of 50 NADH-oxidoreductase (complex I) deficient patients, we carried out the systematic sequence analysis of all mitochondrially encoded complex I subunits (ND1 to ND6 and ND4L) in affected tissues. This approach yielded the unexpectedly high rate of 20% mutation identification in our series. Recurrent heteroplasmic mutations included two hitherto unreported (T10158C and T14487C) and three previously reported mutations (T10191C, T12706C and A13514G) in children with Leigh or Leigh-like encephalopathy. The recurrent mutations consistently involved T-->C transitions (p<10(-4)). This study supports the view that an efficient molecular screening should be based on an accurate identification of respiratory chain enzyme deficiency.

Molecular characterization of 39 de novo sSMC: contribution to prognosis and genetic counselling, a prospective study.

Small supernumerary marker chromosomes (sSMCs) are structurally abnormal chromosomes that cannot be characterized by karyotype. In many prenatal cases of de novo sSMC, the outcome of pregnancy is difficult to predict because the euchromatin content is unclear. This study aimed to determine the presence or absence of euchromatin material of 39 de novo prenatally ascertained sSMC by array-comparative genomic hybridization (array-CGH) or single nucleotide polymorphism (SNP) array. Cases were prospectively ascertained from the study of 65,000 prenatal samples [0.060%; 95% confidence interval (CI), 0.042-0.082]. Array-CGH showed that 22 markers were derived from non-acrocentric markers (56.4%) and 7 from acrocentic markers (18%). The 10 additional cases remained unidentified (25.6%), but 7 of 10 could be further identified using fluorescence in situ hybridization; 69% of de novo sSMC contained euchromatin material, 95.4% of which for non-acrocentric markers. Some sSMC containing euchromatin had a normal phenotype (31% for non-acrocentric and 75% for acrocentric markers). Statistical differences between normal and abnormal phenotypes were shown for the size of the euchromatin material (more or less than 1 Mb, p = 0.0006) and number of genes (more or less than 10, p = 0.0009). This study is the largest to date and shows the utility of array-CGH or SNP array in the detection and characterization of de novo sSMC in a prenatal context.

Comparaison des coûts du diagnostic moléculaire des rétinites pigmentaires entre le séquençage selon Sanger et un séquençage à haut débit

La rétinite pigmentaire (RP) et l'amaurose congénitale de Leber (LCA) sont deux maladies héréditaires classées dans le groupe des rétinopathies pigmentaires. Plus de 100 gènes ou loci ont étés identifiés dans les RP (comptant pour 60% des patients) et 14 gènes pour les LCA (responsables de 70% des cas). A elles-deux et en ne prenant que les plus fréquentes, ces maladies représentent 948 exons à analyser lorsqu'on recherche une mutation chez un patient. La recherche de mutations au moyen du séquençage génomique classique (séquençage selon Sanger) de tous les gènes de ces deux maladies implique des délais d'analyse importants et des coûts très élevés. La méthode de séquençage à haut débit (avec le séquenceur Roche GS Junior) permet grâce à une puce à ADN le séquençage simultané des 948 exons. Le but de mon travail de Master est de comparer ces deux approches afin de déterminer celle qui est la plus économique et la plus efficiente en temps. Pour cela, j'ai d'abord établi la liste de tous les gènes impliqués dans les RP et LCA, puis identifié tous les exons ainsi que les promoteurs et séquences 3' non traduites. J'ai ensuite calculé le coût théorique d'une analyse de tous les gènes avec chacune des méthodes. J'ai également estimé les coûts à facturer à l'assurance concernant le séquençage à haut débit sur la base des coûts facturables à l'assurance de la méthode Sanger et des bénéfices du laboratoire. Le séquençage des 948 exons par le séquençage à haut débit (avec le GS Junior) représente la technique de séquençage la plus économique et la plus efficiente en temps et constitue donc la méthode de choix dans le screening diagnostic des gènes impliqués dans les RP et LCA. Cette méthode est plus rapide, les réactifs et la machine sont moins coûteux et la laborantine peut analyser un nombre plus important d'exons en un temps moindre, donc elle coûtera moins cher au laboratoire. Cette méthode est donc d'un grand intérêt pour les patients, les assurances et le laboratoire. Cette nouvelle technique de séquençage soulève de nouvelles interrogations telles que la décision de savoir quelle information doit être donnée aux médecins, aux assurances et aux patients. Interrogations auxquelles il devient de plus en plus pressant de répondre.

De novo characterization of the Timema cristinae transcriptome facilitates marker discovery and inference of genetic divergence.

Adaptation to different ecological environments can promote speciation. Although numerous examples of such 'ecological speciation' now exist, the genomic basis of the process, and the role of gene flow in it, remains less understood. This is, at least in part, because systems that are well characterized in terms of their ecology often lack genomic resources. In this study, we characterize the transcriptome of Timema cristinae stick insects, a system that has been researched intensively in terms of ecological speciation, but for which genomic resources have not been previously developed. Specifically, we obtained >1 million 454 sequencing reads that assembled into 84,937 contigs representing approximately 18,282 unique genes and tens of thousands of potential molecular markers. Second, as an illustration of their utility, we used these genomic resources to assess multilocus genetic divergence within both an ecotype pair and a species pair of Timema stick insects. The results suggest variable levels of genetic divergence and gene flow among taxon pairs and genes and illustrate a first step towards future genomic work in Timema.

GENETIC AND BIOCHEMICAL ANALYSIS OF THE TESTIS-SPECIFIC ZINC FINGER PROTEIN CTCFL/BORIS

Abstract en FrançaisCTCFL a d'abord été identifié comme un paralogue de la protéine ubiquitaire CTCF en raison de sa forte homologie entre leurs onze « zinc fingers », un domaine de liaison à l'ADN. Parmi ses nombreux rôles, la liaison des zinc fingers de CTCF à la région de contrôle de l'empreinte (ICR) maternelle non-méthylée Igf2/H19, contrôle l'expression empreinte (monoallélique) de H19 et IGF2 dans les cellules somatiques. La méthylation de l'ICR Igf2/H19 paternelle est nécessaire à l'expression empreinte de ces deux gènes. Bien que le mécanisme par lequel l'ICR est méthylé soit mal compris, il est connu que l'établissement de la méthylation se produit pendant le développement des cellules germinales mâles et que les ADN méthyltransférases de novo DNMT3A et DNMT3L sont essentiels. Par conséquent, CTCFL fournit un bon candidat pour un rôle dans la méthylation de l'ICR paternelle Igf2/H19 en raison de son expression restreinte à certains types de cellules où la méthylation de l'ICR a lieu (spermatogonies et spermatocytes) ainsi qu'en raison sa capacité à lier les ICR lgf2/HÎ9 dans ces cellules. Les premiers travaux expérimentaux de cette thèse portent sur le rôle possible des mutations de CTCFL chez les patients atteints du syndrome de Silver-Russell (SRS), où une diminution de la méthylation de l'ICR IGF2/H19 a été observée chez 60% d'entre eux. Admettant que CTCFL pourrait être muté chez ces patients, j'ai examiné les mutations possibles de CTCFL chez 35 d'entre eux par séquençage de l'ADN et analyse du nombre de copies d'exons. N'ayant trouvé aucune mutation chez ces patients, cela suggère que les mutations de CTCFL ne sont pas associées au SRS. Les travaux expérimentaux suivants ont porté sur les modifications post-traductionnelles de CTCFL par la protéine SU MO « small ubiquitin-like modifier » (SUMO). La modification de protéines par SU MO change les interactions avec d'autres molécules (ADN ou protéines). Comme CTCFL régule sans doute l'expression d'un certain nombre de gènes dans le cancer et que plusieurs facteurs de transcription sont régulés par SUMO, j'ai mené des expériences pour déterminer si CTCFL est sumoylé. En effet, j'ai observé que CTCFL est sumoylated in vitro et in vivo et j'ai déterminé les deux résidus d'attachement de SUMO aux lysines 181 et 645. Utilisant les mutants de CTCFL K181R et K645R ne pouvant pas être sumoylated, j'ai évalué les conséquences fonctionnelles de la modification par SUMO. Je n'ai trouvé aucun changement significatif dans la localisation subcellulaire, la demi-vie ou la liaison à l'ADN, mais ai constaté que la sumoylation module à la fois {'activation CTCFL-dépendante et la répression de l'expression génique. Il s'agit de la première modification post-traductionnelle décrite pour CTCFL et les conséquences possibles de cette modification sont discutées pour le cancer et les testicules normaux. Avec cette thèse, j'espère avoir ajouté des résultats importants à l'étude de CTCFL et donné quelques idées pour de futures recherches.AbstractJeremiah Bernier-Latmani, Institute of Pathology, University of Lausanne, CHUVCTCFL was first identified as a paralog of the ubiquitous protein CTCF because of high homology between their respective eleven zinc fingers, a DNA binding domain. Among its many roles, CTCF zinc finger-mediated binding to the unmethylated maternal Igf2/H19 imprinting control region (ICR), controls the imprinted (monoallelic) expression of Igf2 and H19 in somatic cells. Methylation of the paternal Igf2/H19 ICR is necessary for the imprinted expression of the two genes. Although the mechanism by which the ICR is methylated is incompletely understood, it is known that establishment of methylation occurs during male germ cell development and the de novo DNA methyltransferases DNMT3A and DNMT3L are essential. Therefore, CTCFL provided a good candidate to play a role in methylation of the paternal Igf2/H19 ICR because of its restricted expression to cell types where ICR methylation takes place (spermatogonia and spermatocytes) and its ability to bind the Igf2/H19 ICR in these cells. The first experimental work of this thesis investigated the possible role of CTCFL mutations in Silver-Russell syndrome (SRS) patients, where it has been observed that 60% of the patients have reduced methylation of the IGF2/HÎ9 ICR. Reasoning that CTCFL could be mutated in these patients, I screened 35 patients for mutations in CTCFL by DNA sequencing and exon copy number analysis, I did not find any mutations in these patients suggesting that mutations of CTCFL are not associated with SRS. The next experimental work of my thesis focused on posttranslational modification of CTCFL by small ubiquitin-like modifier (SUMO) protein. SUMO modification of proteins changes the interactions with other molecules (DNA or protein). As CTCFL arguably regulates the expression of a number of genes in cancer and many transcription factors are regulated by SUMO, I conducted experiments to assess whether CTCFL is sumoylated. I found that CTCFL is sumoylated in vitro and in vivo and determined the two residues of SUMO attachment to be lysines 181 and 645. Using K181R, K645R mutated CTCFL- which cannot be detected to be sumoylated-1 assessed the functional consequences of SUMO modification. I found no significant changes in subcellular localization, half-life or DNA binding, but found that sumoylation modulates both CTCFL-dependent activation and repression of gene expression. This is the first posttranslational modification described for CTCFL and possible consequences of this modification are discussed in both cancer and normal testis. With this thesis, I hope I have added important findings to the study of CTCFL and provide some ideas for future research.

Visualization and quality assessment of de novo genome assemblies.

Recent technological progress has greatly facilitated de novo genome sequencing. However, de novo assemblies consist in many pieces of contiguous sequence (contigs) arranged in thousands of scaffolds instead of small numbers of chromosomes. Confirming and improving the quality of such assemblies is critical for subsequent analysis. We present a method to evaluate genome scaffolding by aligning independently obtained transcriptome sequences to the genome and visually summarizing the alignments using the Cytoscape software. Applying this method to the genome of the red fire ant Solenopsis invicta allowed us to identify inconsistencies in 7%, confirm contig order in 20% and extend 16% of scaffolds.Scripts that generate tables for visualization in Cytoscape from FASTA sequence and scaffolding information files are publicly available at https://github.com/ksanao/TGNet.