The Impact of the Tree Prior on Molecular Dating of Data Sets Containing a Mixture of Inter- and Intraspecies Sampling

In Bayesian phylogenetic analyses of genetic data, prior probability distributions need to be specified for the model parameters, including the tree. When Bayesian methods are used for molecular dating, available tree priors include those designed for species-level data, such as the pure-birth and birth-death priors, and coalescent-based priors designed for population-level data. However, molecular dating methods are frequently applied to data sets that include multiple individuals across multiple species. Such data sets violate the assumptions of both the speciation and coalescent-based tree priors, making it unclear which should be chosen and whether this choice can affect the estimation of node times. To investigate this problem, we used a simulation approach to produce data sets with different proportions of within- and between-species sampling under the multispecies coalescent model. These data sets were then analyzed under pure-birth, birth-death, constant-size coalescent, and skyline coalescent tree priors. We also explored the ability of Bayesian model testing to select the best-performing priors. We confirmed the applicability of our results to empirical data sets from cetaceans, phocids, and coregonid whitefish. Estimates of node times were generally robust to the choice of tree prior, but some combinations of tree priors and sampling schemes led to large differences in the age estimates. In particular, the pure-birth tree prior frequently led to inaccurate estimates for data sets containing a mixture of inter- and intraspecific sampling, whereas the birth-death and skyline coalescent priors produced stable results across all scenarios. Model testing provided an adequate means of rejecting inappropriate tree priors. Our results suggest that tree priors do not strongly affect Bayesian molecular dating results in most cases, even when severely misspecified. However, the choice of tree prior can be significant for the accuracy of dating results in the case of data sets with mixed inter- and intraspecies sampling. [Bayesian phylogenetic methods; model testing; molecular dating; node time; tree prior.].

Exploring the diversity of Asian Cryptocercus (Blattodea : Cryptocercidae): species delimitation based on chromosome numbers, morphology and molecular analysis

Woodroaches from the genus Cryptocercus Scudder, 1862 are known to display low levels of morphological divergence, yet significant genetic divergence and variability in chromosome number. Compared with Cryptocercus taxa from North America, the diversity of the genus in Asia has received relatively little attention. We performed morphological and karyotypic examinations of multiple taxa from several previously unsampled mountainous areas of central and south-western China, and identified nine candidate species primarily on the basis of chromosome number. We then investigated diversity across all Asian Cryptocercus, through phylogenetic analyses of 135 COI sequences and 74 28S rRNA sequences from individuals of 28 localities, including species delimitation analysis in General Mixed Yule Coalescent (GMYC) and Automatic Barcode Gap Discovery (ABGD). Phylogenetic results indicated that individuals from the same locality constituted well supported clades. The congruence of GMYC and ABGD results were in almost perfect accord, with 28 candidate species described on the basis of karyotypes (including the nine identified in this study). We provide evidence that each valley population in the Hengduan Mountains contains a separate evolving lineage. We conclude that the principal cause of the rich Cryptocercus diversity in China has been the uplift of the Qinghai-Tibet Plateau.

Parallel evolution of mound-building and grass-feeding in Australian nasute termites

Termite mounds built by representatives of the family Termitidae are among the most spectacular constructions in the animal kingdom, reaching 6-8 m in height and housing millions of individuals. Although functional aspects of these structures are well studied, their evolutionary origins remain poorly understood. Australian representatives of the termitid subfamily Nasutitermitinae display a wide variety of nesting habits, making them an ideal group for investigating the evolution of mound building. Because they feed on a variety of substrates, they also provide an opportunity to illuminate the evolution of termite diets. Here, we investigate the evolution of termitid mound building and diet, through a comprehensive molecular phylogenetic analysis of Australian Nasutitermitinae. Molecular dating analysis indicates that the subfamily has colonized Australia on three occasions over the past approximately 20 Myr. Ancestral-state reconstruction showed that mound building arose on multiple occasions and from diverse ancestral nesting habits, including arboreal and wood or soil nesting. Grass feeding appears to have evolved from wood feeding via ancestors that fed on both wood and leaf litter. Our results underscore the adaptability of termites to ancient environmental change, and provide novel examples of parallel evolution of extended phenotypes.

The evolutionary history of Stomatopoda (Crustacea: Malacostraca) inferred from molecular data

The crustacean order Stomatopoda comprises seven superfamilies of mantis shrimps, found in coastal waters of the tropics and subtropics. These marine carnivores bear notable raptorial appendages for smashing or spearing prey. We investigated the evolutionary relationships among stomatopods using phylogenetic analyses of three mitochondrial and two nuclear markers. Our analyses recovered the superfamily Gonodactyloidea as polyphyletic, with Hemisquilla as the sister group to all other extant stomatopods. A relaxed molecular clock, calibrated by seven fossil-based age constraints, was used to date the origin and major diversification events of stomatopods. Our estimates suggest that crown-group stomatopods (Unipeltata) diverged from their closest crustacean relatives about 340 Ma (95% CRI [401–313 Ma]). We found that the specialized smashing appendage arose after the spearing appendage ∼126 Ma (95% CRI [174–87 Ma]). Ancestral state reconstructions revealed that the most recent common ancestor of extant stomatopods had eyes with six midband rows of hexagonal ommatidia. Hexagonal ommatidia are interpreted as plesiomorphic in stomatopods, and this is consistent with the malacostracan ground-plan. Our study provides insight into the evolutionary timescale and systematics of Stomatopoda, although further work is required to resolve with confidence the phylogenetic relationships among its superfamilies.

Discovery of ectosymbiotic Endomicrobium lineages associated with protists in the gut of stolotermitid termites

The genus Endomicrobium is a dominant bacterial group in the gut of lower termites, and most phylotypes are intracellular symbionts of gut protists. Here we report the discovery of Endomicrobium ectosymbionts of termite gut protists. We found that bristle-like Endomicrobium cells attached to the surface of spirotrichosomid protist cells inhabiting the termite Stolotermes victoriensis. Transmission electron microscopy revealed that a putative Endomicrobium cell likely attached to the protist surface via a protrusion from the tip of the bacterium. A phylotype, sharing 98.9% 16S rRNA sequence identity with the Endomicrobium ectosymbionts of the spirotrichosomid protists, was also found on the cell surface of the protist Trichonympha magna in the gut of the termite Porotermes adamsoni. We propose the novel species 'Candidatus Endomicrobium superficiale' for these bacteria. T. magna simultaneously harboured another Endomicrobium ectosymbiont that shared 93.5-94.2% 16S rRNA sequence identities with 'Ca. Endomicrobium superficiale'. Furthermore, Spirotrichonympha-like protists in P. adamsoni guts were associated with an Endomicrobium phylotype that possibly attached to the host flagella. A phylogenetic analysis suggested that these ectosymbiotic lineages have evolved multiple times from free-living Endomicrobium lineages and are relatively distant from the endosymbionts. Our results provide novel insights into the ecology and evolution of the Endomicrobium.

Mitochondrial Phylogenomics Resolves the Global Spread of Higher Termites, Ecosystem Engineers of the Tropics

The higher termites (Termitidae) are keystone species and ecosystem engineers. They have exceptional biomass and play important roles in decomposition of dead plant matter, in soil manipulation, and as the primary food for many animals, especially in the tropics. Higher termites are most diverse in rainforests, with estimated origins in the late Eocene (∼54 Ma), postdating the breakup of Pangaea and Gondwana when most continents became separated. Since termites are poor fliers, their origin and spread across the globe requires alternative explanation. Here, we show that higher termites originated 42-54 Ma in Africa and subsequently underwent at least 24 dispersal events between the continents in two main periods. Using phylogenetic analyses of mitochondrial genomes from 415 species, including all higher termite taxonomic and feeding groups, we inferred 10 dispersal events to South America and Asia 35-23 Ma, coinciding with the sharp decrease in global temperature, sea level, and rainforest cover in the Oligocene. After global temperatures increased, 23-5 Ma, there was only one more dispersal to South America but 11 to Asia and Australia, and one dispersal back to Africa. Most of these dispersal events were transoceanic and might have occurred via floating logs. The spread of higher termites across oceans was helped by the novel ecological opportunities brought about by environmental and ecosystem change, and led termites to become one of the few insect groups with specialized mammal predators. This has parallels with modern invasive species that have been able to thrive in human-impacted ecosystems.

Reconstructing the phylogeny of Blattodea: robust support for interfamilial relationships and major clades

Cockroaches are among the most recognizable of all insects. In addition to their role as pests, they play a key ecological role as decomposers. Despite numerous studies of cockroach phylogeny in recent decades, relationships among most major lineages are yet to be resolved. Here we examine phylogenetic relationships among cockroaches based on five genes (mitochondrial 12S rRNA, 16S rRNA, COII; nuclear 28S rRNA and histone H3), and infer divergence times on the basis of 8 fossils. We included in our analyses sequences from 52 new species collected in China, representing 7 families. These were combined with data from a recent study that examined these same genes from 49 species, resulting in a significant increase in taxa analysed. Three major lineages, Corydioidea, Blaberoidea, and Blattoidea were recovered, the latter comprising Blattidae, Tryonicidae, Lamproblattidae, Anaplectidae, Cryptocercidae and Isoptera. The estimated age of the split between Mantodea and Blattodea ranged from 204.3 Ma to 289.1 Ma. Corydioidea was estimated to have diverged 209.7 Ma (180.5-244.3 Ma 95% confidence interval [CI]) from the remaining Blattodea. The clade Blattoidea diverged from their sister group, Blaberoidea, around 198.3 Ma (173.1-229.1 Ma). The addition of the extra taxa in this study has resulted in significantly higher levels of support for a number of previously recognized groupings.

Purifying selection and concerted evolution of RNA-sensing toll-like receptors in migratory waders

Migratory birds encounter a broad range of pathogens during their journeys, making them ideal models for studying immune gene evolution. Despite the potential value of these species to immunoecology and disease epidemiology, previous studies have typically focused on their adaptive immune gene repertoires. In this study, we examined the evolution of innate immune genes in three long-distance migratory waders (order Charadriiformes). We analysed two parts of the extracellular domains of two Toll-like receptors (TLR3 and TLR7) involved in virus recognition in the Sanderling (Calidris alba), Red-necked Stint (Calidris ruficollis), and Ruddy Turnstone (Arenaria interpres). Our analysis was extended to 50 avian species for which whole-genome sequences were available, including two additional waders. We found that the inferred relationships among avian TLR3 and TLR7 do not match the whole-genome phylogeny of birds. Further analyses showed that although both loci are predominantly under purifying selection, the evolution of the extracellular domain of avian TLR3 has also been driven by episodic diversifying selection. TLR7 was found to be duplicated in all five wader species and in two other orders of birds, Cuculiformes and Passeriformes. The duplication is likely to have occurred in the ancestor of each order, and the duplicated copies appear to be undergoing concerted evolution. The phylogenetic relationships of wader TLR7 matched those of the five wader species, but that of TLR3 did not. Instead, the tree inferred from TLR3 showed potential associations with the species' ecology, including migratory behaviour and exposure to pathogens. Our study demonstrates the importance of combining immunological and ecological knowledge to understand the impact of immune gene polymorphism on the evolutionary ecology of infectious diseases.

The impacts of drift and selection on genomic evolution in insects

Genomes evolve through a combination of mutation, drift, and selection, all of which act heterogeneously across genes and lineages. This leads to differences in branch-length patterns among gene trees. Genes that yield trees with the same branch-length patterns can be grouped together into clusters. Here, we propose a novel phylogenetic approach to explain the factors that influence the number and distribution of these gene-tree clusters. We apply our method to a genomic dataset from insects, an ancient and diverse group of organisms. We find some evidence that when drift is the dominant evolutionary process, each cluster tends to contain a large number of fast-evolving genes. In contrast, strong negative selection leads to many distinct clusters, each of which contains only a few slow-evolving genes. Our work, although preliminary in nature, illustrates the use of phylogenetic methods to shed light on the factors driving rate variation in genomic evolution.

Species Delimitation and Phylogenetic Relationships in Ectobiid Cockroaches (Dictyoptera, Blattodea) from China

We collected Ectobiidae cockroach specimens from 44 locations in the south of the Yangtze valley. We obtained 297 COI sequences specimens and carried out phylogenetic and divergence dating analyses, as well as species delimitation analysis using a General Mixed Yule Coalescent (GMYC) framework. The intraspecific and interspecific sequence divergence in Ectobiidae cockroaches ranged from 0.0 to 7.0% and 4.6 to 30.8%, respectively. GMYC analysis resulted in 53 (confidence interval: 37-65) entities (likelihood ratio = 103.63) including 14 downloaded species. The COI GMYC groups partly corresponded to the ectobiid species and 52 ectobiid species were delimited successfully based on the combination of GMYC result with morphological information. We used the molecular data and 6 cockroach fossil calibrations to obtain a preliminary estimate of the timescale of ectobiid evolution. The major subfamilies in the group were found to have diverged between ~125-110 Ma, and morphospecies pairs were found to have diverged ~10 or more Ma.

A review of the status of Coptotermes (Isoptera : Rhinotermitidae) species in Australia with the description of two new small termite species from northern and eastern Australia

Integrative taxonomy, including molecular, morphological, distributional and biological data, is applied in a review of the taxonomy of the Australian species of the pest termite genus Coptotermes. The validity of the previously described species is discussed, and two new species, Coptotermes nanus, sp. nov. and Coptotermes cooloola, sp. nov., are described from the Kimberley region of Western Australia and south-east Queensland respectively. Their delimitation is based on morphological and distributional data, and the results of generalised mixed Yule-coalescent analysis of mitochondrial DNA sequence data. Images of the external view of the two new species are provided, as well as a key, based on soldier characters, for all Australian species of Coptotermes.

Multiple evolutionary origins of Australian soil-burrowing cockroaches driven by climate change in the Neogene

Parallel evolution is the independent appearance of similar derived phenotypes from similar ancestral forms. It is of key importance in the debate over whether evolution is stochastic and unpredictable, or subject to constraints that limit available phenotypic options. Nevertheless, its occurrence has rarely been demonstrated above the species level. Climate change on the Australian landmass over the last approximately 20 Myr has provided conditions conducive to parallel evolution, as taxa at the edges of shrinking mesic habitats adapted to drier biomes. Here, we investigate the phylogeny and evolution of Australian soil-burrowing and wood-feeding blaberid cockroaches. Soil burrowers (subfamily Geoscapheinae) are found in relatively dry sclerophyllous and scrubland habits, whereas wood feeders (subfamily Panesthiinae) are found in rainforest and wet sclerophyll. We sequenced and analysed mitochondrial and nuclear markers from 142 specimens, and estimated the evolutionary time scale of the two subfamilies. We found evidence for the parallel evolution of soil-burrowing taxa from wood-feeding ancestors on up to nine occasions. These transitions appear to have been driven by periods of aridification during the Miocene and Pliocene across eastern Australia. Our results provide an illuminating example of climate-driven parallel evolution among species.

Comparison of Intracellular "Ca. Endomicrobium Trichonymphae" Genomovars Illuminates the Requirement and Decay of Defense Systems against Foreign DNA

“Candidatus Endomicrobium trichonymphae” (Bacteria; Elusimicrobia) is an obligate intracellular symbiont of the cellulolytic protist genus Trichonympha in the termite gut. A previous genome analysis of “Ca. Endomicrobium trichonymphae” phylotype Rs-D17 (genomovar Ri2008), obtained from a Trichonympha agilis cell in the gut of the termite Reticulitermes speratus, revealed that its genome is small (1.1 Mb) and contains many pseudogenes; it is in the course of reductive genome evolution. Here we report the complete genome sequence of another Rs-D17 genomovar, Ti2015, obtained from a different T. agilis cell present in an R. speratus gut. These two genomovars share most intact protein-coding genes and pseudogenes, showing 98.6% chromosome sequence similarity. However, characteristic differences were found in their defense systems, which comprised restriction-modification and CRISPR/Cas systems. The repertoire of intact restriction-modification systems differed between the genomovars, and two of the three CRISPR/Cas loci in genomovar Ri2008 are pseudogenized or missing in genomovar Ti2015. These results suggest relaxed selection pressure for maintaining these defense systems. Nevertheless, the remaining CRISPR/Cas system in each genomovar appears to be active; none of the “spacer” sequences (112 in Ri2008 and 128 in Ti2015) were shared whereas the “repeat” sequences were identical. Furthermore, we obtained draft genomes of three additional endosymbiotic Endomicrobium phylotypes from different host protist species, and discovered multiple, intact CRISPR/Cas systems in each genome. Collectively, unlike bacteriome endosymbionts in insects, the Endomicrobium endosymbionts of termite-gut protists appear to require defense against foreign DNA, although the required level of defense has likely been reduced during their intracellular lives.

The Genome as an Evolutionary Timepiece

The molecular clock is a valuable and widely used tool for estimating evolutionary rates and timescales in biological research. There has been considerable progress in the theory and practice of molecular clocks over the past five decades. Although the idea of a molecular clock was originally put forward in the context of protein evolution and advanced using various biochemical techniques, it is now primarily applied to analyses of DNA sequences. An interesting but very underappreciated aspect of molecular clocks is that they can be based on genetic data other than DNA or protein sequences. For example, evolutionary timescales can be estimated using microsatellites, protein folds, and even the extent of recombination. These genome features hold great potential for molecular dating, particularly in cases where nucleotide sequences might be uninformative or unreliable. Here we present an outline of the different genetic data types that have been used for molecular dating, and we describe the features that good molecular clocks should possess. We hope that our article inspires further work on the genome as an evolutionary timepiece.

Novel Borrelia species detected in echidna ticks, Bothriocroton concolor, in Australia

BACKGROUND

To date, little has been documented about microorganisms harboured within Australian native ticks or their pathogenic potential. Recently, a Borrelia sp. related to the Relapsing Fever (RF) group was identified in a single tick removed from a wild echidna (Tachyglossus aculeatus). The present study investigated the presence of Borrelia in 97 Bothriocroton concolor ticks parasitizing echidnas in Queensland, New South Wales, and Victoria, Australia, using nested PCR with Borrelia-specific primers targeting the 16S rRNA (16S) and flaB genes.

RESULTS

Borrelia-specific PCR assays confirmed the presence of a novel Borrelia sp. related to the RF and reptile-associated (REP) spirochaetes in 38 (39 %) B. concolor ticks. This novel Borrelia sp. was identified in 41 % of the B. concolor ticks in Queensland and New South Wales, but not in any ticks from Victoria. The resulting flaB sequences (407 bp) were 88 and 86 % similar to the flaB sequences from Borrelia turcica and Borrelia hermsii, respectively. Of the ticks confirmed as Borrelia-positive following the flaB assay, 28 were positive with the 16S assay. Phylogenetic analysis of the 16S sequences (1097 bp) suggests that these sequences belong to a novel Borrelia sp., which forms a unique monophyletic clade that is similar to, but distinct from, RF Borrelia spp. and REP-associated Borrelia spp.

CONCLUSIONS

We conclude that the novel Borrelia sp. identified in this study does not belong to the Borrelia burgdorferi (sensu lato) complex, and that the phylogenetic analysis of the partial 16S gene sequences suggests it forms a unique monophyletic cluster in the genus Borrelia, potentially forming a fourth major group in this genus associated with monotremes in Australia. However, a thorough molecular characterisation will be required to confirm the phylogenetic position of this unique Borrelia sp. The zoonotic potential and pathogenic consequences of this novel Borrelia sp. are unknown at the current time.

Alternative migratory locust phenotypes are associated with differences in the expression of genes encoding the methylation machinery

Despite the importance of locust density-dependent polyphenism as a model system for understanding phenotypic plasticity, there is still much to be learnt about its underlying molecular control. Here we describe the first investigation into the expression of genes encoding the DNA methylation machinery in the migratory locust (Locusta migratoria). We show that the alternative solitarious and gregarious phenotypic states induced by different locust rearing densities are associated with significant differences in the expression of the target genes DNA methyltransferase 1, DNA methyltransferase 2 and methyl-CpG-binding domain protein 2/3. This variation was most pronounced in the embryos of solitarious vs. gregarious mothers. We mapped the embryonic methylation profiles of several intragenic regions and a Long Interspersed Nuclear Element (LINE), each of which is known to be differentially expressed between alternative locust phenotypes or has been directly implicated in phase change. LmI and three genes, adenyl cyclase-associated binding protein 2, choline kinase alpha-like and henna, were methylated. Our results set the stage for future studies investigating the specific role of DNA methylation in the maternal transfer of migratory locust phase polyphenism.

Oceanic dispersal, vicariance and human introduction shaped the modern distribution of the termites Reticulitermes, Heterotermes and Coptotermes

Reticulitermes, Heterotermes and Coptotermes form a small termite clade with partly overlapping distributions. Although native species occur across all continents, the factors influencing their distribution are poorly known. Here, we reconstructed the historical biogeography of these termites using mitochondrial genomes of species collected on six continents. Our analyses showed that Reticulitermes split from Heterotermes + Coptotermesat 59.5 Ma (49.9-69.5 Ma 95% CI), yet the oldest split within Reticulitermes(Eurasia and North America) is 16.1 Ma (13.4-19.5 Ma) and the oldest split within Heterotermes + Coptotermesis 36.0 Ma (33.9-40.5 Ma). We detected 14 disjunctions between biogeographical realms, all of which occurred within the last 34 Ma, not only after the break-up of Pangaea, but also with the continents in similar to current positions. Land dispersal over land bridges explained four disjunctions, oceanic dispersal by wood rafting explained eight disjunctions, and human introduction was the source of two recent disjunctions. These wood-eating termites, therefore, appear to have acquired their modern worldwide distribution through multiple dispersal processes, with oceanic dispersal and human introduction favoured by the ecological traits of nesting in wood and producing replacement reproductives.

Widespread convergence in toxin resistance by predictable molecular evolution

The question about whether evolution is unpredictable and stochastic or intermittently constrained along predictable pathways is the subject of a fundamental debate in biology, in which understanding convergent evolution plays a central role. At the molecular level, documented examples of convergence are rare and limited to occurring within specific taxonomic groups. Here we provide evidence of constrained convergent molecular evolution across the metazoan tree of life. We show that resistance to toxic cardiac glycosides produced by plants and bufonid toads is mediated by similar molecular changes to the sodium-potassium-pump (Na(+)/K(+)-ATPase) in insects, amphibians, reptiles, and mammals. In toad-feeding reptiles, resistance is conferred by two point mutations that have evolved convergently on four occasions, whereas evidence of a molecular reversal back to the susceptible state in varanid lizards migrating to toad-free areas suggests that toxin resistance is maladaptive in the absence of selection. Importantly, resistance in all taxa is mediated by replacements of 2 of the 12 amino acids comprising the Na(+)/K(+)-ATPase H1-H2 extracellular domain that constitutes a core part of the cardiac glycoside binding site. We provide mechanistic insight into the basis of resistance by showing that these alterations perturb the interaction between the cardiac glycoside bufalin and the Na(+)/K(+)-ATPase. Thus, similar selection pressures have resulted in convergent evolution of the same molecular solution across the breadth of the animal kingdom, demonstrating how a scarcity of possible solutions to a selective challenge can lead to highly predictable evolutionary responses.

Biogeographic calibrations for the molecular clock

Molecular estimates of evolutionary timescales have an important role in a range of biological studies. Such estimates can be made using methods based on molecular clocks, including models that are able to account for rate variation across lineages. All clock models share a dependence on calibrations, which enable estimates to be given in absolute time units. There are many available methods for incorporating fossil calibrations, but geological and climatic data can also provide useful calibrations for molecular clocks. However, a number of strong assumptions need to be made when using these biogeographic calibrations, leading to wide variation in their reliability and precision. In this review, we describe the nature of biogeographic calibrations and the assumptions that they involve. We present an overview of the different geological and climatic events that can provide informative calibrations, and explain how such temporal information can be incorporated into dating analyses.

Inhibition of the endosymbiont "Candidatus Midichloria mitochondrii" during 16S rRNA gene profiling reveals potential pathogens in Ixodes ticks from Australia

Background

The Australian paralysis tick (Ixodes holocyclus) is of significant medical and veterinary importance as a cause of dermatological and neurological disease, yet there is currently limited information about the bacterial communities harboured by these ticks and the risk of infectious disease transmission to humans and domestic animals. Ongoing controversy about the presence of Borrelia burgdorferi sensu lato (the aetiological agent of Lyme disease) in Australia increases the need to accurately identify and characterise bacteria harboured by I. holocyclus ticks.

Methods

Universal PCR primers were used to amplify the V1-2 hyper-variable region of bacterial 16S rRNA genes present in DNA samples from I. holocyclus and I. ricinus ticks, collected in Australia and Germany respectively. The 16S amplicons were purified, sequenced on the Ion Torrent platform, and analysed in USEARCH, QIIME, and BLAST to assign genus and species-level taxonomy. Initial analysis of I. holocyclus and I. ricinus identified that > 95 % of the 16S sequences recovered belonged to the tick intracellular endosymbiont “Candidatus Midichloria mitochondrii” (CMM). A CMM-specific blocking primer was designed that decreased CMM sequences by approximately 96 % in both tick species and significantly increased the total detectable bacterial diversity, allowing identification of medically important bacterial pathogens that were previously masked by CMM.

Results

Borrelia burgdorferi sensu lato was identified in German I. ricinus, but not in Australian I. holocyclus ticks. However, bacteria of medical significance were detected in I. holocyclus ticks, including a Borrelia relapsing fever group sp., Bartonella henselae, novel “Candidatus Neoehrlichia” spp., Clostridium histolyticum, Rickettsia spp., and Leptospira inadai.

Conclusions

Abundant bacterial endosymbionts, such as CMM, limit the effectiveness of next-generation 16S bacterial community profiling in arthropods by masking less abundant bacteria, including pathogens. Specific blocking primers that inhibit endosymbiont 16S amplification during PCR are an effective way of reducing this limitation. Here, this strategy provided the first evidence of a relapsing fever Borrelia sp. and of novel “Candidatus Neoehrlichia” spp. in Australia. Our results raise new questions about tick-borne pathogens in I. holocyclus ticks.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-015-0958-3) contains supplementary material, which is available to authorized users.

Metabolomic profiling of 13C-labelled cellulose digestion in a lower termite: insights into gut symbiont function

Termites consume an estimated 3–7 billion tonnes of lignocellulose annually, a role in nature which is unique for a single order of invertebrates. Their food is digested with the help of microbial symbionts, a relationship that has been recognized for 200 years and actively researched for at least a century. Although DNA- and RNA-based approaches have greatly refined the details of the process and the identities of the participants, the allocation of roles in space and time remains unclear. To resolve this issue, a pioneer study is reported using metabolomics to chart the in situ catabolism of ¹³C-cellulose fed to the dampwood species Hodotermopsis sjostedti. The results confirm that the secretion of endogenous cellulases by the host may be significant to the digestive process and indicate that a major contribution by hindgut bacteria is phosphorolysis of cellodextrins or cellobiose. This study provides evidence that essential amino acid acquisition by termites occurs following the lysis of microbial tissue obtained via proctodaeal trophallaxis.

The evolutionary history of termites as inferred from 66 mitochondrial genomes

Termites have colonized many habitats and are among the most abundant animals in tropical ecosystems, which they modify considerably through their actions. The timing of their rise in abundance and of the dispersal events that gave rise to modern termite lineages is not well understood. To shed light on termite origins and diversification, we sequenced the mitochondrial genome of 48 termite species and combined them with 18 previously sequenced termite mitochondrial genomes for phylogenetic and molecular clock analyses using multiple fossil calibrations. The 66 genomes represent most major clades of termites. Unlike previous phylogenetic studies based on fewer molecular data, our phylogenetic tree is fully resolved for the lower termites. The phylogenetic positions of Macrotermitinae and Apicotermitinae are also resolved as the basal groups in the higher termites, but in the crown termitid groups, including Termitinae + Syntermitinae + Nasutitermitinae + Cubitermitinae, the position of some nodes remains uncertain. Our molecular clock tree indicates that the lineages leading to termites and Cryptocercus roaches diverged 170 Ma (153-196 Ma 95% confidence interval [CI]), that modern Termitidae arose 54 Ma (46-66 Ma 95% CI), and that the crown termitid group arose 40 Ma (35-49 Ma 95% CI). This indicates that the distribution of basal termite clades was influenced by the final stages of the breakup of Pangaea. Our inference of ancestral geographic ranges shows that the Termitidae, which includes more than 75% of extant termite species, most likely originated in Africa or Asia, and acquired their pantropical distribution after a series of dispersal and subsequent diversification events.

The origins and radiation of Australian Coptotermes termites: from rainforest to desert dwellers

The termite genus Coptotermes (Rhinotermitidae) is found in Asia, Africa, Central/South America and Australia, with greatest diversity in Asia. Some Coptotermes species are amongst the world's most damaging invasive termites, but the genus is also significant for containing the most sophisticated mound-building termites outside the family Termitidae. These mound-building Coptotermes occur only in Australia. Despite its economic and evolutionary significance, the biogeographic history of the genus has not been well investigated, nor has the evolution of the Australian mound-building species. We present here the first phylogeny of the Australian Coptotermes to include representatives from all described species. We combined our new data with previously generated data to estimate the first phylogeny to include representatives from all continents where the genus is found. We also present the first estimation of divergence dates during the evolution of the genus. We found the Australian Coptotermes to be monophyletic and most closely related to the Asian Coptotermes, with considerable genetic diversity in some Australian taxa possibly representing undescribed species. The Australian mound-building species did not form a monophyletic clade. Our ancestral state reconstruction analysis indicated that the ancestral Australian Coptotermes was likely to have been a tree nester, and that mound-building behaviour has arisen multiple times. The Australian Coptotermes were found to have diversified ∼13million years ago, which plausibly matches with the narrowing of the Arafura Sea allowing Asian taxa to cross into Australia. The first diverging Coptotermes group was found to be African, casting doubt on the previously raised hypothesis that the genus has an Asian origin.

Perioperative complications after noncardiac surgery in patients with insertion of second-generation drug-eluting stents

The perioperative outcomes of noncardiac surgery in patients who have received second-generation drug-eluting stents (DESs) have received limited study. We reviewed the medical records of 1,748 consecutive patients who received DES at our institution (1,789 procedures) from January 1, 2009, to July 1, 2012, to determine the outcomes of subsequent noncardiac surgery. During a median follow-up of 43 months, 221 patients underwent 345 noncardiac surgeries (138 low risk, 130 intermediate risk, and 77 high risk), of which 278 were in patients with previous second-generation DES implantation. The incidence of noncardiac surgery in patients with previous second-generation DES implantation was 4.5% at 1 year, 11.6% at 2 years, and 15.2% at 3 years. The mean time from stent implantation to surgery was 21±12.9 months. Mean age was 66±8 years, 99% were men, and 11% had a perioperative complication, including 5.8% major bleeding, 2.5% acute kidney injury, 2.2% major adverse cardiac event, and 1.4% stroke. Perioperative stent thrombosis occurred in 2 patients (0.7%, 95% confidence interval 0.2% to 2.6%): 1 patient had received a DES 14 months before surgery and had stent thrombosis on the day of surgery and the other had a DES implanted 21 months before surgery and developed stent thrombosis the day after surgery. In conclusion, the incidence of perioperative complications with noncardiac surgery after second-generation DES implantation was 11% and consisted mainly of bleeding (5.8%). The incidence of definite stent thrombosis was 0.7%.

The impact of modelling rate heterogeneity among sites on phylogenetic estimates of intraspecific evolutionary rates and timescales

Phylogenetic analyses of DNA sequence data can provide estimates of evolutionary rates and timescales. Nearly all phylogenetic methods rely on accurate models of nucleotide substitution. A key feature of molecular evolution is the heterogeneity of substitution rates among sites, which is often modelled using a discrete gamma distribution. A widely used derivative of this is the gamma-invariable mixture model, which assumes that a proportion of sites in the sequence are completely resistant to change, while substitution rates at the remaining sites are gamma-distributed. For data sampled at the intraspecific level, however, biological assumptions involved in the invariable-sites model are commonly violated. We examined the use of these models in analyses of five intraspecific data sets. We show that using 6-10 rate categories for the discrete gamma distribution of rates among sites is sufficient to provide a good approximation of the marginal likelihood. Increasing the number of gamma rate categories did not have a substantial effect on estimates of the substitution rate or coalescence time, unless rates varied strongly among sites in a non-gamma-distributed manner. The assumption of a proportion of invariable sites provided a better approximation of the asymptotic marginal likelihood when the number of gamma categories was small, but had minimal impact on estimates of rates and coalescence times. However, the estimated proportion of invariable sites was highly susceptible to changes in the number of gamma rate categories. The concurrent use of gamma and invariable-site models for intraspecific data is not biologically meaningful and has been challenged on statistical grounds; here we have found that the assumption of a proportion of invariable sites has no obvious impact on Bayesian estimates of rates and timescales from intraspecific data.

Periprocedural myocardial injury in chronic total occlusion percutaneous interventions: a systematic cardiac biomarker evaluation study

OBJECTIVES

This study sought to evaluate the incidence, correlates, and clinical implications of periprocedural myocardial injury (PMI) during percutaneous coronary intervention (PCI) of chronic total occlusions (CTO).

BACKGROUND

The risk of PMI during CTO PCI may be underestimated because systematic cardiac biomarker measurement was not performed in published studies.

METHODS

We retrospectively examined PMI among 325 consecutive CTO PCI performed at our institution between 2005 and 2012. Creatine kinase-myocardial band fraction and troponin were measured before PCI and 8 to 12 h and 18 to 24 h after PCI in all patients. PMI was defined as creatine kinase-myocardial band increase ≥ 3 x the upper limit of normal. Major adverse cardiac events during mid-term follow-up were evaluated.

RESULTS

Mean age was 64 ± 8 years. The retrograde approach was used in 26.8% of all procedures. The technical and procedural success was 77.8% and 76.6%, respectively. PMI occurred in 28 patients (8.6%, 95% confidence intervals: 5.8% to 12.2%), with symptomatic ischemia in 7 of those patients. The incidence of PMI was higher in patients treated with the retrograde than the antegrade approach (13.8% vs. 6.7%, p = 0.04). During a median follow-up of 2.3 years, compared with patients without PMI, those with PMI had a higher incidence of major adverse cardiac events (hazard ratio [HR]: 2.25, p = 0.006). Patients with only asymptomatic PMI also had a higher incidence of major adverse cardiac events on follow-up (HR: 2.26, p = 0.013).

CONCLUSIONS

Systematic measurement of cardiac biomarkers post-CTO PCI demonstrates that PMI occurs in 8.6% of patients, is more common with the retrograde approach, and is associated with worse subsequent clinical outcomes during mid-term follow-up.

Prevalence and management of coronary chronic total occlusions in a tertiary Veterans Affairs hospital

OBJECTIVES

We sought to determine the contemporary prevalence and management of coronary chronic total occlusions (CTO) in a veteran population.

BACKGROUND

The prevalence and management of CTOs in various populations has received limited study.

METHODS

We collected clinical and angiographic data in consecutive patients that underwent coronary angiography at our institution between January 2011 and December 2012. Coronary artery disease (CAD) was defined as ≥50% diameter stenosis in ≥1 coronary artery. CTO was defined as total coronary artery occlusion of ≥3 month duration.

RESULTS

Among 1,699 patients who underwent angiography during the study period, 20% did not have CAD, 20% had CAD and prior coronary artery bypass graft surgery (CABG), and 60% had CAD but no prior CABG. The prevalence of CTO among CAD patients with and without prior CABG was 89 and 31%, respectively. Compared to patients without CTO, CTO patients had more co-morbidities, more extensive CAD and were more frequently referred for CABG. Percutaneous coronary intervention (PCI) to any vessel was performed with similar frequency in patients with and without CTO (50% vs. 53%). CTO PCI was performed in 30% of patients without and 15% of patients with prior CABG with high technical (82 and 75%, respectively) and procedural success rates (80 and 73%, respectively).

CONCLUSIONS

In a contemporary veteran population, coronary CTOs are highly prevalent and are associated with more extensive co-morbidities and higher likelihood for CABG referral. PCI was equally likely to be performed in patients with and without CTO.

Maintenance of essential amino acid synthesis pathways in the Blattabacterium cuenoti symbiont of a wood-feeding cockroach

In addition to harbouring intestinal symbionts, some animal species also possess intracellular symbiotic microbes. The relative contributions of gut-resident and intracellular symbionts to host metabolism, and how they coevolve are not well understood. Cockroaches and the termite Mastotermes darwiniensis present a unique opportunity to examine the evolution of spatially separated symbionts, as they harbour gut symbionts and the intracellular symbiont Blattabacterium cuenoti. The genomes of B. cuenoti from M. darwiniensis and the social wood-feeding cockroach Cryptocercus punctulatus are each missing most of the pathways for the synthesis of essential amino acids found in the genomes of relatives from non-wood-feeding hosts. Hypotheses to explain this pathway degradation include: (i) feeding on microbes present in rotting wood by ancestral hosts; (ii) the evolution of high-fidelity transfer of gut microbes via social behaviour. To test these hypotheses, we sequenced the B. cuenoti genome of a third wood-feeding species, the phylogenetically distant and non-social Panesthia angustipennis. We show that host wood-feeding does not necessarily lead to degradation of essential amino acid synthesis pathways in B. cuenoti, and argue that ancestral high-fidelity transfer of gut microbes best explains their loss in strains from M. darwiniensis and C. punctulatus.

Construction and characterization of normalized cDNA libraries by 454 pyrosequencing and estimation of DNA methylation levels in three distantly related termite species

In termites, division of labor among castes, categories of individuals that perform specialized tasks, increases colony-level productivity and is the key to their ecological success. Although molecular studies on caste polymorphism have been performed in termites, we are far from a comprehensive understanding of the molecular basis of this phenomenon. To facilitate future molecular studies, we aimed to construct expressed sequence tag (EST) libraries covering wide ranges of gene repertoires in three representative termite species, Hodotermopsis sjostedti, Reticulitermes speratus and Nasutitermes takasagoensis. We generated normalized cDNA libraries from whole bodies, except for guts containing microbes, of almost all castes, sexes and developmental stages and sequenced them with the 454 GS FLX titanium system. We obtained >1.2 million quality-filtered reads yielding >400 million bases for each of the three species. Isotigs, which are analogous to individual transcripts, and singletons were produced by assembling the reads and annotated using public databases. Genes related to juvenile hormone, which plays crucial roles in caste differentiation of termites, were identified from the EST libraries by BLAST search. To explore the potential for DNA methylation, which plays an important role in caste differentiation of honeybees, tBLASTn searches for DNA methyltransferases (dnmt1, dnmt2 and dnmt3) and methyl-CpG binding domain (mbd) were performed against the EST libraries. All four of these genes were found in the H. sjostedti library, while all except dnmt3 were found in R. speratus and N. takasagoensis. The ratio of the observed to the expected CpG content (CpG O/E), which is a proxy for DNA methylation level, was calculated for the coding sequences predicted from the isotigs and singletons. In all of the three species, the majority of coding sequences showed depletion of CpG O/E (less than 1), and the distributions of CpG O/E were bimodal, suggesting the presence of DNA methylation.

Evidence for Permo-Triassic colonization of the deep sea by isopods

The deep sea is one of the largest ecosystems on Earth and is home to a highly diverse fauna, with polychaetes, molluscs and peracarid crustaceans as dominant groups. A number of studies have proposed that this fauna did not survive the anoxic events that occurred during the Mesozoic Era. Accordingly, the modern fauna is thought to be relatively young, perhaps having colonized the deep sea after the Eocene/Oligocene boundary. To test this hypothesis, we performed phylogenetic analyses of nuclear ribosomal 18S and 28S and mitochondrial cytochrome oxidase I and 16S sequences from isopod crustaceans. Using a molecular clock calibrated with multiple isopod fossils, we estimated the timing of deep-sea colonization events by isopods. Our results show that some groups have an ancient origin in the deep sea, with the earliest estimated dates spanning 232-314 Myr ago. Therefore, anoxic events at the Permian-Triassic boundary and during the Mesozoic did not cause the extinction of all the deep-sea fauna; some species may have gone extinct while others survived and proliferated. The monophyly of the 'munnopsid radiation' within the isopods suggests that the ancestors of this group evolved in the deep sea and did not move to shallow-water refugia during anoxic events.

The aging of biomedical research in the United States

In the past 30 years, the average age of biomedical researchers has steadily increased. The average age of an investigator at the National Institutes of Health (NIH) rose from 39 to 51 between 1980 and 2008. The aging of the biomedical workforce was even more apparent when looking at first-time NIH grantees. The average age of a new investigator was 42 in 2008, compared to 36 in 1980. To determine if the rising barriers at NIH for entry in biomedical research might impact innovative ideas and research, we analyzed the research and publications of Nobel Prize winners from 1980 to 2010 to assess the age at which their pioneering research occurred. We established that in the 30-year period, 96 scientists won the Nobel Prize in medicine or chemistry for work related to biomedicine, and that their groundbreaking research was conducted at an average age of 41—one year younger than the average age of a new investigator at NIH. Furthermore, 78% of the Nobel Prize winners conducted their research before the age of 51, the average age of an NIH principal investigator. This suggested that limited access to NIH might inhibit research potential and novel projects, and could impact biomedicine and the next generation scientists in the United States.

Phylogenomic evidence for the presence of a flagellum and cbb(3) oxidase in the free-living mitochondrial ancestor

The initiation of the intracellular symbiosis that would give rise to mitochondria and eukaryotes was a major event in the history of life on earth. Hypotheses to explain eukaryogenesis fall into two broad and competing categories: those proposing that the host was a phagocytotic proto-eukaryote that preyed upon the free-living mitochondrial ancestor (hereafter FMA), and those proposing that the host was an archaebacterium that engaged in syntrophy with the FMA. Of key importance to these hypotheses are whether the FMA was motile or nonmotile, and the atmospheric conditions under which the FMA thrived. Reconstructions of the FMA based on genome content of Rickettsiales representatives-generally considered to be the closest living relatives of mitochondria-indicate that it was nonmotile and aerobic. We have sequenced the genome of Candidatus Midichloria mitochondrii, a novel and phylogenetically divergent member of the Rickettsiales. We found that it possesses unique gene sets found in no other Rickettsiales, including 26 genes associated with flagellar assembly, and a cbb(3)-type cytochrome oxidase. Phylogenomic analyses show that these genes were inherited in a vertical fashion from an ancestral α-proteobacterium, and indicate that the FMA possessed a flagellum, and could undergo oxidative phosphorylation under both aerobic and microoxic conditions. These results indicate that the FMA played a more active and potentially parasitic role in eukaryogenesis than currently appreciated and provide an explanation for how the symbiosis could have evolved under low levels of oxygen.

Ants and termites increase crop yield in a dry climate

Agricultural intensification has increased crop yields, but at high economic and environmental cost. Harnessing ecosystem services of naturally occurring organisms is a cheaper but under-appreciated approach, because the functional roles of organisms are not linked to crop yields, especially outside the northern temperate zone. Ecosystem services in soil come from earthworms in these cooler and wetter latitudes; what may fulfill their functional role in agriculture in warmer and drier habitats, where they are absent, is unproven. Here we show in a field experiment that ants and termites increase wheat yield by 36% from increased soil water infiltration due to their tunnels and improved soil nitrogen. Our results suggest that ants and termites have similar functional roles to earthworms, and that they may provide valuable ecosystem services in dryland agriculture, which may become increasingly important for agricultural sustainability in arid climates.

The presence of earthworms is known to enhance the quality and moisture of soil in cool and wet climates. Evans et al. show that termites and ants can improve soil quality in warmer and drier climates—their presence results in elevated water infiltration and nitrogen content, leading to increased wheat yields.

Evidence for genetically influenced caste determination in phylogenetically diverse species of the termite genus Reticulitermes

A number of social insect species have recently been shown to have genetically influenced caste determination (GCD), challenging the conventional view that caste determination should be strictly environmental. To date, GCD has been found in phylogenetically isolated species; examples of GCD being present in multiple species of a genus are lacking. Through crossing experiments of neotenic (juvenile) reproductives, we have recently provided the first evidence for a royal versus worker GCD in the termite Reticulitermes speratus. To elucidate whether this system is more widespread, we performed crossing experiments using three additional Reticulitermes species. Offspring caste and sex ratios were found to be highly similar to those found previously in R. speratus, raising the possibility that GCD was present in an ancestral lineage of Reticulitermes, and subsequently maintained throughout several episodes of speciation.

Permanent Genetic Resources added to Molecular Ecology Resources Database 1 June 2010 - 31 July 2010

This article documents the addition of 205 microsatellite marker loci to the Molecular Ecology Resources Database. Loci were developed for the following species: Bagassa guianensis, Bulweria bulwerii, Camelus bactrianus, Chaenogobius annularis, Creontiades dilutus, Diachasmimorpha tryoni, Dioscorea alata, Euhrychiopsis lecontei, Gmelina arborea, Haliotis discus hannai, Hirtella physophora, Melanaphis sacchari, Munida isos, Thaumastocoris peregrinus and Tuberolachnus salignus. These loci were cross-tested on the following species: Halobaena caerulea, Procellaria aequinoctialis, Oceanodroma monteiroi, Camelus ferus, Creontiades pacificus, Dioscorea rotundata, Dioscorea praehensilis, Dioscorea abyssinica, Dioscorea nummularia, Dioscorea transversa, Dioscorea esculenta, Dioscorea pentaphylla, Dioscorea trifida, Hirtella bicornis, Hirtella glandulosa, Licania alba, Licania canescens, Licania membranaceae, Couepia guianensis and 7 undescribed Thaumastocoris species.

Nature versus nurture in social insect caste differentiation

Recent evidence for genetic effects on royal and worker caste differentiation from diverse social insect taxa has put an end to the view that these phenotypes stem solely from a developmental switch controlled by environmental factors. Instead, the relative influences of genotypic and environmental effects on caste vary among species, ranging from largely environmentally controlled phenotypes to almost purely genetic systems. Disentangling the selective forces that generate variation for caste predisposition will require characterizing the genetic mechanisms underlying this variation, and identifying particular life-history strategies and kin structures associated with strong genetic effects on caste.