Cardiac expression of ms1/STARS, a novel gene involved in cardiac development and disease, is regulated by GATA4.

Ms1/STARS is a novel muscle-specific actin-binding protein that specifically modulates the myocardin-related transcription factor (MRTF)-serum response factor (SRF) regulatory axis within striated muscle. This ms1/STARS-dependent regulatory axis is of central importance within the cardiac gene regulatory network and has been implicated in cardiac development and postnatal cardiac function/homeostasis. The dysregulation of ms1/STARS is associated with and causative of pathological cardiac phenotypes, including cardiac hypertrophy and cardiomyopathy. In order to gain an understanding of the mechanisms governing ms1/STARS expression in the heart, we have coupled a comparative genomic in silico analysis with reporter, gain-of-function, and loss-of-function approaches. Through this integrated analysis, we have identified three evolutionarily conserved regions (ECRs), α, SINA, and DINA, that act as cis-regulatory modules and confer differential cardiac cell-specific activity. Two of these ECRs, α and DINA, displayed distinct regulatory sensitivity to the core cardiac transcription factor GATA4. Overall, our results demonstrate that within embryonic, neonatal, and adult hearts, GATA4 represses ms1/STARS expression with the pathologically associated depletion of GATA4 (type 1/type 2 diabetic models), resulting in ms1/STARS upregulation. This GATA4-dependent repression of ms1/STARS expression has major implications for MRTF-SRF signaling in the context of cardiac development and disease.

Structure of a multisubunit complex that promotes DNA branch migration.

The RuvA and RuvB proteins of Escherichia coli, which are induced in response to DNA damage, are important in the formation of heteroduplex DNA during genetic recombination and related recombinational repair processes. In vitro studies show that RuvA binds Holiday junctions and acts as a specificity factor that targets the RuvB ATPase, a hexameric ring protein, to the junction. Together, RuvA and RuvB promote branch migration, an ATP-dependent reaction that increases the length of the heteroduplex DNA. Electron microscopic visualization of RuvAB now provides a new insight into the mechanism of this process. We observe the formation of a tripartite protein complex in which RuvA binds the crossover and is sandwiched between two hexameric rings of RuvB. The Holliday junction within this complex adopts a square-planar structure. We propose a molecular model for branch migration, a unique feature of which is the role played by the two oppositely oriented RuvB ring motors.

Disease-causing mutations improving the branch site and polypyrimidine tract: pseudoexon activation of LINE-2 and antisense Alu lacking the poly(T)-tail.

Cryptic exons or pseudoexons are typically activated by point mutations that create GT or AG dinucleotides of new 5' or 3' splice sites in introns, often in repetitive elements. Here we describe two cases of tetrahydrobiopterin deficiency caused by mutations improving the branch point sequence and polypyrimidine tracts of repeat-containing pseudoexons in the PTS gene. In the first case, we demonstrate a novel pathway of antisense Alu exonization, resulting from an intronic deletion that removed the poly(T)-tail of antisense AluSq. The deletion brought a favorable branch point sequence within proximity of the pseudoexon 3' splice site and removed an upstream AG dinucleotide required for the 3' splice site repression on normal alleles. New Alu exons can thus arise in the absence of poly(T)-tails that facilitated inclusion of most transposed elements in mRNAs by serving as polypyrimidine tracts, highlighting extraordinary flexibility of Alu repeats in shaping intron-exon structure. In the other case, a PTS pseudoexon was activated by an A>T substitution 9 nt upstream of its 3' splice site in a LINE-2 sequence, providing the first example of a disease-causing exonization of the most ancient interspersed repeat. These observations expand the spectrum of mutational mechanisms that introduce repetitive sequences in mature transcripts and illustrate the importance of intronic mutations in alternative splicing and phenotypic variability of hereditary disorders.

Regulation of horizontal gene transfer of the catabolic island ICEclc

ICEclc is a mobile genetic element found in two copies on the chromosome of the bacterium Pseudomonas knackmussii B13. ICEclc harbors genes encoding metabolic pathways for the degradation of chlorocatechols (CLC) and 2-aminophenol (2AP). At low frequencies, ICEclc excises from the chromosome, closes into a circular DNA molecule which can transfer to another bacterium via conjugation. Once in the recipient cell, ICEclc can reintegrate into the chromosome by site-specific recombination. This thesis aimed at identifying the regulatory network underlying the decisions for ICEclc horizontal transfer (HGT). The first chapter is an introduction on integrative and conjugative elements (ICEs) more in general, of which ICEclc is one example. In particular I emphasized the current knowledge of regulation and conjugation machineries of the different classes of ICE. In the second chapter, I describe a transcriptional analysis using microarrays and other experiments to understand expression of ICEclc in exponential and stationary phase. By overlaying transcriptomic profiles with Northern hybridizations and RT- PCR data, we established a transcription map for the entire core region of ICEclc, a region assumed to encode the ICE conjugation process. We also demonstrated how transcription of the ICEclc core is maximal in stationary phase, which correlates to expression of reporter genes fused to key ICEclc promoters. In the third chapter, I present a transcriptome analysis of ICEclc in a variety of different host species, in order to explore whether there are species-specific differences. In the fourth chapter, I focus on the role of a curious ICEclc-encoded TetR-type transcriptional repressor. We find that this gene, which we name mfsR, not only controls its own expression but that of a set of genes for a putative multi-drug efflux pump (mfsABC) as well. By using a combination of biochemical and molecular biology techniques, I could show that MfsR specifically binds to operator boxes in two ICEclc promoters (PmfsR and PmfsA), inhibiting the transcription of both the mfsR and mfsABC-orf38184 operons. Although we could not detect a clear phenotype of an mfsABC deletion, we discuss the implications of pump gene reorganizations in ICEclc and close relatives. In the fifth chapter, we find that mfsR not only controls its own expression and that of the mfsABC operon, but is also indirectly controlling ICEclc transfer. Using gene deletions, microarrays, transfer assays and microscopy-based reporter fusions, we demonstrate that mfsR actually controls a small operon of three regulatory genes. The last gene of this mfsR operon, orf17162, encodes a LysR-type activator that when deleted strongly impairs ICEclc transfer. Interestingly, deletion of mfsR leads to transfer competence in almost all cells, thereby overruling the bistability process in the wild-type. In the final sixth chapter, I discuss the relevance of the present thesis and the resulting perspectives for future studies.

The E.coli RuvAB proteins branch migrate Holliday junctions through heterologous DNA sequences in a reaction facilitated by SSB.

During genetic recombination a heteroduplex joint is formed between two homologous DNA molecules. The heteroduplex joint plays an important role in recombination since it accommodates sequence heterogeneities (mismatches, insertions or deletions) that lead to genetic variation. Two Escherichia coli proteins, RuvA and RuvB, promote the formation of heteroduplex DNA by catalysing the branch migration of crossovers, or Holliday junctions, which link recombining chromosomes. We show that RuvA and RuvB can promote branch migration through 1800 bp of heterologous DNA, in a reaction facilitated by the presence of E.coli single-stranded DNA binding (SSB) protein. Reaction intermediates, containing unpaired heteroduplex regions bound by SSB, were directly visualized by electron microscopy. In the absence of SSB, or when SSB was replaced by a single-strand binding protein from bacteriophage T4 (gene 32 protein), only limited heterologous branch migration was observed. These results show that the RuvAB proteins, which are induced as part of the SOS response to DNA damage, allow genetic recombination and the recombinational repair of DNA to occur in the presence of extensive lengths of heterology.

Remodeling of DNA replication structures by the branch point translocase FANCM

Fanconi anemia (FA) is a genetically heterogeneous chromosome instability syndrome associated with congenital abnormalities, bone marrow failure, and cancer predisposition. Eight FA proteins form a nuclear core complex, which promotes tolerance of DNA lesions in S phase, but the underlying mechanisms are still elusive. We reported recently that the FA core complex protein FANCM can translocate Holliday junctions. Here we show that FANCM promotes reversal of model replication forks via concerted displacement and annealing of the nascent and parental DNA strands. Fork reversal by FANCM also occurs when the lagging strand template is partially single-stranded and bound by RPA. The combined fork reversal and branch migration activities of FANCM lead to extensive regression of model replication forks. These observations provide evidence that FANCM can remodel replication fork structures and suggest a mechanism by which FANCM could promote DNA damage tolerance in S phase

New electrocardiographic criteria for discriminating between Brugada types 2 and 3 patterns and incomplete right bundle branch block.

OBJECTIVES: The aim of this study was to evaluate new electrocardiographic (ECG) criteria for discriminating between incomplete right bundle branch block (RBBB) and the Brugada types 2 and 3 ECG patterns. BACKGROUND: Brugada syndrome can manifest as either type 2 or type 3 pattern. The latter should be distinguished from incomplete RBBB, present in 3% of the population. METHODS: Thirty-eight patients with either type 2 or type 3 Brugada pattern that were referred for an antiarrhythmic drug challenge (AAD) were included. Before AAD, 2 angles were measured from ECG leads V(1) and/or V(2) showing incomplete RBBB: 1) α, the angle between a vertical line and the downslope of the r'-wave, and 2) β, the angle between the upslope of the S-wave and the downslope of the r'-wave. Baseline angle values, alone or combined with QRS duration, were compared between patients with negative and positive results on AAD. Receiver-operating characteristic curves were constructed to identify optimal discriminative cutoff values. RESULTS: The mean β angle was significantly smaller in the 14 patients with negative results on AAD compared to the 24 patients with positive results on AAD (36 ± 20° vs. 62 ± 20°, p < 0.01). Its optimal cutoff value was 58°, which yielded a positive predictive value of 73% and a negative predictive value of 87% for conversion to type 1 pattern on AAD; α was slightly less sensitive and specific compared with β. When the angles were combined with QRS duration, it tended to improve discrimination. CONCLUSIONS: In patients with suspected Brugada syndrome, simple ECG criteria can enable discrimination between incomplete RBBB and types 2 and 3 Brugada patterns.

Horizontal, but not vertical, biotic interactions affect fine-scale plant distribution patterns in a low-energy system

Studies of species range determinants have traditionally focused on abiotic variables (typically climatic conditions), and therefore the recent explicit consideration of biotic interactions represents an important advance in the field. While these studies clearly support the role of biotic interactions in shaping species distributions, most examine only the influence of a single species and/or a single interaction, failing to account for species being subject to multiple concurrent interactions. By fitting species distribution models (SDMs), we examine the influence of multiple vertical (i.e., grazing, trampling, and manuring by mammalian herbivores) and horizontal (i.e., competition and facilitation; estimated from the cover of dominant plant species) interspecific interactions on the occurrence and cover of 41 alpine tundra plant species. Adding plant-plant interactions to baseline SDMs (using five field-quantified abiotic variables) significantly improved models' predictive power for independent data, while herbivore-related variables had only a weak influence. Overall, abiotic variables had the strongest individual contributions to the distribution of alpine tundra plants, with the importance of horizontal interaction variables exceeding that of vertical interaction variables. These results were consistent across three modeling techniques, for both species occurrence and cover, demonstrating the pattern to be robust. Thus, the explicit consideration of multiple biotic interactions reveals that plant-plant interactions exert control over the fine-scale distribution of vascular species that is comparable to abiotic drivers and considerably stronger than herbivores in this low-energy system.

Describing ancient horizontal gene transfers at the nucleotide and gene levels by comparative pathogenicity island genometrics.

MOTIVATION: Lateral gene transfer is a major mechanism contributing to bacterial genome dynamics and pathovar emergence via pathogenicity island (PAI) spreading. However, since few of these genomic exchanges are experimentally reproducible, it is difficult to establish evolutionary scenarios for the successive PAI transmissions between bacterial genera. Methods initially developed at the gene and/or nucleotide level for genomics, i.e. comparisons of concatenated sequences, ortholog frequency, gene order or dinucleotide usage, were combined and applied here to homologous PAIs: we call this approach comparative PAI genometrics. RESULTS: YAPI, a Yersinia PAI, and related islands were compared with measure evolutionary relationships between related modules. Through use of our genometric approach designed for tracking codon usage adaptation and gene phylogeny, an ancient inter-genus PAI transfer was oriented for the first time by characterizing the genomic environment in which the ancestral island emerged and its subsequent transfers to other bacterial genera.

Excision of the Staphylococcal Cassette Chromosome mec (SCCmec) and its crucial role in the horizontal transfer of methicillin resistance in staphylococci.

Staphylococcus aureus est un pathogène humain majeur ayant développé des résistances contre la quasi totalité des antibiotiques disponibles, incluant la très importante famille des β- lactamines. La résistance à cette classe d'antibiotiques est conférée par la « Staphylococcal Cassette Chromosome mec » (SCCmec), qui est un élément génétique mobile capable de s'insérer dans le chromosome bactérien et capable d'être transféré horizontalement chez d'autres staphylocoques. Le mécanisme moléculaire impliqué dans ce transfert horizontal demeure largement inconnu. L'une des premières étapes du transfert est l'excision du SCC mec du chromosome bactérien. Cette excision est promue par des enzymes codées par l'élément SCCmec lui- même et appelées de ce fait « Cassette Chromosome Recombinases » (Ccr). L'un des buts de ce travail de thèse a été de comprendre la régulation de l'expression des gènes codant pour les Ccr recombinases. En utilisant des outils moléculaires originaux, nous avons été en mesure de démontrer en premier lieu que les Ccr recombinases étaient exprimées de façon « bistable », c'est à dire qu'uniquement quelques pourcents de cellules dans une population exprimaient ces gènes à un temps donné. Dans un deuxième temps, nous avons également démontré que l'expression de ces gènes était régulée par des facteurs étrangers au SCC mec. L'expression bistable des recombinases est un concept important. Effectivement, cela permet à la majorité des cellules d'une population de conserver l'élément SCC mec, alors que seulement une petite fraction le perd afin de le rendre disponible pour un transfert. Ainsi, alors que l'élément SCC mec continue de se propager avec la multiplication des bactéries Staphylococcus aureus résistant à la méticilline (SARM), il peut être simultanément transmis à des souches susceptibles (Staphylococcus aureus susceptible à la méticilline, SASM), entraînant l'apparition de nouveaux SARM. De façon très intéressante, le fait que cette bistabilité est contrôlée par les bactéries, et non le SCCmec lui-même, montre que la décision de transférer ou non la cassette SCC mec appartient à la bactérie. En conséquence, il doit exister dans la nature des souches qui sont plus ou moins aptes à effectuer ce transfert. En nous appuyant sur ces observations, nous avons montré que l'excision du SCC mec était effectivement régulée de façon très étroite au cours de la division cellulaire, et ne se passait que pendant un temps limité au début de la croissance. Ce résultat est compatible avec une régulation génétique commandée par la densité cellulaire, qui pourrait être dépendante de la production de signaux extracellulaires, du type que l'on rencontre dans le quorum sensing. Les signaux hypothétiques entraînant l'excision du SCC mec restent inconnus à l'heure actuelle. La connaissance de ces signaux pourrait se révéler très importante afin de développer des stratégies pour interférer avec la dissémination de la résistance au β-lactamines. Deux sujets additionnels ont été logiquement investigués au vu de ces premiers résultats. Premièrement, si certaines souches de SARM sont plus ou moins aptes à déclencher l'excision du SCC mec, de même certaines souches de SASM devraient être plus ou moins aptes à acquérir cet élément. Deuxièmement, afin d'étudier ces mécanismes de transfert au niveau épidémiologique, il nous a été nécessaire de développer des outils nous permettant d'explorer le phénomène à une plus large échelle. Concernant le premier point, il a été postulé que certains SASM seraient réfractaires à l'intégration génomique d'un SCC mec en raison de polymorphismes particuliers à proximité du site d'insertion chromosomique (attB). En étudiant plus de 40 isolais de S. aureus, provenant de porteurs sains, nous avons confirmé ce polymorphisme dans l'environnement à'attB. De plus, nous avons pu montrer que ces régions polymorphiques ont évolué parallèlement à des groupes phylogénétiques bien connus. Ainsi, si des telles régions réfractaires à l'intégration de SCC mec existent, celles-ci devraient ségréger dans des complexes clonaux bien définis qui devraient être facilement identifiables au niveau épidémiologique. Concernant le second point, nous avons été capables de construire un système rapporteur de l'excision du SCCmec, en utilisant un plasmide à faible copie. Ce système consistait en un promoteur fort et un gène codant pour une protéine verte fluorescente (GFP) sous le contrôle d'un promoteur fort séparés à l'aide d'un élément SCC artificiel portant trois terminateurs de transcription. Ainsi, la fluorescence ne s'exprime que si l'élément SCC est excisé du plasmide. Ce système a été testé avec succès dans plusieurs types de staphylocoques, et est actuellement évalué dans d'autres souches et conditions stimulant ou inhibant l'excision. De manière générale, cette dissertation représente parcours scientifique à travers plusieurs aspects d'un problème de santé publique majeur en rapport avec la résistance bactérienne aux antibiotiques. Ce travail s'attaque à des problèmes fondamentaux concernant le transfert horizontal de l'élément SCC mec. De plus, il s'intéresse à des aspects plus généraux de cet élément génétique mobile qui pourraient se révéler très importants en terme de mouvement de gènes au sein des staphylocoques, voir d'autres bactéries gram-positives. Finalement ce travail de thèse met en place le fondamentaux requis pour des recherches futures visant à interférer avec le transfert horizontal de la résistance aux β-lactamines. - Staphylococcus aureus is a major human pathogen. Moreover, S. aureus have developed resistance to almost all available antibiotics, including the important family of β-lactam molecules. Intrinsic resistance to β-lactams is conferred by the Staphylococcal Cassette Chromosome mec (SCCmec), which is a mobile genomic island that inserts into the staphylococcal chromosome and can be horizontally transferred into other staphylococci. However, little is known about the molecular mechanisms involved in this horizontal transfer into naïve strains. One of the first steps in SCC mec horizontal transfer is its excision from the chromosome. Excision is mediated by recombinase enzymes that are encoded by SCC mec itself, and named accordingly Ccr recombinases - for Cassette Chromosome recombinases. One goal of this thesis was to understand the regulation these recombinase genes. By using original molecular tools we could demonstrate first that the Ccr recombinases were expressed in a "bistable" manner, i.e. in only few percentages of the bacterial cells at a given time, and second that they were regulated by determinants that were not encoded on the SCC mec element, but elsewhere on the staphylococcal genome. "Bistable" expression Ccr recombinases is an important concept. It allows SCC mec to be excised and thus available for horizontal transfer, while ensuring that only some cells, but not the whole population, loose their valuable SCC mec genes. Thus, while the SCC mec element expands with the multiplication of the MRSA colony, it can simultaneously be transmitted into methicillin-susceptible S. aureus (MSSA), which convert into new MRSA. Most interestingly, the fact that bistability was regulated by the cells, rather than by SCC mec, indicates that it was the choice of the bacteria to trigger or not SCC mec transfer. As a consequence, there must be, in nature, staphylococcal strains that are more or less prone to sustain SCC mec transfer. Following these seminal observations we found that excision was indeed tightly regulated during bacterial division, and occurred only during a limited period of time at the beginning of bacterial growth. This is compatible with cell-density mediated gene regulation, and may depend on the production of extracellular signal molecules that transmit appropriate orders to neighboring cells, such as in quorum sensing. The potential signal triggering SCCmec excision is as yet unknown. However, it could be critical in promoting the horizontal transfer of methicillin resistance, or for the possible development of means to interfere with it. Two additional hypothesis were logically investigated in the view of these first results. First, if some strains of MRSA might be more prone than others to promote SCC mec excision, then some strains of MS SA might be more or less prone to acquire the element as well. Second, to investigate these multiple mechanisms at an epidemiological level, one would need to develop tools amenable to explore S. aureus strains at a larger scale. Regarding the first issue, it was postulated by others that some MSSA might be refractory to SCC mec integration because they had peculiar DNA polymorphisms in the vicinity of the site-specific chromosomal entry point {attB) of SCC mec. By studying >40 S. aureus isolates from healthy carriers, we confirmed the polymorphism of the attB environment. Moreover, we could show that these polymorphic regions co-evolved with well-known phylogenic clonal clusters. Therefore, if SCCwec-refractory attB environments exist, then they would segregate in well- defined S. aureus clonal clusters that would be easy to identify at the epidemiological level. Regarding the second issue, we were able to construct a new excision reporter system in a low copy number S. aureus plasmid. The reporter system consists in a strong promoter driving a green fluorescent protein {gfp) gene, separated by an artificial SCC-like element carrying three transcriptional terminators. Thus, fluorescence is not expressed unless the SCC-like element is excised. The system has been successfully tested in several aureus and non- aureus staphylococci, and is now being applied to more strains and various excision- triggering or inhibiting conditions. Altogether the dissertation is a scientific journey through various aspects of a salient medical problem with regard to antibiotic resistance and public health threat. The research work tackles fundamental issues about the mechanisms of horizontal transfer of the SCC mec element. Moreover, it also addresses more general features of this mobile element, which could be of larger importance with regard to gene trafficking in staphylococci, and maybe other gram-positive bacteria. Finally, the dissertation sets the fundamentals for future work and possible new ways to interfere with the horizontal transfer of methicillin resistance.

Organization of lin genes and IS6100 among different strains of hexachlorocyclohexane-degrading Sphingomonas paucimobilis: evidence for horizontal gene transfer.

The organization of lin genes and IS6100 was studied in three strains of Sphingomonas paucimobilis (B90A, Sp+, and UT26) which degraded hexachlorocyclohexane (HCH) isomers but which had been isolated at different geographical locations. DNA-DNA hybridization data revealed that most of the lin genes in these strains were associated with IS6100, an insertion sequence classified in the IS6 family and initially found in Mycobacterium fortuitum. Eleven, six, and five copies of IS6100 were detected in B90A, Sp+, and UT26, respectively. IS6100 elements in B90A were sequenced from five, one, and one regions of the genomes of B90A, Sp+, and UT26, respectively, and were found to be identical. DNA-DNA hybridization and DNA sequencing of cosmid clones also revealed that S. paucimobilis B90A contains three and two copies of linX and linA, respectively, compared to only one copy of these genes in strains Sp+ and UT26. Although the copy number and the sequence of the remaining genes of the HCH degradative pathway (linB, linC, linD, and linE) were nearly the same in all strains, there were striking differences in the organization of the linA genes as a result of replacement of portions of DNA sequences by IS6100, which gave them a strange mosaic configuration. Spontaneous deletion of linD and linE from B90A and of linA from Sp+ occurred and was associated either with deletion of a copy of IS6100 or changes in IS6100 profiles. The evidence gathered in this study, coupled with the observation that the G+C contents of the linA genes are lower than that of the remaining DNA sequence of S. paucimobilis, strongly suggests that all these strains acquired the linA gene through horizontal gene transfer mediated by IS6100. The association of IS6100 with the rest of the lin genes further suggests that IS6100 played a role in shaping the current lin gene organization.

Helicase-defective RuvB(D113E) promotes RuvAB-mediated branch migration in vitro.

In Escherichia coli, the RuvA and RuvB proteins interact at Holliday junctions to promote branch migration leading to the formation of heteroduplex DNA. RuvA provides junction-binding specificity and RuvB drives ATP-dependent branch migration. Since RuvB contains sequence motifs characteristic of a DNA helicase and RuvAB exhibit helicase activity in vitro, we have analysed the role of DNA unwinding in relation to branch migration. A mutant RuvB protein, RuvB(D113E), mutated in helicase motif II (the DExx box), has been purified to homogeneity. The mutant protein forms hexameric rings on DNA similar to those formed by wild-type protein and promotes branch migration in the presence of RuvA. However, RuvB(D113E) exhibits reduced ATPase activity and is severely compromised in its DNA helicase activity. Models for RuvAB-mediated branch migration that invoke only limited DNA unwinding activity are proposed.

The Fanconi anemia protein FANCM can promote branch migration of Holliday junctions and replication forks.

Fanconi anemia (FA) is a genetically heterogeneous cancer-prone disorder associated with chromosomal instability and cellular hypersensitivity to DNA crosslinking agents. The FA pathway is suspected to play a crucial role in the cellular response to DNA replication stress. At a molecular level, however, the function of most of the FA proteins is unknown. FANCM displays DNA-dependent ATPase activity and promotes the dissociation of DNA triplexes, but the physiological significance of this activity remains elusive. Here we show that purified FANCM binds to Holliday junctions and replication forks with high specificity and promotes migration of their junction point in an ATPase-dependent manner. Furthermore, we provide evidence that FANCM can dissociate large recombination intermediates, via branch migration of Holliday junctions through 2.6 kb of DNA. Our data suggest a direct role for FANCM in DNA processing, consistent with the current view that FA proteins coordinate DNA repair at stalled replication forks.

SlimCodeML: An Optimized Version of CodeML for the Branch-Site Model

CodeML (part of the PAML package) im- plements a maximum likelihood-based approach to de- tect positive selection on a specific branch of a given phylogenetic tree. While CodeML is widely used, it is very compute-intensive. We present SlimCodeML, an optimized version of CodeML for the branch-site model. Our performance analysis shows that SlimCodeML substantially outperforms CodeML (up to 9.38 times faster), especially for large-scale genomic analyses.