Relaxed Purifying Selection Research Articles

Mobile Introns Shape the Genetic Diversity of Their Host Genes.

Self-splicing introns populate several highly conserved protein-coding genes in fungal and plant mitochondria. In fungi, many of these introns have retained their ability to spread to intron-free target sites, often assisted by intron-encoded endonucleases that initiate the homing process. Here, leveraging population genomic data from Saccharomyces cerevisiae, Schizosaccharomyces pombe, and Lachancea kluyveri, we expose nonrandom patterns of genetic diversity in exons that border self-splicing introns. In particular, we show that, in all three species, the density of single nucleotide polymorphisms increases as one approaches a mobile intron. Through multiple lines of evidence, we rule out relaxed purifying selection as the cause of uneven nucleotide diversity. Instead, our findings implicate intron mobility as a direct driver of host gene diversity. We discuss two mechanistic scenarios that are consistent with the data: either endonuclease activity and subsequent error-prone repair have left a mutational footprint on the insertion environment of mobile introns or nonrandom patterns of genetic diversity are caused by exonic coconversion, which occurs when introns spread to empty target sites via homologous recombination. Importantly, however, we show that exonic coconversion can only explain diversity gradients near intron-exon boundaries if the conversion template comes from outside the population. In other words, there must be pervasive and ongoing horizontal gene transfer of self-splicing introns into extant fungal populations.

Chance and necessity in the genome evolution of endosymbiotic bacteria of insects

An open question in evolutionary biology is how does the selection–drift balance determine the fates of biological interactions. We searched for signatures of selection and drift in genomes of five endosymbiotic bacterial groups known to evolve under strong genetic drift. Although most genes in endosymbiotic bacteria showed evidence of relaxed purifying selection, many genes in these bacteria exhibited stronger selective constraints than their orthologs in free-living bacterial relatives. Remarkably, most of these highly constrained genes had no role in the host–symbiont interactions but were involved in either buffering the deleterious consequences of drift or other host-unrelated functions, suggesting that they have either acquired new roles or their role became more central in endosymbiotic bacteria. Experimental evolution of Escherichia coli under strong genetic drift revealed remarkable similarities in the mutational spectrum, genome reduction patterns and gene losses to endosymbiotic bacteria of insects. Interestingly, the transcriptome of the experimentally evolved lines showed a generalized deregulation of the genome that affected genes encoding proteins involved in mutational buffering, regulation and amino acid biosynthesis, patterns identical to those found in endosymbiotic bacteria. Our results indicate that drift has shaped endosymbiotic associations through a change in the functional landscape of bacterial genes and that the host had only a small role in such a shift.

Positive Selection in Rapidly Evolving Plastid-Nuclear Enzyme Complexes.

Rates of sequence evolution in plastid genomes are generally low, but numerous angiosperm lineages exhibit accelerated evolutionary rates in similar subsets of plastid genes. These genes include clpP1 and accD, which encode components of the caseinolytic protease (CLP) and acetyl-coA carboxylase (ACCase) complexes, respectively. Whether these extreme and repeated accelerations in rates of plastid genome evolution result from adaptive change in proteins (i.e., positive selection) or simply a loss of functional constraint (i.e., relaxed purifying selection) is a source of ongoing controversy. To address this, we have taken advantage of the multiple independent accelerations that have occurred within the genus Silene (Caryophyllaceae) by examining phylogenetic and population genetic variation in the nuclear genes that encode subunits of the CLP and ACCase complexes. We found that, in species with accelerated plastid genome evolution, the nuclear-encoded subunits in the CLP and ACCase complexes are also evolving rapidly, especially those involved in direct physical interactions with plastid-encoded proteins. A massive excess of nonsynonymous substitutions between species relative to levels of intraspecific polymorphism indicated a history of strong positive selection (particularly in CLP genes). Interestingly, however, some species are likely undergoing loss of the native (heteromeric) plastid ACCase and putative functional replacement by a duplicated cytosolic (homomeric) ACCase. Overall, the patterns of molecular evolution in these plastid-nuclear complexes are unusual for anciently conserved enzymes. They instead resemble cases of antagonistic coevolution between pathogens and host immune genes. We discuss a possible role of plastid-nuclear conflict as a novel cause of accelerated evolution.

Positive and strongly relaxed purifying selection drive the evolution of repeats in proteins.

Protein repeats are considered hotspots of protein evolution, associated with acquisition of new functions and novel phenotypic traits, including disease. Paradoxically, however, repeats are often strongly conserved through long spans of evolution. To resolve this conundrum, it is necessary to directly compare paralogous (horizontal) evolution of repeats within proteins with their orthologous (vertical) evolution through speciation. Here we develop a rigorous methodology to identify highly periodic repeats with significant sequence similarity, for which evolutionary rates and selection (dN/dS) can be estimated, and systematically characterize their evolution. We show that horizontal evolution of repeats is markedly accelerated compared with their divergence from orthologues in closely related species. This observation is universal across the diversity of life forms and implies a biphasic evolutionary regime whereby new copies experience rapid functional divergence under combined effects of strongly relaxed purifying selection and positive selection, followed by fixation and conservation of each individual repeat.

The Interplay of Temperature and Genotype on Patterns of Alternative Splicing in Drosophila melanogaster.

Alternative splicing is the highly regulated process of variation in the removal of introns from premessenger-RNA transcripts. The consequences of alternative splicing on the phenotype are well documented, but the impact of the environment on alternative splicing is not yet clear. We studied variation in alternative splicing among four different temperatures, 13, 18, 23, and 29°, in two Drosophila melanogaster genotypes. We show plasticity of alternative splicing with up to 10% of the expressed genes being differentially spliced between the most extreme temperatures for a given genotype. Comparing the two genotypes at different temperatures, we found <1% of the genes being differentially spliced at 18°. At extreme temperatures, however, we detected substantial differences in alternative splicing—with almost 10% of the genes having differential splicing between the genotypes: a magnitude similar to between species differences. Genes with differential alternative splicing between genotypes frequently exhibit dominant inheritance. Remarkably, the pattern of surplus of differences in alternative splicing at extreme temperatures resembled the pattern seen for gene expression intensity. Since different sets of genes were involved for the two phenotypes, we propose that purifying selection results in the reduction of differences at benign temperatures. Relaxed purifying selection at temperature extremes, on the other hand, may cause the divergence in gene expression and alternative splicing between the two strains in rarely encountered environments.

Understanding the evolution of holoparasitic plants: the complete plastid genome of the holoparasite Cytinus hypocistis (Cytinaceae).

Background and Aims Plant plastid genomes are highly conserved in size, gene content and structure; however, parasitic plants are a noticeable exception to this evolutionary stability. Although the evolution of parasites could help to better understand plastome evolution in general, complete plastomes of parasites have been sequenced only for some lineages so far. Here we contribute to filling this gap by providing and analysing the complete plastome sequence of Cytinus hypocistis, the first parasite sequenced for Malvales and a species suspected to have an extremely small genome. Methods We sequenced and assembled de novo the plastid genome of Cytinus hypocistis using a shotgun approach on genomic DNA. Phylogenomic analyses based on coding regions were performed on Malvidae. For each coding region present in Cytinus, we tested for relaxation or intensification of selective pressures in the Cytinus lineage compared with autotrophic Malvales. Key Results Cytinus hypocistis has an extremely divergent genome that is among the smallest sequenced to date (19·4 kb), with only 23 genes and no inverted repeat regions. Phylogenomic analysis confirmed the position of Cytinus within Malvales. All coding regions of Cytinus plastome presented very high substitution rates compared with non-parasitic Malvales. Conclusions Some regions were inferred to be under relaxed negative selection in Cytinus, suggesting that further plastome reduction is occurring due to relaxed purifying selection associated with the loss of photosynthetic activity. On the other hand, increased selection intensity and strong positive selection were detected for rpl22 in the Cytinus lineage, which might indicate an evolutionary role in the host-parasite arms race, a point that needs further research.

Genome Evolution in the Obligate but Environmentally Active Luminous Symbionts of Flashlight Fish.

The luminous bacterial symbionts of anomalopid flashlight fish are thought to be obligately dependent on their hosts for growth and share several aspects of genome evolution with unrelated obligate symbionts, including genome reduction. However, in contrast to most obligate bacteria, anomalopid symbionts have an active environmental phase that may be important for symbiont transmission. Here we investigated patterns of evolution between anomalopid symbionts compared with patterns in free-living relatives and unrelated obligate symbionts to determine if trends common to obligate symbionts are also found in anomalopid symbionts. Two symbionts, “Candidatus Photodesmus katoptron” and “Candidatus Photodesmus blepharus,” have genomes that are highly similar in gene content and order, suggesting genome stasis similar to ancient obligate symbionts present in insect lineages. This genome stasis exists in spite of the symbiont’s inferred ability to recombine, which is frequently lacking in obligate symbionts with stable genomes. Additionally, we used genome comparisons and tests of selection to infer which genes may be particularly important for the symbiont’s ecology compared with relatives. In keeping with obligate dependence, substitution patterns suggest that most symbiont genes are experiencing relaxed purifying selection compared with relatives. However, genes involved in motility and carbon storage, which are likely to be used outside the host, appear to be under increased purifying selection. Two chemoreceptor chemotaxis genes are retained by both species and show high conservation with amino acid sensing genes, suggesting that the bacteria may actively seek out hosts using chemotaxis toward amino acids, which the symbionts are not able to synthesize.

Targeted sequencing for high-resolution evolutionary analyses following genome duplication in salmonid fish: Proof of concept for key components of the insulin-like growth factor axis

High-throughput sequencing has revolutionised comparative and evolutionary genome biology. It has now become relatively commonplace to generate multiple genomes and/or transcriptomes to characterize the evolution of large taxonomic groups of interest. Nevertheless, such efforts may be unsuited to some research questions or remain beyond the scope of some research groups. Here we show that targeted high-throughput sequencing offers a viable alternative to study genome evolution across a vertebrate family of great scientific interest. Specifically, we exploited sequence capture and Illumina sequencing to characterize the evolution of key components from the insulin-like growth (IGF) signalling axis of salmonid fish at unprecedented phylogenetic resolution. The IGF axis represents a central governor of vertebrate growth and its core components were expanded by whole genome duplication in the salmonid ancestor ~95Ma. Using RNA baits synthesised to genes encoding the complete family of IGF binding proteins (IGFBP) and an IGF hormone (IGF2), we captured, sequenced and assembled orthologous and paralogous exons from species representing all ten salmonid genera. This approach generated 299 novel sequences, most as complete or near-complete protein-coding sequences. Phylogenetic analyses confirmed congruent evolutionary histories for all nineteen recognized salmonid IGFBP family members and identified novel salmonid-specific IGF2 paralogues. Moreover, we reconstructed the evolution of duplicated IGF axis paralogues across a replete salmonid phylogeny, revealing complex historic selection regimes - both ancestral to salmonids and lineage-restricted - that frequently involved asymmetric paralogue divergence under positive and/or relaxed purifying selection. Our findings add to an emerging literature highlighting diverse applications for targeted sequencing in comparative-evolutionary genomics. We also set out a viable approach to obtain large sets of nuclear genes for any member of the salmonid family, which should enable insights into the evolutionary role of whole genome duplication before additional nuclear genome sequences become available.

Plastid-Nuclear Interaction and Accelerated Coevolution in Plastid Ribosomal Genes in Geraniaceae.

Plastids and mitochondria have many protein complexes that include subunits encoded by organelle and nuclear genomes. In animal cells, compensatory evolution between mitochondrial and nuclear-encoded subunits was identified and the high mitochondrial mutation rates were hypothesized to drive compensatory evolution in nuclear genomes. In plant cells, compensatory evolution between plastid and nucleus has rarely been investigated in a phylogenetic framework. To investigate plastid–nuclear coevolution, we focused on plastid ribosomal protein genes that are encoded by plastid and nuclear genomes from 27 Geraniales species. Substitution rates were compared for five sets of genes representing plastid- and nuclear-encoded ribosomal subunit proteins targeted to the cytosol or the plastid as well as nonribosomal protein controls. We found that nonsynonymous substitution rates (dN) and the ratios of nonsynonymous to synonymous substitution rates (ω) were accelerated in both plastid- (CpRP) and nuclear-encoded subunits (NuCpRP) of the plastid ribosome relative to control sequences. Our analyses revealed strong signals of cytonuclear coevolution between plastid- and nuclear-encoded subunits, in which nonsynonymous substitutions in CpRP and NuCpRP tend to occur along the same branches in the Geraniaceae phylogeny. This coevolution pattern cannot be explained by physical interaction between amino acid residues. The forces driving accelerated coevolution varied with cellular compartment of the sequence. Increased ω in CpRP was mainly due to intensified positive selection whereas increased ω in NuCpRP was caused by relaxed purifying selection. In addition, the many indels identified in plastid rRNA genes in Geraniaceae may have contributed to changes in plastid subunits.

Molecular evolution analysis of WUSCHEL-related homeobox transcription factor family reveals functional divergence among clades in the homeobox region.

Gene families have been shown to play important roles in plant evolution and are associated with diversification and speciation. Genes of WUSCHEL-related homeobox family of transcription factors have important functions in plant development and are correlated with the appearance of evolutionary novelties. There are several published studies related to this family, but little is known about the relationships among the main clades in the phylogeny and the molecular evolution of the family. In this study, we obtained a well-resolved Bayesian phylogenetic tree establishing the relationships among the main clades and determining the position of Selaginella moellendorffii WOX genes. Moreover, a correlation was identified between the number of genes in the genomes and the events of whole-genome duplications. The intron-exon structure is more consistent across the modern clade, which appeared more recently in the WOX evolutionary history, and coincides with the development of higher complexity in plant species. No positive selection was detected among sites through the branches in the tree. However, with regard to the main clades, functional divergence among certain amino acids in the homeodomain region was found. Relaxed purifying selection could be the main driving force in the evolution of these genes and in agreement with some genes have been demonstrated to be functionally redundant.

The Role of microRNAs in the Repeated Parallel Diversification of Lineages of Midas Cichlid Fish from Nicaragua

Cichlid fishes are an ideal model system for studying biological diversification because they provide textbook examples of rapid speciation. To date, there has been little focus on the role of gene regulation during cichlid speciation. However, in recent years, gene regulation has been recognized as a powerful force linking diversification in gene function to speciation. Here, we investigated the potential role of miRNA regulation in the diversification of six cichlid species of the Midas cichlid lineage (Amphilophus spp.) inhabiting the Nicaraguan crater lakes. Using several genomic resources, we inferred 236 Midas miRNA genes that were used to predict the miRNA target sites on 8,232 Midas 3′-UTRs. Using population genomic calculations of SNP diversity, we found the miRNA genes to be more conserved than protein coding genes. In contrast to what has been observed in other cichlid fish, but similar to what has been typically found in other groups, we observed genomic signatures of purifying selection on the miRNA targets by comparing these sites with the less conserved nontarget portion of the 3′-UTRs. However, in one species pair that has putatively speciated sympatrically in crater Lake Apoyo, we recovered a different pattern of relaxed purifying selection and high genetic divergence at miRNA targets. Our results suggest that sequence evolution at miRNA binding sites could be a critical genomic mechanism contributing to the rapid phenotypic evolution of Midas cichlids.

Evolution of a Sigma Factor: An All-In-One of Gene Duplication, Horizontal Gene Transfer, Purifying Selection, and Promoter Differentiation.

Sigma factors are an essential part of bacterial gene regulation and have been extensively studied as far as their molecular mechanisms and protein structure are concerned. However, their molecular evolution, especially for the alternative sigma factors, is poorly understood. Here, we analyze the evolutionary forces that have shaped the rpoH sigma factors within the alphaproteobacteria. We found that an ancient duplication gave rise to two major groups of rpoH sigma factors and that after this event horizontal gene transfer (HGT) occurred in rpoH1 group. We also noted that purifying selection has differentially affected distinct parts of the gene; singularly, the gene segment that encodes the region 4.2, which interacts with the −35 motif of the RpoH-dependent genes, has been under relaxed purifying selection. Furthermore, these two major groups are clearly differentiated from one another regarding their promoter selectivity, as rpoH1 is under the transcriptional control of σ70 and σ32, whereas rpoH2 is under the transcriptional control of σ24. Our results suggest a scenario in which HGT, gene loss, variable purifying selection and clear promoter specialization occurred after the ancestral duplication event. More generally, our study offers insights into the molecular evolution of alternative sigma factors and highlights the importance of analyzing not only the coding regions but also the promoter regions.

Genome-wide polymorphism and signatures of selection in the symbiotic sea anemone Aiptasia.

BackgroundCoral reef ecosystems are declining in response to global climate change and anthropogenic impacts. Yet patterns of standing genetic variation within cnidarian species, a major determinant of adaptive potential, are virtually unknown at genome-scale resolution. We explore patterns of genome-wide polymorphism and identify candidate loci under selection in the sea anemone Aiptasia, an important laboratory model system for studying the symbiosis between corals and dinoflagellate algae of the genus Symbiodinium.ResultsLow coverage genome sequencing revealed large genetic distances among globally widespread lineages, novel candidate targets of selection, and considerably higher heterozygosity than previously reported for Aiptasia. More than 670,000 single nucleotide polymorphisms were identified among 10 Aiptasia individuals including two pairs of genetic clones. Evolutionary relationships based on genome-wide polymorphism supported the current paradigm of a genetically distinct population from the US South Atlantic that harbors diverse Symbiodinium clades. However, anemones from the US South Atlantic demonstrated a striking lack of shared derived polymorphism. Heterozygosity was an important feature shaping nucleotide diversity patterns: at any given SNP site, more than a third of individuals genotyped were heterozygotes, and heterozygosity within individual genomes ranged from 0.37–0.58 %. Analysis of nonsynonymous and synonymous sites suggested that highly heterozygous regions are evolving under relaxed purifying selection compared to the rest of the Aiptasia genome. Genes previously identified as having elevated evolutionary rates in Aiptasia compared to other cnidarians were found in our study to be under strong purifying selection within Aiptasia. Candidate targets of selection, including lectins and genes involved in Rho GTPase signalling, were identified based on unusual signatures of nucleotide diversity, Tajima’s D, and heterozygosity compared to genome-wide averages.ConclusionsThis study represents the first genome-wide analysis of Tajima’s D in a cnidarian. Our results shed light on patterns of intraspecific genome-wide polymorphism in a model for studies of coral-algae symbiosis and present genetic targets for future research on evolutionary and cellular processes in early-diverging metazoans.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-2488-6) contains supplementary material, which is available to authorized users.

Pleistocene divergence across a mountain range and the influence of selection on mitogenome evolution in threatened Australian freshwater cod species.

Climatic differences across a taxon's range may be associated with specific bioenergetic demands and may result in genetics-based metabolic adaptation, particularly in aquatic ectothermic organisms that rely on heat exchange with the environment to regulate key physiological processes. Extending down the east coast of Australia, the Great Dividing Range (GDR) has a strong influence on climate and the evolutionary history of freshwater fish species. Despite the GDR acting as a strong contemporary barrier to fish movement, many species, and species with shared ancestries, are found on both sides of the GDR, indicative of historical dispersal events. We sequenced complete mitogenomes from the four extant species of the freshwater cod genus Maccullochella, two of which occur on the semi-arid, inland side of the GDR, and two on the mesic coastal side. We constructed a dated phylogeny and explored the relative influences of purifying and positive selection in the evolution of mitogenome divergence among species. Results supported mid- to late-Pleistocene divergence of Maccullochella across the GDR (220-710 thousand years ago), bringing forward previously reported dates. Against a background of pervasive purifying selection, we detected potentially functionally relevant fixed amino acid differences across the GDR. Although many amino acid differences between inland and coastal species may have become fixed under relaxed purifying selection in coastal environments rather than positive selection, there was evidence of episodic positive selection acting on specific codons in the Mary River coastal lineage, which has consistently experienced the warmest and least extreme climate in the genus.

The landscape of somatic mutations in protein coding genes in apparently benign human tissues carries signatures of relaxed purifying selection.

Mutations acquired during development and aging lead to inter- and intra-tissue genetic variations. Evidence linking such mutations to complex traits and diseases is rising. We detected somatic mutations in protein-coding regions in 140 benign tissue samples representing nine tissue-types (bladder, breast, liver, lung, prostate, stomach, thyroid, head and neck) and paired blood from 70 donors. A total of 80% of the samples had 2–39 mutations detectable at tissue-level resolution. Factors such as age and smoking were associated with increased burden of detectable mutations, and tissues carried signatures of distinct mutagenic processes such as oxidative DNA damage and transcription-coupled repair. Using mutational signatures, we predicted that majority of the mutations in blood originated in hematopoietic stem and early progenitor cells. Missense to silent mutations ratio and the persistence of potentially damaging mutations in expressed genes carried signatures of relaxed purifying selection. Our findings have relevance for etiology, diagnosis and treatment of diseases including cancer.

Genes with monoallelic expression contribute disproportionately to genetic diversity in humans.

An unexpectedly large number of human autosomal genes are subject to monoallelic expression (MAE). Our analysis of 4,227 such genes reveals surprisingly high genetic variation across human populations. This increased diversity is unlikely to reflect relaxed purifying selection. Remarkably, MAE genes exhibit elevated recombination rate and increased density of hypermutable sequence contexts. However, these factors do not fully account for the increased diversity. We find that the elevated nucleotide diversity of MAE genes is also associated with greater allelic age: their variants tend to be older and are enriched in polymorphisms shared with Neanderthals and chimpanzees. Both synonymous and nonsynonymous alleles in MAE genes have elevated average population frequencies. We also observed strong enrichment of the MAE signature among genes reported to evolve under balancing selection. We propose that an important biological function of widespread MAE might be generation of cell-to-cell heterogeneity; the increased genetic variation contributes to this heterogeneity.

Comparative analysis of complete mitochondrial genomes suggests that relaxed purifying selection is driving high nonsynonymous evolutionary rate of the NADH2 gene in whitefish (Coregonus ssp.)

Several studies have recently reported evidence for positive selection acting on the mitochondrial genome (mitogenome), emphasizing its potential role in adaptive divergence and speciation. In this study we searched 107 full mitogenomes of recently diverged species and lineages of whitefish (Coregonus ssp.) for signals of positive selection. These salmonids show several distinct morphological and ecological differences that may be associated with energetics and therefore potentially positive selection at the mitogenome level. We found that purifying selection and genetic drift were the predominant evolutionary forces acting on the analyzed mitogenomes. However, the NADH dehydrogenase 2 gene (ND2) showed a highly elevated dN/dS ratio compared to the other mitochondrial genes, which was significantly higher in whitefish compared to other salmonids. We therefore further examined nonsynonymous evolution in ND2 by (i) mapping amino acid changes to a protein model structure which showed that they were located away from key functional residues of the protein, (ii) locating them in the sequences of other species of fish (Salmonidae, Anguillidae, Scombridae and Percidae) only to find pronounced overlap of nonsynonymous regions. We thus conclude that relaxed purifying selection is driving the evolution of ND2 by affecting mostly regions that have lower functional relevance.

Non-canonical structure, function and phylogeny of the Bsister MADS-box gene OsMADS30 of rice (Oryza sativa).

Bsister MADS-box genes play key roles in female reproductive organ and seed development throughout seed plants. This view is supported by their high conservation in terms of sequence, expression and function. In grasses, there are three subclades of Bsister genes: the OsMADS29-, the OsMADS30- and the OsMADS31-like genes. Here, we report on the evolution of the OsMADS30-like genes. Our analyses indicate that these genes evolved under relaxed purifying selection and are rather weakly expressed. OsMADS30, the representative of the OsMADS30-like genes from rice (Oryza sativa), shows strong sequence deviations in its 3' region when compared to orthologues from other grass species. We show that this is due to a 2.4-kbp insertion, possibly of a hitherto unknown helitron, which confers a heterologous C-terminal domain to OsMADS30. This putative helitron is not present in the OsMADS30 orthologues from closely related wild rice species, pointing to a relatively recent insertion event. Unlike other Bsister mutants O. sativa plants carrying a T-DNA insertion in the OsMADS30 gene do not show aberrant seed phenotypes, indicating that OsMADS30 likely does not have a canonical 'Bsister function'. However, imaging-based phenotyping of the T-DNA carrying plants revealed alterations in shoot size and architecture. We hypothesize that sequence deviations that accumulated during a period of relaxed selection in the gene lineage that led to OsMADS30 and the alteration of the C-terminal domain might have been a precondition for a potential neo-functionalization of OsMADS30 in O. sativa.

Evolution and divergence of SBP-box genes in land plants.

BackgroundSquamosa promoter binding protein (SBP)-box family genes encode plant-specific transcription factors that control many important biological functions, including phase transition, inflorescence branching, fruit ripening, and copper homeostasis. Nevertheless, the evolutionary patterns of SBP-box genes and evolutionary forces driving them are still not well understood.Methods104 SBP-box gene candidates of five representative land plants were obtained from Phytozome database (v10.3). Phylogenetic combined with gene structure analyses were used to identify SBP-box gene lineages in land plants. Gene copy number and the sequence and structure features were then compared among these different SBP-box lineages. Selection analysis, relative rate tests and expression divergence were finally used to interpret the evolutionary relationships and divergence of SBP-box genes in land plants.ResultsWe investigated 104 SBP-box genes from moss, Arabidopsis, poplar, rice, and maize. These genes are divided into group I and II, and the latter is further divided into two subgroups (subgroup II-1 and II-2) based on phylogenetic analysis. Interestingly, subgroup II-1 genes have similar sequence and structural features to group I genes, whereas subgroup II-2 genes exhibit intrinsic differences on these features, including high copy numbers and the presence of miR156/miR529 regulation. Further analyses indicate that subgroup II-1 genes are constrained by stronger purifying selection and evolve at a lower substitution rate than II-2 genes, just as group I genes do when compared to II genes. Among subgroup II-2 genes, miR156 targets evolve more rapidly than miR529 targets and experience comparatively relaxed purifying selection. These results suggest that group I and subgroup II-1 genes under strong selective constraint are conserved. By contrast, subgroup II-2 genes evolve under relaxed purifying selection and have diversified through gene copy duplications and changes in miR156/529 regulation, which might contribute to morphological diversifications of land plants.ConclusionsOur results indicate that different evolutionary rates and selection strengths lead to differing evolutionary patterns in SBP-box genes in land plants, providing a guide for future functional diversity analyses of these genes.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-1998-y) contains supplementary material, which is available to authorized users.

Target-Driven Evolution of Scorpion Toxins

It is long known that peptide neurotoxins derived from a diversity of venomous animals evolve by positive selection following gene duplication, yet a force that drives their adaptive evolution remains a mystery. By using maximum-likelihood models of codon substitution, we analyzed molecular adaptation in scorpion sodium channel toxins from a specific species and found ten positively selected sites, six of which are located at the core-domain of scorpion α-toxins, a region known to interact with two adjacent loops in the voltage-sensor domain (DIV) of sodium channels, as validated by our newly constructed computational model of toxin-channel complex. Despite the lack of positive selection signals in these two loops, they accumulated extensive sequence variations by relaxed purifying selection in prey and predators of scorpions. The evolutionary variability in the toxin-bound regions of sodium channels indicates that accelerated substitutions in the multigene family of scorpion toxins is a consequence of dealing with the target diversity. This work presents an example of atypical co-evolution between animal toxins and their molecular targets, in which toxins suffered from more prominent selective pressure from the channels of their competitors. Our discovery helps explain the evolutionary rationality of gene duplication of toxins in a specific venomous species.