Search

BackgroundTrypanosomes are single-celled eukaryotic parasites characterised by the unique biology of their mitochondrial DNA. African livestock trypanosomes impose a major burden on agriculture across sub-Saharan Africa, but are poorly understood compared to those that cause sleeping sickness and Chagas disease in humans. Here we explore the potential of the maxicircle, a component of trypanosome mitochondrial DNA to study the evolutionary history of trypanosomes.ResultsWe used long-read sequencing to completely assemble maxicircle mitochondrial DNA from four previously uncharacterized African trypanosomes, and leveraged these assemblies to scaffold and assemble a further 103 trypanosome maxicircle gene coding regions from published short-read data. While synteny was largely conserved, there were repeated, independent losses of Complex I genes. Comparison of pre-edited and non-edited genes revealed the impact of RNA editing on nucleotide composition, with non-edited genes approaching the limits of GC loss. African tsetse-transmitted trypanosomes showed high levels of RNA editing compared to other trypanosomes. The gene coding regions of maxicircle mitochondrial DNAs were used to construct time-resolved phylogenetic trees, revealing deep divergence events among isolates of the pathogensTrypanosoma bruceiandT. congolense.ConclusionsOur data represents a new resource for experimental and evolutionary analyses of trypanosome phylogeny, molecular evolution and function. Molecular clock analyses yielded a timescale for trypanosome evolution congruent with major biogeographical events in Africa and revealed the recent emergence ofTrypanosoma brucei gambienseandT. equiperdum, major human and animal pathogens.

BackgroundEP300 is a conserved protein in vertebrates, which serves as a key mediator of cellular homeostasis. Mutations and dysregulation of EP300 give rise to severe human developmental disorders and malignancy.Danio reriois a promising model organism to study EP300 related diseases and drugs; however, the effect of EP300 duplicates derived from teleost-specific whole genome duplication should not just be neglected.ResultsIn this study, we obtained EP300 protein sequences of representative teleosts, mammals and sauropsids, with which we inferred a highly supported maximum likelihood tree. We observed that Ep300 duplicates (Ep300a and Ep300b) were widely retained in teleosts and universally expressed in a variety of tissues. Consensus sequences of Ep300a and Ep300b had exactly the same distribution of conserved domains, suggesting that their functions should still be largely overlapped. We analyzed the molecular evolution of Ep300 duplicates in teleosts, using branch-site models, clade models and site models. The results showed that both duplicates were subject to strong positive selection; however, for an extant species, generally at most one copy was under positive selection. At the clade level, there were evident positive correlations between evolutionary rates, the number of positively selected sites and gene expression levels. In Ostariophysi, Ep300a were under stronger positive selection than Ep300b; in Neoteleostei, another species-rich teleost clade, the contrary was the case. We also modeled 3D structures of zf-TAZ domain and its flanking regions of Ep300a and Ep300b ofD. rerioandOryzias latipesand found that in either species the faster evolving copy had more short helixes.ConclusionsCollectively, the two copies of Ep300 have undoubtedly experienced directional divergence in main teleost clades. The divergence of EP300 between teleosts and mammals should be greater than the divergence between different teleost clades. Further studies are needed to clarify to what extent the EP300 involved regulatory network has diverged between teleosts and mammals, which would also help explain the huge success of teleosts.

BackgroundThe aragonite shelled, planktonic gastropod family Atlantidae (shelled heteropods) is likely to be one of the first groups to be impacted by imminent ocean changes, including ocean warming and ocean acidification. With a fossil record spanning at least 100 Ma, atlantids have experienced and survived global-scale ocean changes and extinction events in the past. However, the diversification patterns and tempo of evolution in this family are largely unknown.ResultsBased on a concatenated maximum likelihood phylogeny of three genes (cytochromecoxidase subunit 1 mitochondrial DNA, 28S and 18S ribosomal rRNA) we show that the three extant genera of the family Atlantidae,Atlanta, ProtatlantaandOxygyrus,form monophyletic groups. The genusAtlantais split into two groups, one exhibiting smaller, well ornamented shells, and the other having larger, less ornamented shells. The fossil record, in combination with a fossil-calibrated phylogeny, suggests that large scale atlantid extinction was accompanied by considerable and rapid diversification over the last 25 Ma, potentially driven by vicariance events.ConclusionsNow confronted with a rapidly changing modern ocean, the ability of atlantids to survive past global change crises gives some optimism that they may be able to persist through the Anthropocene.

Background The distribution of genetic diversity and the underlying processes are important for conservation planning but are unknown for most species and have not been well studied in many regions. In East Asia, the Sichuan Basin and surrounding mountains constitute an understudied region that exhibits a “ring” of high species richness overlapping the eastern edge of the global biodiversity hotspot Mountains of Southwest China. We examine the distributional history and genetic diversification of the Emei mustache toad Leptobrachium boringii, a typical “ring” element characterized by disjunct ranges in the mountains, by integrating time-calibrated gene tree, genetic variability, individual-level clustering, inference of population splitting and mixing from allele frequencies, and paleoclimatic suitability modeling. Results The results reveal extensive range dynamics, including secondary contact after long-term isolation via westward dispersal accompanied by variability loss. They allow the proposal of a model that combines recurrent contractions caused by Quaternary climatic changes and some failed expansions under suitable conditions for explaining the shared disjunct distribution pattern. Providing exceptional low-elevation habitats in the hotspot area, the eastern edge harbors both long-term refugial and young immigrant populations. This finding and a synthesis of evidence from other taxa demonstrate that a certain contributor to biodiversity, one that preserves and receives low-elevation elements of the east in this case, can be significant for only a particular part of a hotspot. By clarifying the low variability of these refugial populations, we show that discordant mitochondrial estimates of diversity can be obtained for populations that experienced admixture, which would have unlikely left proportional immigrant alleles for each locus. Conclusions Dispersal after long-term isolation can explain much of the spatial distribution of genetic diversity in this species, while secondary contact and long-term persistence do not guarantee a large variation. The model for the formation of disjunct ranges may apply to many other taxa isolated in the mountains surrounding the Sichuan Basin. Furthermore, this study provides insights into the heterogeneous nature of hotspots and discordant variability obtained from genome-wide and mitochondrial data.

Background Many animals rely heavily on olfaction to navigate their environment. Among rodents, olfaction is crucial for a wide range of social behaviors. The vomeronasal olfactory system in particular plays an important role in mediating social communication, including the detection of pheromones and recognition signals. In this study we examine patterns of vomeronasal type-1 receptor (V1R) evolution in the house mouse and related species within the genus Mus. We report the extent of gene repertoire turnover and conservation among species and clades, as well as the prevalence of positive selection on gene sequences across the V1R tree. By exploring the evolution of these receptors, we provide insight into the functional roles of receptor subtypes as well as the dynamics of gene family evolution. Results We generated transcriptomes from the vomeronasal organs of 5 Mus species, and produced high quality V1R repertoires for each species. We find that V1R clades in the house mouse and relatives exhibit distinct evolutionary trajectories. We identify putative species-specific gene expansions, including a large clade D expansion in the house mouse. While gene gains are abundant, we detect very few gene losses. We describe a novel V1R clade and highlight candidate receptors for future study. We find evidence for distinct evolutionary processes across different clades, from largescale turnover to highly conserved repertoires. Patterns of positive selection are similarly variable, as some clades exhibit abundant positive selection while others display high gene sequence conservation. Based on clade-level evolutionary patterns, we identify receptor families that are strong candidates for detecting social signals and predator cues. Our results reveal clades with receptors detecting female reproductive status are among the most conserved across species, suggesting an important role in V1R chemosensation. Conclusion Analysis of clade-level evolution is critical for understanding species’ chemosensory adaptations. This study provides clear evidence that V1R clades are characterized by distinct evolutionary trajectories. As receptor evolution is shaped by ligand identity, these results provide a framework for examining the functional roles of receptors.

Background Although phylogenomic analyses are increasingly used to reveal evolutionary relationships among ciliates, relatively few nuclear protein-coding gene markers have been tested for their suitability as candidates for inferring phylogenies within this group. In this study, we investigate the utility of the heat-shock protein 90 gene (Hsp90) as a marker for inferring phylogenetic relationships among hypotrich ciliates. Results A total of 87 novel Hsp90 gene sequences of 10 hypotrich species were generated. Of these, 85 were distinct sequences. Phylogenetic analyses based on these data showed that: (1) the Hsp90 gene amino acid trees are comparable to the small subunit rDNA tree for recovering phylogenetic relationships at the rank of class, but lack sufficient phylogenetic signal for inferring evolutionary relationships at the genus level; (2) Hsp90 gene paralogs are recent and therefore unlikely to pose a significant problem for recovering hypotrich clades; (3) definitions of some hypotrich orders and families need to be revised as their monophylies are not supported by various gene markers; (4) The order Sporadotrichida is paraphyletic, but the monophyly of the “core” Urostylida is supported; (5) both the subfamily Oxytrichinae and the genus Urosoma seem to be non-monophyletic, but monophyly of Urosoma is not rejected by AU tests. Conclusions Our results for the first time demonstrate that the Hsp90 gene is comparable to SSU rDNA for recovering phylogenetic relationships at the rank of class, and its paralogs are unlikely to pose a significant problem for recovering hypotrich clades. This study shows the value of careful gene marker selection for phylogenomic analyses of ciliates.

BackgroundThe latest advancements in DNA sequencing technologies have facilitated the resolution of the phylogeny of insects, yet parts of the tree of Holometabola remain unresolved. The phylogeny of Neuropterida has been extensively studied, but no strong consensus exists concerning the phylogenetic relationships within the order Neuroptera. Here, we assembled a novel transcriptomic dataset to address previously unresolved issues in the phylogeny of Neuropterida and to infer divergence times within the group. We tested the robustness of our phylogenetic estimates by comparing summary coalescent and concatenation-based phylogenetic approaches and by employing different quartet-based measures of phylogenomic incongruence, combined with data permutations.ResultsOur results suggest that the order Raphidioptera is sister to Neuroptera + Megaloptera. Coniopterygidae is inferred as sister to all remaining neuropteran families suggesting that larval cryptonephry could be a ground plan feature of Neuroptera. A clade that includes Nevrorthidae, Osmylidae, and Sisyridae (i.e. Osmyloidea) is inferred as sister to all other Neuroptera except Coniopterygidae, and Dilaridae is placed as sister to all remaining neuropteran families. Ithonidae is inferred as the sister group of monophyletic Myrmeleontiformia. The phylogenetic affinities of Chrysopidae and Hemerobiidae were dependent on the data type analyzed, and quartet-based analyses showed only weak support for the placement of Hemerobiidae as sister to Ithonidae + Myrmeleontiformia. Our molecular dating analyses suggest that most families of Neuropterida started to diversify in the Jurassic and our ancestral character state reconstructions suggest a primarily terrestrial environment of the larvae of Neuropterida and Neuroptera.ConclusionOur extensive phylogenomic analyses consolidate several key aspects in the backbone phylogeny of Neuropterida, such as the basal placement of Coniopterygidae within Neuroptera and the monophyly of Osmyloidea. Furthermore, they provide new insights into the timing of diversification of Neuropterida. Despite the vast amount of analyzed molecular data, we found that certain nodes in the tree of Neuroptera are not robustly resolved. Therefore, we emphasize the importance of integrating the results of morphological analyses with those of sequence-based phylogenomics. We also suggest that comparative analyses of genomic meta-characters should be incorporated into future phylogenomic studies of Neuropterida.

Background Tree squirrels (Sciuridae, Sciurini), in particular the highly diverse Neotropical lineages, are amongst the most rapidly diversifying branches of the mammal tree of life but also some of the least known. Negligence of this group by systematists is likely a product of the difficulties in assessing morphological informative traits and of the scarcity or unavailability of fresh tissue samples for DNA sequencing. The highly discrepant taxonomic arrangements are a consequence of the lack of phylogenies and the exclusive phenotypic-based classifications, which can be misleading in a group with conservative morphology. Here we used high-throughput sequencing and an unprecedented sampling of museum specimens to provide the first comprehensive phylogeny of tree squirrels, with a special emphasis on Neotropical taxa. Results We obtained complete or partial mitochondrial genomes from 232 historical and modern samples, representing 40 of the 43 currently recognized species of Sciurini. Our phylogenetic analyses—performed with datasets differing on levels of missing data and taxa under distinct analytical methods—strongly support the monophyly of Sciurini and consistently recovered 12 major clades within the tribe. We found evidence that the diversity of Neotropical tree squirrels is underestimated, with at least six lineages that represent taxa to be named or revalidated. Ancestral state reconstructions of number of upper premolars and number of mammae indicated that alternative conditions of both characters must have evolved multiple times throughout the evolutionary history of tree squirrels. Conclusions Complete mitogenomes were obtained from museum specimens as old as 120 years, reinforcing the potential of historical samples for phylogenetic inferences of elusive lineages of the tree of life. None of the taxonomic arrangements ever proposed for tree squirrels fully corresponded to our phylogenetic reconstruction, with only a few of the currently recognized genera recovered as monophyletic. By investigating the evolution of two morphological traits widely employed in the taxonomy of the group, we revealed that their homoplastic nature can help explain the incongruence between phylogenetic results and the classification schemes presented so far. Based on our phylogenetic results we suggest a tentative supraspecific taxonomic arrangement for Sciurini, employing 13 generic names used in previous taxonomic classifications.

Background Parasites may mediate the success of biological invasions through their effect on host fitness and thus, on host population growth and stability. However, a release from the pressure of parasites is strongly related to the genetic differentiation of the host. In invasive host populations, the number of available genetic variants, allowing them to ‘fight’ the infection, are likely to be influenced by founder events and genetic drift. The level standing genetic variation of invasive populations may be crucial in successfully adapting to new environments and resisting diseases. We studied invasive populations of raccoon that experienced a random reduction in genetic diversity during the establishment and evaluated the relationship between host immune genetic diversity and intestinal parasites infection. Results We distinguished two different genetic clusters that are characterized by different sets of functionally relevant MHC-DRB alleles. Both clusters were characterized by considerably different allele-parasite associations and different levels of parasite infection. The specific resistance MHC-DRB alleles explained the lower prevalence of Digenea parasites. An increased infection intensity was related to the presence of two MHC-DRB alleles. One of these alleles significantly decreased in frequency over time, causing a decrease of Digenea abundance in raccoons in consecutive years. Conclusions Our findings suggest that intestinal parasites can exert selective pressure on an invasive host with lowered levels of immune genetic diversity and contribute to promoting local adaptation over time. The random genetic drift that created the two different genetic clusters in the invasive raccoon range imposed completely different MHC-parasite associations, strongly associated with the infection status of populations. Our findings underline the role of standing genetic variation in shaping host-parasite relationships and provide empirical support that functional genetic variation may be, at least partly, responsible for differences in the success of invasive populations.

Background Sympatric sister species provide an opportunity to investigate the genetic mechanisms and evolutionary forces that maintain species boundaries. The persistence of morphologically and genetically distinct populations in sympatry can only occur if some degree of reproductive isolation exists. A pair of sympatric sister species of Primulina (P. depressa and P. danxiaensis) was used to explore the genetic architecture of hybrid male sterility. Results We mapped one major- and seven minor-effect quantitative trait loci (QTLs) that underlie pollen fertility rate (PFR). These loci jointly explained 55.4% of the phenotypic variation in the F2 population. A Bateson–Dobzhansky–Muller (BDM) model involving three loci was observed in this system. We found genotypic correlations between hybrid male sterility and flower morphology, consistent with the weak but significant phenotypic correlations between PFR and floral traits. Conclusions Hybrid male sterility in Primulina is controlled by a polygenic genetic basis with a complex pattern. The genetic incompatibility involves a three-locus BDM model. Hybrid male sterility is genetically correlated with floral morphology and divergence hitchhiking may occur between them.

Background Gymnosperms represent five of the six lineages of seed plants. However, most sequenced plant mitochondrial genomes (mitogenomes) have been generated for angiosperms, whereas mitogenomic sequences have been generated for only six gymnosperms. In particular, complete mitogenomes are available for all major seed plant lineages except Conifer II (non-Pinaceae conifers or Cupressophyta), an important lineage including six families, which impedes a comprehensive understanding of the mitogenomic diversity and evolution in gymnosperms. Results Here, we report the complete mitogenome of Taxus cuspidata in Conifer II. In comparison with previously released gymnosperm mitogenomes, we found that the mitogenomes of Taxus and Welwitschia have lost many genes individually, whereas all genes were identified in the mitogenomes of Cycas, Ginkgo and Pinaceae. Multiple tRNA genes and introns also have been lost in some lineages of gymnosperms, similar to the pattern observed in angiosperms. In general, gene clusters could be less conserved in gymnosperms than in angiosperms. Moreover, fewer RNA editing sites were identified in the Taxus and Welwitschia mitogenomes than in other mitogenomes, which could be correlated with fewer introns and frequent gene losses in these two species. Conclusions We have sequenced the Taxus cuspidata mitogenome, and compared it with mitogenomes from the other four gymnosperm lineages. The results revealed the diversity in size, structure, gene and intron contents, foreign sequences, and mutation rates of gymnosperm mitogenomes, which are different from angiosperm mitogenomes.

Background Early Cambrian Lagerstätten from China have greatly enriched our perspective on the early evolution of animals, particularly arthropods. However, recent studies have shown that many of these early fossil arthropods were more derived than previously thought, casting uncertainty on the ancestral euarthropod body plan. In addition, evidence from fossilized neural tissues conflicts with external morphology, in particular regarding the homology of the frontalmost appendage. Results Here we redescribe the multisegmented megacheirans Fortiforceps and Jianfengia and describe Sklerolibyon maomima gen. et sp. nov., which we place in Jianfengiidae, fam. nov. (in Megacheira, emended). We find that jianfengiids show high morphological diversity among megacheirans, both in trunk ornamentation and head anatomy, which encompasses from 2 to 4 post-frontal appendage pairs. These taxa are also characterized by elongate podomeres likely forming seven-segmented endopods, which were misinterpreted in their original descriptions. Plesiomorphic traits also clarify their connection with more ancestral taxa. The structure and position of the “great appendages” relative to likely sensory antero-medial protrusions, as well as the presence of optic peduncles and sclerites, point to an overall homology with the anterior head of radiodontans. This is confirmed by our Bayesian phylogeny, which places jianfengiids as the basalmost euarthropods, paraphyletic with other megacheirans, and in contiguity with isoxyids and radiodontans. Conclusions Sklerolibyon and other jianfengiids expand the disparity of megacheirans and suggest that the common euarthropod ancestor possessed a remarkable phenotypic variability associated with the externalized cephalon, as well as endopods that were already heptopodomerous, which differs from previous hypotheses and observations. These animals also demonstrate that the frontalmost pair of arthrodized appendage is homologous between radiodontans and megacheirans, refuting the claim that the radiodontan frontal appendages evolved into the euarthropod labrum, and questioning its protocerebral identity. This evidence based on external anatomy now constitutes a solid benchmark upon which we should address issues of homology, with the help of carefully examined palaeoneurological data.

BackgroundUnderstanding the genomic basis of phenotypic diversity can be greatly facilitated by examining adaptive radiations with hypervariable traits. In this study, we focus on a rapidly diverged species group of mormyrid electric fish in the genusParamormyrops, which are characterized by extensive phenotypic variation in electric organ discharges (EODs). The main components of EOD diversity are waveform duration, complexity and polarity. Using an RNA-sequencing based approach, we sought to identify gene expression correlates for each of these EOD waveform features by comparing 11 specimens ofParamormyropsthat exhibit variation in these features.ResultsPatterns of gene expression amongParamormyropsare highly correlated, and 3274 genes (16%) were differentially expressed. Using our most restrictive criteria, we detected 145–183 differentially expressed genes correlated with each EOD feature, with little overlap between them. The predicted functions of several of these genes are related to extracellular matrix, cation homeostasis, lipid metabolism, and cytoskeletal and sarcomeric proteins. These genes are of significant interest given the known morphological differences between electric organs that underlie differences in the EOD waveform features studied.ConclusionsIn this study, we identified plausible candidate genes that may contribute to phenotypic differences in EOD waveforms among a rapidly diverged group of mormyrid electric fish. These genes may be important targets of selection in the evolution of species-specific differences in mate-recognition signals.

BackgroundSecondary contact between closely related lineages can result in a variety of outcomes, including hybridization, depending upon the strength of reproductive barriers. By examining the extent to which different parts of the genome introgress, it is possible to infer the strength of selection and gain insight into the evolutionary trajectory of lineages. Following secondary contact approximately 8000 years ago in the Pacific Northwest, mule deer (Odocoileus hemionus hemionus) and black-tailed deer (O. h. columbianus) formed a hybrid swarm along the Cascade mountain range despite substantial differences in body size (up to two times) and habitat preference. In this study, we examined genetic population structure, extent of introgression, and selection pressures in freely interbreeding populations of mule deer and black-tailed deer using mitochondrial DNA sequences, 9 microsatellite loci, and 95 SNPs from protein-coding genes.ResultsWe observed bi-directional hybridization and classified approximately one third of the 172 individuals as hybrids, almost all of which were beyond the F1 generation. High genetic differentiation between black-tailed deer and mule deer at protein-coding genes suggests that there is positive divergent selection, though selection on these loci is relatively weak. Contrary to predictions, there was not greater selection on protein-coding genes thought to be associated with immune function and mate choice. Geographic cline analyses were consistent across genetic markers, suggesting long-term stability (over hundreds of generations), and indicated that the center of the hybrid swarm is 20-30 km to the east of the Cascades ridgeline, where there is a steep ecological transition from wet, forested habitat to dry, scrub habitat.ConclusionsOur data are consistent with a genetic boundary between mule deer and black-tailed deer that is porous but maintained by many loci under weak selection having a substantial cumulative effect. The absence of clear reproductive barriers and the consistent centering of geographic clines at a sharp ecotone suggests that ecology is a driver of hybrid swarm dynamics. Adaptive introgression in this study (and others) promotes gene flow and provides valuable insight into selection strength on specific genes and the evolutionary trajectory of hybridizing taxa.

Background Triplophysa stoliczkae is the most widespread species in the genus Triplophysa and may have originated from morphological convergence. To understand the evolutionary history of T. stoliczkae, we employed a multilocus approach to investigate the phylogenetics and the morphological evolution of T. stoliczkae on the Qinghai-Tibetan Plateau. Results All phylogenetic analyses (two mitochondrial and five nuclear loci), a genealogical sorting index and species tree inferences suggested that T. stoliczkae consists of distinct lineages that were not closest relatives. The time estimation indicated that the divergence events between “T. stoliczkae” and other Triplophysa species occurred from approximately 0.10 to 4.51 Ma. The ancestral state analyses supported the independent evolution of T. stoliczkae morphology in distinct lineages. The morphometric analysis and convergence estimates demonstrated significant phenotypic convergence among “T. stoliczkae” lineages. Conclusions Triplophysa stoliczkae includes 4 different lineages with similar morphologies. The increasingly harsh environments that have occurred since the Pliocene have driven the occurrences of scrape-feeding fish in the genus Triplophysa. Morphological adaptations associated with scrape-feeding behavior resulted in convergences and the artificial lumping of four different species in the nominal taxon T. stoliczkae. A taxonomic revision for T. stoliczkae is needed. Electronic supplementary material The online version of this article (10.1186/s12862-019-1503-3) contains supplementary material, which is available to authorized users.

BackgroundThe silkmoths and their relatives constitute the ecologically and taxonomically diverse superfamily Bombycoidea, which includes some of the most charismatic species of Lepidoptera. Despite displaying some of the most spectacular forms and ecological traits among insects, relatively little attention has been given to understanding their evolution and the drivers of their diversity.ResultsTo begin to address this problem, we created a new Bombycoidea-specific Anchored Hybrid Enrichment (AHE) probe set and sampled up to 571 loci for 117 taxa across all major lineages of the Bombycoidea, producing a well-supported phylogeny. The tree was overall consistent with prior morphological and molecular studies, although some taxa (e.g.,ArotrosSchaus) were misplaced in the Bombycidae and here formally transferred to Apatelodidae. We identified important evolutionary patterns (e.g., morphology, biogeography, and differences in speciation and extinction), and our analysis of diversification rates highlights the stark increases that exist within the Sphingidae (hawkmoths) and Saturniidae (wild silkmoths).ConclusionsWe postulate that these rate shifts are due to the well-documented bat-moth “arms race” and differences in selective pressures from insectivorous bats. The study establishes a backbone for future evolutionary, comparative, and taxonomic studies, and presents a modified DNA extraction protocol that allows Lepidoptera specimens to be readily sequenced from pinned natural history collections, succeeding in samples up to 30 years old. Our research highlights the flexibility of AHE to generate genomic data from a wide range of museum specimens, both age and preservation method, and will allow researchers to tap into the wealth of biological data residing in natural history collections around the globe.

Background Marine invertebrates are abundant and diverse on the continental shelf in Antarctica, but little is known about their parasitic counterparts. Endoparasites are especially understudied because they often possess highly modified body plans that pose problems for their identification. Asterophila, a genus of endoparasitic gastropod in the family Eulimidae, forms cysts in the arms and central discs of asteroid sea stars. There are currently four known species in this genus, one of which has been described from the Antarctic Peninsula (A. perknasteri). This study employs molecular and morphological data to investigate the diversity of Asterophila in Antarctica and explore cophylogenetic patterns between host and parasite. Results A maximum-likelihood phylogeny of Asterophila and subsequent species-delimitation analysis uncovered nine well-supported putative species, eight of which are new to science. Most Asterophila species were found on a single host species, but four species were found on multiple hosts from one or two closely related genera, showing phylogenetic conservatism of host use. Both distance-based and event-based cophylogenetic analyses uncovered a strong signal of coevolution in this system, but most associations were explained by non-cospeciation events. Discussion The prevalence of duplication and host-switching events in Asterophila and its asteroid hosts suggests that synchronous evolution may be rare even in obligate endoparasitic systems. The apparent restricted distribution of Asterophila from around the Scotia Arc may be an artefact of concentrated sampling in the area and a low obvious prevalence of infection. Given the richness of parasites on a global scale, their role in promoting host diversification, and the threat of their loss through coextinction, future work should continue to investigate parasite diversity and coevolution in vulnerable ecosystems.

Background Diversity at the Major Histocompatibility Complex (MHC) is critical to health and fitness, such that MHC genotype may predict an individual’s quality or compatibility as a competitor, ally, or mate. Moreover, because MHC products can influence the components of bodily secretions, an individual’s body odors may signal its MHC composition and influence partner identification or mate choice. Here, we investigated MHC-based signaling and recipient sensitivity by testing for odor-gene covariance and behavioral discrimination of MHC diversity and pairwise dissimilarity in a strepsirrhine primate, the ring-tailed lemur (Lemur catta). Methods First, we coupled genotyping of the MHC class II gene, DRB, with gas chromatography-mass spectrometry of genital gland secretions to investigate if functional genetic diversity is signaled by the chemical diversity of lemur scent secretions. We also assessed if the chemical similarity between individuals correlated with their MHC-DRB similarity. Next, we assessed if lemurs discriminated this chemically encoded, genetic information in opposite-sex conspecifics. Results We found that both sexes signaled overall MHC-DRB diversity and pairwise MHC-DRB similarity via genital secretions, but in a sex- and season-dependent manner. Additionally, the sexes discriminated absolute and relative MHC-DRB diversity in the genital odors of opposite-sex conspecifics, suggesting that lemur genital odors function to advertise genetic quality. Conclusions In summary, genital odors of ring-tailed lemurs provide honest information about an individual’s absolute and relative MHC quality. Complementing evidence in humans and Old World monkeys, we suggest that reliance on scent signals to communicate MHC quality may be important across the primate lineage. Electronic supplementary material The online version of this article (10.1186/s12862-019-1486-0) contains supplementary material, which is available to authorized users.

Background Allopatric speciation has played a particularly important role in archipelagic settings where populations evolve in isolation after colonizing different islands. The Indo-Australasian island realm is an unparalleled natural laboratory of biotic diversification. Here we explore how the level of earth-historic isolation has influenced genetic differentiation across the region by investigating phylogeographic patterns in the Pitta sordida species complex. Results We generated a de novo genome and compared population genomics of 29 individuals of Pitta sordida from the entire distributional range and we reconstructed phylogenetic relationship using mitogenomes, a multi-nuclear gene dataset and single nucleotide polymorphisms (SNPs). We found deep divergence between an eastern and a western group of taxa across Indo-Australasia. Within both groups we have identified major lineages that are geographically separated into Philippines, Borneo, western Sundaland, and New Guinea, respectively. Although these lineages are genetically well-differentiated, suggesting a long-term isolation, there are signatures of extensive gene flow within each lineage throughout the Pleistocene, despite the wide geographic range occupied by some of them. We found little evidence of hybridization or introgression among the studied taxa, but forsteni from Sulawesi makes an exception. This individual, belonging to the eastern clade, is genetically admixed between the western and eastern clades. Geographically this makes sense as Sulawesi is not far from Borneo that houses a population of hooded pittas that belongs to the western clade. Conclusions We found that geological vicariance events cannot explain the current genetic differentiation in the Pitta sordida species complex. Instead, the glacial-interglacial cycles may have played a major role therein. During glacials the sea level could be up to 120 m lower than today and land bridges formed within both the Sunda Shelf and the Sahul Shelf permitting dispersal of floral and faunal elements. The geographic distribution of hooded pittas shows the importance of overwater, “stepping-stone” dispersals not only to deep-sea islands, but also from one shelf to the other. The most parsimonious hypothesis is an Asian ancestral home of the Pitta sordida species complex and a colonization from west to east, probably via Wallacea. Electronic supplementary material The online version of this article (10.1186/s12862-019-1481-5) contains supplementary material, which is available to authorized users.

Background Understanding how variation in gene expression contributes to morphological diversity is a major goal in evolutionary biology. Cichlid fishes from the East African Great lakes exhibit striking diversity in trophic adaptations predicated on the functional modularity of their two sets of jaws (oral and pharyngeal). However, the transcriptional basis of this modularity is not so well understood, as no studies thus far have directly compared the expression of genes in the oral and pharyngeal jaws. Nor is it well understood how gene expression may have contributed to the parallel evolution of trophic morphologies across the replicate cichlid adaptive radiations in Lake Tanganyika, Malawi and Victoria. Results We set out to investigate the role of gene expression divergence in cichlid fishes from these three lakes adapted to herbivorous and carnivorous trophic niches. We focused on the development stage prior to the onset of exogenous feeding that is critical for understanding patterns of gene expression after oral and pharyngeal jaw skeletogenesis, anticipating environmental cues. This framework permitted us for the first time to test for signatures of gene expression underlying jaw modularity in convergent eco-morphologies across three independent adaptive radiations. We validated a set of reference genes, with stable expression between the two jaw types and across species, which can be important for future studies of gene expression in cichlid jaws. Next we found evidence of modular and non-modular gene expression between the two jaws, across different trophic niches and lakes. For instance, prdm1a, a skeletogenic gene with modular anterior-posterior expression, displayed higher pharyngeal jaw expression and modular expression pattern only in carnivorous species. Furthermore, we found the expression of genes in cichlids jaws from the youngest Lake Victoria to exhibit low modularity compared to the older lakes. Conclusion Overall, our results provide cross-species transcriptional comparisons of modularly-regulated skeletogenic genes in the two jaw types, implicating expression differences which might contribute to the formation of divergent trophic morphologies at the stage of larval independence prior to foraging. Electronic supplementary material The online version of this article (10.1186/s12862-019-1483-3) contains supplementary material, which is available to authorized users.

BackgroundOur ability to investigate processes shaping the evolutionary diversification of corals (Cnidaria: Anthozoa) is limited by a lack of understanding of species boundaries. Discerning species has been challenging due to a multitude of factors, including homoplasious and plastic morphological characters and the use of molecular markers that are either not informative or have not completely sorted. Hybridization can also blur species boundaries by leading to incongruence between morphology and genetics. We used traditional DNA barcoding and restriction-site associated DNA sequencing combined with coalescence-based and allele-frequency methods to elucidate species boundaries and simultaneously examine the potential role of hybridization in a speciose genus of octocoral,Sinularia.ResultsSpecies delimitations using two widely used DNA barcode markers,mtMutSand 28S rDNA, were incongruent with one another and with the morphospecies identifications, likely due to incomplete lineage sorting. In contrast, 12 of the 15 morphospecies examined formed well-supported monophyletic clades in both concatenated RAxML phylogenies and SNAPP species trees of >6,000 RADSeq loci. DAPC and Structure analyses also supported morphospecies assignments, but indicated the potential for two additional cryptic species. Three morphologically distinct species pairs could not, however, be distinguished genetically. ABBA-BABA tests demonstrated significant admixture between some of those species, suggesting that hybridization may confound species delimitation inSinularia.ConclusionsA genomic approach can help to guide species delimitation while simultaneously elucidating the processes generating diversity in corals. Results support the hypothesis that hybridization is an important mechanism in the evolution of Anthozoa, including octocorals, and future research should examine the contribution of this mechanism in generating diversity across the coral tree of life.

Background Gene duplication has played an important role in the evolution and domestication of flowering plants. Yet little is known about how plant duplicate genes evolve and are retained over long timescales, particularly those arising from small-scale duplication (SSD) rather than whole-genome duplication (WGD) events. Results We address this question in the Poaceae (grass) family by analyzing gene expression data from nine tissues of Brachypodium distachyon, Oryza sativa japonica (rice), and Sorghum bicolor (sorghum). Consistent with theoretical predictions, expression profiles of most grass genes are conserved after SSD, suggesting that functional conservation is the primary outcome of SSD in grasses. However, we also uncover support for widespread functional divergence, much of which occurs asymmetrically via the process of neofunctionalization. Moreover, neofunctionalization preferentially targets younger (child) duplicate gene copies, is associated with RNA-mediated duplication, and occurs quickly after duplication. Further analysis reveals that functional divergence of SSD-derived genes is positively correlated with both sequence divergence and tissue specificity in all three grass species, and particularly with anther expression in B. distachyon. Conclusions Our results suggest that SSD-derived grass genes often undergo rapid functional divergence that may be driven by natural selection on male-specific phenotypes. These observations are consistent with those in several animal species, suggesting that duplicate genes take similar evolutionary trajectories in plants and animals. Electronic supplementary material The online version of this article (10.1186/s12862-019-1415-2) contains supplementary material, which is available to authorized users.

Background A major focus of evolutionary biology is the formation of reproductive barriers leading to divergence and ultimately, speciation. Often, it is not clear whether the separation of populations is complete or if there still is ongoing gene flow in the form of rare cases of admixture, known as isolation with migration. Here, we studied the speciation of two fire ant species, Solenopsis invicta and Solenopsis richteri, both native to South America, both inadvertently introduced to North America in the early twentieth century. While the two species are known to admix in the introduced range, in the native range no hybrids were found. Results We conducted a population genomic survey of native and introduced populations of the two species using reduced representation genomic sequencing of 337 samples. Using maximum likelihood analysis over native range samples, we found no evidence of any gene flow between the species since they diverged. We estimated their time of divergence to 190,000 (100,000–350,000) generations ago. Modelling the demographic history of native and introduced S. invicta populations, we evaluated their divergence times and historic and contemporary population sizes, including the original founder population in North America, which was estimated at 26 (10–93) unrelated singly-mated queens. Conclusions We provide evidence for complete genetic isolation maintained between two invasive species in their natïve range, based, for the first time, on large scale genomic data analysis. The results lay the foundations for further studies into different stages in the formation of genetic barriers in dynamic, invasive populations. Electronic supplementary material The online version of this article (10.1186/s12862-019-1437-9) contains supplementary material, which is available to authorized users.

Background RNA interference (RNAi) related pathways provide defense against viruses and transposable elements, and have been implicated in the suppression of meiotic drive elements. Genes in these pathways often exhibit high levels of adaptive substitution, and over longer timescales show gene duplication and loss—most likely as a consequence of their role in mediating conflict with these parasites. This is particularly striking for Argonaute 2 (Ago2), which is ancestrally the key effector of antiviral RNAi in insects, but has repeatedly formed new testis-specific duplicates in the recent history of the obscura species-group of Drosophila. Results Here we take advantage of publicly available genomic and transcriptomic data to identify six further RNAi-pathway genes that have duplicated in this clade of Drosophila, and examine their evolutionary history. As seen for Ago2, we observe high levels of adaptive amino-acid substitution and changes in sex-biased expression in many of the paralogs. However, our phylogenetic analysis suggests that co-duplications of the RNAi machinery were not synchronous, and our expression analysis fails to identify consistent male-specific expression. Conclusions These results confirm that RNAi genes, including genes of the antiviral and piRNA pathways, have undergone multiple independent duplications and that their history has been particularly labile within the obscura group. However, they also suggest that the selective pressures driving these changes have not been consistent, implying that more than one selective agent may be responsible., BMC Evolutionary Biology, 19, ISSN:1471-2148

Background Recent studies have demonstrated that Bayesian species delimitation based on the multispecies coalescent model can produce inaccurate results by misinterpreting population splits as species divergences. An approach based on the genealogical divergence index (gdi) was shown to be a viable alternative, especially for delimiting allopatric populations where gene flow is low. We implemented these analyses to assess species boundaries in Southeast Asian toads, a group that is understudied and characterized by numerous unresolved species complexes. Results Multilocus phylogenetic analyses showed that deep evolutionary relationships including the genera Sigalegalephrynus, Ghatophryne, Parapelophryne, Leptophryne, Pseudobufo, Rentapia, and Phrynoides remain unresolved. Comparison of genetic divergences revealed that intraspecific divergences among allopatric populations of Pelophyrne signata (Borneo vs. Peninsular Malaysia), Ingerophrynus parvus (Peninsular Malaysia vs. Myanmar), and Leptophryne borbonica (Peninsular Malaysia, Java, Borneo, and Sumatra) are consistent with interspecific divergences of other Southeast Asian bufonid taxa. Conversely, interspecific divergences between Pelophryne guentheri/P. api, Ansonia latiffi/A. leptopus, and I. gollum/I. divergens were low (

Background Animals have a greater diversity of signalling pathways than their unicellular relatives, consistent with the evolution and expansion of these pathways occurring in parallel with the origin of animal multicellularity. However, the genomes of sponges and ctenophores – non-bilaterian basal animals – typically encode no, or far fewer, recognisable signalling ligands compared to bilaterians and cnidarians. For instance, the largest subclass of receptor tyrosine kinases (RTKs) in bilaterians, the Eph receptors (Ephs), are present in sponges and ctenophores, but their cognate ligands, the ephrins, have not yet been detected. Results Here, we use an iterative HMM analysis to identify for the first time membrane-bound ephrins in sponges and ctenophores. We also expand the number of Eph-receptor subtypes identified in these animals and in cnidarians. Both sequence and structural analyses are consistent with the Eph ligand binding domain (LBD) and the ephrin receptor binding domain (RBD) having evolved via the co-option of ancient galactose-binding (discoidin-domain)-like and monodomain cupredoxin domains, respectively. Although we did not detect a complete Eph-ephrin signalling pathway in closely-related unicellular holozoans or in other non-metazoan eukaryotes, truncated proteins with Eph receptor LBDs and ephrin RBDs are present in some choanoflagellates. Together, these results indicate that Eph-ephrin signalling was present in the last common ancestor of extant metazoans, and perhaps even in the last common ancestor of animals and choanoflagellates. Either scenario pushes the origin of Eph-ephrin signalling back much earlier than previously reported. Conclusions We propose that the Eph-LBD and ephrin-RBD, which were ancestrally localised in the cytosol, became linked to the extracellular parts of two cell surface proteins before the divergence of sponges and ctenophores from the rest of the animal kingdom. The ephrin-RBD lost the ancestral capacity to bind copper, and the Eph-LBD became linked to an ancient RTK. The identification of divergent ephrin ligands in sponges and ctenophores suggests that these ligands evolve faster than their cognate receptors. As this may be a general phenomena, we propose that the sequence-structure approach used in this study may be usefully applied to other signalling systems where no, or a small number of, ligands have been identified.

Color and pattern phenotypes have clear implications for survival and reproduction in many species. However, the mechanisms that produce this coloration are still poorly characterized, especially at the genomic level. Here we have taken a transcriptomics-based approach to elucidate the underlying genetic mechanisms affecting color and pattern in a highly polytypic poison frog. We sequenced RNA from the skin from four different color morphs during the final stage of metamorphosis and assembled a de novo transcriptome. We then investigated differential gene expression, with an emphasis on examining candidate color genes from other taxa. Overall, we found differential expression of a suite of genes that control melanogenesis, melanocyte differentiation, and melanocyte proliferation (e.g., tyrpl, lefl, leol, and mitf) as well as several differentially expressed genes involved in purine synthesis and iridophore development (e.g., arfgapl, arfgap2, airc, and gairt). Our results provide evidence that several gene networks known to affect color and pattern in vertebrates play a role in color and pattern variation in this species of poison frog.

Background Sex-biased gene expression is thought to drive the phenotypic differences in males and females in metazoans. Drosophila has served as a primary model for studying male-female differences in gene expression, and its effects on protein sequence divergence. However, the forces shaping evolution of sex-biased expression remain largely unresolved, including the roles of selection and pleiotropy. Research on sex organs in Drosophila, employing original approaches and multiple-species contrasts, provides a means to gain insights into factors shaping the turnover and magnitude (fold-bias) of sex-biased expression. Results Here, using recent RNA-seq data, we studied sex-biased gonadal expression in 10,740 protein coding sequences in four species of Drosophila, D. melanogaster, D. simulans, D. yakuba and D. ananassae (5 to 44 My divergence). Using an approach wherein we identified genes with lineage-specific transitions (LSTs) in sex-biased status (amongst testis-biased, ovary-biased and unbiased; thus, six transition types) standardized to the number of genes with the ancestral state (S-LSTs), and those with clade-wide expression bias status, we reveal several key findings. First, the six categorical types of S-LSTs in sex-bias showed disparate rates of turnover, consistent with differential selection pressures. Second, the turnover in sex-biased status was largely unrelated to cross-tissue expression breadth, suggesting pleiotropy does not restrict evolution of sex-biased expression. Third, the fold-sex-biased expression, for both testis-biased and ovary-biased genes, evolved directionally over time toward higher values, a crucial finding that could be interpreted as a selective advantage of greater sex-bias, and sexual antagonism. Fourth, in terms of protein divergence, genes with LSTs to testis-biased expression exhibited weak signals of elevated rates of evolution (than ovary-biased) in as little as 5 My, which strengthened over time. Moreover, genes with clade-wide testis-specific expression (44 My), a status not observed for any ovary-biased genes, exhibited striking acceleration of protein divergence, which was linked to low pleiotropy. Conclusions By studying LSTs and clade-wide sex-biased gonadal expression in a multi-species clade of Drosophila, we describe evidence that interspecies turnover and magnitude of sex-biased expression have been influenced by selection. Further, whilst pleiotropy was not connected to turnover in sex-biased gonadal expression, it likely explains protein sequence divergence. Electronic supplementary material The online version of this article (10.1186/s12862-019-1377-4) contains supplementary material, which is available to authorized users.

Background The plastid is a semiautonomous organelle with its own genome. Plastid genomes have been widely used as models for studying phylogeny, speciation and adaptive evolution. However, most studies focus on comparisons of plastid genome evolution at high taxonomic levels, and comparative studies of the process of plastome evolution at the infrageneric or intraspecific level remain elusive. Holcoglossum is a small genus of Orchidaceae, consisting of approximately 20 species of recent radiation. This made it an ideal group to explore the plastome mutation mode at the infrageneric or intraspecific level. Results In this paper, we reported 15 complete plastid genomes from 12 species of Holcoglossum and 1 species of Vanda. The plastid genomes of Holcoglossum have a total length range between 145 kb and 148 kb, encoding a set of 102 genes. The whole set of ndh-gene families in Holcoglossum have been truncated or pseudogenized. Hairpin inversion in the coding region of the plastid gene ycf2 has been found. Conclusions Using a comprehensive comparative plastome analysis, we found that all the indels between different individuals of the same species resulted from the copy number variation of the short repeat sequence, which may be caused by replication slippage. Annotation of tandem repeats shows that the variation introduced by tandem repeats is widespread in plastid genomes. The hairpin inversion found in the plastid gene ycf2 occurred randomly in the Orchidaceae.

Background Germ lines are the cell lineages that give rise to the sperm and eggs in animals. The germ lines first arise from primordial germ cells (PGCs) during embryogenesis: these form from either a presumed derived mode of preformed germ plasm (inheritance) or from an ancestral mechanism of inductive cell-cell signalling (induction). Numerous genes involved in germ line specification and development have been identified and functionally studied. However, little is known about the molecular evolutionary dynamics of germ line genes in metazoan model systems. Results Here, we studied the molecular evolution of germ line genes within three metazoan model systems. These include the genus Drosophila (N=34 genes, inheritance), the fellow insect Apis (N=30, induction), and their more distant relative Caenorhabditis (N=23, inheritance). Using multiple species and established phylogenies in each genus, we report that germ line genes exhibited marked variation in the constraint on protein sequence divergence (dN/dS) and codon usage bias (CUB) within each genus. Importantly, we found that de novo lineage-specific inheritance (LSI) genes in Drosophila (osk, pgc) and in Caenorhabditis (pie-1, pgl-1), which are essential to germ plasm functions under the derived inheritance mode, displayed rapid protein sequence divergence relative to the other germ line genes within each respective genus. We show this may reflect the evolution of specialized germ plasm functions and/or low pleiotropy of LSI genes, features not shared with other germ line genes. In addition, we observed that the relative ranking of dN/dS and of CUB between genera were each more strongly correlated between Drosophila and Caenorhabditis, from different phyla, than between Drosophila and its insect relative Apis, suggesting taxonomic differences in how germ line genes have evolved. Conclusions Taken together, the present results advance our understanding of the evolution of animal germ line genes within three well-known metazoan models. Further, the findings provide insights to the molecular evolution of germ line genes with respect to LSI status, pleiotropy, adaptive evolution as well as PGC-specification mode. Electronic supplementary material The online version of this article (10.1186/s12862-019-1363-x) contains supplementary material, which is available to authorized users.

Background The C-C motif chemokine ligand 16 (CCL16) is a potent pro-inflammatory chemokine and a chemoattractant for monocytes and lymphocytes. In normal plasma, it is present at high concentrations and elicits its effects on cells by interacting with cell surface chemokine receptors. In the European rabbit and in rodents such as mouse, rat and guinea pig, CCL16 was identified as a pseudogene, while in the thirteen-lined ground squirrel it appears to be potentially functional. To gain insight into the evolution of this gene in the superorder Glires (rodents and lagomorphs), we amplified the CCL16 gene from eleven Leporidae and seven Ochotonidae species. Results We compared our sequences with CCL16 sequences of twelve rodent species retrieved from public databases. The data show that for all leporid species studied CCL16 is a pseudogene. This is primarily due to mutations at the canonical Cys Cys motif, creating either premature stop codons, or disrupting amino acid replacements. In the Mexican cottontail, CCL16 is pseudogenized due to a frameshift deletion. Additionally, in the exon 1 (signal peptide), there are frameshift deletions present in all leporids studied. In contrast, in Ochotona species, CCL16 is potentially functional, except for an allele in Hoffmann’s pika. In rodents, CCL16 is functional in a number of species, but patterns of pseudogenization similar to those observed in lagomorphs also exist. Conclusions Our results suggest that while functional in the Glires ancestor, CCL16 underwent pseudogenization in some species. This process occurred stochastically or in specific lineages at different moments in the evolution of Glires. These observations suggest that the CCL16 had different evolutionary constrains in the Glires group that could be associated with the CCL16 biological function. Electronic supplementary material The online version of this article (10.1186/s12862-019-1390-7) contains supplementary material, which is available to authorized users.

Background Cilia and flagella are complex cellular structures thought to have first evolved in a last ciliated eukaryotic ancestor due to the conserved 9 + 2 microtubule doublet structure of the axoneme and associated proteins. The Tektin family of coiled-coil domain containing proteins was previously identified in cilia of organisms as diverse as green algae and sea urchin. While studies have shown that some Tektins are necessary for ciliary function, there has been no comprehensive phylogenetic survey of tektin genes. To fill this gap, we sampled tektin sequences broadly among metazoan and unicellular lineages in order to determine how the tektin gene complements evolved in over 100 different extant species. Results Using Bayesian and Maximum Likelihood analyses, we have ascertained with high confidence that all metazoan tektins arose from a single ancestral tektin gene in the last common ancestor of metazoans and choanoflagellates. Gene duplications gave rise to two tektin genes in the metazoan ancestor, and a subsequent expansion to three and four tektin genes in early bilaterian ancestors. While all four tektin genes remained highly conserved in most deuterostome and spiralian species surveyed, most tektin genes in ecdysozoans are highly derived with extensive gene loss in several lineages including nematodes and some crustaceans. In addition, while tektin-1, − 2, and − 4 have remained as single copy genes in most lineages, tektin-3/5 has been duplicated independently several times, notably at the base of the spiralian, vertebrate and hymenopteran (Ecdysozoa) clades. Conclusions We provide a solid description of tektin evolution supporting one, two, three, and four ancestral tektin genes in a holozoan, metazoan, bilaterian, and nephrozoan ancestor, respectively. The isolated presence of tektin in a cryptophyte and a chlorophyte branch invokes events of horizontal gene transfer, and that the last common ciliated eukaryotic ancestor lacked a tektin gene. Reconstructing the evolutionary history of the tektin complement in each extant metazoan species enabled us to pinpoint lineage specific expansions and losses. Our analysis will help to direct future studies on Tektin function, and how gain and loss of tektin genes might have contributed to the evolution of various types of cilia and flagella. Electronic supplementary material The online version of this article (10.1186/s12862-019-1360-0) contains supplementary material, which is available to authorized users.

Background Morphological diversity among closely related animals can be the result of differing growth patterns. The Australian radiation of agamid lizards (Amphibolurinae) exhibits great ecological and morphological diversity, which they have achieved on a continent-wide scale, in a relatively short period of time (30 million years). Amphibolurines therefore make an ideal study group for examining ontogenetic allometry. We used two-dimensional landmark-based geometric morphometric methods to characterise the postnatal growth patterns in cranial shape of 18 species of amphibolurine lizards and investigate the associations between cranial morphology, and life habit and phylogeny. Results For most amphibolurine species, juveniles share a similar cranial phenotype, but by adulthood crania are more disparate in shape and occupy different sub-spaces of the total shape space. To achieve this disparity, crania do not follow a common post-natal growth pattern; there are differences among species in both the direction and magnitude of change in morphospace. We found that these growth patterns among the amphibolurines are significantly associated with ecological life habits. The clade Ctenophorus includes species that undergo small magnitudes of shape change during growth. They have dorsoventrally deep, blunt-snouted skulls (associated with terrestrial lifestyles), and also dorsoventrally shallow skulls (associated with saxicolous lifestyles). The sister clade to Ctenophorus, which includes the bearded dragon (Pogona), frill-neck lizard (Chlamydosaurus), and long-nosed dragon (Gowidon), exhibit broad and robust post-orbital regions and differing snout lengths (mainly associated with scansorial lifestyles). Conclusions Australian agamids show great variability in the timing of development and divergence of growth trajectories which results in a diversity of adult cranial shapes. Phylogenetic signal in cranial morphology appears to be largely overwritten by signals that reflect life habit. This knowledge about growth patterns and skull shape diversity in agamid lizards will be valuable for placing phylogenetic, functional and ecological studies in a morphological context. Electronic supplementary material The online version of this article (10.1186/s12862-018-1335-6) contains supplementary material, which is available to authorized users.

Background Establishing the divergence times of groups of organisms is a major goal of evolutionary biology. This is especially challenging for microbial lineages due to the near-absence of preserved physical evidence (diagnostic body fossils or geochemical biomarkers). Horizontal gene transfer (HGT) can serve as a temporal scaffold between microbial groups and other fossil-calibrated clades, potentially improving these estimates. Specifically, HGT to or from organisms with fossil-calibrated age estimates can propagate these constraints to additional groups that lack fossils. While HGT is common between lineages, only a small subset of HGT events are potentially informative for dating microbial groups. Results Constrained by published fossil-calibrated studies of fungal evolution, molecular clock analyses show that multiple clades of Bacteria likely acquired chitinase homologs via HGT during the very late Neoproterozoic into the early Paleozoic. These results also show that, following these HGT events, recipient terrestrial bacterial clades likely diversified ~ 300–500 million years ago, consistent with established timescales of arthropod and plant terrestrialization. Conclusions We conclude that these age estimates are broadly consistent with the dispersal of chitinase genes throughout the microbial world in direct response to the evolution and ecological expansion of detrital-chitin producing groups. The convergence of multiple lines of evidence demonstrates the utility of HGT-based dating methods in microbial evolution. The pattern of inheritance of chitinase genes in multiple terrestrial bacterial lineages via HGT processes suggests that these genes, and possibly other genes encoding substrate-specific enzymes, can serve as a “standard candle” for dating microbial lineages across the Tree of Life. Electronic supplementary material The online version of this article (10.1186/s12862-019-1357-8) contains supplementary material, which is available to authorized users.

Background Hox genes are key elements in patterning animal development. They are renowned for their, often, clustered organisation in the genome, with supposed mechanistic links between the organisation of the genes and their expression. The widespread distribution and comparable functions of Hox genes across the animals has led to them being a major study system for comparing the molecular bases for construction and divergence of animal morphologies. Echinoderms (including sea urchins, sea stars, sea cucumbers, feather stars and brittle stars) possess one of the most unusual body plans in the animal kingdom with pronounced pentameral symmetry in the adults. Consequently, much interest has focused on their development, evolution and the role of the Hox genes in these processes. In this context, the organisation of echinoderm Hox gene clusters is distinctive. Within the classificatory system of Duboule, echinoderms constitute one of the clearest examples of Disorganized (D) clusters (i.e. intact clusters but with a gene order or orientation rearranged relative to the ancestral state). Results Here we describe two Hox genes (Hox11/13d and e) that have been overlooked in most previous work and have not been considered in reconstructions of echinoderm Hox complements and cluster organisation. The two genes are related to Posterior Hox genes and are present in all classes of echinoderm. Importantly, they do not reside in the Hox cluster of any species for which genomic linkage data is available. Conclusion Incorporating the two neglected Posterior Hox genes into assessments of echinoderm Hox gene complements and organisation shows that these animals in fact have Split (S) Hox clusters rather than simply Disorganized (D) clusters within the Duboule classification scheme. This then has implications for how these genes are likely regulated, with them no longer covered by any potential long-range Hox cluster-wide, or multigenic sub-cluster, regulatory mechanisms. Electronic supplementary material The online version of this article (10.1186/s12862-018-1307-x) contains supplementary material, which is available to authorized users.

Background Variation in the non-coding regions of Y-chromosomes have been shown to influence gene regulation throughout the genome in some systems; a phenomenon termed Y-linked regulatory variation (YRV). This type of sex-specific genetic variance could have important implications for the evolution of male and female traits. If YRV contributes to the additive genetic variation of an autosomally coded trait shared between the sexes (e.g. body size), then selection could facilitate sexually dimorphic evolution via the Y-chromosome. In contrast, if YRV is entirely non-additive (i.e. interacts epistatically with other chromosomes), then Y-chromosomes could constrain trait evolution in both sexes whenever they are selected for the same trait value. The ability for this phenomenon to influence such fundamental evolutionary dynamics remains unexplored. Results Here we address the evolutionary contribution of Y-linked variance by selecting for improved male geotaxis in populations possessing multiple Y-chromosomes (i.e. possessed Y-linked additive and/or epistatic variation) or a single Y-chromosome variant (i.e. possessed no Y-linked variation). We found that males from populations possessing Y-linked variation did not significantly respond to selection; however, males from populations with no Y-linked variation did respond. These patterns suggest the presence of a large quantity of Y-linked epistatic variance in the multi-Y population that dramatically slowed its response. Conclusions Our results imply that YRV is unlikely to facilitate the evolution of sexually dimorphic traits (at least for the trait examined here), but can interfere with the rate of trait evolution in both males and females. This result could have real biological implications as it suggests that YRV can affect how quickly a population responds to new selective pressures (e.g. invasive species, novel pathogens, or climate change). Considering that YRV influences hundreds of genes and is likely typical of other independently-evolved hemizygous chromosomes, YRV-like phenomena may represent common and significant costs to hemizygous sex determination.

Background Transcriptome sequencing has become a method of choice for evolutionary studies in microbial eukaryotes due to low cost and minimal sample requirements. Transcriptome data has been extensively used in phylogenomic studies to infer ancient evolutionary histories. However, its utility in studying cryptic species diversity is not well explored. An empirical investigation was conducted to test the applicability of transcriptome data in resolving two major types of discordances at lower taxonomic levels. These include cases where species have the same morphology but different genetics (cryptic species) and species of different morphologies but have the same genetics. We built a species comparison bioinformatic pipeline that takes into account the nature of transcriptome data in amoeboid microbes exemplifying such discordances. Result Our analyses of known or suspected cryptic species yielded consistent results regardless of the methods of culturing, RNA collection or sequencing. Over 95% of the single copy genes analyzed in samples of the same species sequenced using different methods and cryptic species had intra- and interspecific divergences below 2%. Only a minority of groups (2.91–4.87%) had high distances exceeding 2% in these taxa, which was likely caused by low data quality. This pattern was also observed in suspected genetically similar species with different morphologies. Transcriptome data consistently delineated all taxa above species level, including cryptically diverse species. Using our approach we were able to resolve cryptic species problems, uncover misidentification and discover new species. We also identified several potential barcode markers with varying evolutionary rates that can be used in lineages with different evolutionary histories. Conclusion Our findings demonstrate that transcriptome data is appropriate for understanding cryptic species diversity in microbial eukaryotes. Electronic supplementary material The online version of this article (10.1186/s12862-018-1283-1) contains supplementary material, which is available to authorized users.

Background The males of some sepsid species (Sepsidae: Diptera) have abdominal appendages that are remarkable in several ways. They are sexually dimorphic, have a complex evolutionary history of gain and loss, and can be jointed and thus highly mobile. The sternite brushes are used extensively in complex courtship behaviors that differ considerably between species and during mating. The abdominal appendages have a novel developmental pathway developing from histoblast nests rather than imaginal discs. Results We focus on the evolution of cell number, nest area, and segment length in both sexes to understand how this tissue relates to the formation of novel abdominal appendages. We map histoblast nest size of wandering-phase larvae of 17 species across 10 genera to a phylogenetic tree of Sepsidae and demonstrate that abdominal appendages require significant increases of histoblast nest size and cell number in most species while one species produces small appendages even without such modifications. In species with particularly large appendages, not only the nests on the fourth, but nests in neighboring segments are enlarged (Themira biloba, Themira putris). The loss of abdominal appendages corresponds to the loss of an enlarged fourth histoblast nest, although one species showed an exception to this pattern. One species that constitutes an independent origin of abdominal appendages (Perochaeta dikowi) uses an unusual developmental mechanism in that the histoblast nest sizes are not sexually dimorphic. Conclusions The surprisingly high diversity in histoblast size and degree of sexual dimorphism suggests that the developmental mechanism used for abdominal appendage formation in sepsids is highly adaptable. The presence of appendages usually correlate with increased histoblast cell number and in most cases appendage loss results in a return to ancestral histoblast morphology. However, we also identify several exceptions that indicate the abdominal appendages have a malleable developmental origin that is responsive to selection. Electronic supplementary material The online version of this article (10.1186/s12862-018-1265-3) contains supplementary material, which is available to authorized users.

Background Previous molecular studies on the phylogeny and classification of clupeocephalan fishes revealed numerous new taxonomic entities. For re-analysing these taxa, we perform target gene capturing and subsequent next generation sequencing of putative ortholog exons of major clupeocephalan lineages. Sequence information for the RNA bait design was derived from publicly available genomes of bony fishes. Newly acquired sequence data comprising > 800 exon sequences was subsequently used for phylogenetic reconstructions. Results Our results support monophyletic Otomorpha comprising Alepocephaliformes. Within Ostariophysi, Gonorynchiformes are sister to a clade comprising Cypriniformes, Characiformes, Siluriformes and Gymnotiformes, where the interrelationships of Characiformes, Siluriformes and Gymnotiformes remain enigmatic. Euteleosts comprise four major clades: Lepidogalaxiiformes, Protacanthopterygii, Stomiatii, and Galaxiiformes plus Neoteleostei. The monotypic Lepidogalaxiiformes form the sister-group to all remaining euteleosts. Protacanthopterygii, comprising Argentini-, Esoci- and Salmoniformes, is sister to Stomiatii (Osmeriformes and Stomiatiformes) and Galaxiiformes plus Neoteleostei. Conclusions Several proposed monophyla defined by morphological apomorphies within the Clupeocephalan phylogeny are confirmed by the phylogenetic estimates presented herein. However, other morphologically described groups cannot be reconciled with molecular phylogenies. Thus, numerous morphological apomoprhies of supposed monophyla are called into question. The interpretation of suggested morphological synapomorphies of otomorph fishes is strongly affected by the inclusion of deep-sea inhabiting, and to that effect morphologically adapted Alepocephaliformes. Our revision of these potential synapomorphies, in the context that Alepocephaliformes are otomorph fishes, reveals that only a single character of the total nine characters proposed as synapomorphic for the group is clearly valid for all otomorphs. Three further characters remain possible apomorphies since their status remains unclear in the deep-sea adapted Alepocephaliformes showing developmental lag and lacking a swim bladder. Further, our analysis places Galaxiiformes as sister group to neoteleosts, which contradicts some previous molecular phylogenetic studies. This needs further investigation from a morphological perspective, as suggested synapomophies for several euteleostean lineages are challenged or still lacking. For the verification of results presented herein, a denser phylogenomic-level taxon sampling should be applied. Electronic supplementary material The online version of this article (10.1186/s12862-018-1267-1) contains supplementary material, which is available to authorized users.

BackgroundParasitoidism, a specialized life strategy in which a parasite eventually kills its host, is frequently found within the insect order Hymenoptera (wasps, ants and bees). A parasitoid lifestyle is one of two dominant life strategies within the hymenopteran superfamily Cynipoidea, with the other being an unusual plant-feeding behavior known as galling. Less commonly, cynipoid wasps exhibit inquilinism, a strategy where some species have adapted to usurp other species’ galls instead of inducing their own. Using a phylogenomic data set of ultraconserved elements from nearly all lineages of Cynipoidea, we here generate a robust phylogenetic framework and timescale to understand cynipoid systematics and the evolution of these life histories.ResultsOur reconstructed evolutionary history for Cynipoidea differs considerably from previous hypotheses. Rooting our analyses with non-cynipoid outgroups, the Paraulacini, a group of inquilines, emerged as sister-group to the rest of Cynipoidea, rendering the gall wasp family Cynipidae paraphyletic. The families Ibaliidae and Liopteridae, long considered archaic and early-branching parasitoid lineages, were found nested well within the Cynipoidea as sister-group to the parasitoid Figitidae. Cynipoidea originated in the early Jurassic around 190 Ma. Either inquilinism or parasitoidism is suggested as the ancestral and dominant strategy throughout the early evolution of cynipoids, depending on whether a simple (three states: parasitoidism, inquilinism and galling) or more complex (seven states: parasitoidism, inquilinism and galling split by host use) model is employed.ConclusionsOur study has significant impact on understanding cynipoid evolution and highlights the importance of adequate outgroup sampling. We discuss the evolutionary timescale of the superfamily in relation to their insect hosts and host plants, and outline how phytophagous galling behavior may have evolved from entomophagous, parasitoid cynipoids. Our study has established the framework for further physiological and comparative genomic work between gall-making, inquiline and parasitoid lineages, which could also have significant implications for the evolution of diverse life histories in other Hymenoptera.

BackgroundThe Drosophilidae family is traditionally divided into two subfamilies: Drosophilinae and Steganinae. This division is based on morphological characters, and the two subfamilies have been treated as monophyletic in most of the literature, but some molecular phylogenies have suggested Steganinae to be paraphyletic. To test the paraphyletic-Steganinae hypothesis, here, we used genomic sequences of eight Drosophilidae (three Steganinae and five Drosophilinae) and two Ephydridae (outgroup) species and inferred the phylogeny for the group based on a dataset of 1,028 orthologous genes present in all species (> 1,000,000 bp). This dataset includes three genera that broke the monophyly of the subfamilies in previous works. To investigate possible biases introduced by small sample sizes and automatic gene annotation, we used the same methods to infer species trees from a set of 10 manually annotated genes that are commonly used in phylogenetics.ResultsMost of the 1,028 gene trees depicted Steganinae as paraphyletic with distinct topologies, but the most common topology depicted it as monophyletic (43.7% of the gene trees). Despite the high levels of gene tree heterogeneity observed, species tree inference in ASTRAL, in PhyloNet, and with the concatenation approach strongly supported the monophyly of both subfamilies for the 1,028-gene dataset. However, when using the concatenation approach to infer a species tree from the smaller set of 10 genes, we recovered Steganinae as a paraphyletic group. The pattern of gene tree heterogeneity was asymmetrical and thus could not be explained solely by incomplete lineage sorting (ILS).ConclusionsSteganinae was clearly a monophyletic group in the dataset that we analyzed. In addition to ILS, gene tree discordance was possibly the result of introgression, suggesting complex branching processes during the early evolution of Drosophilidae with short speciation intervals and gene flow. Our study highlights the importance of genomic data in elucidating contentious phylogenetic relationships and suggests that phylogenetic inference for drosophilids based on small molecular datasets should be performed cautiously. Finally, we suggest an approach for the correction and cleaning of BUSCO-derived genomic datasets that will be useful to other researchers planning to use this tool for phylogenomic studies.

Background Ecdysozoa are the moulting protostomes, including arthropods, tardigrades, and nematodes. Both the molecular and fossil records indicate that Ecdysozoa is an ancient group originating in the terminal Proterozoic, and exceptional fossil biotas show their dominance and diversity at the beginning of the Phanerozoic. However, the nature of the ecdysozoan common ancestor has been difficult to ascertain due to the extreme morphological diversity of extant Ecdysozoa, and the lack of early diverging taxa in ancient fossil biotas. Results Here we re-describe Acosmia maotiania from the early Cambrian Chengjiang Biota of Yunnan Province, China and assign it to stem group Ecdysozoa. Acosmia features a two-part body, with an anterior proboscis bearing a terminal mouth and muscular pharynx, and a posterior annulated trunk with a through gut. Morphological phylogenetic analyses of the protostomes using parsimony, maximum likelihood and Bayesian inference, with coding informed by published experimental decay studies, each placed Acosmia as sister taxon to Cycloneuralia + Panarthropoda—i.e. stem group Ecdysozoa. Ancestral state probabilities were calculated for key ecdysozoan nodes, in order to test characters inferred from fossils to be ancestral for Ecdysozoa. Results support an ancestor of crown group ecdysozoans sharing an annulated vermiform body with a terminal mouth like Acosmia, but also possessing the pharyngeal armature and circumoral structures characteristic of Cambrian cycloneuralians and lobopodians. Conclusions Acosmia is the first taxon placed in the ecdysozoan stem group and provides a constraint to test hypotheses on the early evolution of Ecdysozoa. Our study suggests acquisition of pharyngeal armature, and therefore a change in feeding strategy (e.g. predation), may have characterised the origin and radiation of crown group ecdysozoans from Acosmia-like ancestors.

Background Deep mitochondrial divergences were observed in Scandinavian populations of the terrestrial to semi-aquatic annelid Fridericia magna (Clitellata: Enchytraeidae). This raised the need for testing whether the taxon is a single species or a complex of cryptic species. Results A total of 62 specimens from 38 localities were included in the study, 44 of which were used for species delimitation. First, the 44 specimens were divided into clusters using ABGD (Automatic Barcode Gap Discovery) on two datasets, consisting of sequences of the mitochondrial markers COI and 16S. For each dataset, the worms were divided into six not completely congruent clusters. When they were combined, a maximum of seven clusters, or species hypotheses, were obtained, and the seven clusters were used as input in downstream analyses. We tested these hypotheses by constructing haplowebs for two nuclear markers, H3 and ITS, and in both haplowebs the specimens appeared as a single species. Multi-locus species delimitation analyses performed with the Bayesian BPP program also mainly supported a single species. Furthermore, no apparent morphological differences were found between the clusters. Two of the clusters were partially separated from each other and the other clusters, but not strongly enough to consider them as separate species. All 62 specimens were used to visualise the Scandinavian distribution, of the species, and to compare with published COI data from other Fridericia species. Conclusion We show that the morphospecies Fridericia magna is a single species, harbouring several distinct mitochondrial clusters. There is partial genetic separation between some of them, which may be interpreted as incipient speciation. The study shows the importance of rigorous species delimitation using several independent markers when deep mitochondrial divergences might give the false impression of cryptic speciation.

Background For decades, codon usage has been used as a measure of adaptation for translational efficiency and translation accuracy of a gene’s coding sequence. These patterns of codon usage reflect both the selective and mutational environment in which the coding sequences evolved. Over this same period, gene transfer between lineages has become widely recognized as an important biological phenomenon. Nevertheless, most studies of codon usage implicitly assume that all genes within a genome evolved under the same selective and mutational environment, an assumption violated when introgression occurs. In order to better understand the effects of introgression on codon usage patterns and vice versa, we examine the patterns of codon usage in Lachancea kluyveri, a yeast which has experienced a large introgression. We quantify the effects of mutation bias and selection for translation efficiency on the codon usage pattern of the endogenous and introgressed exogenous genes using a Bayesian mixture model, ROC SEMPPR, which is built on mechanistic assumptions about protein synthesis and grounded in population genetics. Results We find substantial differences in codon usage between the endogenous and exogenous genes, and show that these differences can be largely attributed to differences in mutation bias favoring A/T ending codons in the endogenous genes while favoring C/G ending codons in the exogenous genes. Recognizing the two different signatures of mutation bias and selection improves our ability to predict protein synthesis rate by 42% and allowed us to accurately assess the decaying signal of endogenous codon mutation and preferences. In addition, using our estimates of mutation bias and selection, we identify Eremothecium gossypii as the closest relative to the exogenous genes, providing an alternative hypothesis about the origin of the exogenous genes, estimate that the introgression occurred ∼6×108 generation ago, and estimate its historic and current selection against mismatched codon usage. Conclusions Our work illustrates how mechanistic, population genetic models like ROC SEMPPR can separate the effects of mutation and selection on codon usage and provide quantitative estimates from sequence data.

Background The Caribbean offers a unique opportunity to study evolutionary dynamics in insular mammals. However, the recent extinction of most Caribbean non-volant mammals has obstructed evolutionary studies, and poor DNA preservation associated with tropical environments means that very few ancient DNA sequences are available for extinct vertebrates known from the region’s Holocene subfossil record. The endemic Caribbean eulipotyphlan family Nesophontidae (“island-shrews”) became extinct ~ 500 years ago, and the taxonomic validity of many Nesophontes species and their wider evolutionary dynamics remain unclear. Here we use both morphometric and palaeogenomic methods to clarify the status and evolutionary history of Nesophontes species from Hispaniola, the second-largest Caribbean island. Results Principal component analysis of 65 Nesophontes mandibles from late Quaternary fossil sites across Hispaniola identified three non-overlapping morphometric clusters, providing statistical support for the existence of three size-differentiated Hispaniolan Nesophontes species. We were also able to extract and sequence ancient DNA from a ~ 750-year-old specimen of Nesophontes zamicrus, the smallest non-volant Caribbean mammal, including a whole-mitochondrial genome and partial nuclear genes. Nesophontes paramicrus (39-47 g) and N. zamicrus (~ 10 g) diverged recently during the Middle Pleistocene (mean estimated divergence = 0.699 Ma), comparable to the youngest species splits in Eulipotyphla and other mammal groups. Pairwise genetic distance values for N. paramicrus and N. zamicrus based on mitochondrial and nuclear genes are low, but fall within the range of comparative pairwise data for extant eulipotyphlan species-pairs. Conclusions Our combined morphometric and palaeogenomic analyses provide evidence for multiple co-occurring species and rapid body size evolution in Hispaniolan Nesophontes, in contrast to patterns of genetic and morphometric differentiation seen in Hispaniola’s extant non-volant land mammals. Different components of Hispaniola’s mammal fauna have therefore exhibited drastically different rates of morphological evolution. Morphological evolution in Nesophontes is also rapid compared to patterns across the Eulipotyphla, and our study provides an important new example of rapid body size change in a small-bodied insular vertebrate lineage. The Caribbean was a hotspot for evolutionary diversification as well as preserving ancient biodiversity, and studying the surviving representatives of its mammal fauna is insufficient to reveal the evolutionary patterns and processes that generated regional diversity.

Background Hypotrichia are a group with the most complex morphology and morphogenesis within the ciliated protists. The classification of Gastrostyla-like species, a taxonomically difficult group of hypotrichs with a common ventral cirral pattern but various dorsal and ontogenetic patterns, is poorly understood. Hence, systematic relationships within this group and with other taxa in the subclass Hypotrichia remain unresolved. Results 18S rRNA gene sequence of a new Gastrostyla-like taxon was obtained. Phylogenetic analyses based on the 18S rRNA gene sequences indicate that this ciliate represents a new genus that is closely related to Heterourosomoida and Kleinstyla within the oxytrichid clade of the Hypotrichia. However, the position of this cluster remains unresolved. All three genera deviate from the typical oxytrichids by their incomplete (or lack of) dorsal kinety fragmentation during morphogenesis. Morphology and morphogenesis of this newly discovered form, Heterogastrostyla salina nov. gen., nov. spec., are described. Heterogastrostyla nov. gen., is characterised as follows: more than 18 fronto-ventral-transverse cirri, cirral anlagen V and VI develop pretransverse cirri, and dorsal ciliature in Urosomoida-like pattern. Conclusions Similar to the CEUU-hypothesis about convergent evolution of urostylids and uroleptids, we speculate that the shared ventral cirral patterns of Gastrostyla-like taxa might have resulted from convergent evolution.

Background Evolution in one selective environment often latently generates phenotypic change that is manifested only later in different environments, but the complexity of behavior important to fitness in the original environment might influence the character of such latent-phenotype evolution. Using Myxococcus xanthus, a bacterium possessing two motility systems differing in effectiveness on hard vs. soft surfaces, we test (i) whether and how evolution while swarming on one surface—the selective surface—latently alters motility on the alternative surface type and (ii) whether patterns of such latent-phenotype evolution depend on the complexity of ancestral motility, specific ancestral motility genotypes and/or the selective surface of evolution. We analysze an experiment in which populations established from three ancestral genotypes—one with both motility systems intact and two others with one system debilitated—evolved while swarming across either hard or soft agar in six evolutionary treatments. We then compare motility-phenotype patterns across selective vs. alternative surface types. Results Latent motility evolution was pervasive but varied in character as a function of the presence of one or two functional motility systems and, for some individual-treatment comparisons, the specific ancestral genotype and/or selective surface. Swarming rates on alternative vs. selective surfaces were positively correlated generally among populations with one functional motility system but not among those with two. This suggests that opportunities for pleiotropy and epistasis generated by increased genetic complexity underlying behavior can alter the character of latent-phenotype evolution. No tradeoff between motility performance across surface types was detected in the dual-system treatments, even after adaptation on a surface on which one motility system dominates strongly over the other in driving movement, but latent-phenotype evolution was instead idiosyncratic in these treatments. We further find that the magnitude of stochastic diversification at alternative-surface swarming among replicate populations greatly exceeded diversification of selective-surface swarming within some treatments and varied across treatments. Conclusion Collectively, our results suggest that increases in the genetic and mechanistic complexity of behavior can increase the complexity of latent-phenotype evolution outcomes and illustrate that diversification manifested during evolution in one environment can be augmented greatly by diversification of latent phenotypes manifested later., BMC Evolutionary Biology, 20 (1), ISSN:1471-2148

BackgroundAdaptive radiations are characterized by extreme and/or iterative phenotypic divergence; however, such variation does not accumulate evenly across an organism. Instead, it is often partitioned into sub-units, or modules, which can differentially respond to selection. While it is recognized that changing the pattern of modularity or the strength of covariation (integration) can influence the range or rate of morphological evolution, the relationship between shape variation and covariation remains unclear. For example, it is possible that rapid phenotypic change requires concomitant changes to the underlying covariance structure. Alternatively, repeated shifts between phenotypic states may be facilitated by a conserved covariance structure. Distinguishing between these scenarios will contribute to a better understanding of the factors that shape biodiversity. Here, we explore these questions using a diverse Lake Malawi cichlid species complex,Tropheops,that appears to partition habitat by depth.ResultsWe construct a phylogeny ofTropheopspopulations and use 3D geometric morphometrics to assess the shape of four bones involved in feeding (mandible, pharyngeal jaw, maxilla, pre-maxilla) in populations that inhabit deep versus shallow habitats. We next test numerous modularity hypotheses to understand whether fish at different depths are characterized by conserved or divergent patterns of modularity. We further examine rates of morphological evolution and disparity between habitats and among modules. Finally, we raise a singleTropheopsspecies in environments mimicking deep or shallow habitats to discover whether plasticity can replicate the pattern of morphology, disparity, or modularity observed in natural populations.ConclusionsOur data support the hypothesis that conserved patterns of modularity permit the evolution of divergent morphologies and may facilitate the repeated transitions between habitats. In addition, we find the lab-reared populations replicate many trends in the natural populations, which suggests that plasticity may be an important force in initiating depth transitions, priming the feeding apparatus for evolutionary change.

Background Mygalomorph spiders represent a diverse, yet understudied lineage for which genomic level data has only recently become accessible through high-throughput genomic and transcriptomic sequencing methods. The Aptostichus atomarius species complex (family Euctenizidae) includes two coastal dune endemic members, each with inland sister species – affording exploration of dune adaptation associated patterns at the transcriptomic level. We apply an RNAseq approach to examine gene family conservation across the species complex and test for patterns of positive selection along branches leading to dune endemic species. Results An average of ~ 44,000 contigs were assembled for eight spiders representing dune (n = 2), inland (n = 4), and atomarius species complex outgroup taxa (n = 2). Transcriptomes were estimated to be 64% complete on average with 77 spider reference orthologs missing from all taxa. Over 18,000 orthologous gene clusters were identified within the atomarius complex members, > 5000 were detected in all species, and ~ 4700 were shared between species complex members and outgroup Aptostichus species. Gene family analysis with the FUSTr pipeline identified 47 gene families appearing to be under selection in the atomarius ingroup; four of the five top clusters include sequences strongly resembling other arthropod venom peptides. The COATS pipeline identified six gene clusters under positive selection on branches leading to dune species, three of which reflected the preferred species tree. Genes under selection were identified as Cytochrome P450 2c15 (also recovered in the FUSTr analysis), Niemann 2 Pick C1-like, and Kainate 2 isoform X1. Conclusions We have generated eight draft transcriptomes for a closely related and ecologically diverse group of trapdoor spiders, identifying venom gene families potentially under selection across the Aptostichus atomarius complex and chemosensory-associated gene families under selection in dune endemic lineages.

BackgroundAustropotamobius torrentiumis a freshwater crayfish species native to central and south-eastern Europe, with an intricate evolutionary history and the highest genetic diversity recorded in the northern-central Dinarides (NCD). Its populations are facing declines, both in number and size across its entire range. By extanding current knowledge on the genetic diversity of this species, we aim to assist conservation programmes. Multigene phylogenetic analyses were performed using different divergence time estimates based on mitochondrial and, for the first time, nuclear DNA markers on the largest data set analysed so far. In order to reassess taxonomic relationships within this species we applied several species delimitation methods and studied the meristic characters with the intention of finding features that would clearly separate stone crayfish belonging to different phylogroups.ResultsOur results confirmed the existence of high genetic diversity withinA. torrentium, maintained in divergent phylogroups which have their own evolutionary dynamics. A new phylogroup in the Kordun region belonging to NCD has also been discovered. Due to the incongruence between implemented species delimitation approaches and the lack of any morphological characters conserved within lineages, we are of the opinion that phylogroups recovered on mitochondrial and nuclear DNA are cryptic subspecies and distinct evolutionary significant units.ConclusionsGeographically and genetically isolated phylogroups represent the evolutionary legacy ofA. torrentiumand are highly relevant for conservation due to their evolutionary distinctiveness and restricted distribution.