Leonardo de Oliveira Martins

Evolutionary Process of Deep-Sea Bathymodiolus Mussels

Background Since the discovery of deep-sea chemosynthesis-based communities, much work has been done to clarify their organismal and environmental aspects. However, major topics remain to be resolved, including when and how organisms invade and adapt to deep-sea environments; whether strategies for invasion and adaptation are shared by different taxa or unique to each taxon; how organisms extend their distribution and diversity; and how they become isolated to speciate in continuous waters. Deep-sea mussels are one of the dominant organisms in chemosynthesis-based communities, thus investigations of their origin and evolution contribute to resolving questions about life in those communities. Methodology/Principal Finding We investigated worldwide phylogenetic relationships of deep-sea Bathymodiolus mussels and their mytilid relatives by analyzing nucleotide sequences of the mitochondrial cytochrome c oxidase subunit I (COI) and NADH dehydrogenase subunit 4 (ND4) genes. Phylogenetic analysis of the concatenated sequence data showed that mussels of the subfamily Bathymodiolinae from vents and seeps were divided into four groups, and that mussels of the subfamily Modiolinae from sunken wood and whale carcasses assumed the outgroup position and shallow-water modioline mussels were positioned more distantly to the bathymodioline mussels. We provisionally hypothesized the evolutionary history of Bathymodilolus mussels by estimating evolutionary time under a relaxed molecular clock model. Diversification of bathymodioline mussels was initiated in the early Miocene, and subsequently diversification of the groups occurred in the early to middle Miocene. Conclusions/Significance The phylogenetic relationships support the “Evolutionary stepping stone hypothesis,” in which mytilid ancestors exploited sunken wood and whale carcasses in their progressive adaptation to deep-sea environments. This hypothesis is also supported by the evolutionary transition of symbiosis in that nutritional adaptation to the deep sea proceeded from extracellular to intracellular symbiotic states in whale carcasses. The estimated evolutionary time suggests that the mytilid ancestors were able to exploit whales during adaptation to the deep sea.

SimPhy: Phylogenomic Simulation of Gene, Locus, and Species Trees.

We present a fast and flexible software package—SimPhy—for the simulation of multiple gene families evolving under incomplete lineage sorting, gene duplication and loss, horizontal gene transfer—all three potentially leading to species tree/gene tree discordance—and gene conversion. SimPhy implements a hierarchical phylogenetic model in which the evolution of species, locus, and gene trees is governed by global and local parameters (e.g., genome-wide, species-specific, locus-specific), that can be fixed or be sampled from a priori statistical distributions. SimPhy also incorporates comprehensive models of substitution rate variation among lineages (uncorrelated relaxed clocks) and the capability of simulating partitioned nucleotide, codon, and protein multilocus sequence alignments under a plethora of substitution models using the program INDELible. We validate SimPhys output using theoretical expectations and other programs, and show that it scales extremely well with complex models and/or large trees, being an order of magnitude faster than the most similar program (DLCoal-Sim). In addition, we demonstrate how SimPhy can be useful to understand interactions among different evolutionary processes, conducting a simulation study to characterize the systematic overestimation of the duplication time when using standard reconciliation methods. SimPhy is available at https://github.com/adamallo/SimPhy, where users can find the source code, precompiled executables, a detailed manual and example cases.

Phylogenetic Detection of Recombination with a Bayesian Prior on the Distance between Trees

Genomic regions participating in recombination events may support distinct topologies, and phylogenetic analyses should incorporate this heterogeneity. Existing phylogenetic methods for recombination detection are challenged by the enormous number of possible topologies, even for a moderate number of taxa. If, however, the detection analysis is conducted independently between each putative recombinant sequence and a set of reference parentals, potential recombinations between the recombinants are neglected. In this context, a recombination hotspot can be inferred in phylogenetic analyses if we observe several consecutive breakpoints. We developed a distance measure between unrooted topologies that closely resembles the number of recombinations. By introducing a prior distribution on these recombination distances, a Bayesian hierarchical model was devised to detect phylogenetic inconsistencies occurring due to recombinations. This model relaxes the assumption of known parental sequences, still common in HIV analysis, allowing the entire dataset to be analyzed at once. On simulated datasets with up to 16 taxa, our method correctly detected recombination breakpoints and the number of recombination events for each breakpoint. The procedure is robust to rate and transition∶transversion heterogeneities for simulations with and without recombination. This recombination distance is related to recombination hotspots. Applying this procedure to a genomic HIV-1 dataset, we found evidence for hotspots and de novo recombination.

Phylogeography of Loaches of the Genus Lefua (Balitoridae, Cypriniformes) Inferred from Mitochondrial DNA Sequences

Abstract In order to elucidate phylogenetic relationships and intraspecific variations and to infer the evolutionary process of loaches of the genus Lefua, we analyzed nucleotide sequences of the mitochondrial D-loop region of 100 specimens obtained from 97 localities in Japan and Korea. The genus Lefua includes three described species, L. nikkonis, L. echigonia, and L. costata and an undescribed species, Lefua sp. Our results showed that each species of Lefua formed a monophyletic group, indicating clearly that Lefua species can be genetically distinguished from one another. Lefua nikkonis was the most closely related to L. costata, while L. sp. was the most closely related to L. echigonia. Specimens of L. sp. were grouped into two intraspecific populations and specimens of L. echigonia were grouped into six populations. These populations were well separated geographically from one another by mountain ranges and highlands. We estimated the evolutionary time for splitting of the species and intraspecific populations, and speculated on the evolutionary process of the genus Lefua. Species of Lefua are severely threatened. Fundamental genetic information is indispensable for conservation. We presented genetic background in order to protect these threatened loaches.

A Bayesian Supertree Model for Genome-Wide Species Tree Reconstruction

Current phylogenomic data sets highlight the need for species tree methods able to deal with several sources of gene tree/species tree incongruence. At the same time, we need to make most use of all available data. Most species tree methods deal with single processes of phylogenetic discordance, namely, gene duplication and loss, incomplete lineage sorting (ILS) or horizontal gene transfer. In this manuscript, we address the problem of species tree inference from multilocus, genome-wide data sets regardless of the presence of gene duplication and loss and ILS therefore without the need to identify orthologs or to use a single individual per species. We do this by extending the idea of Maximum Likelihood (ML) supertrees to a hierarchical Bayesian model where several sources of gene tree/species tree disagreement can be accounted for in a modular manner. We implemented this model in a computer program called guenomu whose inputs are posterior distributions of unrooted gene tree topologies for multiple gene families, and whose output is the posterior distribution of rooted species tree topologies. We conducted extensive simulations to evaluate the performance of our approach in comparison with other species tree approaches able to deal with more than one leaf from the same species. Our method ranked best under simulated data sets, in spite of ignoring branch lengths, and performed well on empirical data, as well as being fast enough to analyze relatively large data sets. Our Bayesian supertree method was also very successful in obtaining better estimates of gene trees, by reducing the uncertainty in their distributions. In addition, our results show that under complex simulation scenarios, gene tree parsimony is also a competitive approach once we consider its speed, in contrast to more sophisticated models.

Lineage-specific positive selection at the merozoite surface protein 1 (msp1) locus of Plasmodium vivax and related simian malaria parasites

BackgroundThe 200 kDa merozoite surface protein 1 (MSP-1) of malaria parasites, a strong vaccine candidate, plays a key role during erythrocyte invasion and is a target of host protective immune response. Plasmodium vivax, the most widespread human malaria parasite, is closely related to parasites that infect Asian Old World monkeys, and has been considered to have become a parasite of man by host switch from a macaque malaria parasite. Several Asian monkey parasites have a range of natural hosts. The same parasite species shows different disease manifestations among host species. This suggests that host immune responses to P. vivax-related malaria parasites greatly differ among host species (albeit other factors). It is thus tempting to invoke that a major immune target parasite protein such as MSP-1 underwent unique evolution, depending on parasite species that exhibit difference in host range and host specificity.ResultsWe performed comparative phylogenetic and population genetic analyses of the gene encoding MSP-1 (msp1) from P. vivax and nine P. vivax-related simian malaria parasites. The inferred phylogenetic tree of msp1 significantly differed from that of the mitochondrial genome, with a striking displacement of P. vivax from a position close to P. cynomolgi in the mitochondrial genome tree to an outlier of Asian monkey parasites. Importantly, positive selection was inferred for two ancestral branches, one leading to P. inui and P. hylobati and the other leading to P. vivax, P. fieldi and P. cynomolgi. This ancestral positive selection was estimated to have occurred three to six million years ago, coinciding with the period of radiation of Asian macaques. Comparisons of msp1 polymorphisms between P. vivax, P. inui and P. cynomolgi revealed that while some positively selected amino acid sites or regions are shared by these parasites, amino acid changes greatly differ, suggesting that diversifying selection is acting species-specifically on msp1.ConclusionsThe present results indicate that the msp1 locus of P. vivax and related parasite species has lineage-specific unique evolutionary history with positive selection. P. vivax and related simian malaria parasites offer an interesting system toward understanding host species-dependent adaptive evolution of immune-target surface antigen genes such as msp1.

Evolutionary Dynamics of HIV-1 BF and CB Recombinants and Its Parental Counterparts in South America

Analysis of the near full-length genomes revealed that the subtype F appeared in Brazilian HIV-1 epidemics at least 10 years after the subtype B. Notably, the BF recombinant emerged almost simultaneously with the introduction of subtype F in Brazil. Analysis of reverse transcriptase fragments indicated that the C subtype originated in the early 1990s, and the CB recombinant emerged 2 years after the appearance of subtype C. The high growth rate of BF recombinant possibly obscured the prevalence of the pure subtype F. In contrast, subtype C, although appearing 20 years after subtype B, was responsible for a well-defi ned epidemic. Nevertheless, the CB recombinant equally emerged rapidly after the introduction of the second parental (subtype C). Our results suggest that the outcome related to the recombinant profi le are probably infl uenced by the capacity of the newly arriving subtype to establish a critical number of infections before it recombines with the previous circulating subtype.

Codon pairs of the HIV-1 vif gene correlate with CD4+ T cell count.

BackgroundThe human APOBEC3G (A3G) protein activity is associated with innate immunity against HIV-1 by inducing high rates of guanosines to adenosines (G-to-A) mutations (viz., hypermutation) in the viral DNA. If hypermutation is not enough to disrupt the reading frames of viral genes, it may likely increase the HIV-1 diversity. To counteract host innate immunity HIV-1 encodes the Vif protein that binds A3G protein and form complexes to be degraded by cellular proteolysis.MethodsHere we studied the pattern of substitutions in the vif gene and its association with clinical status of HIV-1 infected individuals. To perform the study, unique vif gene sequences were generated from 400 antiretroviral-naïve individuals.ResultsThe codon pairs: 78–154, 85–154, 101–157, 105–157, and 105–176 of vif gene were associated with CD4+ T cell count lower than 500 cells per mm3. Some of these codons were located in the 81LGQGVSIEW89 region and within the BC-Box. We also identified codons under positive selection clustered in the N-terminal region of Vif protein, between 21WKSLVK26 and 40YRHHY44 regions (i.e., 31, 33, 37, 39), within the BC-Box (i.e., 155, 159) and the Cullin5-Box (i.e., 168) of vif gene. All these regions are involved in the Vif-induced degradation of A3G/F complexes and the N-terminal of Vif protein binds to viral and cellular RNA.ConclusionsAdaptive evolution of vif gene was mostly to optimize viral RNA binding and A3G/F recognition. Additionally, since there is not a fully resolved structure of the Vif protein, codon pairs associated with CD4+ T cell count may elucidate key regions that interact with host cell factors. Here we identified and discriminated codons under positive selection and codons under functional constraint in the vif gene of HIV-1.

Unsorted Homology within Locus and Species Trees

The concept of homology lies at the root of evolutionary biology. Since the seminal work of Fitch (1970), three main categories of homology relationships have been defined at the molecular level: orthology, paralogy, and xenology. In brief, if two gene copies arose by duplication they are paralogs, whereas if they arose through speciation they are orthologs. If one of them was transferred from a contemporaneous species, we call them xenologs (Supplementary Fig. S1 in Supplementary Material online, available at http://dx.doi.org/10.5061/dryad.87k57; see Gray and Fitch (1983); Fitch (2000)). Indeed, these terms were coined under a phylogenetic framework in which species were represented by single individuals, and as such they have remained very much intact during the last four decades—although particular cases within these categories have received specific names (Mindell and Meyer 2001). However, advances in sequencing technology have changed the field, and it is now very common to collect data sets containing multiple gene loci and/or multiple individuals per species. In general, such genome-wide data sets not only have unveiled extensive phylogenomic incongruence (Jeffroy et al. 2006; Salichos and Rokas 2013) but have brought back to the spotlight the consideration of how ancestral polymorphisms sort within populations (Edwards 2009). Altogether, phylogenomic data make imperative the explicit distinction between organismal and gene histories. Let us consider phylogenetic relationships at three different levels: species, loci, and gene copies (Fig. 1). The distinction between species/population trees and gene trees has been known for decades (Goodman et al. 1979; Pamilo and Nei 1988; Takahata 1989), whereas the introduction of locus trees into these models is very recent (Rasmussen and Kellis 2012). In brief, a species tree depicts the evolutionary history of the sampled organisms. In this case, the nodes represent speciation events, connected by branches that reflect the population history along these periods, and where their widths represent effective population size (Ne) and their lengths represent time (usually in years or number of generations). Apart from speciations, only evolutionary processes that affect species as a whole are represented at this level, like hybridization. Note that species trees are equivalent to population trees when the organismal units of interest are conspecific populations. In this case, the nodes of the population trees represent isolation events. In general, we will refer to “species” as any diverging, interbreeding group of individuals regardless of its taxonomic rank. On the other hand, a locus tree represents the evolutionary history of the sampled loci for a given gene family (see Rasmussen and Kellis 2012). Since the loci exist inside individuals evolving as part of a population, the locus tree is embedded within the species tree. In a locus tree, the nodes depict either genetic divergence due to speciation in the embedding species tree or locus-level events such as duplication, losses, or horizontal gene transfers, whereas the branch lengths and widths represent time and Ne, respectively. Here, we assume that the locuslevel events get immediately fixed in the population, so these Ne are equivalent to those in the species tree and are the same for every locus. Finally, a gene tree represents the evolutionary history of the sampled gene copies that evolve inside the locus tree. Gene tree nodes indicate coalescent events, which looking forward in time correspond to the process of DNA replication and divergence, and that can occur around the speciation time, well before (deep coalescence) or afterwards (migration in population trees). The branches of the gene tree usually represent amount of substitutions per site, and can also represent number of generations or other measures of time. Importantly, these three historical layers do not necessarily coincide. True species/population trees can differ from true locus trees due to gene duplications, losses, and/or horizontal gene transfers, whereas true gene trees can differ from their embedding locus and species trees if there is incomplete lineage sorting (ILS) (Maddison 1997; Page and Charleston 1997) (and migration in the case of population trees). In this regard,

A Likelihood-based Index of Protein–Protein Binding Affinities with Application to Influenza HA Escape from Antibodies

Abstract In many biological systems, proteins interact with other organic molecules to produce indispensable functions, in which molecular recognition phenomena are essential. Proteins have kept or gained their functions during molecular evolution. Their functions seem to be flexible, and a few amino acid substitutions sometimes cause drastic changes in function. In order to monitor and predict such drastic changes in the early stages in target populations, we need to identify patterns of structural changes during molecular evolution causing decreases or increases in the binding affinity of protein complexes. In previous work, we developed a likelihood-based index to quantify the degree to which a sequence fits a given structure. This index was named the sequence-structure fitness (SSF) and is calculated empirically based on amino acid preferences and pairwise interactions in the structural environment present in template structures. In the present work, we used the SSF to develop an index to measure the binding affinity of protein–protein complexes defined as the log likelihood ratio, contrasting the fitness of the sequences to the structure of the complex and that of the uncomplexed proteins. We applied the developed index to the complexes formed between influenza A hemagglutinin (HA) and four antibodies. The antibody–antigen binding region of HA is under strong selection pressure by the host immune system. Hence, examination of the long-term adaptation of HA to the four antibodies could reveal the strategy of the molecular evolution of HA. Two antibodies cover the HA receptor-binding region, while the other two bind away from the receptor-binding region. By focusing on branches with a significant decline in binding ability, we could detect key amino acid replacements and investigate the mechanism via conditional probabilities. The contrast between the adaptations to the two types of antibodies suggests that the virus adapts to the immune system at the cost of structural change.