Identification and purification of two distinct complexes containing the five RAD51 paralogs.

Cells defective in any of the RAD51 paralogs (RAD51B, RAD51C, RAD51D, XRCC2, and XRCC3) are sensitive to DNA cross-linking agents and to ionizing radiation. Because the paralogs are required for the assembly of DNA damage-induced RAD51 foci, and mutant cell lines are defective in homologous recombination and show genomic instability, their defect is thought to be caused by an inability to promote efficient recombinational repair. Here, we show that the five paralogs exist in two distinct complexes in human cells: one contains RAD51B, RAD51C, RAD51D, and XRCC2 (defined as BCDX2), whereas the other consists of RAD51C with XRCC3. Both protein complexes have been purified to homogeneity and their biochemical properties investigated. BCDX2 binds single-stranded DNA and single-stranded gaps in duplex DNA, in accord with the proposal that the paralogs play an early (pre-RAD51) role in recombinational repair. Moreover, BCDX2 complex binds specifically to nicks in duplex DNA. We suggest that the extreme sensitivity of paralog-defective cell lines to cross-linking agents is owing to defects in the processing of incised cross links and the consequential failure to initiate recombinational repair at these sites.

How confident can we be that orthologs are similar, but paralogs differ?

Homologous genes are classified into orthologs and paralogs, depending on whether they arose by speciation or duplication. It is widely assumed that orthologs share similar functions, whereas paralogs are expected to diverge more from each other. But does this assumption hold up on further examination? We present evidence that orthologs and paralogs are not so different in either their evolutionary rates or their mechanisms of divergence. We emphasize the importance of appropriately designed studies to test models of gene evolution between orthologs and between paralogs. Thus, functional change between orthologs might be as common as between paralogs, and future studies should be designed to test the impact of duplication against this alternative model.

Large-scale analysis of orthologs and paralogs under covarion-like and constant-but-different models of amino acid evolution.

Functional divergence between homologous proteins is expected to affect amino acid sequences in two main ways, which can be considered as proxies of biochemical divergence: a "covarion-like" pattern of correlated changes in evolutionary rates, and switches in conserved residues ("conserved but different"). Although these patterns have been used in case studies, a large-scale analysis is needed to estimate their frequency and distribution. We use a phylogenomic framework of animal genes to answer three questions: 1) What is the prevalence of such patterns? 2) Can we link such patterns at the amino acid level with selection inferred at the codon level? 3) Are patterns different between paralogs and orthologs? We find that covarion-like patterns are more frequently detected than "constant but different," but that only the latter are correlated with signal for positive selection. Finally, there is no obvious difference in patterns between orthologs and paralogs.

Resolving the ortholog conjecture: orthologs tend to be weakly, but significantly, more similar in function than paralogs.

The function of most proteins is not determined experimentally, but is extrapolated from homologs. According to the "ortholog conjecture", or standard model of phylogenomics, protein function changes rapidly after duplication, leading to paralogs with different functions, while orthologs retain the ancestral function. We report here that a comparison of experimentally supported functional annotations among homologs from 13 genomes mostly supports this model. We show that to analyze GO annotation effectively, several confounding factors need to be controlled: authorship bias, variation of GO term frequency among species, variation of background similarity among species pairs, and propagated annotation bias. After controlling for these biases, we observe that orthologs have generally more similar functional annotations than paralogs. This is especially strong for sub-cellular localization. We observe only a weak decrease in functional similarity with increasing sequence divergence. These findings hold over a large diversity of species; notably orthologs from model organisms such as E. coli, yeast or mouse have conserved function with human proteins.

Consequences of lineage-specific gene loss on functional evolution of surviving paralogs: ALDH1A and retinoic acid signaling in vertebrate genomes

Genome duplications increase genetic diversity and may facilitate the evolution of gene subfunctions. Little attention, however, has focused on the evolutionary impact of lineage-specific gene loss. Here, we show that identifying lineage-specific gene loss after genome duplication is important for understanding the evolution of gene subfunctions in surviving paralogs and for improving functional connectivity among human and model organism genomes. We examine the general principles of gene loss following duplication, coupled with expression analysis of the retinaldehyde dehydrogenase Aldh1a gene family during retinoic acid signaling in eye development as a case study. Humans have three ALDH1A genes, but teleosts have just one or two. We used comparative genomics and conserved syntenies to identify loss of ohnologs (paralogs derived from genome duplication) and to clarify uncertain phylogenies. Analysis showed that Aldh1a1 and Aldh1a2 form a clade that is sister to Aldh1a3-related genes. Genome comparisons showed secondarily loss of aldh1a1 in teleosts, revealing that Aldh1a1 is not a tetrapod innovation and that aldh1a3 was recently lost in medaka, making it the first known vertebrate with a single aldh1a gene. Interestingly, results revealed asymmetric distribution of surviving ohnologs between co-orthologous teleost chromosome segments, suggesting that local genome architecture can influence ohnolog survival. We propose a model that reconstructs the chromosomal history of the Aldh1a family in the ancestral vertebrate genome, coupled with the evolution of gene functions in surviving Aldh1a ohnologs after R1, R2, and R3 genome duplications. Results provide evidence for early subfunctionalization and late subfunction-partitioning and suggest a mechanistic model based on altered regulation leading to heterochronic gene expression to explain the acquisition or modification of subfunctions by surviving ohnologs that preserve unaltered ancestral developmental programs in the face of gene loss.

The discovery of Foxl2 paralogs in chondrichthyan, coelacanth and tetrapod genomes reveals an ancient duplication in vertebrates

The Foxl2 (forkhead box L2) gene is an important member of the forkhead domain family, primarily responsible for the development of ovaries during female sex differentiation. The evolutionary studies conducted previously considered the presence of paralog Foxl2 copies only in teleosts. However, to search for possible paralog copies in other groups of vertebrates and ensure that all predicted copies were homolog to the Foxl2 gene, a broad evolutionary analysis was performed, based on the forkhead domain family. A total of 2464 sequences for the forkhead domain were recovered, and subsequently, 64 representative sequences for Foxl2 were used in the evolutionary analysis of this gene. The most important contribution of this study was the discovery of a new subgroup of Foxl2 copies (ortholog to Foxl2B) present in the chondrichthyan Callorhinchus milii, in the coelacanth Latimeria chalumnae, in the avian Taeniopygia guttata and in the marsupial Monodelphis domestica. This new scenario indicates a gene duplication event in an ancestor of gnathostomes. Furthermore, based on the analysis of the syntenic regions of both Foxl2 copies, the duplication event was not exclusive to Foxl2. Moreover, the duplicated copy distribution was shown to be complex across vertebrates, especially in tetrapods, and the results strongly support a loss of this copy in eutherian species. Finally, the scenario observed in this study suggests an update for Foxl2 gene nomenclature, extending the actual suggested teleost naming of Foxl2A and Foxl2B to all vertebrate sequences and contributing to the establishment of a new evolutionary context for the Foxl2 gene. © 2013 Macmillan Publishers Limited All rights reserved.

Ki-1/57 And Cgi-55 Ectopic Expression Impact Cellular Pathways Involved In Proliferation And Stress Response Regulation.

Ki-1/57 (HABP4) and CGI-55 (SERBP1) are regulatory proteins and paralogs with 40.7% amino acid sequence identity and 67.4% similarity. Functionally, they have been implicated in the regulation of gene expression on both the transcriptional and mRNA metabolism levels. A link with tumorigenesis is suggested, since both paralogs show altered expression levels in tumor cells and the Ki-1/57 gene is found in a region of chromosome 9q that represents a haplotype for familiar colon cancer. However, the target genes regulated by Ki-1/57 and CGI-55 are unknown. Here, we analyzed the alterations of the global transcriptome profile after Ki-1/57 or CGI-55 overexpression in HEK293T cells by DNA microchip technology. We were able to identify 363 or 190 down-regulated and 50 or 27 up-regulated genes for Ki-1/57 and CGI-55, respectively, of which 20 were shared between both proteins. Expression levels of selected genes were confirmed by qRT-PCR both after protein overexpression and siRNA knockdown. The majority of the genes with altered expression were associated to proliferation, apoptosis and cell cycle control processes, prompting us to further explore these contexts experimentally. We observed that overexpression of Ki-1/57 or CGI-55 results in reduced cell proliferation, mainly due to a G1 phase arrest, whereas siRNA knockdown of CGI-55 caused an increase in proliferation. In the case of Ki-1/57 overexpression, we found protection from apoptosis after treatment with the ER-stress inducer thapsigargin. Together, our data give important new insights that may help to explain these proteins putative involvement in tumorigenic events.

LGI1 mutations in temporal lobe epilepsies

Background and Objectives: A number of familial temporal lobe epilepsies (TLE) have been recently recognized. Mutations in LGI1 (leucine-rich, glioma-inactivated 1 gene) have been found in a few families with the syndrome of autosomal dominant partial epilepsy with auditory features (ADPEAF). The authors aimed to determine the spectrum of TLE phenotypes with LGI1 mutations, to study the frequency of mutations in ADPEAF, and to examine the role of LGI1 paralogs in ADPEAF without LGI1 mutations. Methods: The authors performed a clinical and molecular analysis on 75 pedigrees comprising 54 with a variety of familial epilepsies associated with TLE and 21 sporadic TLE cases. All were studied for mutations in LGI1. ADPEAF families negative for LGI1 mutations were screened for mutations in LGI2, LGI3, and LGI4. Results: Four families had ADPEAF, 22 had mesial TLE, 11 had TLE with febrile seizures, two had TLE with developmental abnormalities, and 15 had various other TLE syndromes. LGI1 mutations were found in two of four ADPEAF families, but in none of the other 50 families nor in the 21 individuals with sporadic TLE. The mutations were novel missense mutations in exons 1 (c. 124T --> G; C42G) and 8 (c. 1418C --> T; S473L). No mutations in LGI2, LGI3, or LGI4 were found in the other two ADPEAF families. Conclusion: In TLE, mutations in LGI1 are specific for ADPEAF but do not occur in all families. ADPEAF is genetically heterogeneous, but mutations in LGI2, LGI3, or LGI4 did not account for families without LGI1 mutations.

Evolutionary History of Arabidopsis thaliana Aminoacyl-tRNA Synthetase Dual-Targeted Proteins

Aminoacyl-transfer RNA (tRNA) synthetases (aaRS) are key players in translation and act early in protein synthesis by mediating the attachment of amino acids to their cognate tRNA molecules. In plants, protein synthesis may occur in three subcellular compartments (cytosol, mitochondria, and chloroplasts), which requires multiple versions of the protein to be correctly delivered to its proper destination. The organellar aaRS are nuclear encoded and equipped with targeting information at the N-terminal sequence, which enables them to be specifically translocated to their final location. Most of the aaRS families present organellar proteins that are dual targeted to mitochondria and chloroplasts. Here, we examine the dual targeting behavior of aaRS from an evolutionary perspective. Our results show that Arabidopsis thaliana aaRS sequences are a result of a horizontal gene transfer event from bacteria. However, there is no evident bias indicating one single ancestor (Cyanobacteria or Proteobacteria). The dual-targeted aaRS phylogenetic relationship was characterized into two different categories (paralogs and homologs) depending on the state recovered for both dual-targeted and cytosolic proteins. Taken together, our results suggest that the dual-targeted condition is a gain-of-function derived from gene duplication. Selection may have maintained the original function in at least one of the copies as the additional copies diverged.

Drosophila melanogaster lipins are tissue-regulated and developmentally regulated and present specific subcellular distributions

Lipins constitute a novel family of Mg2+-dependent phosphatidate phosphatases that catalyze the dephosphorylation of phosphatidic acid to yield diacylglycerol, an important intermediate in lipid metabolism and cell signaling. Whereas a single lipin is detected in less complex organisms, in mammals there are distinct lipin isoforms and paralogs that are differentially expressed among tissues. Compatible with organism tissue complexity, we show that the single Drosophila Lpin1 ortholog (CG8709, here named DmLpin) expresses at least three isoforms (DmLpinA, DmLpinK and DmLpinJ) in a temporal and spatially regulated manner. The highest levels of lipin in the fat body, where DmLpinA and DmLpinK are expressed, correlate with the highest levels of triacylglycerol (TAG) measured in this tissue. DmLpinK is the most abundant isoform in the central nervous system, where TAG levels are significantly lower than in the fat body. In the testis, where TAG levels are even lower, DmLpinJ is the predominant isoform. Together, these data suggest that DmLpinA might be the isoform that is mainly involved in TAG production, and that DmLpinK and DmLpinJ could perform other cellular functions. In addition, we demonstrate by immunofluorescence that lipins are most strongly labeled in the perinuclear region of the fat body and ventral ganglion cells. In visceral muscles of the larval midgut and adult testis, lipins present a sarcomeric distribution. In the ovary chamber, the lipin signal is concentrated in the internal rim of the ring canal. These specific subcellular localizations of the Drosophila lipins provide the basis for future investigations on putative novel cellular functions of this protein family.

New insights into the long-term evolutionary history of the avian MHC

SummaryGene duplication and neofunctidnalization are important processes in the evolution of phenotypic complexity. They account for important evolutionary novelties that confer ecological adaptation, such as the major histocompatibility complex (MHC), a multigene family with a central role in vertebrates' adaptive immune system. Multigene families, which evolved in large part through duplication, represent promising systems to study the still strongly depbated relative roles of neutral and adaptive processes in the evolution of phenotypic complexity. Detailed knowledge on ecological function and a well-characterized evolutionary history place the mammals' MHC amongst ideal study systems. However mammalian MHCs usually encompass several million base pairs and hold a large number of functional and non-functional duplicate genes, which makes their study complex. Avian MHCs on the other hand are usually way more compact, but the reconstruction of. their evolutionary history has proven notoriously difficult. However, no focused attempt has been undertaken so far to study the avian MHC evolutionary history in a broad phylogenetic context and using adequate gene regions.In the present PhD, we were able to make important contributions to the understanding of the long-term evolution of the avian MHC class II Β (MHCI1B). First, we isolated and characterized MHCIIB genes in barn owl (Tyto alba?, Strigiformes, Tytonidae), a species from an avian lineage in which MHC has not been studied so far. Our results revealed that with only two functional MHCIIB genes the MHC organization of barn owl may be similar to the 'minimal essential' MHC of chicken (Gallus gallus), indicating that simple MHC organization may be ancestral to birds. Taking advantage of the sequence information from barn owl, we studied the evolution of MHCIIB genes in 13 additional species of 'typical' owls (Strigiformes, Strigidae). Phylogenetic analyses revealed that according to their function, in owls the peptide-binding region (PBR) encoding exon 2 and the non-PBR encoding exon 3 evolve by different patterns. Exon 2 exhibited an evolutionary history of positive selection and recombination, while exon 3 traced duplication history and revealed two paralogs evolving divergently from each other in owls, and in a shorebird, the great snipe {Gallinago media). The results from exon 3 were the first ever from birds to demonstrate gene orthology in species that diverged tens of millions of years ago, and strongly questioned whether the taxa studied before provided an adequate picture of avian MHC evolution. In a follow-up study, we aimed at explaining a striking pattern revealed by phylogenetic trees analyzing the owl sequences along with MHCIIB sequences from other birds: One owl paralog (termed DAB1) grouped with sequences of passerines and falcons, while the other (DAB2) grouped with wildfowl, penguins and birds of prey. This could be explained by either a duplication event preceding the evolution of these bird orders, or by convergent evolution of similar sequences in a number of orders. With extensive phylogenetic analyses we were able to show, that indeed a duplication event preceeded the major avian radiation -100 my ago, and that following this duplication, the paralogs evolved under positive selection. Furthermore, we showed that the divergently evolving amino acid residues in the MHCIIB-encoded β-chain potentially interact with the MHCI I α-chain, and that molecular coevolution of the interacting residues may have been involved in the divergent evolution of the MHCIIB paralogs.The findings of this PhD are of particular interest to the understanding of the evolutionary history of the avian MHC and, by providing essential information on long-term gene history in the avian MHC, open promising perspectives for advances in the understanding of the evolution of multigene families in general, and for avian MHC organization in particular. Amongst others I discuss the importance of including protein structure in the phylogenetic study of multigene families, and the roles of ecological versus molecular selection pressures. I conclude by providing a population genomic perspective on avian MHC, which may serve as a basis for future research to investigate the relative roles of neutral processes involving effective population size effects and of adaptation in the evolution of avian MHC diversity and organization.RésuméLa duplication de gènes et leur néo-fonctionnalisation sont des processus importants dans l'évolution de la complexité phénotypique. Ils sont impliqués dans l'apparition d'importantes nouveautés évolutives favorisant l'adaptation écologique, comme c'est le cas pour le complexe majeur d'histocompatibilité

Arbuscular mycorrhizal fungi (Glomeromycota) harbour ancient fungal tubulin genes that resemble those of the chytrids (Chytridiomycota).

The genes encoding alpha- and beta-tubulins have been widely sampled in most major fungal phyla and they are useful tools for fungal phylogeny. Here, we report the first isolation of alpha-tubulin sequences from arbuscular mycorrhizal fungi (AMF). In parallel, AMF beta-tubulins were sampled and analysed to identify the presence of paralogs of this gene. The AMF alpha-tubulin amino acid phylogeny was congruent with the results previously reported for AMF beta-tubulins and showed that AMF tubulins group together at a basal position in the fungal clade and showed high sequence similarities with members of the Chytridiomycota. This is in contrast with phylogenies for other regions of the AMF genome. The amount and nature of substitutions are consistent with an ancient divergence of both orthologs and paralogs of AMF tubulins. At the amino acid level, however, AMF tubulins have hardly evolved from those of the chytrids. This is remarkable given that these two groups are ancient and the monophyletic Glomeromycota probably diverged from basal fungal ancestors at least 500 million years ago. The specific primers we designed for the AMF tubulins, together with the high molecular variation we found among the AMF species we analysed, make AMF tubulin sequences potentially useful for AMF identification purposes.

SwissTargetPrediction: a web server for target prediction of bioactive small molecules.

Bioactive small molecules, such as drugs or metabolites, bind to proteins or other macro-molecular targets to modulate their activity, which in turn results in the observed phenotypic effects. For this reason, mapping the targets of bioactive small molecules is a key step toward unraveling the molecular mechanisms underlying their bioactivity and predicting potential side effects or cross-reactivity. Recently, large datasets of protein-small molecule interactions have become available, providing a unique source of information for the development of knowledge-based approaches to computationally identify new targets for uncharacterized molecules or secondary targets for known molecules. Here, we introduce SwissTargetPrediction, a web server to accurately predict the targets of bioactive molecules based on a combination of 2D and 3D similarity measures with known ligands. Predictions can be carried out in five different organisms, and mapping predictions by homology within and between different species is enabled for close paralogs and orthologs. SwissTargetPrediction is accessible free of charge and without login requirement at http://www.swisstargetprediction.ch.

Specific duplication and dorsoventrally asymmetric expression patterns of Cycloidea-like genes in zygomorphic species of Ranunculaceae.

Floral bilateral symmetry (zygomorphy) has evolved several times independently in angiosperms from radially symmetrical (actinomorphic) ancestral states. Homologs of the Antirrhinum majus Cycloidea gene (Cyc) have been shown to control floral symmetry in diverse groups in core eudicots. In the basal eudicot family Ranunculaceae, there is a single evolutionary transition from actinomorphy to zygomorphy in the stem lineage of the tribe Delphinieae. We characterized Cyc homologs in 18 genera of Ranunculaceae, including the four genera of Delphinieae, in a sampling that represents the floral morphological diversity of this tribe, and reconstructed the evolutionary history of this gene family in Ranunculaceae. Within each of the two RanaCyL (Ranunculaceae Cycloidea-like) lineages previously identified, an additional duplication possibly predating the emergence of the Delphinieae was found, resulting in up to four gene copies in zygomorphic species. Expression analyses indicate that the RanaCyL paralogs are expressed early in floral buds and that the duration of their expression varies between species and paralog class. At most one RanaCyL paralog was expressed during the late stages of floral development in the actinomorphic species studied whereas all paralogs from the zygomorphic species were expressed, composing a species-specific identity code for perianth organs. The contrasted asymmetric patterns of expression observed in the two zygomorphic species is discussed in relation to their distinct perianth architecture.

Gene duplication, population genomics, and species-level differentiation within a tropical mountain shrub.

Gene duplication leads to paralogy, which complicates the de novo assembly of genotyping-by-sequencing (GBS) data. The issue of paralogous genes is exacerbated in plants, because they are particularly prone to gene duplication events. Paralogs are normally filtered from GBS data before undertaking population genomics or phylogenetic analyses. However, gene duplication plays an important role in the functional diversification of genes and it can also lead to the formation of postzygotic barriers. Using populations and closely related species of a tropical mountain shrub, we examine 1) the genomic differentiation produced by putative orthologs, and 2) the distribution of recent gene duplication among lineages and geography. We find high differentiation among populations from isolated mountain peaks and species-level differentiation within what is morphologically described as a single species. The inferred distribution of paralogs among populations is congruent with taxonomy and shows that GBS could be used to examine recent gene duplication as a source of genomic differentiation of nonmodel species.