Evolutionary Comparison of Two Combinatorial Regulators of SBP-Box Genes,MiR156 and MiR529, in Plants

A complete picture of the evolution of miRNA combinatorial regulation requires the synthesis of information on all miRNAs and their targets. MiR156 and miR529 are two combinatorial regulators of squamosa promoter binding protein-like (SBP-box) genes. Previous studies have clarified the evolutionary dynamics of their targets; however, there have been no reports on the evolutionary patterns of two miRNA regulators themselves to date. In this study, we investigated the evolutionary differences between these two miRNA families in extant land plants. Our work found that miR529 precursor, especially of its mature miRNA sequence, has a higher evolutionary rate. Such accelerating evolution of miR529 has significantly effects on its structural stability, and sequence conservation against existence of itself. By contrast, miR156 evolves more rapidly in loop region of the stable secondary structure, which may contribute to its functional diversity. Moreover, miR156 and miR529 genes have distinct rates of loss after identical duplication events. MiR529 genes have a higher average loss rate and asymmetric loss rate in duplicated gene pairs, indicating preferred miR529 gene losses become another predominant mode of inactivation, that are implicated in the contraction of this family. On the contrary, duplicated miR156 genes have a low loss rate, and could serve as another new source for functional diversity. Taken together, these results provide better insight into understanding the evolutionary divergence of miR156 and miR529 family in miRNA combinational regulation network.     

Detecting positive selection within genomes: the problem of biased gene conversion

The identification of loci influenced by positive selection is a major goal of evolutionary genetics. A popular approach is to perform scans of alignments on a genome-wide scale in order to find regions evolving at accelerated rates on a particular branch of a phylogenetic tree. However, positive selection is not the only process that can lead to accelerated evolution. Notably, GC-biased gene conversion (gBGC) is a recombination-associated process that results in the biased fixation of G and C nucleotides. This process can potentially generate bursts of nucleotide substitutions within hotspots of meiotic recombination. Here, we analyse the results of a scan for positive selection on genes on branches across the primate phylogeny. We show that genes identified as targets of positive selection have a significant tendency to exhibit the genomic signature of gBGC. Using a maximum-likelihood framework, we estimate that more than 20 per cent of cases of significantly elevated non-synonymous to synonymous substitution rates ratio (d_N/d_S), particularly in shorter branches, could be due to gBGC. We demonstrate that in some cases, gBGC can lead to very high d_N/d_S (more than 2). Our results indicate that gBGC significantly affects the evolution of coding sequences in primates, often leading to patterns of evolution that can be mistaken for positive selection.     

The Trouble with Sliding Windows and the Selective Pressure in BRCA1

Sliding-window analysis has widely been used to uncover synonymous (silent, d_S) and nonsynonymous (replacement, d_N) rate variation along the protein sequence and to detect regions of a protein under selective constraint (indicated by d_N<d_S) or positive selection (indicated by d_N>d_S). The approach compares two or more protein-coding genes and plots estimates dˆ _S and dˆ _N from each sliding window along the sequence. Here we demonstrate that the approach produces artifactual trends of synonymous and nonsynonymous rate variation, with greater variation in dˆ _S than in dˆ _N. Such trends are generated even if the true d_S and d_N are constant along the whole protein and different codons are evolving independently. Many published tests of negative and positive selection using sliding windows that we have examined appear to be invalid because they fail to correct for multiple testing. Instead, likelihood ratio tests provide a more rigorous framework for detecting signals of natural selection affecting protein evolution. We demonstrate that a previous finding that a particular region of the BRCA1 gene experienced a synonymous rate reduction driven by purifying selection is likely an artifact of the sliding window analysis. We evaluate various sliding-window analyses in molecular evolution, population genetics, and comparative genomics, and argue that the approach is not generally valid if it is not known a priori that a trend exists and if no correction for multiple testing is applied.     

Molecular evolution of <Emphasis Type="Italic">PKD2</Emphasis> gene family in mammals

PKD2 gene encodes a critical cation channel protein that plays important roles in various developmental processes and is usually evolutionarily conserved. In the present study, we analyzed the evolutionary patterns of PKD2 and its homologous genes (PKD2L1, PKD2L2) from nine mammalian species. In this study, we demonstrated the orthologs of PKD2 gene family evolved under a dominant purifying selection force. Our results in combination with the reported evidences from functional researches suggested the entire PKD2 gene family are conserved and perform essential biological roles during mammalian evolution. In rodents, PKD2 gene family members appeared to have evolved more rapidly than other mammalian lineages, probably resulting from relaxation of purifying selection. However, positive selection imposed on synonymous sites also potentially contributed to this case. For the paralogs, our results implied that PKD2L2 genes evolved under a weaker purifying selection constraint than PKD2 and PKD2L1 genes. Interestingly, some loop regions of transmembrane domain of PKD2L2 exhibited higher P _N/P _S ratios than expected, suggesting these regions are more functional divergent in organisms and worthy of special attention. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. Chun Ye, Huan Sun have contributed equally to this work.     

miR156靶基因SPL在植物中的系统分布和进化模式分析

Squamosa promoter binding protein like genes (SPLs)在植物发育过程中具有重要作用。很多SPLs被miR156调节,然而,对于它们在植物中的系统分布和进化模式还知之甚少。本文对9个测序物种（藻类,苔藓,石松,单子叶和双子叶植物）的183个SPLs进行了生物信息学分析。结果表明miR156应答元件（MREs）仅在陆生植物SPLs中发现,藻类中不存在。系统进化分析显示陆生植物SPLs分为两大分支：group I和group II。 MiR156靶基因仅分布于group II,表明它们有着共同的祖先。Group II进一步分为7个亚支(IIa IIg),miR156靶基因分布在除IId外的其余6个亚支的特定SPLs。系统分类与基因结构的相关性反映了SPL靶基因结构上的变化。在进化过程中,它们可能发生外显子的丢失且伴随MRE的丢失。另外,基因重复对SPL靶基因的丰度变化影响很大,尤其是被子植物与低等植物分歧后它们数量明显增加。以拟南芥为模式植物分析发现串联重复和片段重复是SPL靶基因扩张的主要机制。     

Profiling microRNAs in Eucalyptus grandis reveals no mutual relationship between alterations in miR156 and miR172 expression and adventitious root induction during development

Background

The change from juvenile to mature phase in woody plants is often accompanied by a gradual loss of rooting ability, as well as by reduced microRNA (miR) 156 and increased miR172 expression.

Results

We characterized the population of miRNAs of Eucalyptus grandis and compared the gradual reduction in miR156 and increase in miR172 expression during development to the loss of rooting ability. Forty known and eight novel miRNAs were discovered and their predicted targets are listed. The expression pattern of nine miRNAs was determined during adventitious root formation in juvenile and mature cuttings. While the expression levels of miR156 and miR172 were inverse in juvenile and mature tissues, no mutual relationship was found between high miR156 expression and rooting ability, or high miR172 expression and loss of rooting ability. This is shown both in E. grandis and in E. brachyphylla, in which explants that underwent rejuvenation in tissue culture conditions were also examined.

Conclusions

It is suggested that in these Eucalyptus species, there is no correlation between the switch of miR156 with miR172 expression in the stems and the loss of rooting ability.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-524) contains supplementary material, which is available to authorized users.     

The Genomic Selfing Syndrome Accompanies the Evolutionary Breakdown of Heterostyly

The evolutionary transition from outcrossing to selfing can have important genomic consequences. Decreased effective population size and the reduced efficacy of selection are predicted to play an important role in the molecular evolution of the genomes of selfing species. We investigated evidence for molecular signatures of the genomic selfing syndrome using 66 species of Primula including distylous (outcrossing) and derived homostylous (selfing) taxa. We complemented our comparative analysis with a microevolutionary study of P. chungensis, which is polymorphic for mating system and consists of both distylous and homostylous populations. We generated chloroplast and nuclear genomic data sets for distylous, homostylous, and distylous–homostylous species and identified patterns of nonsynonymous to synonymous divergence (d_N/d_S) and polymorphism (π_N/π_S) in species or lineages with contrasting mating systems. Our analysis of coding sequence divergence and polymorphism detected strongly reduced genetic diversity and heterozygosity, decreased efficacy of purifying selection, purging of large-effect deleterious mutations, and lower rates of adaptive evolution in samples from homostylous compared with distylous populations, consistent with theoretical expectations of the genomic selfing syndrome. Our results demonstrate that self-fertilization is a major driver of molecular evolutionary processes with genomic signatures of selfing evident in both old and relatively young homostylous populations.     

Adaptive Evolution of the First Extra Exon in the Murid Rodent Prolactin Gene Family

The prolactin gene family in rodents consists of multiple members that coordinate the processes of reproduction and pregnancy. Some members of this family acquired one or two additional exons between exon 2 and exon 3 of the prototypical 5-exon, 4-intron structure, but the evolutionary importance of this insertion is unclear. Here, we focus on those members and survey this question by molecular evolutionary methods. Phylogenetic analysis shows that those members cluster into two distinct groups. Further analysis shows that the two groups of genes originated before the divergence of mouse and rat but after that of rodents from other mammals. We compared the d _N/d _S values for each branch of the gene tree but found no evidence to support positive selection for any branch. We found strong evidence, however, that one site (11E) of the 13 sites of the first extra exon underwent positive selection by the site-specific models of the maximum-likelihood method. Combining our molecular evolutionary analysis with other known functional evidence, we believe that the insertion of the extra exon implies some functional adaptation.     

Tempo of Degeneration Across Independently Evolved Nonrecombining Regions

Recombination is beneficial over the long term, allowing more effective selection. Despite long-term advantages of recombination, local recombination suppression can evolve and lead to genomic degeneration, in particular on sex chromosomes. Here, we investigated the tempo of degeneration in nonrecombining regions, that is, the function curve for the accumulation of deleterious mutations over time, leveraging on 22 independent events of recombination suppression identified on mating-type chromosomes of anther-smut fungi, including newly identified ones. Using previously available and newly generated high-quality genome assemblies of alternative mating types of 13 Microbotryum species, we estimated degeneration levels in terms of accumulation of nonoptimal codons and nonsynonymous substitutions in nonrecombining regions. We found a reduced frequency of optimal codons in the nonrecombining regions compared with autosomes, that was not due to less frequent GC-biased gene conversion or lower ancestral expression levels compared with recombining regions. The frequency of optimal codons rapidly decreased following recombination suppression and reached an asymptote after ca. 3 Ma. The strength of purifying selection remained virtually constant at d_N/d_S = 0.55, that is, at an intermediate level between purifying selection and neutral evolution. Accordingly, nonsynonymous differences between mating-type chromosomes increased linearly with stratum age, at a rate of 0.015 per My. We thus develop a method for disentangling effects of reduced selection efficacy from GC-biased gene conversion in the evolution of codon usage and we quantify the tempo of degeneration in nonrecombining regions, which is important for our knowledge on genomic evolution and on the maintenance of regions without recombination.     

Unraveling the Distribution and Evolution of miR156 targeted SPLs in Plants by Phylogenetic Analysis

Squamosa promoter binding protein like genes (SPLs) are critical during plant development and mostly regulated by miR156. However, little is known about phylogenetic distribution and evolutionary patterns of miR156 targeted SPLs. In this study, 183 SPLs from nine genome sequenced species representing algae, bryophytes, lycophyte, monocots, and eudicots were computationally analyzed. Our results showed that miR156 responsive elements (MREs) on SPLs were present in land plants but absent from unicellular green algae. Phylogenetic analysis revealed that miR156 targeted SPLs only distributed in group II not group I of land plants, suggesting they originated from a common ancestor. In addition, group II were further divided into seven subgroups (IIa IIg) and miR156 targeted SPLs distributed in some specific members of SPLs from six subgroups except subgroup IId. Such distribution pattern was well elucidated by gene structure evolution of miR156 targeted SPLs based on the correlation of phylogenetic classification and gene structure. They could suffer from the exon loss events combined with MREs loss during evolution. Moreover, gene duplication contributed to the abundance of miR156 targeted SPLs, which had significantly increased after angiosperms and lower plants split. With Arabidopsis as the model species, we found segmental and tandem gene duplications predominated during miR156 targeted SPLs expansion. Taken together, these results provide better insights in understanding the function diversity and evolution of miR156 targeted SPLs in plants.     

miR156靶基因SPL在植物中的系统分布和进化模式分析

Squamosa promoterbinding proteinlike genes （SPLs）在植物发育过程中具有重要作用。很多SPLs被miR156调节,然而,对于它们在植物中的系统分布和进化模式还知之甚少。本文对9个测序物种（藻类,苔藓,石松,单子叶和双子叶植物）的183个SPLs进行了生物信息学分析。结果表明miR156应答元件（MREs）仅在陆生植物SPLs中发现,藻类中不存在。系统进化分析显示陆生植物SPLs分为两大分支：group I和group II。 MiR156靶基因仅分布于group II,表明它们有着共同的祖先。Group II进一步分为7个亚支（IIaIIg）,miR156靶基因分布在除IId外的其余6个亚支的特定SPLs。系统分类与基因结构的相关性反映了SPL靶基因结构上的变化。在进化过程中,它们可能发生外显子的丢失且伴随MRE的丢失。另外,基因重复对SPL靶基因的丰度变化影响很大,尤其是被子植物与低等植物分歧后它们数量明显增加。以拟南芥为模式植物分析发现串联重复和片段重复是SPL靶基因扩张的主要机制。     

The PI/GLO-like locus in orchids: duplication and purifying selection at synonymous sites within Orchidinae (Orchidaceae)

Positive selection and relaxation of purifying constraints after duplication events have driven the functional diversification of gene families involved in development. One example of this occurred within the plant MADS-box genes. The evolution of the orchid flower was driven by duplication events followed by sub- and neo-functionalization of class B DEF-like MADS-box genes, which are present at three to four copies in the orchid genome. In contrast, the orchid PI/GLO-like class B MADS-box genes have been reported thus far as single-copy loci, with the only exception of Habenaria radiata.We isolated a novel PI/GLO-like gene (OrcPI2) in Orchis italica, which is different than the previously characterized OrcPI locus. The presence of two functional paralogs of PI/GLO-like genes in orchids is detectable only within the tribe Orchidinae. Evolutionary analyses revealed an apparent relaxation of purifying selection acting on the two PI/GLO-like paralogs of the Orchidinae when compared to the single-copy PI/GLO-like genes found in other orchid species. Furthermore, by measuring dN/dS (ω) ratios, we show that a high percentage of sites between the two PI/GLO-like paralogs have different evolutionary pressures. Interestingly, the apparent relaxation of selective constraints on the two PI/GLO-like paralogs is due to strong purifying selection at synonymous sites rather than to a high value of nonsynonymous substitution rate. This peculiar evolutionary pattern might be related to molecular processes such as mRNA folding and/or translational efficiency control. These processes could potentially be involved in or predate the functional diversification of the two PI/GLO-like paralogs within Orchidinae.     

Molecular Phylogenies and Evolution of crt Genes in Algae

ABSTRACT

The carotenoids constitute the most widespread class of pigments in nature. Most previous work has concentrated on the identification and characterization of their chemical physical properties and bioavailability. In recent years, significant amounts of research have been conducted in an attempt to analyze the genes and the molecular regulation of the genes involved in the biosynthesis of carotenoids. However, it is important not to lose sight of the early evolution of carotenoid biosynthesis. One of the major obstacles in understanding the evolution of the respective enzymes and their patterns of selection is a lack of a well-supported phylogenic analysis. In the present research, a major long-term objective was to provide a clearer picture of the evolutionary history of genes, together with an evaluation of the patterns of selection in algae. These phylogenies will be important in studies characterizing the evolution of algae. The gene sequences of the enzymes involved in the major steps of the carotenoid biosynthetic pathway in algae (cyanobacteria, rhofophyta, chlorophyta) have been analyzed. Phylogenetic relationships among protein-coding DNA sequences were reconstructed by neighbor-joining (NJ) analysis for the respective carotenoid biosynthetic pathway genes (crt) in algae. The analysis also contains an estimation of the rate of nonsynonymous nucleotide substitutions per nonsynonymous site (d_N), synonymous nucleotide substitution per synonymous site (d_S), and the ratio of nonsynonmous (d_N/d_S) for the test of selection patterns. The phylogenetic trees show that the taxa of some genera have a closer evolutionary relationship with other genera in some gene sequences, which suggests a common ancient origin and that lateral gene transfer has occurred among unrelated genera. The d_N values of crt genes in the early pathway are relatively low, while those of the following steps are slightly higher, while the d_N values of crt genes in chlorophyta are higher than those in cyanobacteria. Most of the d_N/d_S values exceed 1. The phylogenetic analysis revealed that lateral gene transfer may have taken place across algal genomes and the d_N values suggest that most of the early crt genes are well conserved compared to the later crt genes. Furthermore, d_N values also revealed that the crt genes of chlorophyta are more evolutionary than cyanobacteria. The amino acids' changes are mostly adaptive evolution under the influence of positive diversity selection.     

Expansion and Accelerated Evolution of 9-Exon Odorant Receptors in Polistes Paper Wasps

Independent origins of sociality in bees and ants are associated with independent expansions of particular odorant receptor (OR) gene subfamilies. In ants, one clade within the OR gene family, the 9-exon subfamily, has dramatically expanded. These receptors detect cuticular hydrocarbons (CHCs), key social signaling molecules in insects. It is unclear to what extent 9-exon OR subfamily expansion is associated with the independent evolution of sociality across Hymenoptera, warranting studies of taxa with independently derived social behavior. Here, we describe OR gene family evolution in the northern paper wasp, Polistes fuscatus, and compare it to four additional paper wasp species spanning ∼40 million years of evolutionary divergence. We find 200 putatively functional OR genes in P. fuscatus, matching predictions from neuroanatomy, and more than half of these are in the 9-exon subfamily. Most OR gene expansions are tandemly arrayed at orthologous loci in Polistes genomes, and microsynteny analysis shows species-specific gain and loss of 9-exon ORs within tandem arrays. There is evidence of episodic positive diversifying selection shaping ORs in expanded subfamilies. Values of omega (d_N/d_S) are higher among 9-exon ORs compared to other OR subfamilies. Within the Polistes OR gene tree, branches in the 9-exon OR clade experience relaxed negative (relaxed purifying) selection relative to other branches in the tree. Patterns of OR evolution within Polistes are consistent with 9-exon OR function in CHC perception by combinatorial coding, with both natural selection and neutral drift contributing to interspecies differences in gene copy number and sequence.     

Kr/Kc but not dN/dS correlates positively with body mass in birds,raising implications for inferring lineage-specific selection

Background

The ratio of the rates of non-synonymous and synonymous substitution (d_N/d_S) is commonly used to estimate selection in coding sequences. It is often suggested that, all else being equal, d_N/d_S should be lower in populations with large effective size (N_e) due to increased efficacy of purifying selection. As N_e is difficult to measure directly, life history traits such as body mass, which is typically negatively associated with population size, have commonly been used as proxies in empirical tests of this hypothesis. However, evidence of whether the expected positive correlation between body mass and d_N/d_S is consistently observed is conflicting.

Results

Employing whole genome sequence data from 48 avian species, we assess the relationship between rates of molecular evolution and life history in birds. We find a negative correlation between d_N/d_S and body mass, contrary to nearly neutral expectation. This raises the question whether the correlation might be a method artefact. We therefore in turn consider non-stationary base composition, divergence time and saturation as possible explanations, but find no clear patterns. However, in striking contrast to d_N/d_S, the ratio of radical to conservative amino acid substitutions (K_r/K_c) correlates positively with body mass.

Conclusions

Our results in principle accord with the notion that non-synonymous substitutions causing radical amino acid changes are more efficiently removed by selection in large populations, consistent with nearly neutral theory. These findings have implications for the use of d_N/d_S and suggest that caution is warranted when drawing conclusions about lineage-specific modes of protein evolution using this metric.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0542-8) contains supplementary material, which is available to authorized users.