Myc and Max: molecular evolution of a family of proto-oncogene products and their dimerization partner.

The myc gene family encodes a group of transcription factors that regulate cell proliferation and differentiation. These genes are widely studied because of their importance as proto-oncogenes. Phylogenetic analyses are described here for 45 Myc protein sequences representing c-, N-, L-, S-, and B-myc genes. A gene duplication early in vertebrate evolution produced the c-myc lineage and another lineage that later gave rise to the N- and L-myc lineages by another gene duplication. Evolutionary divergence in the myc gene family corresponds closely to the known branching order of the major vertebrate groups. The patterns of sequence evolution are described for five separate highly conserved regions, and these analyses show that differential rates of sequence divergence (= mosaic evolution) have occurred among conserved motifs. Further, the closely related dimerization partner protein Max exhibits significantly less sequence variability than Myc. It is suggested that the reduced variability in max stems from natural selection acting to preserve dimerization capability with products of myc and related genes.

A VLT/FLAMES survey for massive binaries in Westerlund 1. III. The WC9d binary W239 and implications for massive stellar evolution

Context. There is growing evidence that a treatment of binarity amongst OB stars is essential for a full theory of stellar evolution. However the binary properties of massive stars – frequency, mass ratio & orbital separation – are still poorly constrained. Aims. In order to address this shortcoming we have undertaken a multiepoch spectroscopic study of the stellar population of the young massive cluster Westerlund 1. In this paper we present an investigation into the nature of the dusty Wolf-Rayet star and candidate binary W239. Methods. To accomplish this we have utilised our spectroscopic data in conjunction with multi-year optical and near-IR photometric observations in order to search for binary signatures. Comparison of these data to synthetic non-LTE model atmosphere spectra were used to derive the fundamental properties of the WC9 primary. Results. We found W239 to have an orbital period of only ~5.05 days, making it one of the most compact WC binaries yet identified. Analysis of the long term near-IR lightcurve reveals a significant flare between 2004-6. We interpret this as evidence for a third massive stellar component in the system in a long period (>6 yr), eccentric orbit, with dust production occuring at periastron leading to the flare. The presence of a near-IR excess characteristic of hot (~1300 K) dust at every epoch is consistent with the expectation that the subset of persistent dust forming WC stars are short (<1 yr) period binaries, although confirmation will require further observations. Non-LTE model atmosphere analysis of the spectrum reveals the physical properties of the WC9 component to be fully consistent with other Galactic examples. Conclusions. The simultaneous presence of both short period Wolf-Rayet binaries and cool hypergiants within Wd 1 provides compelling evidence for a bifurcation in the post-Main Sequence evolution of massive stars due to binarity. Short period O+OB binaries will evolve directly to the Wolf-Rayet phase, either due to an episode of binary mediated mass loss – likely via case A mass transfer or a contact configuration – or via chemically homogenous evolution. Conversely, long period binaries and single stars will instead undergo a red loop across the HR diagram via a cool hypergiant phase. Future analysis of the full spectroscopic dataset for Wd 1 will constrain the proportion of massive stars experiencing each pathway; hence quantifying the importance of binarity in massive stellar evolution up to and beyond supernova and the resultant production of relativistic remnants.

Rapid evolution of noncoding RNAs: lack of conservation does not mean lack of function

The mammalian transcriptome contains many nonprotein-coding RNAs (ncRNAs), but most of these are of unclear significance and lack strong sequence conservation, prompting suggestions that they might be non-functional. However, certain long functional ncRNAs such as Air and Xist are also poorly conserved. In this article, we systematically analyzed the conservation of several groups of functional ncRNAs, including miRNAs, snoRNAs and longer ncRNAs whose function has been either documented or confidently predicted. As expected, miRNAs and snoRNAs were highly conserved. By contrast, the longer functional non-micro, non-sno ncRNAs were much less conserved with many displaying rapid sequence evolution. Our findings suggest that longer ncRNAs are under the influence of different evolutionary constraints and that the lack of conservation displayed by the thousands of candidate ncRNAs does not necessarily signify an absence of function.

Molecular Systematics and Traits Evolution in Phasmatodea Leach, 1815 (Hexapoda; Insecta)

Phasmatodea Leach, 1815 (Hexapoda; Insecta) is a polyneopteran order which counts approximately 3000 described species, often known for their remarkable forms of mimicry. In this thesis, I provide a comprehensive systematic framework which includes over 180 species never considered in a phylogenetic framework: the latter can facilitate a better understanding of the processes underlying phasmids evolutionary history. The clade represents in fact an incredible testing ground to study trait evolution and its striking disparity of reproductive strategies and wing morphologies have been of great interest to the evolutionary biology community. Phasmids wings represent one of the first and most notable rejection of Dollo’s law and they played a central role in initiating a long- standing debate on the irreversibility of complex traits loss. Macroevolutionary analyses presented here confirm that wings evolution in phasmids is a reversible process even when possible biases - such as systematic uncertainty and trait-dependent diversification rates - are considered. These findings remark how complex traits can evolve in a dynamic, reversible manner and imply that their molecular groundplan can be preserved despite its phenotypical absence. This concept has been further tested with phylogenetic and transcriptomic approaches in two phasmids parthenogenetic lineages and a bisexual congeneric of the European Bacillus species complex. Leveraging a gene co-expression network approach, male gonad associated genes were retrieved in the bisexual species and then their modifications in the parthenogens were charachterized. Pleiotropy appears to constrain gene modifications associated to male reproductive structures after their loss in parthenogens, so that the lost trait molecular groundplan can be largely preserved in both transcription patterns and sequence evolution. Overall, the results presented in this thesis contribute to shape our understanding of the interplay between the phenotypic and molecular levels in trait evolution.

Mitochondrial Inheritance, Mito-nuclear Coevolution, and Sex-associated Genes in Bivalve Molluscs: Transcriptomics and Molecular Evolution

Bivalvia represents an ancient taxon including around 25,000 living species that have adapted to a wide range of environmental conditions, and show a great diversity in body size, shell shapes, and anatomic structure. Bivalves are characterized by highly variable genome sizes and extremely high levels of heterozygosity, which obstacle complete and accurate genome assemblies and hinder further genomic studies. Moreover, some bivalve species presented a stable evolutionary exception to the strictly maternal inheritance of mitochondria, namely doubly uniparental inheritance (DUI), making these species a precious model to study mitochondrial biology. During my PhD, I focused on a DUI species, the Manila clam Ruditapes philippinarum, and my work was two-folded. First, taking advantage of a newly assembled draft genome and a large RNA-seq dataset from different tissues of both sexes, I investigated 1) the role of gene expression and alternative splicing in tissue differentiation; 2) the relationship across tissue specificity, regulatory network connectivity, and sequence evolution; 3) sexual contrasting genetic markers potentially associated with sexual differentiation. The detailed information for this part is in Chapter 2. Second, using the same RNA-seq data, I investigated how nuclear oxidative phosphorylation (OXPHOS) genes coordinate with two divergent mitochondrial genomes in DUI species (mito-nuclear coordination and coevolution). To address this question, I compared transcription, polymorphism, and synonymous codon usage in the mitochondrial and nuclear OXPHOS genes of R. philippinarum in Chapter 3. To my knowledge, this thesis represents the first study exploring the role of alternative splicing in tissue differentiation, and the first study analyzing both transcriptional regulation and sequence evolution to investigate the coordination of OXPHOS genes in bivalves.

Transcriptomic basis for an antiserum against Micrurus corallinus (coral snake) venom

Background: Micrurus corallinus (coral snake) is a tropical forest snake belonging to the family Elapidae. Its venom shows a high neurotoxicity associated with pre- and post-synaptic toxins, causing diaphragm paralysis, which may result in death. In spite of a relatively small incidence of accidents, serum therapy is crucial for those bitten. However, the adequate production of antiserum is hampered by the difficulty in obtaining sufficient amounts of venom from a small snake with demanding breeding conditions. In order to elucidate the molecular basis of this venom and to uncover possible immunogens for an antiserum, we generated expressed sequences tags (ESTs) from its venom glands and analyzed the transcriptomic profile. In addition, their immunogenicity was tested using DNA immunization. Results: A total of 1438 ESTs were generated and grouped into 611 clusters. Toxin transcripts represented 46% of the total ESTs. The two main toxin classes consisted of three-finger toxins (3FTx) (24%) and phospholipases A(2) (PLA(2)s) (15%). However, 8 other classes of toxins were present, including C-type lectins, natriuretic peptide precursors and even high-molecular mass components such as metalloproteases and L-amino acid oxidases. Each class included an assortment of isoforms, some showing evidence of alternative splicing and domain deletions. Five antigenic candidates were selected (four 3FTx and one PLA(2)) and used for a preliminary study of DNA immunization. The immunological response showed that the sera from the immunized animals were able to recognize the recombinant antigens. Conclusion: Besides an improvement in our knowledge of the composition of coral snake venoms, which are very poorly known when compared to Old World elapids, the expression profile suggests abundant and diversified components that may be used in future antiserum formulation. As recombinant production of venom antigens frequently fails due to complex disulfide arrangements, DNA immunization may be a viable alternative. In fact, the selected candidates provided an initial evidence of the feasibility of this approach, which is less costly and not dependent on the availability of the venom.

Intra- and interspecific divergence in the nuclear sequences of the clock gene period in species of the Drosophila buzzatii cluster

We have measured nucleotide variation in the CLOCK/CYCLE heterodimer inhibition domain (CCID) of the clock X-linked gene period in seven species belonging to the Drosophila buzzatii cluster, namely D. buzzatii, Drosophila koepferae, Drosophila antonietae, Drosophila serido, Drosophila gouveai, Drosophila seriema and Drosophila borborema. We detected that the purifying selection is the main force driving the sequence evolution in period, in agreement with the important role of CCID in clock machinery. Our survey revealed that period provides valuable phylogenetic information that allowed to resolve phylogenetic relationships among D. gouveai, D. borborema and D. seriema, which composed a polytomic clade in preliminary studies. The analysis of patterns of intraspecific variation revealed two different lineages of period in D. koepferae, probably reflecting introgressive hybridization from D. buzzatii, in concordance with previous molecular data.

Molecular rates parallel diversification contrasts between carnivorous plant sister lineages

In the carnivorous plant family Lentibulariaceae, the bladderwort lineage (Utricularia and Genlisea) is substantially more species-rich and morphologically divergent than its sister lineage, the butterworts (Pinguicula). Bladderworts have a relaxed body plan that has permitted the evolution of terrestrial, epiphytic, and aquatic forms that capture prey in intricately designed suction bladders or corkscrew-shaped lobster-pot traps. In contrast, the flypaper-trapping butterworts maintain vegetative structures typical of angiosperms. We found that bladderwort genomes evolve significantly faster across seven loci (the trnL intron, the second trnL exon, the trnL-F intergenic spacer, the rps16 intron, rbcL, coxI, and 5.8S rDNA) representing all three genomic compartments. Generation time differences did not show a significant association. We relate these findings to the contested speciation rate hypothesis, which postulates a relationship between increased nucleotide substitution and increased cladogenesis. (C) 2002 The Willi Hennig Society.

Evolutionary genomics and HIV restriction factors.

PURPOSE OF REVIEW: To provide updated insights into innate antiviral immunity and highlight prototypical evolutionary features of well characterized HIV restriction factors. RECENT FINDINGS: Recently, a new HIV restriction factor, Myxovirus resistance 2, has been discovered and the region/residue responsible for its activity identified using an evolutionary approach. Furthermore, IFI16, an innate immunity protein known to sense several viruses, has been shown to contribute to the defense to HIV-1 by causing cell death upon sensing HIV-1 DNA. SUMMARY: Restriction factors against HIV show characteristic signatures of positive selection. Different patterns of accelerated sequence evolution can distinguish antiviral strategies--offense or defence--as well as the level of specificity of the antiviral properties. Sequence analysis of primate orthologs of restriction factors serves to localize functional domains and sites responsible for antiviral action. We use recent discoveries to illustrate how evolutionary genomic analyses help identify new antiviral genes and their mechanisms of action.

Improving the performance of positive selection inference by filtering unreliable alignment regions.

Errors in the inferred multiple sequence alignment may lead to false prediction of positive selection. Recently, methods for detecting unreliable alignment regions were developed and were shown to accurately identify incorrectly aligned regions. While removing unreliable alignment regions is expected to increase the accuracy of positive selection inference, such filtering may also significantly decrease the power of the test, as positively selected regions are fast evolving, and those same regions are often those that are difficult to align. Here, we used realistic simulations that mimic sequence evolution of HIV-1 genes to test the hypothesis that the performance of positive selection inference using codon models can be improved by removing unreliable alignment regions. Our study shows that the benefit of removing unreliable regions exceeds the loss of power due to the removal of some of the true positively selected sites.

A generalized mechanistic codon model.

Models of codon evolution have attracted particular interest because of their unique capabilities to detect selection forces and their high fit when applied to sequence evolution. We described here a novel approach for modeling codon evolution, which is based on Kronecker product of matrices. The 61 × 61 codon substitution rate matrix is created using Kronecker product of three 4 × 4 nucleotide substitution matrices, the equilibrium frequency of codons, and the selection rate parameter. The entities of the nucleotide substitution matrices and selection rate are considered as parameters of the model, which are optimized by maximum likelihood. Our fully mechanistic model allows the instantaneous substitution matrix between codons to be fully estimated with only 19 parameters instead of 3,721, by using the biological interdependence existing between positions within codons. We illustrate the properties of our models using computer simulations and assessed its relevance by comparing the AICc measures of our model and other models of codon evolution on simulations and a large range of empirical data sets. We show that our model fits most biological data better compared with the current codon models. Furthermore, the parameters in our model can be interpreted in a similar way as the exchangeability rates found in empirical codon models.

Differentiated evolutionary rates in alternative exons and the implications for splicing regulation

Background: Alternatively spliced exons play an important role in the diversification of gene function in most metazoans and are highly regulated by conserved motifs in exons and introns. Two contradicting properties have been associated to evolutionary conserved alternative exons: higher sequence conservation and higher rate of non-synonymous substitutions, relative to constitutive exons. In order to clarify this issue, we have performed an analysis of the evolution of alternative and constitutive exons, using a large set of protein coding exons conserved between human and mouse and taking into account the conservation of the transcript exonic structure. Further, we have also defined a measure of the variation of the arrangement of exonic splicing enhancers (ESE-conservation score) to study the evolution of splicing regulatory sequences. We have used this measure to correlate the changes in the arrangement of ESEs with the divergence of exon and intron sequences. Results: We find evidence for a relation between the lack of conservation of the exonic structure and the weakening of the sequence evolutionary constraints in alternative and constitutive exons. Exons in transcripts with non-conserved exonic structures have higher synonymous (dS) and non-synonymous (dN) substitution rates than exons in conserved structures. Moreover, alternative exons in transcripts with non-conserved exonic structure are the least constrained in sequence evolution, and at high EST-inclusion levels they are found to be very similar to constitutive exons, whereas alternative exons in transcripts with conserved exonic structure have a dS significantly lower than average at all EST-inclusion levels. We also find higher conservation in the arrangement of ESEs in constitutive exons compared to alternative ones. Additionally, the sequence conservation at flanking introns remains constant for constitutive exons at all ESE-conservation values, but increases for alternative exons at high ESE-conservation values. Conclusion: We conclude that most of the differences in dN observed between alternative and constitutive exons can be explained by the conservation of the transcript exonic structure. Low dS values are more characteristic of alternative exons with conserved exonic structure, but not of those with non-conserved exonic structure. Additionally, constitutive exons are characterized by a higher conservation in the arrangement of ESEs, and alternative exons with an ESE-conservation similar to that of constitutive exons are characterized by a conservation of the flanking intron sequences higher than average, indicating the presence of more intronic regulatory signals.

The evolutionary conserved BER1 gene is involved in microtubule stability in yeast

In yeast, microtubules are dynamic filaments necessary for spindle and nucleus positioning, as well as for proper chromosome segregation. We identify a function for the yeast gene BER1 (Benomyl REsistant 1) in microtubule stability. BER1 belongs to an evolutionary conserved gene family whose founding member Sensitivity to Red light Reduced is involved in red-light perception and circadian rhythms in Arabidopsis. Here, we present data showing that the ber1Delta mutant is affected in microtubule stability, particularly in presence of microtubule-depolymerising drugs. The pattern of synthetic lethal interactions obtained with the ber1Delta mutant suggests that Ber1 may function in N-terminal protein acetylation. Our work thus suggests that microtubule stability might be regulated through this post-translational modification on yet-to-be determined proteins

Phylogenomics Controlling for Base Compositional Bias Reveals a Single Origin of Eusociality in Corbiculate Bees.

As increasingly large molecular data sets are collected for phylogenomics, the conflicting phylogenetic signal among gene trees poses challenges to resolve some difficult nodes of the Tree of Life. Among these nodes, the phylogenetic position of the honey bees (Apini) within the corbiculate bee group remains controversial, despite its considerable importance for understanding the emergence and maintenance of eusociality. Here, we show that this controversy stems in part from pervasive phylogenetic conflicts among GC-rich gene trees. GC-rich genes typically have a high nucleotidic heterogeneity among species, which can induce topological conflicts among gene trees. When retaining only the most GC-homogeneous genes or using a nonhomogeneous model of sequence evolution, our analyses reveal a monophyletic group of the three lineages with a eusocial lifestyle (honey bees, bumble bees, and stingless bees). These phylogenetic relationships strongly suggest a single origin of eusociality in the corbiculate bees, with no reversal to solitary living in this group. To accurately reconstruct other important evolutionary steps across the Tree of Life, we suggest removing GC-rich and GC-heterogeneous genes from large phylogenomic data sets. Interpreted as a consequence of genome-wide variations in recombination rates, this GC effect can affect all taxa featuring GC-biased gene conversion, which is common in eukaryotes.

Suppression of long-branch attraction artefacts in the animal phylogeny using a site-heterogeneous model

Affiliation: Département de Biochimie, Université de Montréal