Concerted evolution: molecular mechanism and biological implications

Phylogenomic analysis demonstrates a pattern of rare and long-lasting concerted evolution in prokaryotes

Concerted evolution, where paralogs in the same species show higher sequence similarity to each other than to orthologs in other species, is widely found in many species. However, cases of concerted evolution that last for hundreds of millions of years are very rare. By genome-wide analysis of a broad selection of prokaryotes, we provide strong evidence of recurrent concerted evolution in 26 genes, most of which have lasted more than ~500 million years. We find that most concertedly evolving genes are key members of important pathways, and encode proteins from the same complexes and/or pathways, suggesting coevolution of genes via concerted evolution to maintain gene balance. We also present LRCE-DB, a comprehensive online repository of long-lasting concerted evolution. Collectively, our study reveals that although most duplicated genes may diverge in sequence over a long period, on rare occasions this constraint can be breached, leading to unexpected long-lasting concerted evolution in a recurrent manner.

Sishuo Wang and Youhua Chen present an analysis of concerted evolution in prokaryotes using a new computational pipeline, iSeeCE. They find evidence in 26 genes for recurrent concerted evolution, most of which last more than ~500 million years, and provide a database, LRCE-DB, for data exploration.

Phylogenetic Network Analysis Revealed the Occurrence of Horizontal Gene Transfer of 16S rRNA in the Genus Enterobacter

Horizontal gene transfer (HGT) is a ubiquitous genetic event in bacterial evolution, but it seldom occurs for genes involved in highly complex supramolecules (or biosystems), which consist of many gene products. The ribosome is one such supramolecule, but several bacteria harbor dissimilar and/or chimeric 16S rRNAs in their genomes, suggesting the occurrence of HGT of this gene. However, we know little about whether the genes actually experience HGT and, if so, the frequency of such a transfer. This is primarily because the methods currently employed for phylogenetic analysis (e.g., neighbor-joining, maximum likelihood, and maximum parsimony) of 16S rRNA genes assume point mutation-driven tree-shape evolution as an evolutionary model, which is intrinsically inappropriate to decipher the evolutionary history for genes driven by recombination. To address this issue, we applied a phylogenetic network analysis, which has been used previously for detection of genetic recombination in homologous alleles, to the 16S rRNA gene. We focused on the genus Enterobacter, whose phylogenetic relationships inferred by multi-locus sequence alignment analysis and 16S rRNA sequences are incompatible. All 10 complete genomic sequences were retrieved from the NCBI database, in which 71 16S rRNA genes were included. Neighbor-joining analysis demonstrated that the genes residing in the same genomes clustered, indicating the occurrence of intragenomic recombination. However, as suggested by the low bootstrap values, evolutionary relationships between the clusters were uncertain. We then applied phylogenetic network analysis to representative sequences from each cluster. We found three ancestral 16S rRNA groups; the others were likely created through recursive recombination between the ancestors and chimeric descendants. Despite the large sequence changes caused by the recombination events, the RNA secondary structures were conserved. Successive intergenomic and intragenomic recombination thus shaped the evolution of 16S rRNA genes in the genus Enterobacter.

An Unusual Accumulation of Ribosomal Multigene Families and Microsatellite DNAs in the XX/XY Sex Chromosome System in the Trans-Andean Catfish Pimelodella cf. chagresi (Siluriformes:Heptapteridae)

This work constitutes the first cytogenetic characterization of a trans-Andean species of Heptapteridae. The catfish Pimelodella cf. chagresi from the Upper Rio Magdalena was studied, applying standard cytogenetic techniques (Giemsa, C-banding, and argyrophilic nucleolar organizer region [Ag-NOR]) and fluorescence in situ hybridization techniques using repetitive DNA probes: microsatellites (CA₁₅ and GA₁₅) and ribosomal RNA (rRNA) multigene families (18S and 5S recombinant DNA [rDNA] probes). The species showed a unique diploid chromosome number 2n = 50 (32m [metacentrics] +14sm [submetacentrics] +4st [subtelocentrics]) and a XX/XY sex chromosomal system, where the heteromorphic Y-chromosome revealed a conspicuous accumulation of all the assayed domains of repetitive DNA. P. cf. chagresi karyotype shares common features with other Heptapteridae, such as the predominance of metacentric and submetacentric chromosomes, and one pair of subtelomeric nucleolar organizer regions (NORs). These results reflect an independent karyological identity of a trans-Andean species and the relevance of repetitive DNA sequences in the process of sex chromosome differentiation in fish; it is the first case of syntenic accumulation of rRNA multigene families (18S and 5S rDNA) and microsatellite sequences (CA₁₅ and GA₁₅) in a differentiated sex chromosome in Neotropical fish.

SMORE: Synteny Modulator of Repetitive Elements

Several families of multicopy genes, such as transfer ribonucleic acids (tRNAs) and ribosomal RNAs (rRNAs), are subject to concerted evolution, an effect that keeps sequences of paralogous genes effectively identical. Under these circumstances, it is impossible to distinguish orthologs from paralogs on the basis of sequence similarity alone. Synteny, the preservation of relative genomic locations, however, also remains informative for the disambiguation of evolutionary relationships in this situation. In this contribution, we describe an automatic pipeline for the evolutionary analysis of such cases that use genome-wide alignments as a starting point to assign orthology relationships determined by synteny. The evolution of tRNAs in primates as well as the history of the Y RNA family in vertebrates and nematodes are used to showcase the method. The pipeline is freely available.

Peterson, Soreng RJ, Paszko B. Molecular phylogenetics of cool-season grasses in the subtribes Agrostidinae, Anthoxanthinae, Aveninae, Brizinae, Calothecinae, Koeleriinae and Phalaridinae (Poaceae, Pooideae, Poeae, Poeae chloroplast group 1)

Circumscriptions of and relationships among many genera and suprageneric taxa of the diverse grass tribe Poeae remain controversial. In an attempt to clarify these, we conducted phylogenetic analyses of >2400 new DNA sequences from two nuclear ribosomal regions (ITS, including internal transcribed spacers 1 and 2 and the 5.8S gene, and the 3'-end of the external transcribed spacer (ETS)) and five plastid regions (matK, trnL-trnF, atpF-atpH, psbK-psbI, psbA-rps19-trnH), and of more than 1000 new and previously published ITS sequences, focused particularly on Poeae chloroplast group 1 and including broad and increased species sampling compared to previous studies. Deep branches in the combined plastid and combined ITS+ETS trees are generally well resolved, the trees are congruent in most aspects, branch support across the trees is stronger than in trees based on only ITS and fewer plastid regions, and there is evidence of conflict between data partitions in some taxa. In plastid trees, a strongly supported clade corresponds to Poeae chloroplast group 1 and includes Agrostidinae p.p., Anthoxanthinae, Aveninae s.str., Brizinae, Koeleriinae (sometimes included in Aveninae s.l.), Phalaridinae and Torreyochloinae. In the ITS+ETS tree, a supported clade includes these same tribes as well as Sesleriinae and Scolochloinae. Aveninae s.str. and Sesleriinae are sister taxa and form a clade with Koeleriinae in the ITS+ETS tree whereas Aveninae s.str. and Koeleriinae form a clade and Sesleriinae is part of Poeae chloroplast group 2 in the plastid tree. All species of Trisetum are part of Koeleriinae, but the genus is polyphyletic. Koeleriinae is divided into two major subclades: one comprises Avellinia, Gaudinia, Koeleria, Rostraria, Trisetaria and Trisetum subg. Trisetum, and the other Calamagrostis/Deyeuxia p.p. (multiple species from Mexico to South America), Peyritschia, Leptophyllochloa, Sphenopholis, Trisetopsis and Trisetum subg. Deschampsioidea. Graphephorum, Trisetum cernuum, T. irazuense and T. macbridei fall in different clades of Koeleriinae in plastid vs. nuclear ribosomal trees, and are likely of hybrid origin. ITS and matK trees identify a third lineage of Koeleriinae corresponding to Trisetum subsect. Sibirica, and affinities of Lagurus ovatus with respect to Aveninae s.str. and Koeleriinae are incongruent in nuclear ribosomal and plastid trees, supporting recognition of Lagurus in its own subtribe. A large clade comprises taxa of Agrostidinae, Brizinae and Calothecinae, but neither Agrostidinae nor Calothecinae are monophyletic as currently circumscribed and affinities of Brizinae differ in plastid and nuclear ribosomal trees. Within this clade, one newly identified lineage comprises Calamagrostis coarctata, Dichelachne, Echinopogon (Agrostidinae p.p.) and Relchela (Calothecinae p.p.), and another comprises Chascolytrum (Calothecinae p.p.) and Deyeuxia effusa (Agrostidinae p.p.). Within Agrostidinae p.p., the type species of Deyeuxia and Calamagrostis s.str. are closely related, supporting classification of Deyeuxia as a synonym of Calamagrostis s.str. Furthermore, the two species of Ammophila are not sister taxa and are nested among different groups of Calamagrostis s.str., supporting their classification in Calamagrostis. Agrostis, Lachnagrostis and Polypogon form a clade and species of each are variously intermixed in plastid and nuclear ribosomal trees. Additionally, all but one species from South America classified in Deyeuxia sect. Stylagrostis resolve in Holcinae p.p. (Deschampsia). The current phylogenetic results support recognition of the latter species in Deschampsia, and we also demonstrate Scribneria is part of this clade. Moreover, Holcinae is not monophyletic in its current circumscription because Deschampsia does not form a clade with Holcus and Vahlodea, which are sister taxa. The results support recognition of Deschampsia in its own subtribe Aristaveninae. Substantial further changes to the classification of these grasses will be needed to produce generic circumscriptions consistent with phylogenetic evidence. The following 15 new combinations are made: Calamagrostis × calammophila, C. breviligulata, C. breviligulata subsp. champlainensis, C. × don-hensonii, Deschampsia aurea, D. bolanderi, D. chrysantha, D. chrysantha var. phalaroides, D. eminens, D. eminens var. fulva, D. eminens var. inclusa, D. hackelii, D. ovata, and D. ovata var. nivalis. D. podophora; the new name Deschampsia parodiana is proposed; the new subtribe Lagurinae is described; and a second-step lectotype is designated for the name Deyeuxia phalaroides.

Unexpected High Intragenomic Variation in Two of Three Major Pest Thrips Species Does Not Affect Ribosomal Internal Transcribed Spacer 2 (ITS2) Utility for Thrips Identification

The mitochondrial cytochrome oxidase I gene (mtCO1) and the ribosomal internal transcribed spacer 2 region (ITS2) are among the most widely used molecular markers for insect taxonomic characterization. Three economically important species of thrips, Scirtothripsdorsalis, Thripspalmi, and Frankliniellaoccidentalis were selected to examine the extent of intragenomic variation within these two marker regions in the family Thripidae, and determine if this variation would affect the utility of markers in thrips molecular diagnostics. For each species, intragenomic (within individual) variation and intergenomic (among individuals) variation was assessed by cloning and sequencing PCR-amplified copies. Intergenomic variation was generally higher than intragenomic variation except in cases where intergenomic variation was very low, as in mtCO1 from S.dorsalis and F.occidentalis. Intragenomic variation was detected in both markers in all three of the thrips species, however, 2-3 times more intragenomic variation was observed for ITS2 than mtCO1 in both S.dorsalis and T.palmi. Furthermore, levels of intragenomic variation were low for both of the genes in F.occidentalis. In all of the three thrips species, no sex-based clustering of haplotypes was observed in either marker. Unexpected high intragenomic variation in ITS2 for two of three thrips species did not interfere with thrips diagnostics. However, caution should be taken in applying ITS2 to certain studies of S.dorsalis and T.palmi when high levels of intragenomic variation could be problematic or confounding. In such studies, mtCO1 may be a preferable marker. Possible reasons for discrepancies in intragenomic variation among genomic regions are discussed.

Gene Conversion Facilitates Adaptive Evolution on Rugged Fitness Landscapes

Gene conversion is a ubiquitous phenomenon that leads to the exchange of genetic information between homologous DNA regions and maintains coevolving multi-gene families in most prokaryotic and eukaryotic organisms. In this paper, we study its implications for the evolution of a single functional gene with a silenced duplicate, using two different models of evolution on rugged fitness landscapes. Our analytical and numerical results show that, by helping to circumvent valleys of low fitness, gene conversion with a passive duplicate gene can cause a significant speedup of adaptation, which depends nontrivially on the frequency of gene conversion and the structure of the landscape. We find that stochastic effects due to finite population sizes further increase the likelihood of exploiting this evolutionary pathway. A universal feature appearing in both deterministic and stochastic analysis of our models is the existence of an optimal gene conversion rate, which maximizes the speed of adaptation. Our results reveal the potential for duplicate genes to act as a "scratch paper" that frees evolution from being limited to strictly beneficial mutations in strongly selective environments.

Two new species of Ripella (Amoebozoa, Vannellida) and unusual intragenomic variability in the SSU rRNA gene of this genus

Two new species, Ripella decalvata and R. tribonemae (Amoebozoa, Vannellida), are described and the diversity of known strains assigned to the genus analyzed. Ripella spp. are closely similar to each other in the light microscopic characters and sequences of small-subunit (SSU) ribosomal RNA gene, but differences in the cell coat structure and cytochrome oxidase (COI) gene sequences are more prominent. SSU rRNA in R. platypodia CCAP1589/2, R. decalvata and R. tribonemae demonstrates an unusual pattern of intragenomic variation. Sequencing of multiple molecular clones of this gene produced numerous sequence variants in a number of specific sites. These sites were usually terminal parts of several variable helices in all studied strains. Analysis of all known Ripella strains shows that SSU rRNA sites differing between strains of different origin are mainly restricted to these areas of the gene. There are only two sites, which differ between strains, but not within genomes. This intragenomic variability of the SSU rRNA gene, seemingly characteristic of all Ripella spp., was never reported to be so extensive in Amoebozoa. The data obtained show another example of complex organization of rRNA gene cluster in protists and emphasize caution needed when interpreting the metagenomic data based on this marker.

Genomic fossils reveal adaptation of non-autonomous pararetroviruses driven by concerted evolution of noncoding regulatory sequences

The interplay of different virus species in a host cell after infection can affect the adaptation of each virus. Endogenous viral elements, such as endogenous pararetroviruses (PRVs), have arisen from vertical inheritance of viral sequences integrated into host germline genomes. As viral genomic fossils, these sequences can thus serve as valuable paleogenomic data to study the long-term evolutionary dynamics of virus-virus interactions, but they have rarely been applied for this purpose. All extant PRVs have been considered autonomous species in their parasitic life cycle in host cells. Here, we provide evidence for multiple non-autonomous PRV species with structural defects in viral activity that have frequently infected ancient grass hosts and adapted through interplay between viruses. Our paleogenomic analyses using endogenous PRVs in grass genomes revealed that these non-autonomous PRV species have participated in interplay with autonomous PRVs in a possible commensal partnership, or, alternatively, with one another in a possible mutualistic partnership. These partnerships, which have been established by the sharing of noncoding regulatory sequences (NRSs) in intergenic regions between two partner viruses, have been further maintained and altered by the sequence homogenization of NRSs between partners. Strikingly, we found that frequent region-specific recombination, rather than mutation selection, is the main causative mechanism of NRS homogenization. Our results, obtained from ancient DNA records of viruses, suggest that adaptation of PRVs has occurred by concerted evolution of NRSs between different virus species in the same host. Our findings further imply that evaluation of within-host NRS interactions within and between populations of viral pathogens may be important.

Intra- and inter-isolate variation of ribosomal and protein-coding genes in Pleurotus: implications for molecular identification and phylogeny on fungal groups

Background

The internal transcribed spacer (ITS), RNA polymerase II second largest subunit (RPB2), and elongation factor 1-alpha (EF1α) are often used in fungal taxonomy and phylogenetic analysis. As we know, an ideal molecular marker used in molecular identification and phylogenetic studies is homogeneous within species, and interspecific variation exceeds intraspecific variation. However, during our process of performing ITS, RPB2, and EF1α sequencing on the Pleurotus spp., we found that intra-isolate sequence polymorphism might be present in these genes because direct sequencing of PCR products failed in some isolates. Therefore, we detected intra- and inter-isolate variation of the three genes in Pleurotus by polymerase chain reaction amplification and cloning in this study.

Results

Results showed that intra-isolate variation of ITS was not uncommon but the polymorphic level in each isolate was relatively low in Pleurotus; intra-isolate variations of EF1α and RPB2 sequences were present in an unexpectedly high amount. The polymorphism level differed significantly between ITS, RPB2, and EF1α in the same individual, and the intra-isolate heterogeneity level of each gene varied between isolates within the same species. Intra-isolate and intraspecific variation of ITS in the tested isolates was less than interspecific variation, and intra-isolate and intraspecific variation of RPB2 was probably equal with interspecific divergence. Meanwhile, intra-isolate and intraspecific variation of EF1α could exceed interspecific divergence. These findings suggested that RPB2 and EF1α are not desirable barcoding candidates for Pleurotus. We also discussed the reason why rDNA and protein-coding genes showed variants within a single isolate in Pleurotus, but must be addressed in further research.

Conclusions

Our study demonstrated that intra-isolate variation of ribosomal and protein-coding genes are likely widespread in fungi. This has implications for studies on fungal evolution, taxonomy, phylogenetics, and population genetics. More extensive sampling of these genes and other candidates will be required to ensure reliability as phylogenetic markers and DNA barcodes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-017-1046-y) contains supplementary material, which is available to authorized users.

Cadherin genes and evolutionary novelties in the octopus

All animals with large brains must have molecular mechanisms to regulate neuronal process outgrowth and prevent neurite self-entanglement. In vertebrates, two major gene families implicated in these mechanisms are the clustered protocadherins and the atypical cadherins. However, the molecular mechanisms utilized in complex invertebrate brains, such as those of the cephalopods, remain largely unknown. Recently, we identified protocadherins and atypical cadherins in the octopus. The octopus protocadherin expansion shares features with the mammalian clustered protocadherins, including enrichment in neural tissues, clustered head-to-tail orientations in the genome, and a large first exon encoding all cadherin domains. Other octopus cadherins, including a newly-identified cadherin with 77 extracellular cadherin domains, are elevated in the suckers, a striking cephalopod novelty. Future study of these octopus genes may yield insights into the general functions of protocadherins in neural wiring and cadherin-related proteins in complex morphogenesis.

Development of a multilocus-based approach for sponge (phylum Porifera) identification: refinement and limitations

For sponges (phylum Porifera), there is no reliable molecular protocol available for species identification. To address this gap, we developed a multilocus-based Sponge Identification Protocol (SIP) validated by a sample of 37 sponge species belonging to 10 orders from South Australia. The universal barcode COI mtDNA, 28S rRNA gene (D3-D5), and the nuclear ITS1-5.8S-ITS2 region were evaluated for their suitability and capacity for sponge identification. The highest Bit Score was applied to infer the identity. The reliability of SIP was validated by phylogenetic analysis. The 28S rRNA gene and COI mtDNA performed better than the ITS region in classifying sponges at various taxonomic levels. A major limitation is that the databases are not well populated and possess low diversity, making it difficult to conduct the molecular identification protocol. The identification is also impacted by the accuracy of the morphological classification of the sponges whose sequences have been submitted to the database. Re-examination of the morphological identification further demonstrated and improved the reliability of sponge identification by SIP. Integrated with morphological identification, the multilocus-based SIP offers an improved protocol for more reliable and effective sponge identification, by coupling the accuracy of different DNA markers.

Evolution of Fungal U3 snoRNAs: Structural Variation and Introns

The U3 small nucleolar RNA (snoRNA) is an essential player in the initial steps of ribosomal RNA biogenesis which is ubiquitously present in Eukarya. It is exceptional among the small nucleolar RNAs in its size, the presence of multiple conserved sequence boxes, a highly conserved secondary structure core, its biogenesis as an independent gene transcribed by polymerase III, and its involvement in pre-rRNA cleavage rather than chemical modification. Fungal U3 snoRNAs share many features with their sisters from other eukaryotic kingdoms but differ from them in particular in their 5’ regions, which in fungi has a distinctive consensus structure and often harbours introns. Here we report on a comprehensive homology search and detailed analysis of the evolution of sequence and secondary structure features covering the entire kingdom Fungi.

ITS non-concerted evolution and rampant hybridization in the legume genus Lespedeza (Fabaceae)

The internal transcribed spacer (ITS) as one part of nuclear ribosomal DNA is one of the most extensively sequenced molecular markers in plant systematics. The ITS repeats generally exhibit high-level within-individual homogeneity, while relatively small-scale polymorphism of ITS copies within individuals has often been reported in literature. Here, we identified large-scale polymorphism of ITS copies within individuals in the legume genus Lespedeza (Fabaceae). Divergent paralogs of ITS sequences, including putative pseudogenes, recombinants, and multiple functional ITS copies were sometimes detected in the same individual. Thirty-seven ITS pseudogenes could be easily detected according to nucleotide changes in conserved 5.8S motives, the significantly lower GC contents in at least one of three regions, and the lost ability of 5.8S rDNA sequence to fold into a conserved secondary structure. The distribution patterns of the putative functional clones were highly different between the traditionally recognized two subgenera, suggesting different rates of concerted evolution in two subgenera which could be attributable to their different extents/frequencies of hybridization, confirmed by our analysis of the single-copy nuclear gene PGK. These findings have significant implications in using ITS marker for reconstructing phylogeny and studying hybridization.

Selection Signatures in the First Exon of Paralogous Receptor Kinase Genes from the Sym2 Region of the Pisum sativum L. Genome

During the initial step of the symbiosis between legumes (Fabaceae) and nitrogen-fixing bacteria (rhizobia), the bacterial signal molecule known as the Nod factor (nodulation factor) is recognized by plant LysM motif-containing receptor-like kinases (LysM-RLKs). The fifth chromosome of barrel medic (Medicago truncatula Gaertn.) contains a cluster of paralogous LysM-RLK genes, one of which is known to participate in symbiosis. In the syntenic region of the pea (Pisum sativum L.) genome, three genes have been identified: PsK1 and PsSym37, two symbiosis-related LysM-RLK genes with known sequences, and the unsequenced PsSym2 gene which presumably encodes a LysM-RLK and is associated with increased selectivity to certain Nod factors. In this work, we identified a new gene encoding a LysM-RLK, designated as PsLykX, within the Sym2 genomic region. We sequenced the first exons (corresponding to the protein receptor domain) of PsSym37, PsK1, and PsLykX from a large set of pea genotypes of diverse origin. The nucleotide diversity of these fragments was estimated and groups of haplotypes for each gene were revealed. Footprints of selection pressure were detected via comparative analyses of SNP distribution across the first exons of these genes and their homologs MtLYK2, MtLYK3, and MtLYK4 from M. truncatula retrieved from the Medicago Hapmap project. Despite the remarkable similarity among all the studied genes, they exhibited contrasting selection signatures, possibly pointing to diversification of their functions. Signatures of balancing selection were found in LysM1-encoding parts of PsSym37 and PsK1, suggesting that the diversity of these parts may be important for pea LysM-RLKs. The first exons of PsSym37 and PsK1 displayed signatures of purifying selection, as well as MtLYK2 of M. truncatula. Evidence of positive selection affecting primarily LysM domains was found in all three investigated M. truncatula genes, as well as in the pea gene PsLykX. The data suggested that PsLykX is a promising candidate for PsSym2, which has remained elusive for more than 30 years.

The D3-D5 region of large subunit ribosomal DNA provides good resolution of German limnic and terrestrial nematode communities

Reliable and well-developed DNA barcode databases are indispensable for the identification of microscopic life. However, effectiveness of molecular barcoding in identifying terrestrial specimens, and nematodes in particular, has received little attention. In this study, ca 600 ribosomal large subunit DNA fragments (D3-D5 region) were successfully amplified for 79 limnic and terrestrial nematode species sampled at 147 locations across Germany. Distinctive DNA motifs in the LSU region were identified in 80% of all species examined. For 13 supposedly single morphospecies, 2-7 LSU barcode groups were detected with a wide range of intraspecific variations (0.09-7.9%). This region seems to be more suitable for the assessment of limno-terrestrial nematode diversity than the frequently used mitochondrial gene COI, as amplification success of the latter fragment is low for several nematode species. Our reference database for nematodes may serve as a starting point for applied and fundamental studies for these ubiquitous, ecologically highly relevant, organisms.

Geographic Population Structure in Epstein-Barr Virus Revealed by Comparative Genomics

Epstein-Barr virus (EBV) latently infects the majority of the human population and is implicated as a causal or contributory factor in numerous diseases. We sequenced 27 complete EBV genomes from a cohort of Multiple Sclerosis (MS) patients and healthy controls from Italy, although no variants showed a statistically significant association with MS. Taking advantage of the availability of ∼130 EBV genomes with known geographical origins, we reveal a striking geographic distribution of EBV sub-populations with distinct allele frequency distributions. We discuss mechanisms that potentially explain these observations, and their implications for understanding the association of EBV with human disease.

Evolution of Protein Domain Repeats in Metazoa

Repeats are ubiquitous elements of proteins and they play important roles for cellular function and during evolution. Repeats are, however, also notoriously difficult to capture computationally and large scale studies so far had difficulties in linking genetic causes, structural properties and evolutionary trajectories of protein repeats. Here we apply recently developed methods for repeat detection and analysis to a large dataset comprising over hundred metazoan genomes. We find that repeats in larger protein families experience generally very few insertions or deletions (indels) of repeat units but there is also a significant fraction of noteworthy volatile outliers with very high indel rates. Analysis of structural data indicates that repeats with an open structure and independently folding units are more volatile and more likely to be intrinsically disordered. Such disordered repeats are also significantly enriched in sites with a high functional potential such as linear motifs. Furthermore, the most volatile repeats have a high sequence similarity between their units. Since many volatile repeats also show signs of recombination, we conclude they are often shaped by concerted evolution. Intriguingly, many of these conserved yet volatile repeats are involved in host-pathogen interactions where they might foster fast but subtle adaptation in biological arms races. KEY WORDS: protein evolution, domain rearrangements, protein repeats, concerted evolution.

Orthologs, turn-over, and remolding of tRNAs in primates and fruit flies

Background

Transfer RNAs (tRNAs) are ubiquitous in all living organism. They implement the genetic code so that most genomes contain distinct tRNAs for almost all 61 codons. They behave similar to mobile elements and proliferate in genomes spawning both local and non-local copies. Most tRNA families are therefore typically present as multicopy genes. The members of the individual tRNA families evolve under concerted or rapid birth-death evolution, so that paralogous copies maintain almost identical sequences over long evolutionary time-scales. To a good approximation these are functionally equivalent. Individual tRNA copies thus are evolutionary unstable and easily turn into pseudogenes and disappear. This leads to a rapid turnover of tRNAs and often large differences in the tRNA complements of closely related species. Since tRNA paralogs are not distinguished by sequence, common methods cannot not be used to establish orthology between tRNA genes.

Results

In this contribution we introduce a general framework to distinguish orthologs and paralogs in gene families that are subject to concerted evolution. It is based on the use of uniquely aligned adjacent sequence elements as anchors to establish syntenic conservation of sequence intervals. In practice, anchors and intervals can be extracted from genome-wide multiple sequence alignments. Syntenic clusters of concertedly evolving genes of different families can then be subdivided by list alignments, leading to usually small clusters of candidate co-orthologs. On the basis of recent advances in phylogenetic combinatorics, these candidate clusters can be further processed by cograph editing to recover their duplication histories. We developed a workflow that can be conceptualized as stepwise refinement of a graph of homologous genes. We apply this analysis strategy with different types of synteny anchors to investigate the evolution of tRNAs in primates and fruit flies. We identified a large number of tRNA remolding events concentrated at the tips of the phylogeny. With one notable exception all phylogenetically old tRNA remoldings do not change the isoacceptor class.

Conclusions

Gene families evolving under concerted evolution are not amenable to classical phylogenetic analyses since paralogs maintain identical, species-specific sequences, precluding the estimation of correct gene trees from sequence differences. This leaves conservation of syntenic arrangements with respect to “anchor elements” that are not subject to concerted evolution as the only viable source of phylogenetic information. We have demonstrated here that a purely synteny-based analysis of tRNA gene histories is indeed feasible. Although the choice of synteny anchors influences the resolution in particular when tight gene clusters are present, and the quality of sequence alignments, genome assemblies, and genome rearrangements limits the scope of the analysis, largely coherent results can be obtained for tRNAs. In particular, we conclude that a large fraction of the tRNAs are recent copies. This proliferation is compensated by rapid pseudogenization as exemplified by many very recent alloacceptor remoldings.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2927-4) contains supplementary material, which is available to authorized users.

Intragenomic polymorphism and intergenomic recombination in the ribosomal RNA genes of strains belonging to a yeast species Pichia membranifaciens

A concerted evolution model has been proposed to explain the observed lack of sequence variation among the multiple ribosomal RNA (rRNA) gene copies in many different eukaryotic species. Recent studies on the level of intragenomic variations in the rRNA gene repeats of fungi resulted in controversial conclusions. In this study, we clearly showed that significant polymorphisms of the internal transcribed spacers (ITS1 and ITS2) of ribosomal DNA (rDNA) exist within the genome of a strain of the yeast species Pichia membranifaciens. More interestingly, we showed that the intragenomic ITS sequence polymorphisms were formed by intergenomic rDNA recombination among different P. membranifaciens strains with significantly different ITS sequences. Intergenomic rDNA recombination was also responsible for the diversification of rDNA sequences in different strains of the species. After the events bring together different rDNA types in individual genomes of the P. membranifaciens strains compared, rDNA sequence heterogeneity has remained in the genome of one but eliminated by homogenisation in the genomes of other strains. Our findings show new clue for further investigation on the mechanism of concerted evolution of rRNA genes in eukaryotes.

Contrasting soil bacterial community structure between the phyla Acidobacteria and Proteobacteria in tropical Southeast Asian and temperate Japanese forests

Soil bacterial community structures of six dominant phyla (Acidobacteria, Proteobacteria, Verrucomicrobia, Planctomycetes, Bacteroidetes and Actinobacteria) and unclassified bacteria detected in tropical Sarawakian and temperate Japanese forests were compared based on 16S rRNA gene sequence variation. The class composition in each phylum was similar among the studied forests; however, significant heterogeneities of class frequencies were detected. Acidobacteria and Proteobacteria were the most dominant phyla in all six forests, but differed in the level of bacterial species diversity, pattern of species occurrence and association pattern of species composition with physicochemical properties in soil. Species diversity among Acidobacteria was approximately half that among Proteobacteria, based on the number of clusters and the Chao1 index, even though a similar number of sequence reads were obtained for these two phyla. In contrast, species diversity within Planctomycetes and Bacteroidetes was nearly as high as within Acidobacteria, despite many fewer sequence reads. The density of species (the number of sequence reads per cluster) correlated negatively with species diversity, and species density within Acidobacteria was approximately twice that within Proteobacteria. Although the percentage of forest-specific species was high for all bacterial groups, sampling site-specific species varied among bacterial groups, indicating limited inter-forest migration and differential movement of bacteria in forest soil. For five of the seven bacterial groups, including Acidobacteria, soil pH appeared to strongly influence species composition, but this association was not observed for Proteobacterial species. Topology of UPGMA trees and pattern of NMDS plots among the forests differed among the bacterial groups, suggesting that each bacterial group has adapted and evolved independently in each forest.

The confounding complexity of innate immune receptors within and between teleost species

Teleost genomes encode multiple multigene families of immunoglobulin domain-containing innate immune receptors (IIIRs) with unknown function and no clear mammalian orthologs. However, the genomic organization of IIIR gene clusters and the structure and signaling motifs of the proteins they encode are similar to those of mammalian innate immune receptor families such as the killer cell immunoglobulin-like receptors (KIRs), leukocyte immunoglobulin-like receptors (LILRs), Fc receptors, triggering receptors expressed on myeloid cells (TREMs) and CD300s. Teleost IIIRs include novel immune-type receptors (NITRs); diverse immunoglobulin domain containing proteins (DICPs); polymeric immunoglobulin receptor-like proteins (PIGRLs); novel immunoglobulin-like transcripts (NILTs) and leukocyte immune-type receptors (LITRs). The accumulation of genomic sequence data has revealed that IIIR gene clusters in zebrafish display haplotypic and gene content variation. This intraspecific genetic variation, as well as significant interspecific variation, frequently confounds the identification of definitive orthologous IIIR sequences between teleost species. Nevertheless, by defining which teleost lineages encode (and do not encode) different IIIR families, predictions can be made about the presence (or absence) of specific IIIR families in each teleost lineage. It is anticipated that further investigations into available genomic resources and the sequencing of a variety of multiple teleost genomes will identify additional IIIR families and permit the modeling of the evolutionary origins of IIIRs.

Rare or rarely detected? Ceraceosorus guamensis sp. nov.: a second described species of Ceraceosorales and the potential for underdetection of rare lineages with common sampling techniques

Ceraceosorales is a monotypic order in Ustilaginomycotina. Its namesake, Ceraceosorus bombacis, was described as a phytopathogen of Bombax ceiba in India. In this study, we describe Ceraceosorus guamensis sp. nov., collected on the South Pacific island of Guam, which appears to represent the second isolation of any member of this order in over 40 years. Ceraceosorus species are monokaryotic and filamentous in culture, producing conidia on potato dextrose agar. However, both species behave yeast-like when cultured on corn meal agar. The internal transcribed spacer (ITS) region (spanning the ITS1-5.8S-ITS2) in both species of Ceraceosorus is highly heterogeneous containing multiple disparate copies that can vary intragenomically by up to 3.5 %. Moreover, this region could not be amplified using the fungal ITS primers most frequently used for culture-independent methods of assessing fungal biodiversity. This fact, combined with the extremely slow growth rates on commonly employed media, may indicate that members of this lineage are potentially underdetected by current sampling methods.

Horizontal transfer - imperative mission of acellular life forms, Acytota

Acytota is a kingdom of life covering satellites, plasmids, transposable elements, viroids and viruses, all outside the conventional tree of life but satisfying most life definitions. This review focuses on some aspects of Acytota, their "genomes" and life styles, the dominance of transposable elements and their evolutionary influence on other life forms in order to vindicate the Acytota as a life kingdom no more polyphyletic than other kingdoms and its members no more parasitic than other life forms.

Daphnia parasite dynamics across multiple Caullerya epidemics indicate selection against common parasite genotypes

Studies of parasite population dynamics in natural systems are crucial for our understanding of host-parasite coevolutionary processes. Some field studies have reported that host genotype frequencies in natural populations change over time according to parasite-driven negative frequency-dependent selection. However, the temporal patterns of parasite genotypes have rarely been investigated. Moreover, parasite-driven negative frequency-dependent selection is contingent on the existence of genetic specificity between hosts and parasites. In the present study, the population dynamics and host-genotype specificity of the ichthyosporean Caullerya mesnili, a common endoparasite of Daphnia water fleas, were analysed based on the observed sequence variation in the first internal transcribed spacer (ITS1) of the ribosomal DNA. The Daphnia population of lake Greifensee (Switzerland) was sampled and subjected to parasite screening and host genotyping during C. mesnili epidemics of four consecutive years. The ITS1 of wild-caught C. mesnili-infected Daphnia was sequenced using the 454 pyrosequencing platform. The relative frequencies of C. mesnili ITS1 sequences differed significantly among years: the most abundant C. mesnili ITS1 sequence decreased and rare sequences increased over the course of the study, a pattern consistent with negative frequency-dependent selection. However, only a weak signal of host-genotype specificity between C. mesnili and Daphnia genotypes was detected. Use of cutting edge genomic techniques will allow further investigation of the underlying micro-evolutionary relationships within the Daphnia-C. mesnili system.

Polymorphic Mucin-Like Proteins in Schistosoma mansoni, a Variable Antigen and a Key Component of the Compatibility Between the Schistosome and Its Snail Host

The arms race between vertebrate hosts and parasites has led to diversification systems able to generate huge repertoires of immune recognition receptors and antigenic variants. Until recently, the invertebrate immunity was considered to be poorly specific, and consequently, antigenic variability was not expected to be high for their respective parasites. In the present chapter, we show how the study of the interaction between the snail Biomphalaria glabrata and its parasite Schistosome mansoni has shaken this paradigm. We show that the fate of the interaction between the snail and its parasite is at least partly the result of the concordance of highly variable repertoires of immune recognition receptors in the snail and corresponding antigenic variants in the parasite. We call these antigenic variants of the schistosome Schistosoma mansoni polymorphic mucins (SmPoMucs). We show that their high level of diversification is the result of a complex cascade of mechanisms, thus presenting evidence for antigenic variation in a parasite infecting an invertebrate species.

Balancing selection and recombination as evolutionary forces caused population genetic variations in golden pheasant MHC class I genes

BACKGROUND

The major histocompatibility complex (MHC) genes are vital partners in the acquired immune processes of vertebrates. MHC diversity may be directly associated with population resistance to infectious pathogens. Here, we screened for polymorphisms in exons 2 and 3 of the IA1 and IA2 genes in 12 golden pheasant populations across the Chinese mainland to characterize their genetic variation levels, to understand the effects of historical positive selection and recombination in shaping class I diversity, and to investigate the genetic structure of wild golden pheasant populations.

RESULTS

Among 339 individual pheasants, we identified 14 IA1 alleles in exon 2 (IA1-E2), 11 IA1-E3 alleles, 27 IA2-E2 alleles, and 28 IA2-E3 alleles. The non-synonymous substitution rate was significantly greater than the synonymous substitution rate at sequences in the IA2 gene encoding putative peptide-binding sites but not in the IA1 gene; we also found more positively selected sites in IA2 than in IA1. Frequent recombination events resulted in at least 9 recombinant IA2 alleles, in accordance with the intermingling pattern of the phylogenetic tree. Although some IA alleles are widely shared among studied populations, large variation occurs in the number of IA alleles across these populations. Allele frequency analysis across 2 IA loci showed low levels of genetic differentiation among populations on small geographic scales; however, significant genetic differentiation was observed between pheasants from the northern and southern regions of the Yangtze River. Both STRUCTURE analysis and F-statistic (F ST ) value comparison classified those populations into 2 major groups: the northern region of the Yangtze River (NYR) and the southern region of the Yangtze River (SYR).

CONCLUSIONS

More extensive polymorphisms in IA2 than IA1 indicate that IA2 has undergone much stronger positive-selection pressure during evolution. Moreover, the recombination events detected between the genes and the intermingled phylogenetic pattern indicate that interlocus recombination accounts for much of the allelic variation in IA2. Analysis of the population differentiation implied that homogenous balancing selection plays an important part in maintaining an even distribution of MHC variations. The natural barrier of the Yangtze River and heterogeneous balancing selection might help shape the NYR-SYR genetic structure in golden pheasants.

Tubulanus tamias sp. nov. (Nemertea: Palaeonemertea) with Two Different Types of Epidermal Eyes

Based on specimens collected subtidally (∼10 m in depth) in Tomioka Bay, Japan, we describe the palaeonemertean Tubulanus tamias sp. nov., which differs from all its congeners in body coloration. In molecular phylogenetic analyses based on partial sequences of the nuclear 18S and 28S rRNA genes and histone H3, as well as the mitochondrial 16S rRNA and cytochrome c oxidase subunit I genes, among selected palaeonemerteans, T. tamias nested with part of the congeners in Tubulanus, while the genus as currently diagnosed appears to be non-monophyletic. Molecular cloning detected polymorphism in 28S rDNA sequences in a single individual of T. tamias, indicating incomplete concerted evolution of multiple copies. Tubulanus tamias is peculiar among tubulanids in having 9-10 pigment-cup eyes in the epidermis on either side of the head anterior to the cerebral sensory organs, and remarkably there are two types of eyes. The anterior 8-9 pairs of eyes, becoming larger from anterior to posterior, are completely embedded in the epidermis and proximally abutting the basement membrane; each pigment cup contains bundle of up to seven, rod-shaped structure that resemble a rhabdomeric photoreceptor cell. In contrast, the posterior-most pair of eyes, larger than most of the anterior ones, have an optical cavity filled with long cilia and opening to the exterior, thus appearing to have ciliary-type photoreceptor cells. The size and arrangement of the eyes indicate that the posterior-most pair of eyes are the remnant of the larval (or juvenile) eyes.

Molecular Differentiation of the African Yellow Fever Vector Aedes bromeliae (Diptera: Culicidae) from Its Sympatric Non-vector Sister Species, Aedes lilii

INTRODUCTION

Yellow fever continues to be a problem in sub-Saharan Africa with repeated epidemics occurring. The mosquito Aedes bromeliae is a major vector of yellow fever, but it cannot be readily differentiated from its non-vector zoophilic sister species Ae. lilii using morphological characters. Genetic differences have been reported between anthropophilic Ae. bromeliae and zoophilic Ae. lilii and between forest and domestic populations. However, due to the application of different molecular markers and non-overlapping populations employed in previous studies, interpretation of species delimitation is unclear.

METHODOLOGY/PRINCIPLE FINDINGS

DNA sequences were generated from specimens of Ae. simpsoni s.l. from the Republic of Benin, Tanzania and Uganda for two nuclear genes apolipophorin 2 (apoLp2) and cytochrome p450 (CYPJ92), the ribosomal internal transcribed spacer region (ITS) and the mitochondrial cytochrome c oxidase (COI) barcoding region. Nuclear genes apoLp2 and CYPJ92 were unable to differentiate between species Ae. bromeliae and Ae. lilii due to ancestral lineage sorting, while ITS sequence data provided clear topological separation on a phylogeny. The standard COI barcoding region was shown to be subject to species introgression and unable to clearly distinguish the two taxa. Here we present a reliable direct PCR-based method for differentiation of the vector species Ae. bromeliae from its isomorphic, sympatric and non-biomedically important sister taxon, Ae. lilii, based on the ITS region. Using molecular species verification, we describe novel immature habitats for Ae. lilii and report both sympatric and allopatric populations. Whereas only Ae. lilii is found in the Republic of Benin and only Ae. bromeliae in Tanzania, both species are sympatric in Uganda.

CONCLUSIONS/SIGNIFICANCE

Our accurate identification method will allow informed distribution and detailed ecological studies that will facilitate assessment of arboviral disease risk and development of future targeted vector control.

Patterns of Gene Conversion in Duplicated Yeast Histones Suggest Strong Selection on a Coadapted Macromolecular Complex

We find evidence for interlocus gene conversion in five duplicated histone genes from six yeast species. The sequences of these duplicated genes, surviving from the ancient genome duplication, show phylogenetic patterns inconsistent with the well-resolved orthology relationships inferred from a likelihood model of gene loss after the genome duplication. Instead, these paralogous genes are more closely related to each other than any is to its nearest ortholog. In addition to simulations supporting gene conversion, we also present evidence for elevated rates of radical amino acid substitutions along the branches implicated in the conversion events. As these patterns are similar to those seen in ribosomal proteins that have undergone gene conversion, we speculate that in cases where duplicated genes code for proteins that are a part of tightly interacting complexes, selection may favor the fixation of gene conversion events in order to maintain high protein identities between duplicated copies.

High-Throughput Sequencing-The Key to Rapid Biodiversity Assessment of Marine Metazoa?

The applications of traditional morphological and molecular methods for species identification are greatly restricted by processing speed and on a regional or greater scale are generally considered unfeasible. In this context, high-throughput sequencing, or metagenetics, has been proposed as an efficient tool to document biodiversity. Here we evaluated the effectiveness of 454 pyrosequencing in marine metazoan community analysis using the 18S rDNA: V1-V2 region. Multiplex pyrosequencing of the V1-V2 region was used to analyze two pooled samples of DNA, one comprising 118 and the other 37 morphologically identified species, and one natural sample taken directly from a North Sea zooplankton community. A DNA reference library comprising all species represented in the pooled samples was created by Sanger sequencing, and this was then used to determine the optimal similarity threshold for species delineation. The optimal threshold was found at 99% species similarity, with 85% identification success. Pyrosequencing was able to identify between fewer species: 67% and 78% of the species in the two pooled samples. Also, a large number of sequences for three species that were not included in the pooled samples were amplified by pyrosequencing, suggesting preferential amplification of some genotypes and the sensitivity of this approach to even low levels of contamination. Conversely, metagenetic analysis of the natural zooplankton sample identified many more species (particularly gelatinous zooplankton and meroplankton) than morphological analysis of a formalin-fixed sample from the same sampling site, suggesting an increased level of taxonomic resolution with pyrosequencing. The study demonstrated that, based on the V1-V2 region, 454 sequencing does not provide accurate species differentiation and reliable taxonomic classification, as it is required in most biodiversity monitoring. The analysis of artificially prepared samples indicated that species detection in pyrosequencing datasets is complicated by potential PCR-based biases and that the V1-V2 marker is poorly resolved for some taxa.

Long-Lasting Gene Conversion Shapes the Convergent Evolution of the Critical Methanogenesis Genes

Methanogenesis and its key small-molecule methyltransferase Mtr complex are poorly understood despite their pivotal role in Earth’s global carbon cycle. Mtr complex is encoded by a conserved mtrEDCBAFGH operon in most methanogens. Here we report that two discrete lineages, Methanococcales and Methanomicrobiales, have a noncanonical mtr operon carrying two copies of mtrA resulting from an ancient duplication. Compared to mtrA-1, mtrA-2 acquires a distinct transmembrane domain through domain shuffling and gene fusion. However, the nontransmembrane domains (MtrA domain) of mtrA-1 and mtrA-2 are homogenized by gene conversion events lasting throughout the long history of these extant methanogens (over 2410 million years). Furthermore, we identified a possible recruitment of ancient nonmethanogenic methyltransferase genes to establish the methanogenesis pathway. These results not only provide novel evolutionary insight into the methanogenesis pathway and methyltransferase superfamily but also suggest an unanticipated long-lasting effect of gene conversion on gene evolution in a convergent pattern.

Intragenomic variations of multicopy ITS2 marker in Agrodiaetus blue butterflies (Lepidoptera, Lycaenidae)

The eukaryotic ribosomal DNA cluster consists of multiple copies of three genes, 18S, 5. 8S and 28S rRNAs, separated by multiple copies of two internal transcribed spacers, ITS1 and ITS2. It is an important, frequently used marker in both molecular cytogenetic and molecular phylogenetic studies. Despite this, little is known about intragenomic variations within the copies of eukaryotic ribosomal DNA genes and spacers. Here we present data on intraindividual variations of ITS2 spacer in three species of Agrodiaetus Hübner, 1822 blue butterflies revealed by cloning technique. We demonstrate that a distinctly different intragenomic ITS2 pattern exists for every individual analysed. ITS2 sequences of these species show significant intragenomic variation (up to 3.68% divergence), setting them apart from each other on inferred phylogenetic tree. This variation is enough to obscure phylogenetic relationships at the species level.

Concerted Evolution of Duplicate Control Regions in the Mitochondria of Species of the Flatfish Family Bothidae (Teleostei: Pleuronectiformes)

Mitogenomes of flatfishes (Pleuronectiformes) exhibit the greatest diversity of gene rear-rangements in teleostean fishes. Duplicate control regions (CRs) have been found in the mito-genomes of two flatfishes, Samariscus latus (Samaridae) and Laeops lanceolata (Bothidae), which is rare in teleosts. It has been reported that duplicate CRs have evolved in a concerted fashion in fishes and other animals, however, whether concerted evo-lution exists in flatfishes remains unknown. In this study, based on five newly sequenced and six previously reported mitogenomes of lefteye flounders in the Bothidae, we explored whether duplicate CRs and concerted evolution exist in these species. Results based on the present study and previous reports show that four out of eleven bothid species examined have duplicate CRs of their mitogenomes. The core regions of the duplicate CRs of mitogenomes in the same species have identical, or nearly identical, sequences when compared to each other. This pattern fits the typical characteristics of concerted evolution. Additionally, phylogenetic and ancestral state reconstruction analysis also provided evidence to support the hypothesis that duplicate CRs evolved concertedly. The core region of concerted evolution is situated at the conserved domains of the CR of the mitogenome from the termination associated sequences (TASs) to the conserved sequence blocks (CSBs). Commonly, this region is con-sidered to regulate mitochondrial replication and transcription. Thus, we hypothesize that the cause of concerted evolution of the duplicate CRs in the mtDNAs of these four bothids may be related to some function of the conserved sequences of the CRs during mitochondrial rep-lication and transcription. We hope our results will provide fresh insight into the molecular mechanisms related to replication and evolution of mitogenomes.

Trans-species polymorphism at antimicrobial innate immunity cathelicidin genes of Atlantic cod and related species

Natural selection, the most important force in evolution, comes in three forms. Negative purifying selection removes deleterious variation and maintains adaptations. Positive directional selection fixes beneficial variants, producing new adaptations. Balancing selection maintains variation in a population. Important mechanisms of balancing selection include heterozygote advantage, frequency-dependent advantage of rarity, and local and fluctuating episodic selection. A rare pathogen gains an advantage because host defenses are predominantly effective against prevalent types. Similarly, a rare immune variant gives its host an advantage because the prevalent pathogens cannot escape the host's apostatic defense. Due to the stochastic nature of evolution, neutral variation may accumulate on genealogical branches, but trans-species polymorphisms are rare under neutrality and are strong evidence for balancing selection. Balanced polymorphism maintains diversity at the major histocompatibility complex (MHC) in vertebrates. The Atlantic cod is missing genes for both MHC-II and CD4, vital parts of the adaptive immune system. Nevertheless, cod are healthy in their ecological niche, maintaining large populations that support major commercial fisheries. Innate immunity is of interest from an evolutionary perspective, particularly in taxa lacking adaptive immunity. Here, we analyze extensive amino acid and nucleotide polymorphisms of the cathelicidin gene family in Atlantic cod and closely related taxa. There are three major clusters, Cath1, Cath2, and Cath3, that we consider to be paralogous genes. There is extensive nucleotide and amino acid allelic variation between and within clusters. The major feature of the results is that the variation clusters by alleles and not by species in phylogenetic trees and discriminant analysis of principal components. Variation within the three groups shows trans-species polymorphism that is older than speciation and that is suggestive of balancing selection maintaining the variation. Using Bayesian and likelihood methods positive and negative selection is evident at sites in the conserved part of the genes and, to a larger extent, in the active part which also shows episodic diversifying selection, further supporting the argument for balancing selection.

Molecular taxonomy provides new insights into anopheles species of the neotropical arribalzagia series

Phylogenetic analysis of partial mitochondrial cytochrome oxidase c subunit I (COI) and nuclear internal transcribed spacer 2 (ITS2) sequences were used to evaluate initial identification and to investigate phylogenetic relationships of seven Anopheles morphospecies of the Arribalzagia Series from Colombia. Phylogenetic trees recovered highly supported clades for An. punctimaculas.s., An. calderoni, An. malefactor s.l., An. neomaculipalpus, An. apicimacula s.l., An. mattogrossensis and An. peryassui. This study provides the first molecular confirmation of An. malefactorfrom Colombia and discovered conflicting patterns of divergence for the molecular markers among specimens from northeast and northern Colombia suggesting the presence of two previously unrecognized Molecular Operational Taxonomic Units (MOTUs). Furthermore, two highly differentiated An. apicimacula MOTUs previously found in Panama were detected. Overall, the combined molecular dataset facilitated the detection of known and new Colombian evolutionary lineages, and constitutes the baseline for future research on their bionomics, ecology and potential role as malaria vectors.

The link between independent acquisition of intracellular gamma-endosymbionts and concerted evolution in Tremblaya princeps

Many insect species establish mutualistic symbiosis with intracellular bacteria that complement their unbalanced diets. The betaproteobacterium "Candidatus Tremblaya" maintains an ancient symbiosis with mealybugs (Hemiptera: Pseudococcidae), which are classified in subfamilies Phenacoccinae and Pseudococcinae. Most Phenacoccinae mealybugs have "Candidatus Tremblaya phenacola" as their unique endosymbiont, while most Pseudococcinae mealybugs show a nested symbiosis (a bacterial symbiont placed inside another one) where every "Candidatus Tremblaya princeps" cell harbors several cells of a gammaproteobacterium. Genomic characterization of the endosymbiotic consortium from Planococcus citri, composed by "Ca. Tremblaya princeps" and "Candidatus Moranella endobia," unveiled several atypical features of the former's genome, including the concerted evolution of paralogous loci. Its comparison with the genome of "Ca. Tremblaya phenacola" PAVE, single endosymbiont of Phenacoccus avenae, suggests that the atypical reductive evolution of "Ca. Tremblaya princeps" could be linked to the acquisition of "Ca. Moranella endobia," which possess an almost complete set of genes encoding proteins involved in homologous recombination. In order to test this hypothesis, we performed comparative genomics between "Ca. Tremblaya phenacola" and "Ca. Tremblaya princeps" and searched for the co-occurrence of concerted evolution and homologous recombination genes in endosymbiotic consortia from four unexplored mealybug species, Dysmicoccus boninsis, Planococcus ficus, Pseudococcus longispinus, and Pseudococcus viburni. Our results support a link between concerted evolution and nested endosymbiosis.

Detecting regular sound changes in linguistics as events of concerted evolution

Background

Concerted evolution is normally used to describe parallel changes at different sites in a genome, but it is also observed in languages where a specific phoneme changes to the same other phoneme in many words in the lexicon—a phenomenon known as regular sound change. We develop a general statistical model that can detect concerted changes in aligned sequence data and apply it to study regular sound changes in the Turkic language family.

Results

Linguistic evolution, unlike the genetic substitutional process, is dominated by events of concerted evolutionary change. Our model identified more than 70 historical events of regular sound change that occurred throughout the evolution of the Turkic language family, while simultaneously inferring a dated phylogenetic tree. Including regular sound changes yielded an approximately 4-fold improvement in the characterization of linguistic change over a simpler model of sporadic change, improved phylogenetic inference, and returned more reliable and plausible dates for events on the phylogenies. The historical timings of the concerted changes closely follow a Poisson process model, and the sound transition networks derived from our model mirror linguistic expectations.

Conclusions

We demonstrate that a model with no prior knowledge of complex concerted or regular changes can nevertheless infer the historical timings and genealogical placements of events of concerted change from the signals left in contemporary data. Our model can be applied wherever discrete elements—such as genes, words, cultural trends, technologies, or morphological traits—can change in parallel within an organism or other evolving group.

•

Linguistic evolution is dominated by events of concerted evolutionary change

•

Modeling concerted evolution improves phylogenetic inference and dating

•

Events of concerted change conform closely to a Poisson process

•

Our model can be applied to genes, languages, cultures, and technological change

A systematic computational analysis of biosynthetic gene cluster evolution: lessons for engineering biosynthesis

Bacterial secondary metabolites are widely used as antibiotics, anticancer drugs, insecticides and food additives. Attempts to engineer their biosynthetic gene clusters (BGCs) to produce unnatural metabolites with improved properties are often frustrated by the unpredictability and complexity of the enzymes that synthesize these molecules, suggesting that genetic changes within BGCs are limited by specific constraints. Here, by performing a systematic computational analysis of BGC evolution, we derive evidence for three findings that shed light on the ways in which, despite these constraints, nature successfully invents new molecules: 1) BGCs for complex molecules often evolve through the successive merger of smaller sub-clusters, which function as independent evolutionary entities. 2) An important subset of polyketide synthases and nonribosomal peptide synthetases evolve by concerted evolution, which generates sets of sequence-homogenized domains that may hold promise for engineering efforts since they exhibit a high degree of functional interoperability, 3) Individual BGC families evolve in distinct ways, suggesting that design strategies should take into account family-specific functional constraints. These findings suggest novel strategies for using synthetic biology to rationally engineer biosynthetic pathways.

Bacterial secondary metabolites mediate a broad range of microbe-microbe and microbe-host interactions, and are widely used in human medicine, agriculture and manufacturing. Despite recent advances in synthetic biology, efforts to engineer their biosynthetic genes for the production of unnatural variants are frustrated by a high failure rate. In an effort to better understand what types of genetic changes are most likely to lead to successful improvements, we systematically analyzed the ways in which biosynthetic genes naturally evolve to generate new compounds. We show that large gene clusters appear to evolve through the merger of sub-clusters, which function independently, and are promising units for cluster engineering. Moreover, a subset of gene clusters evolve by concerted evolution, which generates sets of interoperable domains that may enable predictable domain swapping. Finally, many biosynthetic gene clusters evolve in family-specific modes that differ greatly from each other. Overall, this quantitative perspective on the ways in which gene clusters naturally evolve suggests novel strategies for using synthetic biology to engineer the production of unnatural metabolites.

Differential amplification of satellite PaB6 in chromosomally hypervariable Prospero autumnale complex (Hyacinthaceae)

BACKGROUND AND AIMS

Chromosomal evolution, including numerical and structural changes, is a major force in plant diversification and speciation. This study addresses genomic changes associated with the extensive chromosomal variation of the Mediterranean Prospero autumnale complex (Hyacinthaceae), which includes four diploid cytotypes each with a unique combination of chromosome number (x = 5, 6, 7), rDNA loci and genome size.

METHODS

A new satellite repeat PaB6 has previously been identified, and monomers were reconstructed from next-generation sequencing (NGS) data of P. autumnale cytotype B(6)B(6) (2n = 12). Monomers of all other Prospero cytotypes and species were sequenced to check for lineage-specific mutations. Copy number, restriction patterns and methylation levels of PaB6 were analysed using Southern blotting. PaB6 was localized on chromosomes using fluorescence in situ hybridization (FISH).

KEY RESULTS

The monomer of PaB6 is 249 bp long, contains several intact and truncated vertebrate-type telomeric repeats and is highly methylated. PaB6 is exceptional because of its high copy number and unprecedented variation among diploid cytotypes, ranging from 10(4) to 10(6) copies per 1C. PaB6 is always located in pericentromeric regions of several to all chromosomes. Additionally, two lineages of cytotype B(7)B(7) (x = 7), possessing either a single or duplicated 5S rDNA locus, differ in PaB6 copy number; the ancestral condition of a single locus is associated with higher PaB6 copy numbers.

CONCLUSIONS

Although present in all Prospero species, PaB6 has undergone differential amplification only in chromosomally variable P. autumnale, particularly in cytotypes B(6)B(6) and B(5)B(5). These arose via independent chromosomal fusions from x = 7 to x = 6 and 5, respectively, accompanied by genome size increases. The copy numbers of satellite DNA PaB6 are among the highest in angiosperms, and changes of PaB6 are exceptionally dynamic in this group of closely related cytotypes of a single species. The evolution of the PaB6 copy numbers is discussed, and it is suggested that PaB6 represents a recent and highly dynamic system originating from a small pool of ancestral repeats.

The organization and evolution of the Responder satellite in species of the Drosophila melanogaster group: dynamic evolution of a target of meiotic drive

Background

Satellite DNA can make up a substantial fraction of eukaryotic genomes and has roles in genome structure and chromosome segregation. The rapid evolution of satellite DNA can contribute to genomic instability and genetic incompatibilities between species. Despite its ubiquity and its contribution to genome evolution, we currently know little about the dynamics of satellite DNA evolution. The Responder (Rsp) satellite DNA family is found in the pericentric heterochromatin of chromosome 2 of Drosophila melanogaster. Rsp is well-known for being the target of Segregation Distorter (SD)— an autosomal meiotic drive system in D. melanogaster. I present an evolutionary genetic analysis of the Rsp family of repeats in D. melanogaster and its closely-related species in the melanogaster group (D. simulans, D. sechellia, D. mauritiana, D. erecta, and D. yakuba) using a combination of available BAC sequences, whole genome shotgun Sanger reads, Illumina short read deep sequencing, and fluorescence in situ hybridization.

Results

I show that Rsp repeats have euchromatic locations throughout the D. melanogaster genome, that Rsp arrays show evidence for concerted evolution, and that Rsp repeats exist outside of D. melanogaster, in the melanogaster group. The repeats in these species are considerably diverged at the sequence level compared to D. melanogaster, and have a strikingly different genomic distribution, even between closely-related sister taxa.

Conclusions

The genomic organization of the Rsp repeat in the D. melanogaster genome is complex—it exists of large blocks of tandem repeats in the heterochromatin and small blocks of tandem repeats in the euchromatin. My discovery of heterochromatic Rsp-like sequences outside of D. melanogaster suggests that SD evolved after its target satellite and that the evolution of the Rsp satellite family is highly dynamic over a short evolutionary time scale (<240,000 years).

Electronic supplementary material

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

Diversity and evolution of myxozoan minicollagens and nematogalectins

Background

Myxozoa are a diverse group of metazoan parasites with a very simple organization, which has for decades eluded their evolutionary origin. Their most prominent and characteristic feature is the polar capsule: a complex intracellular structure of the myxozoan spore, which plays a role in host infection. Striking morphological similarities have been found between myxozoan polar capsules and nematocysts, the stinging structures of cnidarians (corals, sea anemones and jellyfish) leading to the suggestion that Myxozoa and Cnidaria share a more recent common ancestry. This hypothesis has recently been supported by phylogenomic evidence and by the identification of a nematocyst specific minicollagen gene in the myxozoan Tetracapsuloides bryosalmonae. Here we searched genomes and transcriptomes of several myxozoan taxa for the presence of additional cnidarian specific genes and characterized these genes within a phylogenetic context.

Results

Illumina assemblies of transcriptome or genome data of three myxozoan species (Enteromyxum leei, Kudoa iwatai, and Sphaeromyxa zaharoni) and of the enigmatic cnidarian parasite Polypodium hydriforme (Polypodiozoa) were mined using tBlastn searches with nematocyst-specific proteins as queries. Several orthologs of nematogalectins and minicollagens were identified. Our phylogenetic analyses indicate that myxozoans possess three distinct minicollagens. We found that the cnidarian repertoire of nematogalectins is more complex than previously thought and we identified additional members of the nematogalectin family. Cnidarians were found to possess four nematogalectin/ nematogalectin-related genes, while in myxozoans only three genes could be identified.

Conclusions

Our results demonstrate that myxozoans possess a diverse array of genes that are taxonomically restricted to Cnidaria. Characterization of these genes provide compelling evidence that polar capsules and nematocysts are homologous structures and that myxozoans are highly degenerate cnidarians. The diversity of minicollagens was higher than previously thought, with the presence of three minicollagen genes in myxozoans. Our phylogenetic results suggest that the different myxozoan sequences are the results of ancient divergences within Cnidaria and not of recent specializations of the polar capsule. For both minicollagen and nematogalectin, our results show that myxozoans possess less gene copies than their cnidarian counter parts, suggesting that the polar capsule gene repertoire was simplified with their reduced body plan.

Electronic supplementary material

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

Trypsin isozymes in the lobster Panulirus argus (Latreille, 1804): from molecules to physiology

Trypsin enzymes have been studied in a wide variety of animal taxa due to their central role in protein digestion as well as in other important physiological and biotechnological processes. Crustacean trypsins exhibit a high number of isoforms. However, while differences in properties of isoenzymes are known to play important roles in regulating different physiological processes, there is little information on this aspect for decapod trypsins. The aim of this review is to integrate recent findings at the molecular level on trypsin enzymes of the spiny lobster Panulirus argus, into higher levels of organization (biochemical, organism) and to interpret those findings in relation to the feeding ecology of these crustaceans. Trypsin in lobster is a polymorphic enzyme, showing isoforms that differ in their biochemical features and catalytic efficiencies. Molecular studies suggest that polymorphism in lobster trypsins may be non-neutral. Trypsin isoenzymes are differentially regulated by dietary proteins, and it seems that some isoenzymes have undergone adaptive evolution coupled with a divergence in expression rate to increase fitness. This review highlights important but poorly studied issues in crustaceans in general, such as the relation among trypsin polymorphism, phenotypic (digestive) flexibility, digestion efficiency, and feeding ecology.