Comparative Evolutionary Approach Research Articles

Identification of potential auxin response candidate genes for soybean rapid canopy coverage through comparative evolution and expression analysis.

Throughout domestication, crop plants have gone through strong genetic bottlenecks, dramatically reducing the genetic diversity in today's available germplasm. This has also reduced the diversity in traits necessary for breeders to develop improved varieties. Many strategies have been developed to improve both genetic and trait diversity in crops, from backcrossing with wild relatives, to chemical/radiation mutagenesis, to genetic engineering. However, even with recent advances in genetic engineering we still face the rate limiting step of identifying which genes and mutations we should target to generate diversity in specific traits. Here, we apply a comparative evolutionary approach, pairing phylogenetic and expression analyses to identify potential candidate genes for diversifying soybean (Glycine max) canopy cover development via the nuclear auxin signaling gene families, while minimizing pleiotropic effects in other tissues. In soybean, rapid canopy cover development is correlated with yield and also suppresses weeds in organic cultivation. We identified genes most specifically expressed during early canopy development from the TIR1/AFB auxin receptor, Aux/IAA auxin co-receptor, and ARF auxin response factor gene families in soybean, using principal component analysis. We defined Arabidopsis thaliana and model legume species orthologs for each soybean gene in these families allowing us to speculate potential soybean phenotypes based on well-characterized mutants in these model species. In future work, we aim to connect genetic and functional diversity in these candidate genes with phenotypic diversity in planta allowing for improvements in soybean rapid canopy cover, yield, and weed suppression. Further development of this and similar algorithms for defining and quantifying tissue- and phenotype-specificity in gene expression may allow expansion of diversity in valuable phenotypes in important crops.

Evolution and Diversity of TGF-β Pathways are Linked with Novel Developmental and Behavioral Traits.

Transforming growth factor-β (TGF-β) signaling is essential for numerous biologic functions. It is a highly conserved pathway found in all metazoans including the nematode Caenorhabditis elegans, which has also been pivotal in identifying many components. Utilizing a comparative evolutionary approach, we explored TGF-β signaling in nine nematode species and revealed striking variability in TGF-β gene frequency across the lineage. Of the species analyzed, gene duplications in the DAF-7 pathway appear common with the greatest disparity observed in Pristionchus pacificus. Specifically, multiple paralogues of daf-3, daf-4 and daf-7 were detected. To investigate this additional diversity, we induced mutations in 22 TGF-β components and generated corresponding double, triple, and quadruple mutants revealing both conservation and diversification in function. Although the DBL-1 pathway regulating body morphology appears highly conserved, the DAF-7 pathway exhibits functional divergence, notably in some aspects of dauer formation. Furthermore, the formation of the phenotypically plastic mouth in P. pacificus is partially influenced through TGF-β with the strongest effect in Ppa-tag-68. This appears important for numerous processes in P. pacificus but has no known function in C. elegans. Finally, we observe behavioral differences in TGF-β mutants including in chemosensation and the establishment of the P. pacificus kin-recognition signal. Thus, TGF-β signaling in nematodes represents a stochastic genetic network capable of generating novel functions through the duplication and deletion of associated genes.

Genome-wide analyses reveal genes subject to positive selection in Toxoplasma gondii

Toxoplasma gondii is an important protozoan pathogen that infects many wild and domestic animals and causes infections in immunocompromised humans. However, there has been little investigation of the molecular evolutionary trajectories of this pathogenic protozoa using comparative genomics data. Here, we employed a comparative evolutionary genomics approach to identify genes that are under site- and lineage-specific positive selection in nine strains of T. gondii, including two closely related species, Neospora caninum and Hammondia hammondi. Based on the analyses of five coccidian core genomes, 4.5% of the 5788 core genome genes showed strong signals for positive selection in the site model. In addition, the branch-site model analyses in the nine T. gondii core genomes indicated that 2 to 20 genes underwent significant positive selection along each lineage leading to T. gondii strains. Many of the protein products encoded by the positively selected genes are secretory or surface proteins that have previously been implicated in host pathogenesis. The adaptive changes in these positively selected genes might be related to dynamic interactions between the host immune systems and might play a crucial role in the infection and pathogenic processes of T. gondii.

Challenges Facing the Study of the Evolutionary Origins of Human Right-Handedness and Language

A growing consensus favors the predominance of the human left hemisphere in manipulation and language (speech and gestures). However, the mechanisms underlying brain lateralization for noncommunication and communication functions are still unclear. Many studies emphasize the ambiguous nature of the relationship between the directions of brain lateralization for manipulation and for language. A comparative evolutionary approach investigating lateralization and communication mechanisms in our closest living relatives can improve our understanding of human hemispheric specialisation. We review theories concerning the evolutionary origins of human right-handedness and language and studies, highlighting the relevance of a comparative evolutionary approach. We discuss four methodological issues related to the study of handedness and communication: 1) categorisation of signaling based on criteria of intentionality; 2) use of a comprehensive multimodal and multifactorial approach; 3) investigation of laterality in intraspecific gestures; and 4) comparison of manual laterality between gestures and noncommunication actions. Deeper investigations of the multifaceted nature of intentionality in communication and the multidimensionality of manual laterality and communication in humans and other primates are needed. Finally, we make four recommendations to enhance our understanding of the evolutionary origins of human right-handedness and language: 1) implementing a comprehensive and integrated investigation of laterality and communication; 2) considering socioecologically relevant contexts; 3) investigating intraspecific aspects of communication in humans and other primates that vary in their degree of sociality; and 4) using appropriate, standardized methods of data collection and analysis.

Comparative evolutionary approach to pain perception in fishes

Comparative evolutionary approach to pain perception in fishes

Comparative sequence analyses of rhodopsin and RPE65 reveal patterns of selective constraint across hereditary retinal disease mutations.

Retinitis pigmentosa (RP) comprises several heritable diseases that involve photoreceptor, and ultimately retinal, degeneration. Currently, mutations in over 50 genes have known links to RP. Despite advances in clinical characterization, molecular characterization of RP remains challenging due to the heterogeneous nature of causal genes, mutations, and clinical phenotypes. In this study, we compiled large datasets of two important visual genes associated with RP: rhodopsin, which initiates the phototransduction cascade, and the retinoid isomerase RPE65, which regenerates the visual cycle. We used a comparative evolutionary approach to investigate the relationship between interspecific sequence variation and pathogenic mutations that lead to degenerative retinal disease. Using codon-based likelihood methods, we estimated evolutionary rates (d N/d S) across both genes in a phylogenetic context to investigate differences between pathogenic and nonpathogenic amino acid sites. In both genes, disease-associated sites showed significantly lower evolutionary rates compared to nondisease sites, and were more likely to occur in functionally critical areas of the proteins. The nature of the dataset (e.g., vertebrate or mammalian sequences), as well as selection of pathogenic sites, affected the differences observed between pathogenic and nonpathogenic sites. Our results illustrate that these methods can serve as an intermediate step in understanding protein structure and function in a clinical context, particularly in predicting the relative pathogenicity (i.e., functional impact) of point mutations and their downstream phenotypic effects. Extensions of this approach may also contribute to current methods for predicting the deleterious effects of candidate mutations and to the identification of protein regions under strong constraint where we expect pathogenic mutations to occur.

The evolution of bat nucleic acid-sensing Toll-like receptors.

We characterized the nucleic acid-sensing Toll-like receptors (TLR) of a New World bat species, the common vampire bat (Desmodus rotundus), and through a comparative molecular evolutionary approach searched for general adaptation patterns among the nucleic acid-sensing TLRs of eight different bats species belonging to three families (Pteropodidae, Vespertilionidae and Phyllostomidae). We found that the bat TLRs are evolving slowly and mostly under purifying selection and that the divergence pattern of such receptors is overall congruent with the species tree, consistent with the evolution of many other mammalian nuclear genes. However, the chiropteran TLRs exhibited unique mutations fixed in ligand-binding sites, some of which involved nonconservative amino acid changes and/or targets of positive selection. Such changes could potentially modify protein function and ligand-binding properties, as some changes were predicted to alter nucleic acid binding motifs in TLR 9. Moreover, evidence for episodic diversifying selection acting specifically upon the bat lineage and sublineages was detected. Thus, the long-term adaptation of chiropterans to a wide variety of environments and ecological niches with different pathogen profiles is likely to have shaped the evolution of the bat TLRs in an order-specific manner. The observed evolutionary patterns provide evidence for potential functional differences between bat and other mammalian TLRs in terms of resistance to specific pathogens or recognition of nucleic acids in general.

Gene body methylation is conserved between plant orthologs and is of evolutionary consequence

DNA methylation is a common feature of eukaryotic genomes and is especially common in noncoding regions of plants. Protein coding regions of plants are often methylated also, but the extent, function, and evolutionary consequences of gene body methylation remain unclear. Here we investigate gene body methylation using an explicit comparative evolutionary approach. We generated bisulfite sequencing data from two tissues of Brachypodium distachyon and compared genic methylation patterns to those of rice (Oryza sativa ssp. japonica). Gene body methylation was strongly conserved between orthologs of the two species and affected a biased subset of long, slowly evolving genes. Because gene body methylation is conserved over evolutionary time, it shapes important features of plant genome evolution, such as the bimodality of G+C content among grass genes. Our results superficially contradict previous observations of high cytosine methylation polymorphism within Arabidopsis thaliana genes, but reanalyses of these data are consistent with conservation of methylation within gene regions. Overall, our results indicate that the methylation level is a long-term property of individual genes and therefore of evolutionary consequence.

Gene Repertoire Evolution of Streptococcus pyogenes Inferred from Phylogenomic Analysis with Streptococcus canis and Streptococcus dysgalactiae

Streptococcus pyogenes, is an important human pathogen classified within the pyogenic group of streptococci, exclusively adapted to the human host. Our goal was to employ a comparative evolutionary approach to better understand the genomic events concomitant with S. pyogenes human adaptation. As part of ascertaining these events, we sequenced the genome of one of the potential sister species, the agricultural pathogen S. canis, and combined it in a comparative genomics reconciliation analysis with two other closely related species, Streptococcus dysgalactiae and Streptococcus equi, to determine the genes that were gained and lost during S. pyogenes evolution. Genome wide phylogenetic analyses involving 15 Streptococcus species provided convincing support for a clade of S. equi, S. pyogenes, S. dysgalactiae, and S. canis and suggested that the most likely S. pyogenes sister species was S. dysgalactiae. The reconciliation analysis identified 113 genes that were gained on the lineage leading to S. pyogenes. Almost half (46%) of these gained genes were phage associated and 14 showed significant matches to experimentally verified bacteria virulence factors. Subsequent to the origin of S. pyogenes, over half of the phage associated genes were involved in 90 different LGT events, mostly involving different strains of S. pyogenes, but with a high proportion involving the horse specific pathogen S. equi subsp. equi, with the directionality almost exclusively (86%) in the S. pyogenes to S. equi direction. Streptococcus agalactiae appears to have played an important role in the evolution of S. pyogenes with a high proportion of LGTs originating from this species. Overall the analysis suggests that S. pyogenes adaptation to the human host was achieved in part by (i) the integration of new virulence factors (e.g. speB, and the sal locus) and (ii) the construction of new regulation networks (e.g. rgg, and to some extent speB).

Dissecting Arabidopsis Gβ Signal Transduction on the Protein Surface

The heterotrimeric G-protein complex provides signal amplification and target specificity. The Arabidopsis (Arabidopsis thaliana) Gβ-subunit of this complex (AGB1) interacts with and modulates the activity of target cytoplasmic proteins. This specificity resides in the structure of the interface between AGB1 and its targets. Important surface residues of AGB1, which were deduced from a comparative evolutionary approach, were mutated to dissect AGB1-dependent physiological functions. Analysis of the capacity of these mutants to complement well-established phenotypes of Gβ-null mutants revealed AGB1 residues critical for specific AGB1-mediated biological processes, including growth architecture, pathogen resistance, stomata-mediated leaf-air gas exchange, and possibly photosynthesis. These findings provide promising new avenues to direct the finely tuned engineering of crop yield and traits.

Acireductone Dioxygenase 1 (ARD1) Is an Effector of the Heterotrimeric G Protein β Subunit in Arabidopsis

Heterotrimeric G protein complexes are conserved from plants to mammals, but the complexity of each system varies. Arabidopsis thaliana contains one Gα, one Gβ (AGB1), and at least three Gγ subunits, allowing it to form three versions of the heterotrimer. This plant model is ideal for genetic studies because mammalian systems contain hundreds of unique heterotrimers. The activation of these complexes promotes interactions between both the Gα subunit and the Gβγ dimer with enzymes and scaffolds to propagate signaling to the cytoplasm. However, although effectors of Gα and Gβ are known in mammals, no Gβ effectors were previously known in plants. Toward identifying AGB1 effectors, we genetically screened for dominant mutations that suppress Gβ-null mutant (agb1-2) phenotypes. We found that overexpression of acireductone dioxygenase 1 (ARD1) suppresses the 2-day-old etiolated phenotype of agb1-2. ARD1 is homologous to prokaryotic and eukaryotic ARD proteins; one function of ARDs is to operate in the methionine salvage pathway. We show here that ARD1 is an active metalloenzyme, and AGB1 and ARD1 both control embryonic hypocotyl length by modulating cell division; they also may contribute to the production of ethylene, a product of the methionine salvage pathway. ARD1 physically interacts with AGB1, and ARD enzymatic activity is stimulated by AGB1 in vitro. The binding interface on AGB1 was deduced using a comparative evolutionary approach and tested using recombinant AGB1 mutants. A possible mechanism for AGB1 activation of ARD1 activity was tested using directed mutations in a loop near the substrate-binding site.

Comparative analysis of the primate X-inactivation center region and reconstruction of the ancestral primate XIST locus

Here we provide a detailed comparative analysis across the candidate X-Inactivation Center (XIC) region and the XIST locus in the genomes of six primates and three mammalian outgroup species. Since lemurs and other strepsirrhine primates represent the sister lineage to all other primates, this analysis focuses on lemurs to reconstruct the ancestral primate sequences and to gain insight into the evolution of this region and the genes within it. This comparative evolutionary genomics approach reveals significant expansion in genomic size across the XIC region in higher primates, with minimal size alterations across the XIST locus itself. Reconstructed primate ancestral XIC sequences show that the most dramatic changes during the past 80 million years occurred between the ancestral primate and the lineage leading to Old World monkeys. In contrast, the XIST locus compared between human and the primate ancestor does not indicate any dramatic changes to exons or XIST-specific repeats; rather, evolution of this locus reflects small incremental changes in overall sequence identity and short repeat insertions. While this comparative analysis reinforces that the region around XIST has been subject to significant genomic change, even among primates, our data suggest that evolution of the XIST sequences themselves represents only small lineage-specific changes across the past 80 million years.

The evolutionary origins of cognitive control.

The question of domain-specific versus domain-general processing is an ongoing source of inquiry surrounding cognitive control. Using a comparative evolutionary approach, Stout (2010) proposed two components of cognitive control: coordinating hierarchical action plans and social cognition. This article reports additional molecular and experimental evidence supporting a domain-general attentional process coordinating hierarchical action plans, with the earliest such control processing originating in the capacity of dynamic foraging behaviors-predating the vertebrate-invertebrate divergence (c. 700 million years ago). Further discussion addresses evidence required for additional, domain-specific, cognitive control processes, noting that proposed social processes may simply provide emotionally valenced representational information to the above hierarchical process.

The cell biology of touch

The sense of touch detects forces that bombard the body’s surface. In metazoans, an assortment of morphologically and functionally distinct mechanosensory cell types are tuned to selectively respond to diverse mechanical stimuli, such as vibration, stretch, and pressure. A comparative evolutionary approach across mechanosensory cell types and genetically tractable species is beginning to uncover the cellular logic of touch reception.

Comparative Vertebrate Evolutionary Analyses of Type I Collagen: Potential of COL1a1 Gene Structure and Intron Variation for Common Bone-Related Diseases

Collagen type I alpha 1 (COL1a1), which encodes the primary subunit of type I collagen, the main structural and most abundant protein in vertebrates, harbors hundreds of mutations linked to human diseases like osteoporosis and osteogenesis imperfecta. Previous studies have attempted to predict the phenotypic severity associated with type I collagen mutations, yet an evolutionary analysis that compares historical and recent selective pressures, including across noncoding regions, has never been conducted. Here, we use a comparative genomic and species evolutionary analysis representing ∼450 My of vertebrate history to investigate functional constraints associated with both exons and introns of the >17-kb COL1a1 gene. We find that although the COL1a1 amino acid sequence is highly conserved, there are both spatial and temporal signatures of varying selective constraint across protein domains. Furthermore, sites of high evolutionary constraint significantly correlate with the location of disease-associated mutations, the latter of which also cluster with respect to specific severity classes typically categorized in clinical studies. Finally, we find that COL1a1 introns are significantly short in length with high GC content, patterns that are shared across highly diverged vertebrates, and which may be a signature of strong stabilizing selection for high COL1a1 gene expression. In conclusion, although previous studies focused on COL1a1 coding regions, the current results implicate introns as areas of high selective constraint and targets of bone-related phenotypic variation. From a broader perspective, our comparative evolutionary approach provides further resolution to models predicting mutations associated with bone-related function and disease severity.

Prediction of Protein–Protein Interfaces on G-Protein β Subunits Reveals a Novel Phospholipase C β2 Binding Domain

Gβ subunits from heterotrimeric G-proteins (guanine nucleotide-binding proteins) directly bind diverse proteins, including effectors and regulators, to modulate a wide array of signaling cascades. These numerous interactions constrained the evolution of the molecular surface of Gβ. Although mammals contain five Gβ genes comprising two classes (Gβ1-like and Gβ5-like), plants and fungi have a single ortholog, and organisms such as Caenorhabditis elegans and Drosophila melanogaster contain one copy from each class. A limited number of crystal structures of complexes containing Gβ subunits and complementary biochemical data highlight specific sites within Gβs needed for protein interactions. It is difficult to determine from these interaction sites what, if any, additional regions of the Gβ molecular surface comprise interaction interfaces essential to Gβ's role as a nexus in numerous signaling cascades. We used a comparative evolutionary approach to identify five known and eight previously unknown putative interfaces on the surface of Gβ. We show that one such novel interface occurs between Gβ and phospholipase C β2 (PLC-β2), a mammalian Gβ interacting protein. Substitutions of residues within this Gβ–PLC-β2 interface reduce the activation of PLC-β2 by Gβ1, confirming that our de novo comparative evolutionary approach predicts previously unknown Gβ–protein interfaces. Similarly, we hypothesize that the seven remaining untested novel regions contribute to putative interfaces for other Gβ interacting proteins. Finally, this comparative evolutionary approach is suitable for application to any protein involved in a significant number of protein–protein interactions.

Temporal pattern of loss/persistence of duplicate genes involved in signal transduction and metabolic pathways after teleost-specific genome duplication

BackgroundRecent genomic studies have revealed a teleost-specific third-round whole genome duplication (3R-WGD) event occurred in a common ancestor of teleost fishes. However, it is unclear how the genes duplicated in this event were lost or persisted during the diversification of teleosts, and therefore, how many of the duplicated genes contribute to the genetic differences among teleosts. This subject is also important for understanding the process of vertebrate evolution through WGD events. We applied a comparative evolutionary approach to this question by focusing on the genes involved in long-term potentiation, taste and olfactory transduction, and the tricarboxylic acid cycle, based on the whole genome sequences of four teleosts; zebrafish, medaka, stickleback, and green spotted puffer fish.ResultsWe applied a state-of-the-art method of maximum-likelihood phylogenetic inference and conserved synteny analyses to each of 130 genes involved in the above biological systems of human. These analyses identified 116 orthologous gene groups between teleosts and tetrapods, and 45 pairs of 3R-WGD-derived duplicate genes among them. This suggests that more than half [(45×2)/(116+45)] = 56.5%) of the loci, probably more than ten thousand genes, present in a common ancestor of the four teleosts were still duplicated after the 3R-WGD. The estimated temporal pattern of gene loss suggested that, after the 3R-WGD, many (71/116) of the duplicated genes were rapidly lost during the initial 75 million years (MY), whereas on average more than half (27.3/45) of the duplicated genes remaining in the ancestor of the four teleosts (45/116) have persisted for about 275 MY. The 3R-WGD-derived duplicates that have persisted for a long evolutionary periods of time had significantly larger number of interacting partners and longer length of protein coding sequence, implying that they tend to be more multifunctional than the singletons after the 3R-WGD.ConclusionWe have shown firstly the temporal pattern of gene loss process after 3R-WGD on the basis of teleost phylogeny and divergence time frameworks. The 3R-WGD-derived duplicates have not undergone constant exponential decay, suggesting that selection favoured the long-term persistence of a subset of duplicates that tend to be multi-functional. On the basis of these results obtained from the analysis of 116 orthologous gene groups, we propose that more than ten thousand of 3R-WGD-derived duplicates have experienced lineage-specific evolution, that is, the differential sub-/neo-functionalization or secondary loss between lineages, and contributed to teleost diversity.

Post-duplication charge evolution of phosphoglucose isomerases in teleost fishes through weak selection on many amino acid sites

BackgroundThe partitioning of ancestral functions among duplicated genes by neutral evolution, or subfunctionalization, has been considered the primary process for the evolution of novel proteins (neofunctionalization). Nonetheless, how a subfunctionalized protein can evolve into a more adaptive protein is poorly understood, mainly due to the limitations of current analytical methods, which can detect only strong selection for amino acid substitutions involved in adaptive molecular evolution. In this study, we employed a comparative evolutionary approach to this question, focusing on differences in the structural properties of a protein, specifically the electric charge, encoded by fish-specific duplicated phosphoglucose isomerase (Pgi) genes.ResultsFull-length cDNA cloning, RT-PCR based gene expression analyses, and comparative sequence analyses showed that after subfunctionalization with respect to the expression organ of duplicate Pgi genes, the net electric charge of the PGI-1 protein expressed mainly in internal tissues became more negative, and that of PGI-2 expressed mainly in muscular tissues became more positive. The difference in net protein charge was attributable not to specific amino acid sites but to the sum of various amino acid sites located on the surface of the PGI molecule.ConclusionThis finding suggests that the surface charge evolution of PGI proteins was not driven by strong selection on individual amino acid sites leading to permanent fixation of a particular residue, but rather was driven by weak selection on a large number of amino acid sites and consequently by steady directional and/or purifying selection on the overall structural properties of the protein, which is derived from many modifiable sites. The mode of molecular evolution presented here may be relevant to various cases of adaptive modification in proteins, such as hydrophobic properties, molecular size, and electric charge.

Comparative evolutionary psychology of sperm competition.

A comparative evolutionary psychological perspective predicts that species that recurrently faced similar adaptive problems may have evolved similar psychological mechanisms to solve these problems. Sperm competition provides an arena in which to assess the heuristic value of such a comparative evolutionary perspective. The sperm competition that results from female infidelity and polyandry presents a similar class of adaptive problems for individuals across many species. The authors first describe mechanisms of sperm competition in insects and in birds. They suggest that the adaptive problems and evolved solutions in these species provide insight into human anatomy, physiology, psychology, and behavior. The authors then review recent theoretical and empirical arguments for the existence of sperm competition in humans and discuss proposed adaptations in humans that have analogs in insects or birds. The authors conclude by highlighting the heuristic value of a comparative evolutionary psychological approach in this field.

Natural selection and mammalian BRCA1 sequences: elucidating functionally important sites relevant to breast cancer susceptibility in humans

Comparison of orthologous gene sequences is emerging as a powerful approach to elucidating functionally important positions in human disease genes. Using a diverse array of 132 mammalian BRCA1 (exon 11) sequences, we evaluated the functional significance of specific sites in the context of selection information (purifying, neutral, or diversifying) as well as the ability to extract such information from alignments that index varying degrees of mammalian diversity. Small data sets of either closely related taxa (Primates) or divergent placental taxa were unable to distinguish sites conserved due to purifying selection from sites conserved due to chance (false-positive rate = 65%-99%). Increasing the number of placental taxa to 57 greatly reduced the potential false-positive rate (0%-1.5%). Using the larger data set, we ranked the oncogenic risk of human missense mutations using a novel method that incorporates site-specific selection level and severity of the amino acid change evaluated against the amino acids present in other mammalian taxa. In addition to sites undergoing positive selection in Marsupialia, Laurasiatheria, Euarchontoglires, and Primates, we identified sites most likely to be undergoing divergent selection pressure in different lineages and six pairs of potentially interacting sites. Our results demonstrate the necessity of including large numbers of sequences to elucidate functionally important sites of a protein when using a comparative evolutionary approach.