Selection In Lineages Research Articles

Adaptive Evolution of the Myo6 Gene in Old World Fruit Bats (Family: Pteropodidae)

Myosin VI (encoded by the Myo6 gene) is highly expressed in the inner and outer hair cells of the ear, retina, and polarized epithelial cells such as kidney proximal tubule cells and intestinal enterocytes. The Myo6 gene is thought to be involved in a wide range of physiological functions such as hearing, vision, and clathrin-mediated endocytosis. Bats (Chiroptera) represent one of the most fascinating mammal groups for molecular evolutionary studies of the Myo6 gene. A diversity of specialized adaptations occur among different bat lineages, such as echolocation and associated high-frequency hearing in laryngeal echolocating bats, large eyes and a strong dependence on vision in Old World fruit bats (Pteropodidae), and specialized high-carbohydrate but low-nitrogen diets in both Old World and New World fruit bats (Phyllostomidae). To investigate what role(s) the Myo6 gene might fulfill in bats, we sequenced the coding region of the Myo6 gene in 15 bat species and used molecular evolutionary analyses to detect evidence of positive selection in different bat lineages. We also conducted real-time PCR assays to explore the expression levels of Myo6 in a range of tissues from three representative bat species. Molecular evolutionary analyses revealed that the Myo6 gene, which was widely considered as a hearing gene, has undergone adaptive evolution in the Old World fruit bats which lack laryngeal echolocation and associated high-frequency hearing. Real-time PCR showed the highest expression level of the Myo6 gene in the kidney among ten tissues examined in three bat species, indicating an important role for this gene in kidney function. We suggest that Myo6 has undergone adaptive evolution in Old World fruit bats in relation to receptor-mediated endocytosis for the preservation of protein and essential nutrients.

New insights into ligand-receptor pairing and coevolution of relaxin family peptides and their receptors in teleosts.

Relaxin-like peptides (RLN/INSL) play diverse roles in reproductive and neuroendocrine processes in placental mammals and are functionally associated with two distinct types of receptors (RXFP) for each respective function. The diversification of RLN/INSL and RXFP gene families in vertebrates was predominantly driven by whole genome duplications (2R and 3R). Teleosts preferentially retained duplicates of genes putatively involved in neuroendocrine regulation, harboring a total of 10-11 receptors and 6 ligand genes, while most mammals have equal numbers of ligands and receptors. To date, the ligand-receptor relationships of teleost Rln/Insl peptides and their receptors have largely remained unexplored. Here, we use selection analyses based on sequence data from 5 teleosts and qPCR expression data from zebrafish to explore possible ligand-receptor pairings in teleosts. We find support for the hypothesis that, with the exception of RLN, which has undergone strong positive selection in mammalian lineages, the ligand and receptor genes shared between mammals and teleosts appear to have similar pairings. On the other hand, the teleost-specific receptors show evidence of subfunctionalization. Overall, this study underscores the complexity of RLN/INSL and RXFP ligand-receptor interactions in teleosts and establishes theoretical background for further experimental work in nonmammals.

Differential strengths of selection on S-RNases from Physalis and Solanum(Solanaceae)

BackgroundThe S-RNases of the Solanaceae are highly polymorphic self-incompatibility (S-) alleles subject to strong balancing selection. Relatively recent diversification of S-alleles has occurred in the genus Physalis following a historical restriction of S-allele diversity. In contrast, the genus Solanum did not undergo a restriction of S-locus diversity and its S-alleles are generally much older. Because recovery from reduced S-locus diversity should involve increased selection, we employ a statistical framework to ask whether S-locus selection intensities are higher in Physalis than Solanum. Because different S-RNase lineages diversify in Physalis and Solanum, we also ask whether different sites are under selection in different lineages.ResultsMaximum-likelihood and Bayesian coalescent methods found higher intensities of selection and more sites under significant positive selection in the 48 Physalis S-RNase alleles than the 49 from Solanum. Highest posterior densities of dN/dS (ω) estimates show that the strength of selection is greater for Physalis at 36 codons. A nested maximum likelihood method was more conservative, but still found 16 sites with greater selection in Physalis. Neither method found any codons under significantly greater selection in Solanum. A random effects likelihood method that examines data from both taxa jointly confirmed higher selection intensities in Physalis, but did not find different proportions of sites under selection in the two datasets. The greatest differences in strengths of selection were found in the most variable regions of the S-RNases, as expected if these regions encode self-recognition specificities. Clade-specific likelihood models indicated some codons were under greater selection in background Solanum lineages than in specific lineages of Physalis implying that selection on sites may differ among lineages.ConclusionsLikelihood and Bayesian methods provide a statistical approach to testing differential selection across populations or species. These tests appear robust to the levels of polymorphism found in diverse S-allele collections subject to strong balancing selection. As predicted, the intensity of selection at the S-locus was higher in the taxon with more recent S-locus diversification. This is the first confirmation by statistical test of differing selection intensities among self-incompatibility alleles from different populations or species.

Signature of selection on the rhodopsin gene in the marine radiation of American seven‐spined gobies (Gobiidae, Gobiosomatini)

In comparison with terrestrial and freshwater ecosystems, information about speciation modes and the role of selection in marine environments is scarce. Recent studies have indicated that spectral adaptation could play an important role in the diversification of marine species flocks. Natural selection influences specific amino acids (AAs) that are involved in the spectral tuning mechanism of visual pigment genes. To study the wider occurrence and the characteristics of spectral adaptation in marine radiations, a reinterpretation of the rhodopsin (RH1) data of American seven-spined gobies (genus Elacatinus; Gobiidae; Teleostei) was carried out. Reanalysis revealed that some AAs, which are well known in the literature as spectral tuning sites, are variable in Elacatinus. Those crucial AA substitutions originated polyphyletically, indicating convergent evolution within the genus Elacatinus. Moreover, statistical tests based on the d(N)/d(S) ratio detected selection in several phylogenetic lineages and at specific AAs. Many of these AAs were previously shown to be under selection in other marine radiations. Therefore, the current phylogenetic approach provided an extended list of AAs that are probably involved in spectral tuning, and which should be validated by mutagenic experiments.

The phylogeny of monkey beetles based on mitochondrial and ribosomal RNA genes (Coleoptera: Scarabaeidae: Hopliini)

Monkey beetles (Hopliini) are a large clade of flower and leaf feeding species within the Scarabaeidae (chafers) with greatest diversity in southern Africa. Their internal relationships and sister group affinities have not been studied with DNA methods. We used partial gene sequences for 28S rRNA, cytochrome oxidase I (cox1) and 16S rRNA (rrnL) for 158 species, representing most recognized subfamilies of Scarabaeidae, including 46 species of Hopliini. Combined analyses using maximum likelihood and Bayesian inference under the two preferred alignment parameters recovered the Hopliini as monophyletic. Hopliines were inserted at the base of a clade of Cetoniinae+Rutelinae+Dynastinae, being either recovered as their immediate sister group, or as part of an expanded set of basal branches that also includes the tribe Macrodactylini which has been classified as part of the Melolonthinae (may chafers). At the level of subtribes, we found Hopliina paraphyletic with respect to Pachycnemina which also includes the monophyletic clade of Heterochelina and Gymnolomina. Trait mapping under parsimony on the preferred tree resulted in inferences of three independent origins of sexual dimorphism, which coincided with shifts to ‘flower-embedding’ pollination. In contrast, night active taxa, which are general phyllophages as other pleurostict chafers, never show clear sexual dimorphism. South African lineages include several deep-branching lineages. The exceptional morphological and phylogenetic diversity of the South African fauna may therefore be due to their antiquity, in addition to sexual selection in the day-active lineages. Phylogenetic studies of the endemic South African plant radiations have demonstrated the repeated evolutionary shift to beetle pollination, but it remains to be investigated if this is driven by the hopliine pollinators present in the bioregion or by a propensity of the local plant lineages favoring this pollination syndrome.

Population genetics of the cytoplasm and the units of selection on mitochondrial DNA in Drosophila melanogaster

Biological variation exists across a nested set of hierarchical levels from nucleotides within genes to populations within species to lineages within the tree of life. How selection acts across this hierarchy is a long-standing question in evolutionary biology. Recent studies have suggested that genome size is influenced largely by the balance of selection, mutation and drift in lineages with different population sizes. Here we use population cage and maternal transmission experiments to identify the relative strength of selection at an individual and cytoplasmic level. No significant trends were observed in the frequency of large (L) and small (S) mtDNAs across 14 generations in population cages. In all replicate cages, new length variants were observed in heteroplasmic states indicating that spontaneous length mutations occurred in these experimental populations. Heteroplasmic flies carrying L genomes were more frequent than those carrying S genomes suggesting an asymmetric mutation dynamic from larger to smaller mtDNAs. Mother-offspring transmission of heteroplasmy showed that the L mtDNA increased in frequency within flies both between and within generations despite sampling drift of the same intensity as occurred in population cages. These results suggest that selection for mtDNA size is stronger at the cytoplasmic than at the organismal level. The fixation of novel mtDNAs within and between species requires a transient intracellular heteroplasmic stage. The balance of population genetic forces at the cytoplasmic and individual levels governs the units of selection on mtDNA, and has implications for evolutionary inference as well as for the effects of mtDNA mutations on fitness, disease and aging.

Unlocking the Barley Genome by Chromosomal and Comparative Genomics

We used a novel approach that incorporated chromosome sorting, next-generation sequencing, array hybridization, and systematic exploitation of conserved synteny with model grasses to assign ~86% of the estimated ~32,000 barley (Hordeum vulgare) genes to individual chromosome arms. Using a series of bioinformatically constructed genome zippers that integrate gene indices of rice (Oryza sativa), sorghum (Sorghum bicolor), and Brachypodium distachyon in a conserved synteny model, we were able to assemble 21,766 barley genes in a putative linear order. We show that the barley (H) genome displays a mosaic of structural similarity to hexaploid bread wheat (Triticum aestivum) A, B, and D subgenomes and that orthologous genes in different grasses exhibit signatures of positive selection in different lineages. We present an ordered, information-rich scaffold of the barley genome that provides a valuable and robust framework for the development of novel strategies in cereal breeding.

Base composition, selection, and phylogenetic significance of indels in the recombination activating gene-1 in vertebrates

BackgroundThe Recombination Activating Proteins, RAG1 and RAG2, play a crucial role in the immune response in vertebrates. Among the nuclear markers currently used for phylogenetic purposes, Rag1 has especially enjoyed enormous popularity, since it successfully contributed to elucidating the relationships among and within a large variety of vertebrate lineages. We here report on a comparative investigation of the genetic variation, base composition, presence of indels, and selection in Rag1 in different vertebrate lineages (Actinopterygii, Amphibia, Aves, Chondrichthyes, Crocodylia, Lepidosauria, Mammalia, and Testudines) through the analysis of 582 sequences obtained from Genbank. We also analyze possible differences between distinct parts of the gene with different type of protein functions.ResultsIn the vertebrate lineages studied, Rag1 is over 3 kb long. We observed a high level of heterogeneity in base composition at the 3rd codon position in some of the studied vertebrate lineages and in some specific taxa. This result is also paralleled by taxonomic differences in the GC content at the same codon position. Moreover, positive selection occurs at some sites in Aves, Lepidosauria and Testudines. Indels, which are often used as phylogenetic characters, are more informative across vertebrates in the 5' than in the 3'-end of the gene. When the entire gene is considered, the use of indels as phylogenetic character only recovers one major vertebrate clade, the Actinopterygii. However, in numerous cases insertions or deletions are specific to a monophyletic group.ConclusionsRag1 is a phylogenetic marker of undoubted quality. Our study points to the need of carrying out a preliminary investigation on the base composition and the possible existence of sites under selection of this gene within the groups studied to avoid misleading resolution. The gene shows highly heterogeneous base composition, which affects some taxa in particular and contains sites under positive selection in some vertebrate lineages in the 5'-end. The first part of the gene (5'-end) is more variable than the second (3'-end), and less affected by a heterogeneous base composition. However, in some vertebrate lineages the 5'-end of the gene is not yet widely used for phylogenetic studies.

Pervasive, genome-wide positive selection leading to functional divergence in the bacterial genus Campylobacter

An open question in bacterial genomics is the role that adaptive evolution of the core genome plays in diversification and adaptation of bacterial species, and how this might differ between groups of bacteria occupying different environmental circumstances. The genus Campylobacter encompasses several important human and animal enteric pathogens, with genome sequence data available for eight species. We estimate the Campylobacter core genome at 647 genes, with 92.5% of the nonrecombinant core genome loci under positive selection in at least one lineage and the same gene frequently under positive selection in multiple lineages. Tests are provided that reject recombination, saturation, and variation in codon usage bias as factors contributing to this high level of selection. We suggest this genome-wide adaptive evolution may result from a Red Queen macroevolutionary dynamic, in which species are involved in competition for resources within the mammalian and/or vertebrate gastrointestinal tract. Much reduced levels of positive selection evident in Streptococcus, as reported by the authors in an earlier work, may be a consequence of these taxa inhabiting less species-rich habitats, and more unique niches. Despite many common loci under positive selection in multiple Campylobacter lineages, we found no evidence for molecular adaptive convergence at the level of the same or adjacent codons, or even protein domains. Taken collectively, these results describe the diversification of a bacterial genus that involves pervasive natural selection pressure across virtually the entire genome, with this adaptation occurring in different ways in different lineages, despite the species tendency toward a common gastrointestinal habitat.

Molecular adaptation of rbcL in the heterophyllous aquatic plant Potamogeton.

BackgroundHeterophyllous aquatic plants show marked phenotypic plasticity. They adapt to environmental changes by producing different leaf types: submerged, floating and terrestrial leaves. By contrast, homophyllous plants produce only submerged leaves and grow entirely underwater. Heterophylly and submerged homophylly evolved under selective pressure modifying the species-specific optima for photosynthesis, but little is known about the evolutionary outcome of habit. Recent evolutionary analyses suggested that rbcL, a chloroplast gene that encodes a catalytic subunit of RuBisCO, evolves under positive selection in most land plant lineages. To examine the adaptive evolutionary process linked to heterophylly or homophylly, we analyzed positive selection in the rbcL sequences of ecologically diverse aquatic plants, Japanese Potamogeton.Principal FindingsPhylogenetic and maximum likelihood analyses of codon substitution models indicated that Potamogeton rbcL has evolved under positive Darwinian selection. The positive selection has operated specifically in heterophyllous lineages but not in homophyllous ones in the branch-site models. This suggests that the selective pressure on this chloroplast gene was higher for heterophyllous lineages than for homophyllous lineages. The replacement of 12 amino acids occurred at structurally important sites in the quaternary structure of RbcL, two of which (residue 225 and 281) were identified as potentially under positive selection.Conclusions/SignificanceOur analysis did not show an exact relationship between the amino acid replacements and heterophylly or homophylly but revealed that lineage-specific positive selection acted on the Potamogeton rbcL. The contrasting ecological conditions between heterophyllous and homophyllous plants have imposed different selective pressures on the photosynthetic system. The increased amino acid replacement in RbcL may reflect the continuous fine-tuning of RuBisCO under varying ecological conditions.

Possible diversifying selection in the imprinted gene, MEDEA, in Arabidopsis.

Coevolutionary conflict among imprinted genes that influence traits such as offspring growth may arise when maternal and paternal genomes have different evolutionary optima. This conflict is expected in outcrossing taxa with multiple paternity, but not self-fertilizing taxa. MEDEA (MEA) is an imprinted plant gene that influences seed growth. Disagreement exists regarding the type of selection acting on this gene. We present new data and analyses of sequence diversity of MEA in self-fertilizing and outcrossing Arabidopsis and its relatives, to help clarify the form of selection acting on this gene. Codon-based branch analysis among taxa (PAML) suggests that selection on the coding region is changing over time, and nonsynonymous substitution is elevated in at least one outcrossing branch. Codon-based analysis of diversity within outcrossing Arabidopsis lyrata ssp. petraea (OmegaMap) suggests that diversifying selection is acting on a portion of the gene, to cause elevated nonsynonymous polymorphism. Providing further support for balancing selection in A. lyrata, Hudson, Kreitman and Aguadé analysis indicates that diversity/divergence at silent sites in the MEA promoter and genic region is elevated relative to reference genes, and there are deviations from the neutral frequency spectrum. This combination of positive selection as well as balancing and diversifying selection in outcrossing lineages is consistent with other genes influence by evolutionary conflict, such as disease resistance genes. Consistent with predictions that conflict would be eliminated in self-fertilizing taxa, we found no evidence of positive, balancing, or diversifying selection in A. thaliana promoter or genic region.

Evolutionary Switch and Genetic Convergence on rbcL following the Evolution of C4 Photosynthesis

Rubisco is responsible for the fixation of CO2 into organic compounds through photosynthesis and thus has a great agronomic importance. It is well established that this enzyme suffers from a slow catalysis, and its low specificity results into photorespiration, which is considered as an energy waste for the plant. However, natural variations exist, and some Rubisco lineages, such as in C4 plants, exhibit higher catalytic efficiencies coupled to lower specificities. These C4 kinetics could have evolved as an adaptation to the higher CO2 concentration present in C4 photosynthetic cells. In this study, using phylogenetic analyses on a large data set of C3 and C4 monocots, we showed that the rbcL gene, which encodes the large subunit of Rubisco, evolved under positive selection in independent C4 lineages. This confirms that selective pressures on Rubisco have been switched in C4 plants by the high CO2 environment prevailing in their photosynthetic cells. Eight rbcL codons evolving under positive selection in C4 clades were involved in parallel changes among the 23 independent monocot C4 lineages included in this study. These amino acids are potentially responsible for the C4 kinetics, and their identification opens new roads for human-directed Rubisco engineering. The introgression of C4-like high-efficiency Rubisco would strongly enhance C3 crop yields in the future CO2-enriched atmosphere.

Recurrent events of positive selection in independent Drosophila lineages at the spermatogenesis gene roughex.

Our understanding of the role of positive selection in the evolution of genes with male-biased expression can be hindered by two observations. First, male-biased genes tend to be overrepresented among lineage-specific genes. Second, novel genes are prone to experience bursts of adaptive evolution shortly after their formation. A thorough study of the forces acting on male-biased genes therefore would benefit from phylogenywide analyses that could distinguish evolutionary trends associated with gene formation and later events, while at the same time tackling the interesting question of whether adaptive evolution is indeed idiosyncratic. Here we investigate the roughex (rux) gene, a dose-dependent regulator of Drosophila spermatogenesis with a C-terminal domain responsible for nuclear localization that shows a distinct amino acid sequence in the melanogaster subgroup. We collected polymorphism and divergence data in eight populations of six Drosophila species, for a total of 99 rux sequences, to study rates and patterns of evolution at this male-biased gene. Our results from two phylogeny-based methods (PAML and HyPhy) as well as from population genetics analyses (McDonald-Kreitman-based tests) indicate that amino acid replacements have contributed disproportionately to divergence, consistent with adaptive evolution at the Rux protein. Analyses based on extant variation show also the signature of recent selective sweeps in several of the populations surveyed. Most important, we detect the significant and consistent signature of positive selection in several independent Drosophila lineages, which evidences recurrent and concurrent events of adaptive evolution after rux formation.

Exploring Variation in the d N /d S Ratio Among Sites and Lineages Using Mutational Mappings: Applications to the Influenza Virus

We use a likelihood-based method for mapping mutations on a phylogeny in a way that allows for both site-specific and lineage-specific variation in selection intensity. The method accounts for many of the potential sources of bias encountered in mapping of mutations on trees while still being computationally efficient. We apply the method to a previously published influenza data set to investigate hypotheses about changes in selection intensity in influenza strains. Influenza virus is sometimes propagated in chicken cells for several generations before sequencing, a process that has been hypothesized to induce mutations adapting the virus to the lab medium. Our analysis suggests that there are approximately twice as many replacement substitutions in lineages propagated in chicken eggs as in lineages that are not. Previous studies have attempted to predict which viral strains future epidemics may arise from using inferences regarding positive selection. The assumption is that future epidemics are more likely to arise from the strains in which positive selection on the so-called "trunk lineages" of the evolutionary tree is most pervasive. However, we find no difference in the strength of selection in the trunk lineages versus other evolutionary lineages. Our results suggest that it may be more difficult to use inferences regarding the strength of selection on mutations to make predictions regarding viral epidemics than previously thought.

Purifying selection and positive selection on the myxovirus resistance gene in mammals and chickens

The myxovirus resistance gene (Mx) expresses antiviral activity in many species, e.g. mouse, human and chicken. It is not clear if the antiviral activity of Mx has evolved in these species to inhibit a set of species-specific pathogens, nor what factors drive Mx evolution in different animal lineages. Therefore, it is important to determine the evolutionary pattern of Mx and positively selected sites which affect the antiviral activity of the Mx gene in mammals and birds. We used sequence comparisons among species to detect positively selected sites by conducting phylogenetic analysis. The two-ratio model was significantly better than the one-ratio model in four species (mouse, rat, chicken and duck, p < 0.05). Although selection pressure varied among different lineages, Mx had strong purifying selection in mammals and positive selection in chicken and duck lineages. Relative rate test revealed that Mx evolved faster in chickens than in ducks (Tajima's relative rate test, χ 2 = 7.17, p < 0.01). In the further analysis using a branch-site model A test, 8 sites were positively selected in the chicken lineage while no positive selection signals were observed for any site in the other lineages. The branch-site model A test had a ω value of 4.374 for the chicken lineage (2Δ ℓ = 14.20, d.f. = 1, p < 0.001). Comparisons of all currently available Mx mRNA sequences showed that these predicted positively selected sites had been fixed in the chicken lineage, suggesting that the chicken Mx gene evolved within the species to resist newly challenging environments. There is an increased selection constraint leading to mammals, while positive selection has acted on the chicken Mx.

Detecting positive darwinian selection in brain-expressed genes during human evolution

To understand the genetic basis that underlies the phenotypic divergence between human and non-human primates, we screened a total of 7176 protein-coding genes expressed in the human brain and compared them with the chimpanzee orthologs to identify genes that show evidence of rapid evolution in the human lineage. Our results showed that the nonsynonymous/synonymous substitution (Ka/Ks) ratio for genes expressed in the brain of human and chimpanzee is 0.3854, suggesting that the brain-expressed genes are under functional constraint. The X-linked human brain-expressed genes evolved more rapidly than autosomal ones. We further dissected the molecular evolutionary patterns of 34 candidate genes by sequencing representative primate species to identify lineage-specific adaptive evolution. Fifteen out of the 34 candidate genes showed evidence of positive Darwinian selection in human and/or chimpanzee lineages. These genes are predicted to play diverse functional roles in embryonic development, spermatogenesis and male fertility, signal transduction, sensory nociception, and neural function. This study together with others demonstrated the usefulness and power of phylogenetic comparison of multiple closely related species in detecting lineage-specific adaptive evolution, and the identification of the positively selected brain-expressed genes may add new knowledge to the understanding of molecular mechanism of human origin.

Widespread positive selection in the photosynthetic Rubisco enzyme.

BackgroundRubisco enzyme catalyzes the first step in net photosynthetic CO2 assimilation and photorespiratory carbon oxidation and is responsible for almost all carbon fixation on Earth. The large subunit of Rubisco is encoded by the chloroplast rbcL gene, which is widely used for reconstruction of plant phylogenies due to its conservative nature. Plant systematicists have mainly used rbcL paying little attention to its function, and the question whether it evolves under Darwinian selection has received little attention. The purpose of our study was to evaluate how common is positive selection in Rubisco among the phototrophs and where in the Rubisco structure does positive selection occur.ResultsWe searched for positive selection in rbcL sequences from over 3000 species representing all lineages of green plants and some lineages of other phototrophs, such as brown and red algae, diatoms, euglenids and cyanobacteria. Our molecular phylogenetic analysis found the presence of positive selection in rbcL of most analyzed land plants, but not in algae and cyanobacteria. The mapping of the positively selected residues on the Rubisco tertiary structure revealed that they are located in regions important for dimer-dimer, intradimer, large subunit-small subunit and Rubisco-Rubisco activase interactions, and that some of the positively selected residues are close to the active site.ConclusionOur results demonstrate that despite its conservative nature, Rubisco evolves under positive selection in most lineages of land plants, and after billions of years of evolution Darwinian selection still fine-tunes its performance. Widespread positive selection in rbcL has to be taken into account when this gene is used for phylogenetic reconstructions.

SPIDER DRAGLINE SILK: CORRELATED AND MOSAIC EVOLUTION IN HIGH-PERFORMANCE BIOLOGICAL MATERIALS

. –The evolution of biological materials is a critical, yet poorly understood, component in the generation of biodiversity. For example, the diversification of spiders is correlated with evolutionary changes in the way they use silk, and the material properties of these fibers, such as strength, toughness, extensibility, and stiffness, have profound effects on ecological function. Here, we examine the evolution of the material properties of dragline silk across a phylogenetically diverse sample of species in the Araneomorphae (true spiders). The silks we studied are generally stronger than other biological materials and tougher than most biological or man‐made fibers, but their material properties are highly variable; for example, strength and toughness vary more than fourfold among the 21 species we investigated. Furthermore, associations between different properties are complex. Some traits, such as strength and extensibility, seem to evolve independently and show no evidence of correlation or trade‐off across species, even though trade‐offs between these properties are observed within species. Material properties retain different levels of phylogenetic signal, suggesting that traits such as extensibility and toughness may be subject to different types or intensities of selection in several spider lineages. The picture that emerges is complex, with a mosaic pattern of trait evolution producing a diverse set of materials across spider species. These results show that the properties of biological materials are the target of selection, and that these changes can produce evolutionarily and ecologically important diversity.

Functional Diversification of B MADS-Box Homeotic Regulators of Flower Development: Adaptive Evolution in Protein–Protein Interaction Domains after Major Gene Duplication Events

B-class MADS-box genes have been shown to be the key regulators of petal and stamen specification in several eudicot model species such as Arabidopsis thaliana, Antirrhinum majus, and Petunia hybrida. Orthologs of these genes have been found across angiosperms and gymnosperms, and it is thought that the basic regulatory function of B proteins is conserved in seed plant lineages. The evolution of B genes is characterized by numerous duplications that might represent key elements fostering the functional diversification of duplicates with a deep impact on their role in the evolution of the floral developmental program. To evaluate this, we performed a rigorous statistical analysis with B gene sequences. Using maximum likelihood and Bayesian methods, we estimated molecular substitution rates and determined the selective regimes operating at each residue of B proteins. We implemented tests that rely on phylogenetic hypotheses and codon substitution models to detect significant differences in substitution rates (DSRs) and sites under positive adaptive selection (PS) in specific lineages before and after duplication events. With these methods, we identified several protein residues fixed by PS shortly after the origin of PISTILLATA-like and APETALA3-like lineages in angiosperms and shortly after the origin of the euAP3-like lineage in core eudicots, the 2 main B gene duplications. The residues inferred to have been fixed by positive selection lie mostly within the K domain of the protein, which is key to promote heterodimerization. Additionally, we used a likelihood method that accommodates DSRs among lineages to estimate duplication dates for AP3-PI and euAP3-TM6, calibrating with data from the fossil record. The dates obtained are consistent with angiosperm origins and diversification of core eudicots. Our results strongly suggest that novel multimer formation with other MADS proteins could have been crucial for the functional divergence of B MADS-box genes. We thus propose a mechanism of functional diversification and persistence of gene duplicates by the appearance of novel multimerization capabilities after duplications. Multimer formation in different combinations of regulatory proteins can be a mechanistic basis for the origin of novel regulatory functions and a gene regulatory mechanism for the appearance of morphological innovations.

Genetic conflicts during meiosis and the evolutionary origins of centromere complexity

Centromeric DNA evolves rapidly, ranging in size and complexity over several orders of magnitude. Traditional attempts at studying centromeres have left unexplained the causes underlying this complexity and rapid evolution. Instead of directly studying centromeric DNA sequence, our approach has been to study the proteins that epigenetically determine centromere identity. We have discovered that centromeric histones (CenH3s) have evolved under positive selection in multiple lineages, suggesting an involvement in recurrent genetic conflict. Our hypothesis is that 'centromere-drive' is the source of this conflict. Under this model, centromeres compete via microtubule attachments for preferential transmission in female meioses occurring in animals and plants. Since only one of four meiotic products will become the egg, this competition confers a selfish advantage to chromosomes that can make more microtubule attachments, resulting in runaway expansions of centromeric satellites. While beneficial to the 'driving' chromosome, these expansions can have deleterious effects on the fitness of an organism and of the species. CenH3s as well as other heterochromatin proteins have evolved under positive selection to suppress the deleterious consequences of 'centromere-drive' by restoring meiotic parity.