Protein-Protein Interactions in Clathrin Vesicular Assembly: Radial Distribution of Evolutionary Constraints in Interfaces

In eukaryotic organisms clathrin-coated vesicles are instrumental in the processes of endocytosis as well as intracellular protein trafficking. Hence, it is important to understand how these vesicles have evolved across eukaryotes, to carry cargo molecules of varied shapes and sizes. The intricate nature and functional diversity of the vesicles are maintained by numerous interacting protein partners of the vesicle system. However, to delineate functionally important residues participating in protein-protein interactions of the assembly is a daunting task as there are no high-resolution structures of the intact assembly available. The two cryoEM structures closely representing intact assembly were determined at very low resolution and provide positions of C alpha atoms alone. In the present study, using the method developed by us earlier, we predict the protein-protein interface residues in clathrin assembly, taking guidance from the available low-resolution structures. The conservation status of these interfaces when investigated across eukaryotes, revealed a radial distribution of evolutionary constraints, i.e., if the members of the clathrin vesicular assembly can be imagined to be arranged in spherical manner, the cargo being at the center and clathrins being at the periphery, the detailed phylogenetic analysis of these members of the assembly indicated high-residue variation in the members of the assembly closer to the cargo while high conservation was noted in clathrins and in other proteins at the periphery of the vesicle. This points to the strategy adopted by the nature to package diverse proteins but transport them through a highly conserved mechanism.

Differential Evolutionary Constraints in the Evolution of Chemoreceptors: A Murine and Human Case Study

Chemoreception is among the most important sensory modalities in animals. Organisms use the ability to perceive chemical compounds in all major ecological activities. Recent studies have allowed the characterization of chemoreceptor gene families. These genes present strikingly high variability in copy numbers and pseudogenization degrees among different species, but the mechanisms underlying their evolution are not fully understood. We have analyzed the functional networks of these genes, their orthologs distribution, and performed phylogenetic analyses in order to investigate their evolutionary dynamics. We have modeled the chemosensory networks and compared the evolutionary constraints of their genes in Mus musculus, Homo sapiens, and Rattus norvegicus. We have observed significant differences regarding the constraints on the orthologous groups and network topologies of chemoreceptors and signal transduction machinery. Our findings suggest that chemosensory receptor genes are less constrained than their signal transducing machinery, resulting in greater receptor diversity and conservation of information processing pathways. More importantly, we have observed significant differences among the receptors themselves, suggesting that olfactory and bitter taste receptors are more conserved than vomeronasal receptors.

Modularity and evolutionary constraints in a baculovirus gene regulatory network

Abstract Background The structure of regulatory networks remains an open question in our understanding of complex biological systems. Interactions during complete viral life cycles present unique opportunities to understand how host-parasite network take shape and behave. The Anticarsia gemmatalis multiple nucleopolyhedrovirus (AgMNPV) is a large double-stranded DNA virus, whose genome may encode for 152 open reading frames (ORFs). Here we present the analysis of the ordered cascade of the AgMNPV gene expression. Results We observed an earlier onset of the expression than previously reported for other baculoviruses, especially for genes involved in DNA replication. Most ORFs were expressed at higher levels in a more permissive host cell line. Genes with more than one copy in the genome had distinct expression profiles, which could indicate the acquisition of new functionalities. The transcription gene regulatory network (GRN) for 149 ORFs had a modular topology comprising five communities of highly interconnected nodes that separated key genes that are functionally related on different communities, possibly maximizing redundancy and GRN robustness by compartmentalization of important functions. Core conserved functions showed expression synchronicity, distinct GRN features and significantly less genetic diversity, consistent with evolutionary constraints imposed in key elements of biological systems. This reduced genetic diversity also had a positive correlation with the importance of the gene in our estimated GRN, supporting a relationship between phylogenetic data of baculovirus genes and network features inferred from expression data. We also observed that gene arrangement in overlapping transcripts was conserved among related baculoviruses, suggesting a principle of genome organization. Conclusions Albeit with a reduced number of nodes (149), the AgMNPV GRN had a topology and key characteristics similar to those observed in complex cellular organisms, which indicates that modularity may be a general feature of biological gene regulatory networks.

Inference of evolutionary forces acting on human biological pathways

Because natural selection is likely to act on multiple genes underlying a given phenotypic trait, we study here the potential effect of ongoing and past selection on the genetic diversity of human biological pathways. We first show that genes included in gene sets are generally under stronger selective constraints than other genes and that their evolutionary response is correlated. We then introduce a new procedure to detect selection at the pathway level based on a decomposition of the classical McDonald–Kreitman test extended to multiple genes. This new test, called 2DNS, detects outlier gene sets and takes into account past demographic effects and evolutionary constraints specific to gene sets. Selective forces acting on gene sets can be easily identified by a mere visual inspection of the position of the gene sets relative to their two-dimensional null distribution. We thus find several outlier gene sets that show signals of positive, balancing, or purifying selection but also others showing an ancient relaxation of selective constraints. The principle of the 2DNS test can also be applied to other genomic contrasts. For instance, the comparison of patterns of polymorphisms private to African and non-African populations reveals that most pathways show a higher proportion of nonsynonymous mutations in non-Africans than in Africans, potentially due to different demographic histories and selective pressures.

Evolutionary analyses of hedgehog and Hoxd-10 genes in fish species closely related to the zebrafish

The study of development has relied primarily on the isolation of mutations in genes with specific functions in development and on the comparison of their expression patterns in normal and mutant phenotypes. Comparative evolutionary analyses can complement these approaches. Phylogenetic analyses of Sonic hedgehog (Shh) and Hoxd-10 genes from 18 cyprinid fish species closely related to the zebrafish provide novel insights into the functional constraints acting on Shh. Our results confirm and extend those gained from expression and crystalline structure analyses of this gene. Unexpectedly, exon 1 of Shh is found to be almost invariant even in third codon positions among these morphologically divergent species suggesting that this exon encodes for a functionally important domain of the hedgehog protein. This is surprising because the main functional domain of Shh had been thought to be that encoded by exon 2. Comparisons of Shh and Hoxd-10 gene sequences and of resulting gene trees document higher evolutionary constraints on the former than on the latter. This might be indicative of more general evolutionary patterns in networks of developmental regulatory genes interacting in a hierarchical fashion. The presence of four members of the hedgehog gene family in cyprinid fishes was documented and their homologies to known hedgehog genes in other vertebrates were established.

Phylogeny versus body size as determinants of food web structure

Food webs are the complex networks of trophic interactions that stoke the metabolic fires of life. To understand what structures these interactions in natural communities, ecologists have developed simple models to capture their main architectural features. However, apparently realistic food webs can be generated by models invoking either predator-prey body-size hierarchies or evolutionary constraints as structuring mechanisms. As a result, this approach has not conclusively revealed which factors are the most important. Here we cut to the heart of this debate by directly comparing the influence of phylogeny and body size on food web architecture. Using data from 13 food webs compiled by direct observation, we confirm the importance of both factors. Nevertheless, phylogeny dominates in most networks. Moreover, path analysis reveals that the size-independent direct effect of phylogeny on trophic structure typically outweighs the indirect effect that could be captured by considering body size alone. Furthermore, the phylogenetic signal is asymmetric: closely related species overlap in their set of consumers far more than in their set of resources. This is at odds with several food web models, which take only the view-point of consumers when assigning interactions. The echo of evolutionary history clearly resonates through current food webs, with implications for our theoretical models and conservation priorities.

Convergent structure of multitrophic communities over three continents

Ecological theory predicts that communities using the same resources should have similar structure, but evolutionary constraints on colonization and niche shifts may hamper such convergence. Multitrophic communities of wasps exploiting fig fruits, which first evolved about 75MYA, do not show long-term “inheritance” of taxonomic (lineage) composition or species diversity. However, communities on three continents have converged ecologically in the presence and relative abundance of five insect guilds that we define. Some taxa fill the same niches in each community (phylogenetic niche conservatism). However, we show that overall convergence in ecological community structure depends also on a combination of niche shifts by resident lineages and local colonizations of figs by other insect lineages. Our study explores new ground, and develops new heuristic tools, in combining ecology and phylogeny to address patterns in the complex multitrophic communities of insect on plants, which comprise a large part of terrestrial biodiversity.

On the Evolution of Hexose Transporters in Kinetoplastid Potozoans

Glucose, an almost universally used energy and carbon source, is processed through several well-known metabolic pathways, primarily glycolysis. Glucose uptake is considered to be the first step in glycolysis. In kinetoplastids, a protozoan group that includes relevant human pathogens, the importance of glucose uptake in different phases of the life cycles is well established, and hexose transporters have been proposed as targets for therapeutic drugs. However, little is known about the evolutionary history of these hexose transporters. Hexose transporters contain an intracellular N- and C-termini, and 12 transmembrane spans connected by alternate intracellular and extracellular loops. In the present work we tested the hypothesis that the evolutionary rate of the transmembrane span is different from that of the whole sequence and that it is possible to define evolutionary units inside the sequence. The phylogeny of whole molecules was compared to that of their transmembrane spans and the loops connecting the transmembrane spans. We show that the evolutionary units in these proteins primarily consist of clustered rather than individual transmembrane spans. These analyses demonstrate that there are evolutionary constraints on the organization of these proteins; more specifically, the order of the transmembrane spans along the protein is highly conserved. Finally, we defined a signature sequence for the identification of kinetoplastid hexose transporters.

The evolution of antennal courtship in diplazontine parasitoid wasps (Hymenoptera, Ichneumonidae, Diplazontinae)

Background As predicted by theory, traits associated with reproduction often evolve at a comparatively high speed. This is especially the case for courtship behaviour which plays a central role in reproductive isolation. On the other hand, courtship behavioural traits often involve morphological and behavioural adaptations in both sexes; this suggests that their evolution might be under severe constraints, for instance irreversibility of character loss. Here, we use a recently proposed method to retrieve data on a peculiar courtship behavioural trait, i.e. antennal coiling, for 56 species of diplazontine parasitoid wasps. On the basis of a well-resolved phylogeny, we reconstruct the evolutionary history of antennal coiling and associated morphological modifications to study the mode of evolution of this complex character system. Results Our study reveals a large variation in shape, location and ultra-structure of male-specific modifications on the antennae. As for antennal coiling, we find either single-coiling, double-coiling or the absence of coiling; each state is present in multiple genera. Using a model comparison approach, we show that the possession of antennal modifications is highly correlated with antennal coiling behaviour. Ancestral state reconstruction shows that both antennal modifications and antennal coiling are highly congruent with the molecular phylogeny, implying low levels of homoplasy and a comparatively low speed of evolution. Antennal coiling is lost on two independent occasions, and never reacquired. A zero rate of regaining antennal coiling is supported by maximum parsimony, maximum likelihood and Bayesian approaches. Conclusions Our study provides the first comparative evidence for a tight correlation between male-specific antennal modifications and the use of the antennae during courtship. Antennal coiling in Diplazontinae evolved at a comparatively low rate, and was never reacquired in any of the studied taxa. This suggests that the loss of antennal coiling is irreversible on the timescale examined here, and therefore that evolutionary constraints have greatly influenced the evolution of antennal courtship in this group of parasitoid wasps. Further studies are needed to ascertain whether the loss of antennal coiling is irreversible on larger timescales, and whether evolutionary constraints have influenced courtship behavioural traits in a similar way in other groups.

Stable virulence levels in the HIV epidemic of Switzerland over two decades

OBJECTIVE: To determine whether the virulence of HIV-1 has been changing since its introduction into Switzerland. DESIGN: A prospective cohort study of HIV-1 infected individuals with well-characterized pre-therapy disease history. METHODS: To minimize the effect of recently imported viruses and ethnicity-associated host factors, the analysis was restricted to the white, north-west-European majority population of the cohort. Virulence was characterized by the decline slope of the CD4 cell count (n = 817 patients), the decline slope of the CD4:CD8 ratio (n = 815 patients) and the viral setpoint (n = 549 patients) in untreated patients with sufficient data points. Linear regression models were used to detect correlations between the date of diagnosis (ranging between 1984 and 2003) and the virulence markers, controlling for gender, exposure category, age and CD4 cell count at entry. RESULTS: We found no correlation between any of the virulence markers and the date of diagnosis. Inspection of short-term trends confirmed that virulence has fluctuated around a stable level over time. CONCLUSIONS: The lack of long-term time trends in the virulence markers indicates that HIV-1 is not evolving towards increasing or decreasing virulence at a perceptible rate. Both highly virulent and attenuated strains have apparently been unable to spread at the population level. This result suggests that either the evolution of virulence may be slow or inhibited due to evolutionary constraints, or HIV-1 may have already evolved to optimal virulence in the human host.

IDENTIFICATION AND ANALYSIS OF A NOVEL ROLE FOR THE TOUSLED-LIKE KINASE IN REGULATING MITOTIC SPINDLE DYNAMICS

Deregulation of kinase activity is one example of how cells become cancerous by evading evolutionary constraints. The Tousled kinase (Tsl) was initially identified in Arabidopsis thaliana as a developmentally important kinase. There are two mammalian orthologues of Tsl and one orthologue in C. elegans, TLK-1, which is essential for embryonic viability and germ cell development. Depletion of TLK-1 leads to embryonic arrest large, distended nuclei, and ultimately embryonic lethality. Prior to terminal arrest, TLK-1-depleted embryos undergo aberrant mitoses characterized by poor metaphase chromosome alignment, delayed mitotic progression, lagging chromosomes, and supernumerary centrosomes. I discovered an unanticipated requirement for TLK-1 in mitotic spindle assembly and positioning. Normally, in the newly-fertilized zygote (P0) the maternal pronucleus migrates toward the paternal pronucleus at the posterior end of the embryo. After pronuclear meeting, the pronuclear-centrosome complex rotates 90° during centration to align on the anteroposterior axis followed by nuclear envelope breakdown (NEBD). However, in TLK-1-depleted P0 embryos, the centrosome-pronuclear complex rotation is significantly delayed with respect to NEBD and chromosome congression, Additionally, centrosome positions over time in tlk-1(RNAi) early embryos revealed a defect in posterior centrosome positioning during spindle-pronuclear centration, and 4D analysis of centrosome positions and movement in newly fertilized embryos showed aberrant centrosome dynamics in TLK-1-depleted embryos. Several mechanisms contribute to spindle rotation, one of which is the anchoring of astral microtubules to the cell cortex. Attachment of these microtubules to the cortices is thought to confer the necessary stability and forces in order to rotate the centrosome-pronuclear complex in a timely fashion. Analysis of a microtubule end-binding protein revealed that TLK-1-depleted embryos exhibit a more stochastic distribution of microtubule growth toward the cell cortices, and the types of microtubule attachments appear to differ from wild-type embryos. Additionally, fewer astral microtubules are in the vicinity of the cell cortex, thus suggesting that the delayed spindle rotation could be in part due to a lack of appropriate microtubule attachments to the cell cortex. Together with recently published biochemical data revealing the Tousled-like kinases associate with components of the dynein microtubule motor complex in humans, these data suggest that Tousled-like kinases play an important role in mitotic spindle assembly and positioning.

Mitochondrial carbonic anhydrase CA VB: Differences in tissue distribution and pattern of evolution from those of CA VA suggest distinct physiological roles

A cDNA for a second mouse mitochondrial carbonic anhydrase (CA) called CA VB was identified by homology to the previously characterized murine CA V, now called CA VA. The full-length cDNA encodes a 317-aa precursor that contains a 33-aa classical mitochondrial leader sequence. Comparison of products expressed from cDNAs for murine CA VB and CA VA in COS cells revealed that both expressed active CAs that localized in mitochondria, and showed comparable activities in crude extracts and in mitochondria isolated from transfected COS cells. Northern blot analyses of total RNAs from mouse tissues and Western blot analyses of mouse tissue homogenates showed differences in tissue-specific expression between CA VB and CA VA. CA VB was readily detected in most tissues, while CA VA expression was limited to liver, skeletal muscle, and kidney. The human orthologue of murine CA VB was recently reported also. Comparison of the CA domain sequence of human CA VB with that reported here shows that the CA domains of CA VB are much more highly conserved between mouse and human (95% identity) than the CA domains of mouse and human CA VAs (78% identity). Analysis of phylogenetic relationships between these and other available human and mouse CA isozyme sequences revealed that mammalian CA VB evolved much more slowly than CA VA, accepting amino acid substitutions at least 4.5 times more slowly since each evolved from its respective human–mouse ancestral gene around 90 million years ago. Both the differences in tissue distribution and the much greater evolutionary constraints on CA VB sequences suggest that CA VB and CA VA have evolved to assume different physiological roles.

Computer simulations of enzyme catalysis: Finding out what has been optimized by evolution

The origin of the catalytic power of enzymes is discussed, paying attention to evolutionary constraints. It is pointed out that enzyme catalysis reflects energy contributions that cannot be determined uniquely by current experimental approaches without augmenting the analysis by computer simulation studies. The use of energy considerations and computer simulations allows one to exclude many of the popular proposals for the way enzymes work. It appears that the standard approaches used by organic chemists to catalyze reactions in solutions are not used by enzymes. This point is illustrated by considering the desolvation hypothesis and showing that it cannot account for a large increase in kcat relative to the corresponding kcage for the reference reaction in a solvent cage. The problems associated with other frequently invoked mechanisms also are outlined. Furthermore, it is pointed out that mutation studies are inconsistent with ground state destabilization mechanisms. After considering factors that were not optimized by evolution, we review computer simulation studies that reproduced the overall catalytic effect of different enzymes. These studies pointed toward electrostatic effects as the most important catalytic contributions. The nature of this electrostatic stabilization mechanism is far from being obvious because the electrostatic interaction between the reacting system and the surrounding area is similar in enzymes and in solution. However, the difference is that enzymes have a preorganized dipolar environment that does not have to pay the reorganization energy for stabilizing the relevant transition states. Apparently, the catalytic power of enzymes is stored in their folding energy in the form of the preorganized polar environment.

Integração morfológica e modularidade em crânios das espécies do grupo Rhinella granulosa

Os conceitos e métodos provindos das teorias de integração morfológica e de genética quantitativa formam o arcabouço teórico para o estudo da evolução de estruturas complexas, compostas de múltiplos caracteres que interagem entre si. Nesse trabalho, utilizamos o crânio como modelo de estrutura complexa e estudamos sua diversificação nas espécies de sapo do grupo Rhinella granulosa. As perguntas do trabalho foram: (1) A organização da (co)variação é similar entre as espécies?; (2) A organização da (co)variação é modular nas espécies, conforme expectativas baseadas em desenvolvimento ou função?; (3) Fatores externos, como filogenia e clima, estruturam a similaridade no padrão de covariação entre as espécies?; (4) A diversificação da morfologia média do crânio se deu por deriva ou seleção natural?; (5) A divergência na morfologia média do crânio está associada à variação climática entre as espécies?; e finalmente (6) Restrições evolutivas atuaram na divergência entre as espécies? Os espécimes foram escaneados e validamos o uso de imagens 3D para a mensuração de 21 distâncias lineares. Os crânios das espécies foram representados como matrizes fenotípicas (P) de covariância e de correlação entre as distâncias. A similaridade entre as P das espécies é alta. As P de todas as espécies se conformam a um padrão modular compatível com interações funcionais entre ossos. As diferenças entre as P concentram-se no rostro e são associadas a diferenças no clima entre as espécies. Detectamos sinal de seleção natural nos nós mais basais da filogenia e variação local no crânio está associada à variação na sazonalidade da chuva entre as espécies. Restrições evolutivas atuaram na diversificação do crânio das espécies, defletindo as respostas evolutivas para tamanho. Concluímos que tanto seleção estabilizadora e direcional, conectadas à variação climática, quanto restrições evolutivas atuaram na diversificação do crânio das espécies